Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги / Обыкновенные дифференциальные уравнения

..pdf
Скачиваний:
1
Добавлен:
12.11.2023
Размер:
2.7 Mб
Скачать

Пусть соответствующие формулы будут:

x1 = ω

1 ψ(

 

 

ψ1,

 

 

 

2 ,ψ...,

 

 

n1 ),

 

 

 

 

 

 

 

 

x2 = ω

2 ψ(

 

 

ψ1,

 

 

 

 

2 ,ψ...,

 

 

 

n1 ),

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

…………………………

 

 

x

= ω

n1

ψ(

 

 

 

ψ

,

 

 

,ψ...,

 

n1

).

 

 

 

2

 

n1

 

 

 

 

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 , ...ψ ,

 

n1 )

равны

величинам

 

 

xi (i =1, 2, ..., n 1); поэтому из

формулы

(7.16) получаем:

при xn = xn0

f = ϕ ( x1, x2 , ..., xn1 ), т.е. условие (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

, ...,

 

 

n1

 

,

 

 

 

 

 

 

 

 

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)

………………………

 

 

 

 

 

 

ψ n1 (z, x1,..., xn )= Cn1.

 

173

Общее решение уравнения (7.19) имеет вид

V (ψ 0ψ, 1,…ψ , n1),

где Φ – произвольная дифференцируемая функция. Из предыдущего вывода следует, что уравнение

Φ (ψ 0ψ, 1,…ψ , n1) =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

,,

 

 

n1

,

 

 

 

 

= 0.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

xn

 

 

 

 

xn

 

 

 

m

 

 

 

 

 

 

 

 

 

 

 

xn

 

 

 

 

 

 

 

xn

 

 

 

 

 

 

 

 

Разрешая относительно последнего аргумента и затем относительно f, получим:

x

x

x

 

f = xnmΨ

1

,

2

,,

n1

,

 

 

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,, xn01 ; эта система имеет вид:

174

 

 

 

ψ 0 (z, x1,, xn1, xn )= z0 ,

 

 

 

 

 

ψ

i

(z, x ,, x

n1

, x )=

x0

(i= 1, 2, , n1).

(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,, n1).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

∂ ψ k

=

∂ ψ k

+

 

∂ ψ k

 

∂ φ

(k = 0,1, 2,, n 1; j =1, 2,, n).

 

 

 

 

 

 

 

 

xj

xj

z

xj

 

 

 

 

 

Замечая, что функции ψ

0ψ,

1,…ψ ,

n1 удовлетворяют уравнению (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

, …,

 

 

 

 

 

 

 

 

∂ ψ n1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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,, n1).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Таким

образом,

 

ψ 0ψ,

1,…ψ , n1

являются

системой

n

решений

уравнения

в частных

производных

 

 

с

n

независимыми переменными

 

x1, x2,, xn :

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

P1 (ϕ , x1

,, xn )

u

+ … + Pn (ϕ , x1 ,, xn ) u = 0.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

x

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Отсюда, в силу замечания, существует тождественная

по

 

x1, x2,, xn

зависимость между ψ 0ψ,

1,…ψ ,

n1 :

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Φ Ψ( Ψ0, 1,…Ψ ,

n=1)

0

 

 

 

 

 

 

 

 

 

 

 

 

(7.26)

175

Следовательно, существует такая функция Φ (ψ 0ψ, 1,…ψ , n1), что по подстановке в нее вместо 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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

соотношения

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Φ (z2 y, y+ 2 zx− − y )= 0.

 

 

 

 

 

 

Но данное уравнение имеет еще решение z = x + y.

Если выражение

V =z xy подставить

в левую

часть уравнения

для

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 ,, xn01 изменяются

в окрестности значений

 

 

x10 ,,

 

n01 , а xn

 

 

 

 

 

 

x

 

z0 ,

изменяется в окрестности xn0 .

Решим для уравнения (7.17) задачу Коши: найти решение этого уравнения, которое при xn = xn0 обращается в данную дифференцируемую

функцию [определенную в окрестности значений (x10, x20,, xn01)]:

ϕ (x1,x2,,xn1).

Подставим в интегралы (7.24) вместо xn начальное значение xn0 , результат обозначим через ψ i , имеем:

ψ

 

(z, x ,, x

,

 

0 ) = ψ

 

(i =0,1, 2,

, n1).

i

x

i

 

 

1

 

n1

 

n

 

Разрешим уравнения (7.27) относительно z, x1, x2,, xn1 :

 

 

 

 

 

z = ϕ 0 (xn , z0 , x10,, xn01),

 

 

xi = ϕ i (xn , z0 , x10 ,, xn01) (i =1, 2,, n1).

Подставим в эти формулы вместо

xn числовое значение

z0, x10 ,, xn01 через ψ

0ψ,

 

1,…ψ ,

 

n1. Получим выражение вида:

 

 

(7.27)

xn0 и заменим

z = ω 0ψ(

 

ψ0 ,

 

 

 

1,ψ…,

 

 

 

 

 

 

 

 

 

 

(7.28)

 

 

 

 

 

 

n1),

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

xi = ω iψ(

ψ0 ,

 

 

1,ψ…,

 

n1)=

(i 1, 2,…−, n 1).

 

Легко видеть, что формулы (7.28) есть результат решения уравнений

(7.27) относительно z,x1,,xn1.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Решение задачи Коши дается формулой

 

V (z, x1,, xn ) ≡ ω 0 (ψ 0ψ, 1,…ψ , n1)

(7.29)

−ϕ [ω 1(ψ 0ψ, 1,…ψ , n1),,ω

n1(ψ

0ψ, 1,…ψ , n1)] = 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

 

n1

∂φ

∂ω

∂ψ

∂ω

∂ψ

i

 

1

∂ω

∂ψ

i

 

n1

 

 

 

 

 

i

 

0

+

 

 

+…+

 

 

.

∂ω

 

 

 

z∂ψ

 

 

 

 

 

z

i=1

∂ψ

0

 

 

1

 

z ∂ψ

n1

 

 

 

 

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,, n1), в силу равенств (7.27); функции ω 0ω, 1,…ω , n1 от этих аргументов ввиду равенств (7.28) дадут

соответственно z,x1,,xn1, и мы получим: z (x1,x2,,xn1) при xn = xn0 . Таким образом, формула (7.22) (при надлежащем выборе первых

интегралов) дает все частные решения, определяемые начальными данными Коши.

Само построение показывает, что решение задачи Коши – единственное в классе решений, удовлетворяющих условиям непрерывной дифференцируемости и одновременного необращения в нуль коэффициентов уравнения в точках решения.

Пример 5. Для уравнения примера 4 § 7.2 найти решение, удовлетворяющее данным Коши:

z = 2x при y = 0.

178

Исходим из первых интегралов ψ 0 z2y= 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

Соседние файлы в папке книги