книги / Обыкновенные дифференциальные уравнения
..pdfПусть соответствующие формулы будут:
x1 = ω |
1 ψ( |
|
|
ψ1, |
|
|
|
2 ,ψ..., |
|
|
n−1 ), |
|
|||||||||||
|
|
|
|
|
|
|
|||||||||||||||||
x2 = ω |
2 ψ( |
|
|
ψ1, |
|
|
|
|
2 ,ψ..., |
|
|
|
n−1 ), |
|
|||||||||
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
||||||||||||||||
………………………… |
|
||||||||||||||||||||||
|
|||||||||||||||||||||||
x |
= ω |
n−1 |
ψ( |
|
|
|
ψ |
, |
|
|
,ψ..., |
|
n−1 |
). |
|||||||||
|
|
|
2 |
|
|||||||||||||||||||
n−1 |
|
|
|
|
1 |
|
|
|
|
|
|
|
|
При этом, когда ψ i принимают значения
ψ |
0= ψ |
i |
( x , x , ..., x |
) , |
|
|
i |
1 2 |
n |
|
соответствующие функции ω i принимают значения
(7.15)
xi0 (i =1, 2, ..., n −1) . При
наличии частных производных у функций ψ i функции ω i также являются
дифференцируемыми. Покажем, что искомое решение уравнения (7.12), удовлетворяющее начальному условию (7.13), есть
f = ϕ (ω |
1ψ( ψ1 , |
2 ,ψ..., |
|
ωn −1 )ψ, ψ2 ( |
1 ,ψ |
2 , ..., |
n −1 ), ..., |
||||
|
ω |
n −1 |
(ψ |
ψ |
, |
2 |
, ...ψ , |
n −1 |
)). |
(7.16) |
|
|
|
|
1 |
|
|
|
|
|
|||
Действительно, |
во-первых, |
|
|
выражение |
(7.16), |
являясь функцией от |
частных решений ψ i , само является решением уравнения (7.12); далее, если
положить |
xn = xn0 , то величины |
ψ i обращаются в ψ |
i , в силу формул (7.14). |
||||||||
Но по |
формулам (7.15) |
ω i (ψ |
1ψ, |
|
|
2 , ...ψ , |
|
n−1 ) |
равны |
величинам |
|
|
|
||||||||||
xi (i =1, 2, ..., n −1); поэтому из |
формулы |
(7.16) получаем: |
при xn = xn0 |
f = ϕ ( x1, x2 , ..., xn−1 ), т.е. условие (7.13) выполнено.
Из построения решения очевидно, что оно однозначно определено начальными данными (7.13).
Пример 1. Найти общее решение уравнения
|
|
x |
∂ f |
|
+ x |
∂ |
|
f |
+ ... + x |
|
∂ |
|
f |
|
= 0. |
|
||||||||||
|
|
|
|
|
2 ∂ x2 |
n∂ |
|
|
|
|
|
|||||||||||||||
|
|
1 |
∂ x1 |
|
|
|
|
|
xn |
|
|
|
||||||||||||||
Запишем |
соответствующую |
|
|
|
систему |
обыкновенных |
||||||||||||||||||||
дифференциальных уравнений: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
dx1 |
= |
|
dx2 |
= ... = |
dxn |
. |
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
|
|
|
|
|
|
|
x1 |
|
|
x2 |
|
|
|
|
|
|
xn |
|
|
|
|
|
|
|||
Независимая система ее первых интегралов есть |
|
|||||||||||||||||||||||||
|
x1 |
= C1 |
, |
x2 |
= C2 , ..., |
|
|
xn −1 |
= Cn −1 |
|
(xn ≠ 0 ). |
|
||||||||||||||
|
xn |
|
|
|
|
|
|
|||||||||||||||||||
|
|
|
xn |
|
|
|
|
|
|
xn |
|
|
|
|
|
|
|
|
|
|||||||
Общее решение данного уравнения есть |
|
|
|
|
|
|
||||||||||||||||||||
|
|
|
|
f = Ψ |
|
|
x |
|
|
x |
|
x |
|
|
|
|
|
|||||||||
|
|
|
|
|
|
1 |
, |
|
2 |
, ..., |
|
|
n−1 |
|
, |
|
||||||||||
|
|
|
|
|
|
|
xn |
|
||||||||||||||||||
|
|
|
|
|
|
|
|
|
xn |
|
|
xn |
|
|
|
171
т.е. наиболее общая однородная функция n переменных нулевого измерения (обращение теоремы Эйлера об однородных функциях для функций нулевого измерения).
Пример 2. Проинтегрировать уравнение
y dz − x dz = 0 dx dy
(искомая функция обозначена буквой z).
Соответствующая система обыкновенных уравнений сводится к одному уравнению:
|
|
|
|
dx |
= − |
dy |
|
|
|
|
|
|
|
|
x |
|
|||
|
|
|
|
y |
|
|
|||
и имеет |
место |
интеграл |
x2 + y2 = C. |
Общее решение данного |
уравнения |
||||
z = ϕ (x 2+ |
y 2 ), |
где ϕ |
есть произвольная функция, геометрически |
||||||
представляет любую поверхность вращения с осью вращения |
0 z. Задача |
Коши: z = f (x) при y = 0, где |
f |
есть данная функция; ее решение: функция |
||
ψ есть x2 + y2 , |
следовательно, функция ψ |
есть x2 , откуда x = ψ ; искомое |
||
решение: z = f ( |
ψ )= f ( x 2+ |
y 2 |
). |
|
7.2.Линейное неоднородное уравнение
вчастных производных первого порядка
1.Обозначим искомую функцию через z, независимые переменные
через x1, x2,..., xn .
Названное в заголовке уравнение имеет вид
|
|
|
P1 |
∂ z |
+ P2 |
|
∂ |
z |
+ ... + Pn ∂ |
|
z |
= R , |
(7.17) |
||
|
|
|
|
∂ |
x |
|
|
|
|||||||
|
|
|
|
∂ x |
2 |
∂ |
x |
n |
|
||||||
где 1 2 |
n |
– |
1 |
|
|
|
|
|
|
|
|
||||
функции (непрерывные и непрерывно дифференциру- |
|||||||||||||||
P, P , ..., P , R |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
емые) от |
x1, x2,..., xn , |
z. Линейные однородные уравнения (7.12) |
являются |
||||||||||||
частным случаем уравнений типа (7.17), когда правая часть R ≡ |
0 и когда |
||||||||||||||
коэффициенты 1 |
2 |
n при производных не зависят от искомой функции. |
|||||||||||||
|
P, P,..., P |
|
|
|
|
|
|
|
|
|
|
Уравнение вида (7.17) может быть приведено к однородному линейному уравнению следующим образом. Мы будем искать удовлетворяющую уравнению (7.17) функцию z от независимых переменных
x1, x2,..., xn в неявном виде: |
|
V ( z, x1, x2 , ..., xn ) = 0, |
(7.18) |
так что искомой функцией будет V.
172
Из формулы (7.18) получаются значения частных производных |
∂ z |
: |
|||||||
|
|
|
|
|
|
|
|
∂ xi |
|
|
∂ z |
= −∂ |
|
V |
/∂ |
V |
(i = 1, 2, ..., n ) . |
|
|
|
∂ xi |
|
|
|
|
|
|||
|
∂ |
xi ∂ |
z |
|
|
|
Внося эти выражения в данное уравнение (7.17), умножив обе части на
∂V
∂z (мы предполагаем, очевидно, этот множитель не равным тождественно
нулю и рассматриваем окрестность точки, где он не обращается в нуль) и перенося все члены в левую часть, получаем соотношение
|
∂ V |
∂ |
V |
∂ |
V |
∂ |
V |
|
(7.19) |
||||
|
|
|
|
|
|
|
|
|
|
|
|
||
P1 ∂ x |
+ P2 ∂ |
x |
2 |
+ ... + Pn∂ |
x |
n |
+ R∂ |
z |
= 0. |
||||
|
|||||||||||||
1 |
|
|
|
|
|
|
|
|
|
Из вывода соотношения (7.19) следует: в предположении, что (7.18) определяет z как решение уравнения (7.17), соотношение (7.19) должно
удовлетворяться тождественно по x1, x2,..., xn , при условии, что вместо z
подставлено его значение, определенное формулой (7.18). Если мы потребуем больше, а именно – чтобы искомая функция V удовлетворяла
соотношению (7.19) тождественно относительно x1, x2,..., xn и z (т.е. рассматривая в уравнении (7.19), наряду с x1, x2,..., xn , также и z как
независимое переменное), то (7.19) окажется линейным однородным дифференциальным уравнением первого порядка с искомой функцией V и
n +1 независимыми переменными x1, x2,..., xn , z. Очевидно, каждое решение этого уравнения (7.19), содержащее z, будучи приравнено к нулю, дает соотношение (7.18), которое определит функцию z от x1, x2,..., xn ,
удовлетворяющую данному уравнению (7.17). Заметим, что мы наложили дополнительное ограничение, потребовав, чтобы соотношение (7.19)
выполнялось тождественно относительно x1, x2,..., xn , z; поэтому мы не
можем утверждать заранее, что указанным приемом получим все решения уравнения (7.17). Мы еще вернемся к этому вопросу.
Напишем систему обыкновенных дифференциальных уравнений, соответствующую линейному уравнению в частных производных (7.19):
dx1 |
= |
dx2 |
= … = |
dxn |
= |
dz |
. |
(7.20) |
|
|
|
|
|||||
P1 |
|
P2 |
|
Pn R |
|
Эта система n уравнений имеет n независимых первых интегралов; пусть это будут:
ψ 0 (z, x1,..., xn ) = C0 , |
|
|
ψ 1 (z, x1,..., xn )= C1, |
|
|
|
(7.21) |
|
……………………… |
|
|
|
|
|
|
|
|
ψ n−1 (z, x1,..., xn )= Cn−1. |
|
173
Общее решение уравнения (7.19) имеет вид
V =Φ (ψ 0ψ, 1,…ψ , n−1),
где Φ – произвольная дифференцируемая функция. Из предыдущего вывода следует, что уравнение
Φ (ψ 0ψ, 1,…ψ , n−1) =0 |
(7.22) |
определяет (если удовлетворяются условия теоремы существования для неявных функций) z как функцию от x1,x2,...,xn, причем эта функция удовлетворяет данному уравнению (7.17).
Пример 3. Найти решение уравнения
|
|
x1 |
∂ f |
+ x2 ∂ f + … + xn∂ |
|
|
f |
= mf , |
|
|
|||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||
|
|
|
∂ |
x |
|
|
|
|
|
∂ x |
2 |
|
|
|
|
|
|
|
∂ |
x |
n |
|
|
|
|
|
|
||||||||||||
где m – постоянное. |
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
Система |
обыкновенных |
|
|
|
|
дифференциальных |
уравнений, |
||||||||||||||||||||||||||||||||
соответствующая этому уравнению в частных производных, есть |
|
||||||||||||||||||||||||||||||||||||||
|
|
|
|
dx1 |
= |
|
dx2 |
= … = |
dxn |
= |
df |
; |
|
|
|
||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||
|
|
|
|
|
x1 |
|
|
|
|
x2 |
|
|
|
|
|
|
|
xn |
|
|
|
mf |
|
|
|||||||||||||||
система первых интегралов |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
x |
|
|
|
|
x |
2 |
|
|
|
|
|
|
|
|
|
x |
n −1 |
|
|
|
|
|
|
|
|
|
|
|
|
f |
|
|
||||||
|
1 |
= C1 |
, |
|
|
|
|
|
|
|
|
= C2 ,… |
|
|
|
|
= Cn −1 , |
|
|
|
|
= Cn . |
|
||||||||||||||||
|
|
|
|
|
xn |
|
|
|
|
xn |
|
|
m |
|
|||||||||||||||||||||||||
|
xn |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
xn |
|
|
||||||||||
Решение, содержащее произвольную функцию Φ |
, будет: |
|
|||||||||||||||||||||||||||||||||||||
|
|
|
|
|
x |
|
|
x |
|
|
|
x |
|
f |
|
|
|
|
|
|
|
|
|||||||||||||||||
|
|
|
Φ |
|
|
|
|
|
|
1 |
, |
2 |
,…, |
|
|
n−1 |
, |
|
|
|
|
= 0. |
|
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
xn |
|
|
|
|
xn |
|
|
|
m |
|
|
|
|
|
|
||||||||||||||
|
|
|
|
|
xn |
|
|
|
|
|
|
|
xn |
|
|
|
|
|
|
|
|
Разрешая относительно последнего аргумента и затем относительно f, получим:
x |
x |
x |
|
|||
f = xnmΨ |
1 |
, |
2 |
,…, |
n−1 |
, |
|
|
xn |
||||
xn |
xn |
|
где Ψ – произвольная функция. Это – полное обращение теоремы Эйлера об однородных функциях (см. пример 4).
2. Формула (7.22) решения дифференциального уравнения (7.17) выведена при дополнительном требовании, чтобы уравнение (7.19) удовлет-
ворялось функцией V тождественно по x1, x2,…, xn, z . Посмотрим, насколько
общей является эта формула, т.е. какие частные решения в ней заключаются. Пусть
z =ϕ (x1, x2,…, xn ) |
(7.23) |
есть какое-либо решение уравнения (7.17). Возьмем в системе первых интегралов уравнений (7.20) в качестве произвольных постоянных начальные
значения z0, x10,…, xn0−1 ; эта система имеет вид:
174
|
|
|
ψ 0 (z, x1,…, xn−1, xn )= z0 , |
|
|
|||||||||||
|
|
|
ψ |
i |
(z, x ,…, x |
n−1 |
, x )= |
x0 |
(i= 1, 2, …, n− 1). |
(7.24) |
||||||
|
|
|
|
|
|
1 |
|
|
|
n |
i |
|
|
|||
Подставим в левые части равенств (7.24) вместо z его выражение (7.23); |
||||||||||||||||
для получившихся функций от x1,x2,…,xn введем обозначения: |
|
|||||||||||||||
|
|
ψ kϕ( (x1,x2,…,xn),x1,x2,…,xn) ≡Ψ k (x1,x2,…,xn) |
|
|||||||||||||
Мы имеем: |
|
|
|
|
|
|
|
|
(k =0,1 2,…, n−1). |
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
∂ ψ k |
= |
∂ ψ k |
+ |
|
∂ ψ k |
|
∂ φ |
(k = 0,1, 2,…, n −1; j =1, 2,…, n). |
|
||||||
|
|
|
|
|
|
|||||||||||
|
∂ xj |
∂ xj |
∂ z∂ |
xj |
|
|
|
|
|
|||||||
Замечая, что функции ψ |
0ψ, |
1,…ψ , |
n−1 удовлетворяют уравнению (7.19), |
мы получаем после подстановки в соответствующие тождества на место z его
выражения (7.23) систему тождеств: |
∂ψ |
|
|
|
∂ψ |
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
n |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
∑Pj (ϕ (x1,…, xn ), x1,…, xn ) |
|
|
|
k + R(φ, x1,…, xn ) |
|
|
|
k |
= 0 |
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||
|
|
|
j=1 |
|
|
|
|
|
|
|
∂ xj |
|
|
|
|
∂ z |
|
|
|
|
|
|||||||
|
|
|
|
|
(k = 0,1, …, n −1). |
|
|
|
|
|
|
|
|
|
|
|
|
(7.25) |
||||||||||
|
|
Поскольку, по предположению, (7.23) есть решение уравнения (7.17), |
||||||||||||||||||||||||||
то имеем тождество |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
∑Pj (ϕ , x1,…, xn ) ∂ϕ − R(ϕ , x1,…, xn ) = 0. |
|
|
|
|
|
|
|
|
|
||||||||||||||||
|
|
|
n |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
j=1 |
|
|
|
|
|
∂ xj |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
Умножаем последнее равенство последовательно |
на |
|
|
∂ ψ |
0 |
, |
∂ ψ |
1 |
, …, |
|||||||||||||||||
|
|
|
|
|
|
|
||||||||||||||||||||||
|
∂ ψ n−1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
∂ z |
∂ z |
||||
|
и прибавляем |
к |
соответствующим |
равенствам |
(7.25). |
В |
|
силу |
||||||||||||||||||||
|
|
|
||||||||||||||||||||||||||
|
∂ z |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
определения функций Ψ k |
и их производных по xj получаем: |
|
|
|
|
|
||||||||||||||||||||||
|
|
|
P1 (φ( x1 ,…, xn ), x1 ,… |
, xn ) |
∂ ψ |
|
+ … + Pn |
(φ, x1 |
,…, xn ) |
|
∂ ψ |
|
|
|
= 0 |
|
|
|
||||||||||
|
|
|
k |
|
|
k |
|
|
|
|||||||||||||||||||
|
|
|
∂ x |
|
∂ x |
|
|
|
||||||||||||||||||||
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
n |
|
|
|
|
|
||||||
|
|
|
|
|
(k =0,1 2,…, n−1). |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
Таким |
образом, |
|
ψ 0ψ, |
1,…ψ , n−1 |
являются |
системой |
n |
решений |
||||||||||||||||||
уравнения |
в частных |
производных |
|
|
с |
n |
независимыми переменными |
|||||||||||||||||||||
|
x1, x2,…, xn : |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
P1 (ϕ , x1 |
,…, xn ) |
∂ u |
+ … + Pn (ϕ , x1 ,…, xn ) ∂ u = 0. |
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
|
|
|
|
|
|
∂ x |
|
|
|
|
|
∂ x |
n |
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
Отсюда, в силу замечания, существует тождественная |
по |
|
x1, x2,…, xn |
|||||||||||||||||||||||||
зависимость между ψ 0ψ, |
1,…ψ , |
n−1 : |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
Φ Ψ( Ψ0, 1,…Ψ , |
n=−1) |
0 |
|
|
|
|
|
|
|
|
|
|
|
|
(7.26) |
175
Следовательно, существует такая функция Φ (ψ 0ψ, 1,…ψ , n−1), что по подстановке в нее вместо z выражения ϕ (x1,…,xn) мы получаем тождественно
нуль. Таким образом, любое решение z при указанных условиях относительно коэффициентов Pi и R удовлетворяет соотношению вида (7.22).
В этом смысле (7.22) определяет общее решение.
Заметим, что в противоположность однородному линейному уравнению в частных производных мы могли здесь доказать представимость при помощи формулы (7.22) частных решений лишь при выполнении коэффициентами уравнения добавочных условий непрерывности частных производных от коэффициентов и необращения в нуль одновременно всех коэффициентов Pi. Если эти условия не выполнены, уравнение (7.17) может иметь решения, не входящие в формулу (7.22), они соответствуют обращению в нуль левой части (7.19) не тождественно, а лишь в силу соотношения V = 0 . Такие решения называются специальными.
Пример 4.
∂ z |
(1 |
+ z − x − y ) + |
∂ |
z |
= 2. |
|
|
|
|
|
|||
|
∂ |
|
||||
∂ x |
|
y |
Уравнение (7.19) имеет здесь вид
|
|
|
|
|
|
(1 + z − x − y ) |
∂ V |
+ |
∂ |
|
V |
|
+ 2∂ |
V |
= 0. |
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
∂ x |
|
∂ y |
|
|
|
|
∂ z |
|
|
|
|
||||||||
|
|
|
Соответствующая |
|
система |
|
обыкновенных |
дифференциальных |
|||||||||||||||||||||||
уравнений есть |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
dx |
|
|
|
= |
dy |
= |
dz |
. |
|
|
|
|
||||||||||
|
|
|
|
|
|
|
1 + |
z − x − y |
1 |
2 |
|
|
|
|
|
||||||||||||||||
|
|
|
Первые интегралы: |
1) |
z −2y =C1; |
2) |
|
из интегрируемой комбинации |
|||||||||||||||||||||||
|
dy |
= |
dz − dy − dx |
|
имеем y + 2 |
z − x − y = C2 . |
|
Общее решение получается из |
|||||||||||||||||||||||
1 |
− z − x − y |
|
|||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
соотношения |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
Φ (z− 2 y, y+ 2 z− x− − y )= 0. |
|
|
|
||||||||||||||||||||||
|
|
|
Но данное уравнение имеет еще решение z = x + y. |
Если выражение |
|||||||||||||||||||||||||||
V =z −x− y подставить |
в левую |
часть уравнения |
для |
V, мы |
получим |
||||||||||||||||||||||||||
− |
z − x − y = − V ; |
это |
|
выражение обращается |
в |
нуль |
только |
в силу |
|||||||||||||||||||||||
равенства V = 0. |
В |
данном |
примере |
|
|
в |
точках |
специального |
решения |
||||||||||||||||||||||
производные от коэффициентов перестают быть ограниченными. |
|
||||||||||||||||||||||||||||||
|
|
|
3. Пусть в начальной точке |
|
x10 , |
|
20 ,…, |
|
n0 мы имеем |
|
|
||||||||||||||||||||
|
|
|
|
|
|
|
|||||||||||||||||||||||||
|
|
|
|
|
|
|
|||||||||||||||||||||||||
|
|
|
x |
x |
|
|
|||||||||||||||||||||||||
|
|
|
z0 , |
|
|
Pn (z0 , x10 ,…, xn0 ) ≠ 0.
176
Тогда в системе (7.20) можно взять xn за независимое переменное. Систему первых интегралов обыкновенных дифференциальных уравнений (7.20)
возьмем в виде, данном |
формулами (7.24), причем |
начальные значения |
|||||
z0 , x10 ,…, xn0−1 изменяются |
в окрестности значений |
|
|
x10 ,…, |
|
n0−1 , а xn |
|
|
|
||||||
|
|
|
|||||
|
x |
||||||
|
z0 , |
изменяется в окрестности xn0 .
Решим для уравнения (7.17) задачу Коши: найти решение этого уравнения, которое при xn = xn0 обращается в данную дифференцируемую
функцию [определенную в окрестности значений (x10, x20,…, xn0−1)]:
ϕ (x1,x2,…,xn−1).
Подставим в интегралы (7.24) вместо xn начальное значение xn0 , результат обозначим через ψ i , имеем:
ψ |
|
(z, x ,…, x |
, |
|
0 ) = ψ |
|
(i =0,1, 2, |
, n−1). |
||||
i |
x |
i |
||||||||||
|
|
1 |
|
n−1 |
|
n |
|
… |
||||
Разрешим уравнения (7.27) относительно z, x1, x2,…, xn−1 : |
||||||||||||
|
|
|
|
|
z = ϕ 0 (xn , z0 , x10,…, xn0−1), |
|
||||||
|
xi = ϕ i (xn , z0 , x10 ,…, xn0−1) (i =1, 2,…, n−1). |
|||||||||||
Подставим в эти формулы вместо |
xn числовое значение |
|||||||||||
z0, x10 ,…, xn0−1 через ψ |
0ψ, |
|
1,…ψ , |
|
n−1. Получим выражение вида: |
|||||||
|
|
(7.27)
xn0 и заменим
z = ω 0ψ( |
|
ψ0 , |
|
|
|
1,ψ…, |
|
|
|
|
|
|
|
|
|
|
(7.28) |
||||||||
|
|
|
|
|
|
n−1), |
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
xi = ω iψ( |
ψ0 , |
|
|
1,ψ…, |
|
n−1)= |
(i 1, 2,…−, n 1). |
|
|||||||||||||||||
Легко видеть, что формулы (7.28) есть результат решения уравнений |
|||||||||||||||||||||||||
(7.27) относительно z,x1,…,xn−1. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
Решение задачи Коши дается формулой |
|
||||||||||||||||||||||||
V (z, x1,…, xn ) ≡ ω 0 (ψ 0ψ, 1,…ψ , n−1) − |
(7.29) |
||||||||||||||||||||||||
−ϕ [ω 1(ψ 0ψ, 1,…ψ , n−1),…,ω |
n−1(ψ |
0ψ, 1,…ψ , n−1)] = 0. |
|||||||||||||||||||||||
|
|||||||||||||||||||||||||
Докажем прежде всего, что уравнение (7.29) определяет z как |
|||||||||||||||||||||||||
однозначную, непрерывную и дифференцируемую функцию от |
x1,x2,…,xn |
||||||||||||||||||||||||
в окрестности значений |
x10 , |
|
20 ,…, |
|
n0 . Для этого достаточно показать, что |
||||||||||||||||||||
x |
x |
||||||||||||||||||||||||
|
|
|
|
∂ V |
|
|
∂ V |
|
|
|
|||||||||||||||
|
|
|
|
|
|
≡ |
|
|
|
|
≠ 0. |
|
|||||||||||||
|
|
|
|
|
|
|
|
||||||||||||||||||
|
|
|
|
∂ z |
0 |
|
∂ z |
z= |
z |
0 |
0 |
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
x = |
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
x |
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
i i |
|
Вычисляем эту производную:
177
∂ V |
= |
∂ω ∂ ψ |
0 |
+ |
∂ω ∂ |
ψ |
1 |
+ … + |
∂ω |
0 |
∂ |
|
ψ |
n −1 |
− |
|
|
0 |
|
0 |
|
|
|
|
|
||||||||
∂ z |
∂ψ 0∂ |
|
∂ 1 |
|
|
∂ n −1 |
|
|
|
|||||||
|
z ∂ψ |
|
z |
∂ψ |
|
|
z |
|
n−1 |
∂φ |
∂ω |
∂ψ |
∂ω |
∂ψ |
i |
|
1 |
∂ω |
∂ψ |
i |
|
n−1 |
|
||||||
−∑ |
|
|
|
|
i |
|
0 |
+ |
|
|
+…+ |
|
|
. |
||||||
∂ω |
|
∂ |
|
|
z∂ψ |
|
|
|
|
|
z |
|||||||||
i=1 |
∂ψ |
0 |
|
|
∂ |
1 |
|
z ∂ψ |
∂ |
n−1 |
|
|
||||||||
|
|
i |
|
|
|
|
|
|
|
|
|
|
|
|
|
Согласно формулам (7.28) мы имеем:
|
∂ω |
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
0 ≡ |
∂ω |
0 |
|
|
|
= 1; |
|||
|
|
|
|
|
|
|
|
||||||
∂ψ |
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|||||||
|
|
|
xn |
=xn |
∂ψ |
|
|
0 |
z |
=z0 |
|||
|
0 z=z0 |
|
|
|
|
|
0 |
||||||
|
|
|
|
xi =xi0 |
|
|
|
|
|
xi |
=xi |
|
|
∂ω |
|
|
|
|
z |
|
|
аналогично |
|
0 |
|
≡ |
∂ |
|
= 0, |
||
|
|
|
|
|
|
||||
|
|
|
|||||||
|
|
∂ψ |
i 0 |
|
∂ |
xi 0 |
|
того, имеем ли i = j или i ≠ |
|
j . |
||
Далее получим: |
|
|
|
|
|
∂ψ |
|
|
= 1, |
|
∂ z |
0 |
|
|
|
|
|
0 |
Итак,
∂ω |
|
|
|
или |
|
в зависимости от |
|
а также |
|
i |
|
|
|
||
|
∂ψ |
|
|
=1 |
|
0, |
|
|
j |
|
|
|
|
∂ψ |
|
= 0 (i =1, 2,…, n −1). |
|
|
∂ z |
i |
|
|
|
0 |
∂ V |
|
= 1; |
||
|
|
|
||
∂ z |
||||
|
|
0 |
следовательно, эта производная не обращается в нуль в окрестности точки (z0 , x10 ,…, xn0 ) , и формула (7.29) определяет z как функцию от x1,x2,…,xn .
Поскольку формула (7.29) является специальным видом уравнения (7.22), то полученная из нее функция z есть решение уравнения (7.17). Наконец, легко видеть, что она решает поставленную задачу Коши. Действительно, при
xn = xn0 функции ψ i обратятся в ψ i (i =0,1 2,…, n−1), в силу равенств (7.27); функции ω 0ω, 1,…ω , n−1 от этих аргументов ввиду равенств (7.28) дадут
соответственно z,x1,…,xn−1, и мы получим: z =ϕ (x1,x2,…,xn−1) при xn = xn0 . Таким образом, формула (7.22) (при надлежащем выборе первых
интегралов) дает все частные решения, определяемые начальными данными Коши.
Само построение показывает, что решение задачи Коши – единственное в классе решений, удовлетворяющих условиям непрерывной дифференцируемости и одновременного необращения в нуль коэффициентов уравнения в точках решения.
Пример 5. Для уравнения примера 4 § 7.2 найти решение, удовлетворяющее данным Коши:
z = 2x при y = 0.
178
Исходим из первых интегралов ψ 0 ≡ z− 2y= C1, ψ 1≡ |
+y |
|
2 −z −x |
|||||||||
Подставляя значение y = 0, получаем: |
|
|
|
|
|
|
|
|||||
z = ψ |
0 |
, |
2 z − x = ψ |
1 |
, откуда z = ψ |
0 |
, x= ψ |
|
− |
ψ |
12 |
. |
|
|
|||||||||||
|
|
|
||||||||||
|
|
|
|
|
|
0 |
4 |
|
||||
|
|
|
|
|
|
|
|
|
|
|
Подставляя эти значения в данное начальное уравнение и заменяя ψ i , получим искомое решение:
ψ 0 − 2ψ 0 |
+ |
ψ 12 |
= 0, |
|
|||
|
2 |
|
или 2ψ 0 − ψ 12 = 0. Иначе:
2z − 4 y − y2 − 4 y z − x − y − 4z + 4x + 4 y = 0,
или
4 y z − x − y = 4x − 2z − y2 ,
откуда
y= C2 .
ψ i через
z = 2x + 3 y2 − 2y x − y + y2
2 2
(знак минус перед радикалом, как это следует из проверки, соответствует знаку плюс перед радикалом в данном уравнении).
4. Геометрический смысл линейного неоднородного уравнения с двумя независимыми переменными.
Рассмотрим уравнение
P |
∂ z |
+ Q |
∂ |
z |
= R. |
|
|
|
|
|
|||
|
∂ |
|
||||
|
∂ x |
y |
В этом случае система трех переменных x, y, z допускает простое
истолкование – как координаты точки трехмерного пространства. Эта интерпретация поможет нам глубже проникнуть в факты, связанные с линейным уравнением в частных производных. Введем обозначения Монжа для частных производных:
|
∂ z |
≡ p, |
∂ |
z |
≡ q |
|
|
|
|
|
|
|
|
||
|
|
∂ |
|
|
|||
|
∂ x |
y |
|
||||
и напишем уравнение в виде |
|
|
|
|
|
||
|
Pp + Qq = R, |
(7.30) |
|||||
где P, Q, R – заданные (дифференцируемые) функции от x, y, z. |
|
||||||
Искомое решение z = f (x, y) |
|
|
представляет собой |
уравнение |
поверхности (интегральная поверхность), p и q – угловые коэффициенты ее касательной плоскости в точке с координатами (x, y, z).
Уравнение этой плоскости есть
179
Z −z = p(X −x)+q(Y − y). |
(Т) |
||||||
Уравнение (7.30) относит каждой точке пространства вектор |
|
||||||
|
|
(P,Q,R). |
|
|
(V) |
||
Уравнение прямой, проходящей через точку (x, y, z) в направлении, |
|||||||
даваемом векторным полем, есть |
|
|
|
|
|
|
|
|
X − x |
= |
Y − y |
= |
Z − z |
|
|
|
|
|
R . |
|
|||
|
P |
|
Q |
(D) |
Уравнение в частных производных (7.30) показывает, что прямая (D) лежит в плоскости (T).
Кривые в пространстве, которые в каждой точке касаются соответствующего вектора (P,Q,R) и вообще называются векторными
кривыми, для уравнения (7.30) называются характеристическими кривыми,
или характеристиками. Их дифференциальные уравнения есть
dx |
= |
dy |
= |
dz |
. |
(7.31) |
|
|
|
||||
P Q R |
|
Мы покажем, что каждая интегральная поверхность составляется из характеристик. Действительно, поставим себе задачу: найти на данной интегральной поверхности z = f (x, y) такие кривые, которые в каждой точке касаются вектора (V); очевидно, для проекций этих кривых на плоскость
(x, y) мы получаем уравнение dy = Q , где в правой части надо заменить z dx P
через функцию f (x, y). Но мы можем добавить еще одно дифференциальное
уравнение; действительно, при перемещении по поверхности дифференциалы связаны соотношением
dz = pdx + qdy,
или
dz = p + q dy = p + q Q = Pp + Qq .
dx |
dx |
P |
P |
|||
В силу уравнения |
(7.30) это |
дает |
dz |
= |
R |
, т.е. мы пришли |
|
|
|||||
|
|
|
dx P |
кдифференциальным уравнениям характеристик (7.31). Заметим, что система
(29) может быть проинтегрирована без знания интегральной поверхности, и мы получим семейство характеристик от двух параметров, обладающее тем
свойством, что через каждую точку (x0, y0, z0) пространства (точнее, той
области, где выполнены условия существования и единственности решения) проходит одна характеристическая кривая. Если взять два первых интеграла системы(7.31):
180