- •Лекция №1 Линейные ду в частных производных второго порядка.
- •Классификация уравнений в частных производных 2го порядка.
- •Тип {m,0,0} {0,m,0} – эллиптический тип. Все собственные значения отличны от нуля и одного знака.
- •Приведение к каноническому виду.
- •Характеристики.
- •Приведение к каноническому виду уравнения 2го порядка с двумя независимыми переменными.
- •Лекция №2 Основные уравнения математической физики. Уравнения колебаний
- •Лекция №3 Уравнения Максвелла.
- •Задача Коши для одномерного волнового уравнения.
- •Физическая интерпретация решения одномерного волнового уравнения.
- •Решение задачи Коши для волнового уравнения.
- •Лекция №4 Устойчивость задачи Коши к входным данным.
- •Устойчивость неоднородного волнового уравнения задачи Коши к малым изменениям неоднородности (к правой части).
- •Методы решения краевых задач. Метод Фурье.
- •Гиперболическая краевая задача.
- •Решение неоднородной задачи методом Фурье.
- •Свободные колебания круглой мембраны.
- •Уравнения параболического типа.
- •Метод разделения переменных для решения задач уравнений теплопроводности. Метод Фурье.
- •Задача Коши для параболического уравнения.
- •Физический смысл фундаментального решения.
- •Эллиптические уравнения.
- •Сингулярное решение уравнения Лапласа.
- •Частные случаи сингулярных решений.
- •Функция Грина оператора Лапласа.
- •Сущность метода функции Грина.
- •Теория потенциала.
- •Интеграл Гаусса.
- •Интегральные уравнения.
- •Интегральное уравнение с вырожденным ядром.
- •Интегральное уравнение с непрерывным ядром.
- •Решение интегральных уравнений методом последовательных приближений.
- •Решение уравнений Фредгольма заменой интеграла суммой.
- •Классические и обобщённые решения.
- •Метод взвешенных невязок.
- •Связь метода Галёркина и Релея-Ритца.
Лекция №3 Уравнения Максвелла.
Е, Н – электрическая и магнитная напряжённости. Е=(Е1,Е2,Е3), Н=(Н1,Н2,Н3). В качестве переменных могут быть х=(х1,х2,х3), t – время.
div(E)=
div(H)=0
rot E = –(H)/t
rot H = (E)/t + J + j
D=E – индукция электрического поля, В=Н – индукция магнитного поля, j=E – ток, – проводимость, – изолятор пространственный. Для решения уравнений первого порядка не существует метода решения. Rot rot A = grad div A – div grad A — основная формула для решений уравнений Максвелла. С помощью неё можно решить уравнения 3), 4). Rot rot H = /t( rot E) + rot J;
Grad div H – div grad H = /t((–/t[H])) + rot J +. – гиперболический тип. Аналогично с вектором Е. Если в среде присутствует проводящий элемент (медь,), то добавляется дополнительное слагаемое. + 2H/t2.
Задача Коши для одномерного волнового уравнения.
Волновое уравнение – это уравнение колебания струны.
(1) 2U/t2 – a2 2U/x2 = 0; –<x<, t>0 AUxx+2BUxy+CUyy+ (…) = 0. =В2 – АС. В=0, А=1, С = –а2, = а2 => уравнение гиперболического типа. Существует и вторая каноническая форма для канонического вида: dx/dt = a, dx/dt = –a – характеристики. dx/dt = (ВВ2 – АС)/А.
x – at=C1, x + at = C2 => = x – at, = x + at – новые переменные (4). Выполняя замену, получим: 2U/ = 0 (5). Начальные условия: –<x<, t>0 => U(x,t)|t=0=(x), Ut(x,t)|t=0 = (x) (2). Пусть решение этой задачи (1) есть некоторая функция: (3) U(x,t) = f(x – at) + g(x + at). Если f и g С2(R), то и решение U(x,t)C2(R). Докажем, что это правильно:
Из уравнения (5) => /(U/) = 0;
U/(, ) = V(, ) – новое обозначение. Тогда (5) примет вид: (6) V/ = d/d(V|=const) = 0, V|=const – не зависит от , т.е. V(, ) = C() (7).
U/(, ) = V(, ) => d/d(U| = const) = C(). Проинтегрируем обыкновенное дифференциальное уравнение и получим: U| = const = C()d + C1()
U = f() + g()
U = f(x – at) + g(x + at). Решение волнового уравнения для задачи Коши есть две волны.
Физическая интерпретация решения одномерного волнового уравнения.
f(x – at).
Функция U(x, t) – есть отклонение в точке х в момент времени t. Рассмотрим некоторую фиксированную точку х0. Пусть из точки х0 в положительном направлении оси х в момент времени t=0 начинает двигаться наблюдатель со скоростью, равной а. Тогда, в момент времени t1 он окажется в точке х1 = х0 + аt1. Величина отклонения, которую наблюдатель будет видеть в точке х1 в момент времени t1 будет U(x1, t1) = f(x1 – at1) = f(x0), т.е. будет той же самой, что и в t0. Начальный профиль будет двигаться со скоростью а в положительном направлении оси Ох. U(x,t)|t=0 = (x) – начальный профиль. f(x – at) называется прямой бегущей волной. Аналогичные рассуждения относительно g.
g(x + at) – называется обратная бегущая волна. В отрицательном направлении оси Ох.
Вывод: любое решение U(x,t) волнового уравнения (1) может быть представлено суперпозицией 2х волн, прямой и обратной.
Решение задачи Коши для волнового уравнения.
2U/t2 – a22U/x2 = f(x,t), –<x< t>0
U(x,t)|t=0 = (x), Ut(x,t)|t=0 = (x) (2)
I. Задача Коши для однородного уравнения с заданными начальными условиями.
2V/t2 – a2 2V/x2 =0; V|t=0 = (x), Vt|t=0 = (x)
II. Задача Коши для неоднородного уравнения, но с однородными начальными условиями.
2W/t2 – a2 2W/x2 = f(x,t); W|t=0 = 0, Wt|t=0 = 0
I: Решение этой задачи будет таким: V(x,t) = F(x+at) + (x–at) (3). Среди всех таких решений будем искать такое, которое бы удовлетворяло начальным условиям. V(x,0) = (x) = (x)+F(x), Vt(x,0) = (x) = –a(x) + aF(x) (4). Второе уравнение из (4) проинтегрируем и получим два уравнения для определения 2х функций F и :
(х) + F(x) =(x), –Ф(x) + F(x) = 1/a{x0,x}(z)dz + C (5).
Из (5) => F(x)=(x)/2+1/(2a){x0,x}(z)dz+C/2, Ф(x)=(x)/2 – 1/(2a){x0,x}(z)dz – C/2 (6). Из (6)(3) с учётом, что х – аt и х + аt: U(x,t) = 1/2[(x–at) + (x+at)] + 1/(2a){x–at, x+at}(z)dz (7) – формула Даламбера – описывает свободные колебания бесконечной струны при заданных начальных возмущениях. Получили явный вид решения I.
Вывод: Предположив существование решения задачи Коши получили формулу (7), значит, это решение единственно. Если функция (х) обладает производной 1го и 2го порядка, а (х) производную 1го порядка, то формула (7) даёт искомое решение задачи Коши и доказать это можно подстановкой U(x,t) в уравнение I. Построив решение задачи Коши, тем самым доказали его существование. Такой метод построения решения называется методом характеристик.
II. Будем строить вспомогательную функцию (x,t,), которая удовлетворяет однородному уравнению: tt = a2xx и |t= = 0; t|t= = f(x,) – это есть задача I, которая только что решена. w(x,t) = {0,t}(x,t,)d (8) – решение задачи II. Проверим: по правилам дифференцирования интеграла с переменным верхним пределом по параметру, найдём: wt(x,t) = |t= (0, начал. усл.)+ {0,t}t(x,t,)d (9). Из (8), (9) => w(x,t) удовлетворяет нулевым начальным условиям: w(x,t)|t=0 = 0, wt(x,t)|t=0 = 0. Продифференцируем (9) ещё раз и используем начальные условия: wtt = t|t= + {0,t}tt(x,t,)d = f(x,t) + {0,t}tt(x,t,)d (10). wxx – ? wxx = {0,t}xx(x,t,)d (11). a2 wxx = {0,t}a2 tt(x,t,)d. Из (10), (11) => wtt – a2 wxx = f(x,t) + {0,t}[ tt – a2xx](0)d. Значит, функция w(x,t) представленная в виде (8) есть решение задачи II. Найдём её явное выражение: для нахождения функции (x,t,) можно воспользоваться формулой Даламбера: (x,t,) = 1/(2a){x–a(t–), x+a(t–)}f(z,)dz (12). Из (12)(8) => w(x,t) = 1/(2a){0,t}{x–a(t–), x+a(t–)}f(z,)dzd (13). Решение задачи: U(x,t) = V(x,t) + W(x,t).