Дифференциальные уравнения 2
.pdfДля построения решений однородного линейного дифференциального уравнения (34.2) при произвольном n рассматривают линии, которые являются характеристиками и задаются дифференциальной системой в симметрической форме
|
|
dx1 |
|
|
|
dxn |
|
|
|
du |
||
|
|
|
|
|
= . . . = |
|
|
|
|
= |
|
. |
f |
(x |
, . . . , x |
n |
) |
f (x |
, . . . , x |
n |
) |
0 |
|||
1 |
1 |
|
|
|
n 1 |
|
|
|
|
|
Интегрирование квазилинейного дифференциального уравнения (34.1) может быть сведено к интегрированию вспомогательного однородного линейного уравнения с частными про-
изводными первого порядка вида |
|
|
|
|
|||
f1(x1, . . . , xn, u) |
∂v |
+ . . . + fn(x1, . . . , xn, u) |
∂v |
+ |
|
||
|
|
|
|||||
|
∂x1 |
|
∂xn |
(34.11) |
|||
+g(x1, . . . , xn, u) |
∂v |
= 0, |
|
|
|
||
∂u |
|
|
|
||||
|
|
|
|
|
|
|
в котором искомой функцией является функция v, зависящая от x1, . . . , xn, u. Если функция
v = Ψ(x1, . . . , xn, u)
есть решение дифференциального уравнения (34.11), то функция
u = u(x1, . . . , xn),
удовлетворяющая функциональному уравнению
Ψ(x1, . . . , xn, u) = 0, |
(34.12) |
является решением исходного дифференциального уравнения (34.1). В самом деле, дифференцируя тождество
Ψ(x1, . . . , xn, u(x1, . . . , xn)) = 0
171
по переменным x1, . . . , xn, находим |
|
|
|
|
|
|
|
|
|
||||||||||||||||
|
|
|
|
|
|
∂u(x1, . . . , xn) |
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
= |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
∂Ψ(x1,∂x. . .i |
|
|
∂x |
|
|
|
|
|
|
|
|
|
|||||||||
= |
|
|
|
, xn, u)i |
, |
i = 1, n. |
|
(34.13) |
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
− |
1 |
|
|
|
n |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
∂Ψ(x , . . . , x , u) |
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
∂u |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
u=u(x1,...,xn) |
|
|
|
|
|
|
|
|
|
|||||
Поскольку при этом |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
∂Ψ(x1, . . . , xn, u) |
|
|
||||||||||||
|
|
f1(x1, . . . , xn, u) |
|
|
|
|
|
+ . . . + |
|
||||||||||||||||
|
|
|
|
∂x1 |
|
|
|
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
∂Ψ(x1, . . . , xn, u) |
|
|
||||||||
|
|
|
|
+fn(x1, . . . , xn, u) |
|
|
|
|
+ |
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
∂xn |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
∂Ψ(x1, . . . , xn, u) |
|
|
|||||||||
|
|
|
|
+g(x1, . . . , xn, u) |
|
|
|
= 0, |
|
||||||||||||||||
|
|
|
|
|
|
|
|
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
∂u |
|
|
|
|
|
|
|
|
|
то, используя соотношения (34.13), получаем |
|
|
|||||||||||||||||||||||
f1(x1, . . . , xn, u) |
∂u |
+ . . . + fn(x1, . . . , xn, u) |
∂u |
− |
|||||||||||||||||||||
|
|
||||||||||||||||||||||||
∂x1 |
∂xn |
−g(x1, . . . , xn, u) = 0,
т.е. решение u = u(x1, . . . , xn) функционального уравнения (34.12) есть решение квазилинейного дифференциального уравнения (34.1).
Отметим, что поскольку интегрирование однородного дифференциального уравнения (34.11) сводится к отысканию первых интегралов соответствующей дифференциальной систе-
мы в симметрической форме |
|
|
||||
|
dx1 |
dxn |
|
|||
|
|
|
= . . . = |
|
= |
|
|
f1(x1, . . . , xn, u) |
fn(x1, . . . , xn, u) |
||||
|
|
|
du |
|
(34.14) |
|
|
= |
|
, |
|
||
|
|
|
||||
|
g(x1, . . . , xn, u) |
|
то при построении решений квазилинейного дифференциального уравнения (34.1) его можно не приводить к виду (34.11).
172
Пусть Ψ1(x1, . . . , xn, u), . . . , Ψn(x1, . . . , xn, u) есть базис первых интегралов дифференциальной системы (34.14). Тогда при произвольной дифференцируемой функции H функция
H(Ψ1, . . . , Ψn)
также является первым интегралом дифференциальной системы (34.14). Поэтому решение u = u(x1, . . . , xn) функционального уравнения
H(Ψ(x1, . . . , xn, u), . . . , Ψn(x1, . . . , xn, u)) = 0
является решением дифференциального уравнения (34.1). Это решение зависит от произвольной функции H и поэтому является общим решением квазилинейного дифференциального уравнения (34.1).
Аналогичным образом выполняется построение решений линейного неоднородного уравнения с частными производными.
Пример 34.4. Дифференциальному уравнению
x1 |
∂u |
+ x2 |
∂u |
= u |
|
|
|||
|
∂x1 |
∂x2 |
соответствует дифференциальная система в симметрической
форме |
|
|
|
|
|
|
|
|
||
|
dx1 |
|
|
dx2 |
|
du |
||||
|
|
= |
|
|
= |
|
. |
|
||
|
x1 |
|
x2 |
|
||||||
|
|
|
|
|
u |
|||||
Функции |
|
x2 |
|
|
|
|
u |
|||
Ψ1(x1, x2, u) = |
, Ψ2(x1, x2, u) = |
|||||||||
x1 |
x1 |
|||||||||
|
|
|
|
|
|
|
образуют базис первых интегралов этой дифференциальной системы. Поэтому общее решение исходного дифференциального уравнения можно получить из соотношения
|
x2 |
|
u |
|
H( |
|
, |
|
) = 0, |
x1 |
x1 |
|||
|
|
173 |
|
где H есть произвольная дифференцируемая функция. Рассмотрим теперь квазилинейное дифференциальное урав-
нение (34.1) с начальным условием
u(x1, . . . , xn)|x1=ξ = φ(x2, . . . , xn). |
(34.15) |
Решением указанной начальной задачи является функция u = u(x1, . . . , xn), определяемая как решение функционального уравнения
U(Ψ1(x1, . . . , xn, u), . . . , Ψn(x1, . . . , xn, u))− −φ(X2(Ψ1(x1, . . . , xn, u), . . . , Ψn(x1, . . . , xn, u)), . . . , (34.16)
Xn(Ψ1(x1, . . . , xn, u), . . . , Ψn(x1, . . . , xn, u))) = 0,
где функции Ψ1, . . . , Ψn, образуют базис первых интегралов дифференциальной системы (34.14), а функции u = U(C1,
. . . , Cn), xi = Xi(C1, . . . , Cn), i = 2, n, определяются из си-
стемы функциональных уравнений
Ψ1(ξ, x2, . . . , xn) = C1,
. . . . . . . . . . . . . . . . . . . . . . . . . . . (34.17) Ψn(ξ, x2, . . . , xn) = Cn.
Всамом деле, переменные x2, . . . , xn, u, входят в левую часть соотношения (34.16) через интегралы Ψ1, . . . , Ψn дифференциальной системы (34.14). Поэтому функция u = u(x1, . . . , xn), определяемая соотношением (34.16), является решением дифференциального уравнения (34.1). Подставляя это решение в соотношение (34.16), получаем тождество, имеющее при
x1 = ξ вид
U(Ψ1(ξ, x2, . . . , xn, u(ξ, x2, . . . , xn)), . . . ,
Ψn(ξ, x2, . . . , xn, u(ξ, x2, . . . , xn)))− −φ(X2(Ψ1(ξ, x2, . . . , xn, u(ξ, x2, . . . , xn)), . . . ,
174
Ψn(ξ, x2, . . . , xn, u(ξ, x2, . . . , xn))), . . . , Xn(Ψ1(ξ, x2, . . . , xn, u(ξ, x2, . . . , xn)), . . . , Ψn(ξ, x2, . . . , xn, u(ξ, x2, . . . , xn)))) = 0.
Отсюда на основании определения функций U, X2, . . . , Xn, из системы (34.17) получаем
U(C1, . . . , Cn) − φ(X2(C1, . . . , Cn), . . . , Xn(C1, . . . , Cn)) = 0,
то есть
u(ξ, x2, . . . , xn) = φ(ξ, x2, . . . , xn).
Поэтому функция u = u(x1, . . . , xn), являясь решением дифференциального уравнения (34.1), удовлетворяет также и начальному условию (34.15).
Таким образом, построение решения дифференциального уравнения (34.1), удовлетворяющего начальному условию (34.15), может быть проведено по следующей схеме:
1)находим базис первых интегралов дифференциальной системы в симметрической форме (34.14), соответствующей исходному квазилинейному дифференциальному уравнению (34.1);
2)составляем систему функциональных уравнений (34.17), которую разрешаем относительно переменных x2, . . . , xn, u;
3)составляем функциональное уравнение (34.16), решение u = u(x1, . . . , xn) которого является решением исходной начальной задачи;
4)выписываем (если возможно) уравнение (34.16) относительно функции u для получения аналитического выражения для решения.
175
§ 35. Уравнения Пфаффа.
Уравнение Пфаффа представляет собой обобщение дифференциального уравнения первого порядка в симметрической форме и в трехмерном случае имеет вид
P (x, y, z)dx + Q(x, y, z)dy + R(x, y, z)dz = 0, |
(35.1) |
где функции P, Q и R дважды непрерывно дифференцируемы по своим аргументам на области G R3.
Интегралом уравнения Пфаффа (35.1) будем называть такую зависимость переменных x, y, z, при которой дифференциалы dx, dy, dz обращают уравнение (35.1) в тождество на области G. Если указанная зависимость представима в виде u(x, y, z) = 0 (или в параметрической форме x = x(u, v), y = y(u, v), z = z(u, v)), то будем ее называть
двумерным интегралом или интегральной поверхностью уравнения Пфаффа (35.1). Если же интеграл уравнения Пфаффа (35.1) представим в виде системы соотношений
u(x, y, z) = 0, v(x, y, z) = 0
(или в параметрической форме x = x(t), y = y(t), z = z(t)), то будем ее называть одномерным интегралом или инте-
гральной кривой (линией) этого уравнения.
−→
Рассмотрим в пространстве Oxyz векторное поле F , определяемое коэффициентами уравнения Пфаффа (35.1):
−→ −→ −→ −→
F (x, y, z) = P (x, y, z) i + Q(x, y, z) j + R(x, y, z) k .
Если (двумерный или одномерный) интеграл уравнения Пфаффа (35.1) задан векторным уравнением
−→ −→ −→ −→ −→ r = r (x, y, z) = x i + y j + z k ,
то |
|
−→ |
d r = dx−→i + dy−→ |
||
−→ |
j |
+ dz k , |
|
176
где dx, dy, dz в силу определения интеграла уравнения Пфаффа (35.1) вычислены с учетом существующей зависимости между переменными x, y, z. Поэтому уравнение Пфаффа (35.1) равносильно уравнению в векторной форме
−→ |
• |
−→ |
(35.2) |
F (x, y, z) |
|
d r (x, y, z) = 0. |
Если Π есть интегральная поверхность уравнения (35.2),
−→
то векторное поле F в каждой точке этой поверхности орто-
гонально плоскости, касательной к Π в данной точке. Други-
ми словами, поверхность Π есть нормальная трансверсаль
−→
характеристик поля F .
Аналогично, если L есть интегральная линия уравнения
−→
(35.2), то в каждой точке этой линии векторное поле F орто-
гонально касательной линии к L в данной точке, то есть ли-
ния L есть нормальная трансверсаль характеристик поля
−→
F .
Верно и обратное: каждая нормальная трансверсаль харак-
−→
теристик поля F является интегральной поверхностью (линией) уравнения (35.2).
Для векторного поля
−→ −→ −→ −→
F = P i + Q j + R k
ротор (вихрь) |
|
−→i |
|
−→ |
|
|
|||
|
|
|
j |
|
rot−→ |
|
P |
|
|
|
Q |
|||
|
∂ |
|
∂ |
|
|
|
|
||
F = |
|
|
|
|
|
|
∂x |
|
∂y |
|
|
|
||
|
|
|
|
|
|
|
|
|
|
∂. ∂z
−→R
Теорема 35.1. Условие |
|
|
|
|
|
||
rot−→ |
• |
−→ |
|
(x, y, z) |
|
G |
(35.3) |
F (x, y, z) |
F (x, y, z) = 0, |
||||||
|
|
|
|
|
является необходимым для существования двумерных интегралов уравнения Пфаффа (35.1) на области G.
177
Доказательство. Пусть соотношение u(x, y, z) = 0 задает интегральную поверхность Π на области G. Поскольку вектор
(∂u∂x, ∂u∂y , ∂u∂z )
нормален к поверхности Π, то из соотношения (35.2) следует, что в любой точке (x, y, z) Π векторы
(∂u∂x, ∂u∂y , ∂u∂z )
и(P, Q, R) параллельны, т.е. существует такая скалярная функция µ(x, y, z), что
P = µ |
∂u |
, Q = µ |
∂u |
|
, R = µ |
∂u |
. |
||||||||||||||||||
|
∂y |
|
|
|
|||||||||||||||||||||
|
|
∂x |
|
|
|
|
|
|
|
|
|
|
|
|
∂z |
||||||||||
Поэтому скалярное произведение rot−→ |
• −→ |
|
может быть за- |
||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
F |
|
|
F |
|
|
||||
писано в виде |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
∂ |
|
|
|
|
∂ |
|
∂ |
||||||||||||
|
|
|
|
|
|
µ |
∂u |
µ |
∂u |
µ |
∂u |
|
|||||||||||||
|
|
|
|
|
∂x |
∂y |
∂z |
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
rot−→ |
|
|
−→ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
. |
||
|
|
µ |
|
|
|
µ |
|
|
|
µ |
|
|
|
|
|||||||||||
F |
|
|
F = |
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
• |
|
|
∂x |
|
|
|
∂y |
|
∂z |
|
|||||||||||||||
|
|
|
|
|
|
|
|
∂u |
|
|
|
|
∂u |
|
|
∂u |
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
∂x |
|
|
|
|
∂y |
|
|
|
∂z |
|
|
||||||
Раскрывая данный определитель, |
|
|
с учетом равенства сме- |
шанных производных функции u, приходим к соотношению (35.3). Теорема 35.1 доказана.
Условие (35.3) будем называть условием интегрируемости уравнения Пфаффа (35.1). В координатной форме оно
записывается в виде |
|
|
|
|
|
|
|
|
|
|
||||
|
∂R |
|
∂Q |
∂P |
|
∂R |
|
∂Q |
|
∂P |
||||
P ( |
|
− |
|
) + Q( |
|
− |
|
) |
+ R( |
|
|
− |
|
) = 0. (35.4) |
∂y |
∂z |
∂z |
∂x |
∂x |
∂y |
178
−→ −→
Если rotF = 0 , т.е.
∂R |
= |
∂Q |
, |
∂P |
= |
∂R |
, |
∂Q |
= |
∂P |
, |
(35.5) |
|
∂y |
∂z |
∂z |
∂x |
∂x |
∂y |
||||||||
|
|
|
|
|
|
|
то условие (35.4) заведомо выполнено. В этом случае будем говорить, что условие интегрируемости выполнено в усиленной форме. При этом, как известно из математического анализа, выражение
P dx + Qdy + Rdz
является полным дифференциалом некоторой функции и поэтому в односвязной области в силу (35.5) не зависящий от пути интегрирования криволинейный интеграл второго рода
(x,y,z) |
|
u(x, y, z) = (x0,y∫0,z0) P (t, τ, θ)dt+ |
(35.6) |
+Q(t, τ, θ)dτ + R(t, τ, θ)dθ |
|
задает интегральную поверхность u(x, y, z) = |
0 уравнения |
Пфаффа (35.1). |
|
В общем случае из выполнения условия интегрируемости (35.3) следует существование такой функции µ(x, y, z), называемой интегрирующим множителем уравнения Пфаффа (35.1), что выражение
µP dx + µQdy + µRdz
является полным дифференциалом некоторой функции. В этом случае в односвязной области не зависящий от пути интегрирования криволинейный интеграл второго рода
(x,y,z) |
|
|
u(x, y, z) = (x0,y∫0,z0) |
µ(t, τ, θ)P (t, τ, θ)dt+ |
(35.7) |
+µ(t, τ, θ)Q(t, τ, θ)dτ + µ(t, τ, θ)R(t, τ, θ)dθ
179
задает интегральную поверхность u(x, y, z) = 0 уравнения Пфаффа (35.1).
Пусть условие (35.3) выполнено и Π есть интегральная поверхность уравнения Пфаффа (35.1). Тогда любая линия L
на поверхности Π является интегральной линией уравнения
−→
Пфаффа (35.1), т.к. характеристики векторного поля F ортогональны поверхности Π , а поэтому и любой линии, расположенной на поверхности Π. Поэтому, если u(x, y, z) = 0 есть двумерный интеграл уравнения Пфаффа (35.1), то при любой функции v(x, y, z) система соотношений
u(x, y, z) = 0, v(x, y, z) = 0
(если только она определяет линии в пространстве R3) задает одномерный интеграл уравнения Пфаффа (35.1).
Рассмотрим теперь задачу нахождения одномерных интегралов уравнения Пфаффа (35.1) без предположения о выполнении условия (35.3). Пусть z = z(x, y) есть уравнение произвольной поверхности Π. Если интегральная линия L уравнения Пфаффа (35.1) расположена на Π, то одновременно выполняются соотношения
z = z(x, y), P dx + Qdy + Rdz = 0, dz = ∂x∂z dx + ∂y∂zdy.
Поэтому |
|
|
|
|
|
|
P (x, y, z(x, y)) + R(x, y, z(x, y)) |
∂z(x, y) |
|
dx+ |
|
||
∂x |
|
|
|
|||
( |
|
|
) |
|
(35.8) |
|
∂z(x, y) |
|
|
||||
+(Q(x, y, z(x, y)) + R(x, y, z(x, y)) |
|
|
)dy = |
0. |
||
|
∂y |
Уравнение (35.8) представляет собой дифференциальное уравнение первого порядка в симметрической форме. Если φ(x, y) = 0 есть частный интеграл дифференциального уравнения (35.8), то линия L, определяемая соотношением φ(x, y) =
180