л.1 ур-я в ч. производных 1-го порядка (Лекции по УрМатФизу)
Описание файла
Документ из архива "Лекции по УрМатФизу", который расположен в категории "". Всё это находится в предмете "уравнения математической физики (умф)" из 6 семестр, которые можно найти в файловом архиве МГТУ им. Н.Э.Баумана. Не смотря на прямую связь этого архива с МГТУ им. Н.Э.Баумана, его также можно найти и в других разделах. .
Онлайн просмотр документа "л.1 ур-я в ч. производных 1-го порядка"
Текст из документа "л.1 ур-я в ч. производных 1-го порядка"
Линейные и квазилинейные уравнения в частных производных первого порядка.
I. Линейным однородным уравнением в частных производных первого порядка в R3 называют равенство
P(x,y,z) + Q(x,y,z) + R(x,y,z) = 0 . (1)
Здесь u = u(x,y,z) – искомая функция, P(x,y,z) , Q(x,y,z) , R(x,y,z) - достаточно гладкие функции своих аргументов. Если рассматривать их как компоненты векторного поля в R3 , то равенство (1) геометрически означает ортогональность этого векторного поля градиенту u искомой функции. С каждым гладким векторным полем связано понятие векторных линий (или интегральных кривых) этого поля, определяемых как решения системы О.Д.У вида
которая, введением параметра независимого t , может быть переписана в обычном виде
Хорошо известно, что решения системы (2) могут быть выписаны в виде первых интегралов
1 (x, y, z) = C1 ; 2 (x, y, z) = C2 ; 3 (x, y, z) = C3 .
Взяв два из этих интегралов, 1 (x, y, z) = C1 ; 2 (x, y, z) = C2 , получим двухпараметрическое семейство линий, называемых характеристиками . Покажем, прежде всего, что левая часть любого интеграла (x, y, z) = C системы (2) удовлетворяет уравнению (1). Действительно, вдоль любой интегральной кривой d 0, так как =C.
Следовательно, вдоль такой кривой
d = dx + dy + dz 0 . При этом, dx = Pdt ; dy = Qdt ; dz = Rdt, если x = x(t); y = y(t) ; z = z(t) – интегральная кривая в параметрической форме. Подставляя эти выражения в тождество, получим
( P + Q + R)dt 0 , или, так как dt 0 ,
P(x,y,z) + Q(x,y,z) + R(x,y,z) 0 . (3)
Тождество (3) означает, что интеграл удовлетворяет уравнению (1). Если теперь взять два любых первых интеграла 1 (x, y, z) = C1 ; 2 (x, y, z) = C2 системы (2), то
(1 (x, y, z); 2 (x, y, z)), где - произвольная дифференцируемая функция, также будет ее интегралом, поскольку вдоль интегральных кривых (1 (x, y, z); 2 (x, y, z)) = ( C1; C2 ) = C . Получаем, что
u(x ,y, z) = (1 (x, y, z); 2 (x, y, z)) (4)
является решением уравнения (1). Покажем, что такое решение является единственным, то есть, если (x,y,z) - решение уравнения (1), то существует такая дифференцируемая функция (p, q) двух переменных, что (x,y,z) = (1 (x, y, z); 2 (x, y, z)).
Действительно, так как все три функции (x,y,z) ; 1 (x, y, z); 2 (x, y, z) являются решениями уравнения (1), то
Равенства (5) – однородная линейная система трех алгебраических уравнений относительно функций P(x,y,z) , Q(x,y,z) , R(x,y,z), которые, очевидно, не равны нулю тождественно. Следовательно, определитель этой системы должен быть равен нулю:
Но этот определитель является якобианом функций (x,y,z) ; 1 (x, y, z); 2 (x, y, z), причем, ввиду независимости первых интегралов1 (x, y, z); 2 (x, y, z) , по крайней мере, один из миноров второго порядка ; ; отличен от нуля. В курсе математического анализа доказывается теорема о том, что в этой ситуации функция (x,y,z) является функционально зависимой от функций 1 (x, y, z); 2 (x, y, z). Что выражается в виде равенства
=(1 (x, y, z); 2 (x, y, z). Это и доказывает наше утверждение.
II. Нелинейным (квазилинейным) неоднородным уравнением в частных производных первого порядка называют уравнение вида
P(x,y,z) + Q(x,y,z) = R(x,y,z) . (7)
Здесь z = z(x, y) – искомое решение, от которого зависят также коэффициенты уравнения (7). Это уравнение можно привести к линейному однородному, если искать решение в виде неявной функции
u(x,y,z(x, y)) = 0. (8)
Действительно, если z(x, y)) – решение уравнения (7), превращающее равенство (8) в тождество, то + 0; + 0; и подставляя значения производных = / ; = / , в уравнение (7), получим уравнение (1):
P(x,y,z) + Q(x,y,z) + R(x,y,z) = 0 .
Следовательно, процедура построения решения уравнения (7) сводится к предыдущей, но решение z = z(x, y) находится как неявная функция
(1 (x, y, z); 2 (x, y, z)= 0. (9)
III. Предыдущие решения не были единственными, в силу произвольности функции . Задачи на построение единственного решения ставится следующим образом.
Найти решение уравнения (1), или (7), проходящее через кривую, заданную уравнениями
1 ( x, y, z)= 0, 2 ( x, y, z)= 0.
Задача решается следующим образом. Если два первых интеграла 1 (x, y, z) = C1 ;
2 (x, y, z) = C2 системы (2) построены, составляется система четырех уравнений
1 (x, y, z) = C1 ;
2 (x, y, z) = C2 ; (10)
1 ( x, y, z)= 0 ;
2 ( x, y, z)= 0,
из которой последовательно исключаются переменные x, y, z. Остается одно равенство( C1; C2 ) = 0 , или C1 = (C2 ), связывающее константы C1; C2 , причем функция , или явно определена. После этого единственное решение задачи выписывается в виде (4), если решается задача с уравнением (1), или в виде (9), если решается задача с уравнением (7).
Замечание. Задача может иметь неединственное решение, если окажется, что два последних равенства системы (10) являются первыми интегралами системы (2), то есть, линия, ими определяемая является характеристикой. В этом случае через такую линию проходит бесконечное множество интегральных поверхностей.
Примеры.
1. Найти общий интеграл уравнения
Составляя характеристическую систему уравнений dx = dy = dz , находим два первых интеграла этой системы в виде x – y = C1 ; z – x = C2 . Общее решение выписывается при этом в виде ( x – y ; z – x) = 0, где - произвольная дифференцируемая функция. Можно выписать решение и в виде, разрешенном относительно z :
z = x + ( x – y ), где - также произвольная дифференцируемая функция.
2. Построить интегральную поверхность уравнения
x – y = 0, проходящую через кривую x = 0 ; z = y2 .
Интегрируя систему уравнений = = , получаем два первых интеграла:
z = C1 ; x2 + y2 = C2 . Исключая x , y , z из системы уравнений
z = C1 ; x2 + y2 = C2 ; x = 0 ; z = y2 , получаем связь между C1 и C2 в виде C1 = C2 , откуда следует z = x2 + y2 .