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

Дифференциальные уравнения 2

.pdf
Скачиваний:
79
Добавлен:
18.02.2016
Размер:
515.6 Кб
Скачать

Для построения решений однородного линейного дифференциального уравнения (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 функция

H1, . . . , Ψ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), определяемая как решение функционального уравнения

U1(x1, . . . , xn, u), . . . , Ψn(x1, . . . , xn, u))− −φ(X21(x1, . . . , xn, u), . . . , Ψn(x1, . . . , xn, u)), . . . , (34.16)

Xn1(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 = ξ вид

U1(ξ, x2, . . . , xn, u(ξ, x2, . . . , xn)), . . . ,

Ψn(ξ, x2, . . . , xn, u(ξ, x2, . . . , xn)))− −φ(X21(ξ, x2, . . . , xn, u(ξ, x2, . . . , xn)), . . . ,

174

Ψn(ξ, x2, . . . , xn, u(ξ, x2, . . . , xn))), . . . , Xn1(ξ, 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,y0,z0) P (t, τ, θ)dt+

(35.6)

+Q(t, τ, θ)+ R(t, τ, θ)

 

задает интегральную поверхность 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,y0,z0)

µ(t, τ, θ)P (t, τ, θ)dt+

(35.7)

+µ(t, τ, θ)Q(t, τ, θ)+ µ(t, τ, θ)R(t, τ, θ)

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