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

А. Р. Лакерник. Высшая математика.Краткий курс

.pdf
Скачиваний:
5
Добавлен:
10.10.2023
Размер:
3.65 Mб
Скачать

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

шение получается из (14.5) при некотором c0 C . Отсюда ϕ(x0,c0) = y0 , т.е. функция (14.5) удовлетворяет условию 2.

Пусть, наоборот, функция (14.5) удовлетворяет условию 2. Пусть y (x), x U (x0) – произвольное решение (14.3), график которого лежит в области D, и пусть в точке x0 это решение принимает некоторое значе-

íèå y0: y(x0) = y0 ( (x0,y0) D ). Согласно условию 2, существует зна- чение постоянной c0 C , при котором функция (14.5) удовлетворяет

условию ϕ(x0,c0) = y0 . Но так как по теореме 14.1 данному начальному условию может удовлетворять только одно решение уравнен ия (14.3), то функция ϕ(x,c0) , x U (x0), и является решением (14.3), т.е. функция (14.5) удовлетворяет условию 2. x

Определение 14.5. Равенство вида Ф (x,y,c) = 0, неявно задающее общее решение, называется общим интегралом дифференциального уравнения первого порядка (14.3).

14.3. Некоторые типы дифференциальных уравнений первого порядка и методы их решений

В этом разделе будут рассмотрены некоторые уравнения вид а y¢ = f (x,y) и указаны методы решения таких уравнений (функция f (x,y) будет предполагаться удовлетворяющей условиям теоремы 14.1).

1. Уравнения с разделяющимися переменными

Такими уравнениями называются уравнения вида

y¢ = f (x)g(y) èëè

dy

= f (x)g(y).

(14.6)

dx

 

 

 

Ð å ø å í è å

 

Предполагая, что g(y) ¹ 0 , запишем последнее равенство в виде

dy

 

= f (x)dx

(таким образом, мы сумели «разделить переменные» в

g(y)

 

 

уравнении (14.6)). Считая, что y = y(x) есть решение исходного уравнения, мы видим, что последнее равенство есть равенство ди фференциалов двух функций от x, которое может выполняться тогда, и только тогда, когда сами эти функции, или интегралы от их дифферен циалов, отличаются на произвольную постоянную:

14. Дифференциальные уравнения первого порядка...

ò

dy(x)

= ò f (x)dx + c , èëè

ò

dy

= ò f (x)dx + c .

(14.7)

 

 

g(y(x))

g(y)

Равенство (14.7), имеющее вид Ф (x,y,c) = 0, на самом деле и является общим интегралом исходного дифференциального урав нения (14.6).

В качестве примеров рассмотрим уравнения с разделяющими ся переменными (см. разд. 14.1).

Пример 1. Решить дифференциальное уравнение y¢ = - yx . Ð å ø å í è å

dydx = - yx ; dyy = - dxx ; ò dyy = -ò dxx + c; ln | y |= - ln | x | +c.

Из полученного равенства вида (14.7) в данном примере можно в ыразить y. Заменяя в формуле c на ln | c | (то и другое – произвольные постоянные), имеем:

ln | y |= ln | c | -ln | x |; ln | y |= ln

| c |

;

| y |=

c

;

y = ±

c

, èëè ( ±c можно

| x |

x

x

 

 

 

 

 

 

заменить на c) y = cx . В этой формуле c – произвольное число (мы решали

уравнение при условии y ¹ 0 , т.е. c ¹ 0 , но получаемая при c = 0 функция y = 0 тоже, очевидно, является решением исходного уравнения).

Пример 2. Решить дифференциальное уравнение y¢ = ky . Р е ш е н и е

dy

= ky;

dy

= kdx; ò

dy

= k òdx + c; ln | y |= kx + c.

dx

y

y

Из последнего равенства вида (14.7) выражаем y: | y |= ekx+c = ekx × ec ; y = ±ecekx .

Òàê êàê ec – произвольное положительное число, то, заменяя в последнем равенстве ±ec на c, имеем y = cekx. В этом равенстве c – произвольное число (получаемая при c = 0 функция y = 0 тоже является решением исходного уравнения).

Таким образом, в обоих примерах получили те самые решения, которые уже обсуждались в разд. 14.1. Вернувшись к определени ю 14.4 разд. 14.2 (условия 1 и 2), мы видим, что в разд. 14.1 уже было ус-

262

263

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

тановлено, что функции y = c è y = cekx являются общими решения-

x

 

ми рассмотренных дифференциальных уравнений.

 

Уравнениями с разделяющимися переменными являются такж е

уравнения вида

 

P(x)Q(y)dx + R(x)T(y)dy = 0,

(14.8)

в которых переменные разделяются после деления на R(x) и Q(y) (предполагается, что такое деление возможно, т.е. знаменате ли отличны от 0):

T (y) dy = - P(x) dx. Q(y) R(x)

Уравнения, сводящиеся к уравнениям с разделяющимися переменными

К уравнениям с разделяющимися переменными приводятся ур авнения вида

y¢ = f (ax + by + c),

(14.9)

где a ¹ 0, b ¹ 0 , c – некоторые постоянные (при a или b, равном 0, правая часть уравнения зависит только от одной переменной, т .е. (14.9) уже является уравнением с разделяющимися переменными).

 

 

 

Решение

 

 

 

 

 

 

 

 

 

Сделаем в уравнении (14.9) замену u = ax + by + c,

ãäå u = u(x) –

новая неизвестная функция. Тогда

y = 1(u - ax - c),

y¢ = 1(u¢

- a),

è

(14.9) принимает вид 1

 

 

 

b

 

 

 

 

 

b

 

 

(u′ − a) = f (u) ; u¢ = bf (u) + a, что является урав-

 

b

 

 

 

 

 

du

 

 

 

 

 

 

 

нением с разделяющимися переменными:

 

 

= dx .

 

 

 

bf (u) + a

 

 

 

Пример. Решить дифференциальное уравнение y¢ = 3x - 2y + 5 .

 

 

 

 

 

Ð å ø å í è å

 

 

 

 

 

 

 

 

 

u = 3x - 2y + 5; y = 1(3x - u + 5); y¢ =

1(3- u¢);

1(3 - u¢) = u; 3 - 2u = u¢;

 

 

2

 

 

 

 

2

 

2

 

 

 

 

 

 

 

du

du

 

du

 

 

1

d(3 - 2u)

 

 

 

dx = 3 - 2u;

 

= dx; ò

 

=

òdx + c; -

2 ò

 

 

 

= x + c;

 

 

3 - 2u

 

3 - 2u

3 - 2u

 

 

 

 

264

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

14. Дифференциальные уравнения первого порядка...

-

1

- 2u |= x + c; ln | 3 - 2u |= -2x -2c; | 3 -2u|= e2x2c; 3 -2u = ±e2xe2c;

2 ln | 3

 

 

 

 

u = ± 1 e2ce2x +

3 .

 

 

 

 

 

2

2

 

Заменяя здесь ±

1 e2c

на c, имеем: u = ce2x

+ 3

, откуда

 

 

 

2

 

2

 

y= 12 æç3x - ce2x - 23 + 5ö÷ = - 2c e2x + 23 x + 47 ,

èø

или, заменяя - 2c íà c,

y = ce2x + 32 x + 47.

2. Однородные уравнения первого порядка

Такими уравнениями называются уравнения вида

 

 

æ y ö

 

 

 

 

y¢ = f ç

 

÷

 

 

(14.10)

 

 

 

 

 

è x ø

 

 

y

 

 

 

 

 

 

 

(в этих уравнениях правая часть зависит только от отношения

 

).

x

Решение

 

 

 

Сделаем в этом уравнении замену

y

= u , где u = u(x) – новая не-

x

известная функция: y = ux; y

 

 

 

 

 

 

 

 

= u x + u. Тогда уравнение примет вид

u¢x + u = f (u) , или u¢x = f (u) - u . Но последнее уравнение является

уравнением

с разделяющимися

переменными:

du x

= f (u) - u ;

 

du

= dx

 

 

 

 

 

 

 

 

 

 

dx

 

 

и решается так же, как все такие уравнения:

 

 

f (u) - u

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

du

dx

 

 

 

 

 

 

 

 

ò

 

 

 

= ò x

 

+ c.

 

 

 

 

 

 

f (u) - u

 

 

 

 

Пример. Решить дифференциальное уравнение y¢ =

x2 + y2

.

 

xy

 

 

 

Ð å ø å í è å

 

 

 

 

 

 

 

 

 

 

 

 

 

Переписав это уравнение в виде y¢ =

x

+

y

, видим, что оно является

 

y

 

 

 

 

 

 

y

 

 

 

 

x

 

 

однородным. После замены

= u , y = ux ;

y¢ = u¢x + u уравнение принима-

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

265

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

åò âèä u¢x + u = u1 + u ; u¢x = u1 . Разделяем переменные в последнем уравне-

 

du

1

 

 

dx

. Далее имеем: òudu =

ò

dx

 

u2

 

 

íèè:

dx x = u

; udu =

x

x

+ c;

2

= ln | x |

+ ln | c |;

u2 = 2ln | cx |;

u2 = lnc2x2; u2 = lncx2; u = ±

 

lncx2 ;

y = ux = ±x

lncx2 .

Нетрудно проверить, что последняя функция дает общее реше ние исходного дифференциального уравнения: непосредственной подстановкой ее в уравнение проверяем, что при любом c > 0 она является решением этого уравнения. Задавая произвольное начальное условие y (x0) = y0 для чисел x0 è y0 одного или разных знаков, определяем знак «+» или «–» и зна-

 

y2

 

y2

/ x2

 

1

y2

/ x2

 

 

 

 

 

чение c: lncx2 =

0

; cx2

= e 0

0

; c =

 

e 0

0 .

 

 

 

 

x2

x2

 

 

 

 

0

0

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

0

 

 

 

 

 

 

 

 

Уравнения, сводящиеся к однородным

 

 

 

 

 

К однородным уравнениям приводятся уравнения вида

 

 

 

 

æ a x + b y + c

ö

a

 

b

 

 

 

y¢ = f ç

 

1

1

 

1

 

÷,

1

¹

1

.

(14.11)

 

 

 

 

 

 

 

 

a

 

 

 

 

ç a x + b y + c

÷

 

b

 

 

 

è

 

2

2

 

2

 

ø

2

 

2

 

 

Заметим, что если

a1

=

b1

 

a

b

 

 

 

 

 

ìåò âèä

 

 

2

 

 

 

2

 

 

 

 

 

 

 

 

 

æ ka x + kb y + c

ö

 

 

y¢ = f ç

2

2

 

1

÷

=

a x + b y + c

 

ç

 

÷

 

 

è

2

2

2

 

ø

 

 

= k , òî a1 = ka2 , b1 = kb2 , è (14.11) ïðè-

æ k(a x + b y) + c

ö

 

f ç

2

2

1

÷

= g(a2x + b2y),

a x + b y + c

ç

÷

 

è

2

2

2

ø

 

где g – некоторая функция, т.е. примет вид (14.9) и тем самым будет приводиться к уравнению с разделяющимися переменными.

Если бы в уравнении (14.11) c1 = c2 = 0, то это уравнение имело бы

âèä

æ a x + b y ö

 

y¢ = f ç

1

1

÷

=

 

 

ç a x + b y ÷

 

è 2

2

ø

 

æ a + b

y

ö

 

x

÷

 

ç 1 1

,

f ç

 

 

 

 

÷

 

y

 

ç a2 + b2

÷

 

 

 

 

è

 

x ø

 

т.е. вид (14.10), и являлось бы однородным. Поэтому будем пытаться путем некоторой замены (аргумента x и искомой функции y) обратить эти коэффициенты в 0. Положим, что

266

14. Дифференциальные уравнения первого порядка...

ì ~ a x = x+ ;

í ~

î y = y+b,

ãäå α и b – некоторые числа. Тогда

 

 

~

~

 

 

dy

 

 

~

 

~

 

 

 

 

 

 

 

 

 

 

d y

 

 

 

 

 

 

dx = d x, dy = d y è y¢ =

 

 

=

 

 

= y',

 

 

 

 

 

 

 

 

~

 

 

 

 

 

 

 

 

 

 

dx

d x

 

 

 

 

 

~ ~ ~

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где y = y(x). Уравнение (14.11) теперь принимает вид

 

 

 

 

 

~

~

ö

æ

 

~

~

 

 

ö

æ a1(x+a)+ b1(y+b)+ c1

a1x + b1y + a1a+ b1b+ c1

y¢ = f ç

 

 

 

÷

= f ç

 

 

 

 

 

 

 

÷.

 

~

~

 

~

~

 

 

 

 

è a2

(x+a)+ b2(y+b)+ c2

ø

è a2x+ b2 y + a2a+ b2b+ c2 ø

Теперь подберем α и b так, чтобы

ìa a + b b + c

= 0;

Эта система

íï

1

1

1

= 0.

 

 

 

 

 

 

ïa a + b b + c

 

 

 

 

 

 

 

 

 

î

2

2

2

 

 

 

 

имеет единственное решение, так как ее определитель a1

b1

¹ 0 , èáî

 

 

 

 

 

 

 

 

 

 

 

 

a

b

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

2

 

 

по условию строки этого определителя не пропорциональны. При таких α и b наше уравнение, как было показано выше, становится однородным.

Пример. Решить дифференциальное уравнение y¢ = x + y - 3 . x - y -1

Ð å ø å í è å

 

~

~

~

~

 

~ ~

 

 

ìx = x+a

x+a+ y+b–3

 

x+ y+a+b–3

í

 

Þy¢ =

 

 

 

=

 

 

;

~

~

~

 

~ ~

 

 

 

x+a – y –b–1

 

x – y

+a –b–1

î y = y+b

 

 

 

 

 

 

 

 

α и b должны удовлетворять системе

ìa + b - 3 = 0;

í

 

 

 

 

 

 

 

 

 

îa - b -1 = 0.

 

Cкладывая и вычитая уравнения этой системы, имеем

 

 

 

 

2a - 4 = 0; a = 2; 2b - 2 = 0; b = 1

, ò.å. x = x+2, y = y+1.

 

 

 

 

 

~

~

 

 

 

 

При такой замене

 

 

~

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

~ ~

 

1+

 

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x+y

 

~

 

 

 

 

 

 

 

~

 

~

 

 

x

 

 

 

 

 

y¢ =

~ ~

; y¢ =

~

.

 

 

 

 

 

 

x– y

 

1–

 

y

 

 

 

 

 

 

 

 

~

 

 

 

 

 

 

 

 

 

 

 

 

x

 

~

 

 

 

В последнем однородном уравнении сделаем замену

y

~ ~

 

~

= u, y = ux,

 

 

 

 

 

 

 

 

 

 

x

 

~

~

 

 

 

 

 

 

 

 

 

 

 

y' = u¢ x+ u. Тогда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

267

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

 

 

~

 

1+ u

~

 

 

1+ u

 

 

 

 

1+ u - u + u2

1+ u2

 

 

 

 

 

 

 

 

u¢x+ u =

1– u

; u¢x =

 

1- u

- u =

 

 

 

 

 

1- u

 

 

=

1- u .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

du

 

 

 

 

1+ u

2

 

 

 

1– u

 

 

 

 

 

 

~

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

~

 

 

 

 

 

 

du =

d x

 

 

 

 

 

 

 

 

 

 

 

 

 

Разделяем переменные:

 

 

~

x =

1– u

,

 

 

 

 

 

2

 

 

~

 

 

и интегрируем:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d x

1+ u

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1– u

 

 

 

~

 

 

 

 

 

 

 

 

d u

1 d(1+ u

2

)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d x

 

 

 

 

 

 

 

 

 

 

 

 

 

~

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

du =

 

~

 

+ c;

 

 

 

 

 

 

2

 

ò

 

 

1+ u2

 

= ln

 

x

 

+ ln

 

c

 

;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ò1+ u2

 

 

 

ò x

 

 

 

 

 

ò1+ u2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

2

 

 

 

 

 

 

~

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

~

~

 

 

2

 

 

 

 

arctg u – 2 ln(1+ u

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

) = ln cx ; arctg u = ln

 

1+ u

 

 

+ ln cx = ln cx

1+ u

 

.

 

 

Теперь вернемся к переменным x и y; подставляя в эту формулу

 

 

 

 

 

 

 

 

 

 

~

 

 

 

 

 

 

 

~

 

 

 

 

 

 

 

~

 

 

 

y –1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

 

 

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x = x –

2, y = y –1, u =

 

 

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

~

x –2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

имеем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

arctg

y - 1

= ln | c(x - 2) |

1 +

æ y

- 1 ö2

èëè arctg

 

 

y -1

 

= lnc (x - 2)

2

+(y

-1)

2

.

 

ç

 

 

 

 

 

÷

 

 

 

 

 

 

 

 

 

 

 

x - 2

 

 

 

 

 

 

 

 

 

è x

- 2 ø

 

 

 

 

 

 

 

 

 

 

 

x - 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Последнее равенство есть общий интеграл исходного диффе ренциального уравнения.

3. Линейные уравнения первого порядка

Уравнение y¢ = f (x,y) называется линейным, если его правая часть f (x,y) является линейной функцией y. Обычно линейные уравнения записывают в виде

y¢ + P(x)y =Q(x).

(14.12)

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

Первый способ (метод Бернулли). Будем искать решение уравнения (14.12) в виде y = uv , где u = u(x), v = v(x) – некоторые функции. Тогда y¢ = u¢v + uv¢ и (14.12) принимает вид u¢v + uv¢ + P(x)uv =Q(x). Перепишем последнюю формулу следующим образом:

u¢v + u(v¢ + P(x)v) =Q(x).

(14.13)

Теперь выберем функцию v (x) такой, чтобы

v¢ + P(x)v = 0,

(14.14)

268

14.Дифференциальные уравнения первого порядка...

àзатем найдем все функции u (x), при которых справедливо равенство (14.13) (т.е. при нахождении решения в виде произведения двух со множителей выбираем один из этих сомножителей так, как нам уд обно, а затем находим второй сомножитель так, чтобы их произведен ие было решением).

Разделяя переменные, имеем

dv

= -P(x)v;

dv

= -P(x)dx;

ò

dv

= -òP(x)dx + c; ln | v |= -òP(x)dx + c;

dx

v

v

 

 

| v |= e

òP(x)dx+c

= e

òP(x)dx

e

c

;

 

c òP(x)dx

èëè v = ce

òP(x)dx

.

 

 

 

 

 

 

v = ±e e

 

Так как нам достаточно найти лишь одно решение уравнения (14.14), то возьмем в последней формуле c = 1. Тогда

 

v = eòP(x)dx.

 

 

 

 

(14.15)

Далее из (14.13) и (14.15) имеем: u¢v =Q(x); u¢

Q(x)

 

=

 

( v(x) ¹ 0 )

v(x)

Q(x)

 

 

P(x)dx

 

 

 

 

u = ò v(x)dx + c = òQ(x)eò

 

 

 

dx + c;

 

é

òP(x)dx

dx + c

ù

òP(x)dx

.

(14.16)

y = êòQ(x)e

 

úe

 

 

 

ë

 

 

 

û

 

 

 

 

 

Ниже будет проверено, что функция (14.16) и есть общее решение исходного линейного уравнения (14.12).

Пример. Решить дифференциальное уравнение y¢ - x3 y + x = 0 . Ð å ø å í è å

Обычно в примерах используют не готовую формулу (14.16), а пров о- дят для каждого конкретного уравнения те действия, которы е к ней привели.

Будем искать решение уравнения в виде y = uv . Тогда y¢ = u¢v + uv¢ и

уравнение принимает вид

u¢v + uv¢ -

3

 

 

 

 

 

æ

 

3

 

ö

+ x = 0 . Ïî-

x

uv + x = 0 ; u¢v + uçv¢ -

x

v ÷

 

 

 

 

 

 

 

 

è

 

 

ø

 

 

 

 

 

3

 

 

dv

 

3

dv

 

3dx

 

 

ò

dv

= 3ò

dx

 

требуем, чтобы

v¢ - x v = 0 . Тогда

dx

=

x v ;

v

=

x

;

 

v

 

x

+ c ;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

269

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

ln | v |= 3ln | x | + ln | c | ; ln | v |= ln | cx3 | ; | v |=| cx3 | ; v = ±cx3 , èëè v = cx3 . Принимая здесь c = 1 , получаем, что v = x3.

Теперь u¢x3 + x = 0 , u¢ = - x12 , откуда u = -ò dxx2 + c = 1x + c è

y= æç 1x + cö÷ x3 = cx3 + x2.

èø

Второй способ (метод вариации произвольной постоянной). Рассмотрим так называемое линейное однородное уравнение, соотве тствующее данному уравнению (14.12):

 

 

 

y¢ + P(x)y = 0.

(14.17)

 

Решаем это уравнение (с разделяющимися переменными):

dy

 

dy

dy

 

dx

= -P(x)y;

y = -P(x)dx;

ò y = -òP(x)dx + c; ln | y |= -òP(x)dx + c;

 

| y |= eòP(x)dx+c = eòP(x)dxec; y = ±eceòP(x)dx

èëè

 

 

 

y = ceòP(x)dx.

(14.18)

 

Теперь будем искать решение уравнения (14.12) по той же форму-

ле (14.18), считая, что в ней c = c(x) (отсюда и название метода). Тогда

y¢ = c¢eòP(x)dx + c(eòP(x)dx )¢ = c¢eòP(x)dx + ceòP(x)dx (-òP(x)dx )¢ = = c¢eòP(x)dx - ceòP(x)dxP(x).

Подставляя эту производную в формулу (14.12), имеем

c¢eòP(x)dx - cP(x)eòP(x)dx + P(x)ce òP(x)dx =Q(x ), èëè c¢eòP(x)dx =Q(x)

(т.е. члены, содержащие с, всегда сокращаются, остается только член, содержащий c¢).

 

Отсюда c¢ =Q(x)e

òP(x)dx

; c = òQ(x)e

òP(x)dx

~ ~

заменяем наc)

 

 

 

 

 

dx + c è ( c

y =

é

òQ(x)e

òP(x)dx

dx + c

ù

òP(x)dx

, т.е. мы опять получили формулу

ê

 

úe

 

 

 

ë

 

 

 

 

û

 

 

 

 

 

 

(14.16).

14. Дифференциальные уравнения первого порядка...

Пример. Решить дифференциальное уравнение

(1+ x2 )y¢ - 2xy = (1+ x2 )2 .

Ð å ø å í è å

Переписав это уравнение в виде y¢ - 1+2xx2 y = 1+ x2 , видим, что оно действительно является линейным. Соответствующее однородное уравнение

имеет вид

y¢ -

 

2x

 

 

 

y = 0

. Решаем это уравнение:

 

 

 

1+ x2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dy

 

 

 

2x

 

dy

2x

dy

 

2x

 

 

 

 

 

 

 

dx =

 

 

y;

 

y =

 

dx; ò

 

= ò

 

dx + c;

 

 

 

 

 

 

 

 

1+ x2

 

1+ x2

y

1+ x2

 

 

 

 

 

 

ln | y |= ò

d(1+ x2 )

+ ln | c |= ln |1+ x2 | + ln | c |= ln | c(1+ x2 ) |

 

 

 

 

 

 

 

1+ x2

 

 

 

 

 

| y |=| c(1+ x2 ) |; y = ±c(1+ x2 )

èëè y = c(1+ x2 ).

 

 

 

Теперь ищем решение исходного уравнения по этой же формул е, счи-

òàÿ, ÷òî â íåé c = c(x) :

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

2x

 

 

 

2

 

2

 

 

 

~

~

~

 

c¢(1+ x

 

) + c× 2x -

 

 

c(1

+ x

 

) = 1

+ x

; c¢ = 1; c = òdx + c = x +

c

è ( c

çàìå-

 

1+ x2

 

íÿåì íà c) y = (x + c)(1+ x2 ) .

Теперь проверим, что функция (14.16) действительно является общим решением уравнения y¢ + P(x)y =Q(x).

Во-первых, при каждом значении c эта функция будет решением уравнения:

 

é

òQ(x)eò

P(x )dx

ù'

P(x )dx

é

òQ(x)eò

P(x )dx

ù

´

y¢+ P(x)y= ê

dx+cú e

 

ò

 

+ ê

 

dx+ cú

 

ë

 

 

û

 

 

 

ë

 

 

 

 

û

 

æ

òP(x )dx ö'

é

òQ(x)e

òP(x )dx

 

 

ù

òP(x )dx

 

 

´ç e

 

÷

+ P(x)ê

 

 

dx+cú e

=

 

 

è

 

ø

ë

 

 

 

 

 

 

û

 

 

 

 

=Q(x)e

òP(x)dx – òP(x )dx

é

 

òP(x )dx

 

 

ù

òP(x )dx

 

 

 

e

– Q(x)e

 

dx+c

ú

e

P(x)+

 

 

 

 

 

ê

 

 

 

 

 

 

 

 

 

 

 

 

 

ë

 

 

 

 

 

û

 

 

 

 

+ P(x)éQ(x)e

òP(x )dx dx+ cù eòP(x)dx

=Q(x).

ê

ú

 

ë

û

 

Во-вторых, зададим произвольное начальное условие y(x0) = y0 и покажем, что существует значение c, при котором функция (14.16) удовлетворяет этому начальному условию:

270

271

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

é

òQ(x)e

òP(x)dx

 

 

 

+ c

ù

òP(x)dx

= y0

;

 

 

 

ê

 

dx

 

x=x0

úe

 

 

 

 

 

 

ë

 

 

 

 

 

 

û

 

x=x0

 

c= y0 eòP(x)dx x=x0 - òQ(x)eòP(x)dxdx x=x0 .

4.Уравнения Бернулли

Такими уравнениями называются уравнения вида

y¢ + P(x)y =Q(x)yα ,

(14.19)

где a ¹ 0, a ¹ 1 (при a = 0 получаем линейное уравнение, при a = 1 – уравнение с разделяющимися переменными).

Легко убедиться, что к уравнениям Бернулли применим любой из описанных выше методов решения линейных уравнений.

Пример. Решить дифференциальное уравнение y¢ - 4x y = x y . Ð å ø å í è å

 

1-й способ.

 

 

 

 

 

 

 

 

 

 

 

 

 

y = uv; u¢v + uv¢ -

4

uv

 

 

 

 

æ

4

ö

 

 

x

= x uv; u¢v + uçv¢ -

x

v ÷ = x uv.

 

 

 

 

 

 

 

 

 

 

è

ø

 

Потребуем, чтобы v¢ -

4 v = 0 , тогда

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

dv

dx

dv

dx

 

 

 

 

 

 

 

 

 

 

v

= 4 x ; ò v = 4ò

x + c; ln | v |= 4ln | x | + c; ln | x |= ln | cx4 |; v = ±cx4; v = cx4 .

Ïðè c = 1

v = x4. Тогда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

du

dx

 

du

dx

 

 

 

u¢x4 = x ux2; u¢x =

u;

u = x ;

ò

 

= ò

x

+ c; 2 u = ln | x | + ln | c |= ln | cx |;

u

u = 14 ln2 | cx |; y = uv = 41 x4 ln2 | cx | .

2-й способ.

Решаем соответствующее однородное уравнение y¢ = 4x y. Имеем

dy

=

4

dy

=

4

ò

dy

= 4ò

dx

+ c; ln | y |= 4ln | x | + ln | c |;

dx

x y;

y

x dx;

y

x

272

14. Дифференциальные уравнения первого порядка...

ln | y |= ln | cx4 |; y = ±cx4; y = cx4.

Теперь ищем решение исходного уравнения по этой формуле, считая, что в ней c = c(x) :

4

 

3

4

4

 

4

 

 

c x

 

+ c4x

 

x cx

 

= x cx

 

; c x =

 

 

 

 

 

~

 

~

 

 

1

2

 

 

 

 

 

 

 

4 ln

2 c = ln x + ln c = ln c x ;

c =

 

 

dc

 

dx

 

ò

 

dc

= ò

dx

 

~

;

 

c;

 

=

x

;

 

 

 

x

+ c

 

c

 

 

c

 

 

 

 

 

 

 

 

 

 

 

 

 

~

 

 

 

 

 

 

4

 

1

 

 

4

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

c x ; y = c(x )x

 

=

4

x

 

ln

 

cx

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5. Уравнения в полных дифференциалах

Уравнение

P(x,y)dx +Q(x,y)dy = 0

(14.20)

называется уравнением в полных дифференциалах, если его левая часть является (полным) дифференциалом некоторой функции двух переменных, т.е. существует функция u(x,y), такая, что

du(x,y) = P(x,y)dx +Q(x,y)dy.

В этом случае уравнение (14.20) имеет вид du(x,y) = 0 , что выполняется в том, и только в том случае, когда u(x,y) = c , где с – некоторая (произвольная) постоянная. Последнее равенство является общим интегралом уравнения (14.20).

Учитывая изложенное выше, найдем условия, при которых ура в- нение (14.20) будет уравнением в полных дифференциалах, и при в ы- полнении этих условий укажем способ нахождения функции u(x,y).

Будем предполагать, что (x,y)ОD , где D – некоторая область на плоскости 0xy. Напомним некоторые определения, приведенные в разд. 12.1. Под словом «область» понимается открытое связное множе ство. Множество называется открытым, если вместе с каждой своей точкой оно содержит и некоторую окрестность этой точки. Множеств о называется связным, если любые две его точки можно соединить н епрерывной кривой, целиком принадлежащей множеству.

Область D называется односвязной, если вместе с каждым замкнутым самонепересекающимся контуром она содержит и област ь, ограниченную этим контуром (т.е. односвязная область – это обла сть без «дырок»). Область, изображенная на рис. 78, не является односв язной.

273

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

Ðèñ. 78

Теорема 14.2. Пусть функции P(x,y), Q(x,y), P(x,y) è Q(x,y)

y x

непрерывны в области D. Тогда, для того чтобы в этой области выражение P(x,y)dx +Q(x,y)dy являлось полным дифференциалом некоторой функции двух переменных, необходимо, а в предположении од носвязности области и достаточно, чтобы при (x,y) D

P(x,y)

=

Q(x,y).

(14.21)

y

 

 

x

 

¡ Необходимость. Пусть существует функция u(x,y) , такая, что для (x,y) D du(x,y) = P(x,y)dx +Q(x,y)dy. Тогда согласно свойствам дифференциала функции двух переменных,

 

 

 

 

 

 

P(x,y) =

 

 

u(x,y)

, Q(x,y) =

u(x,y)

,

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

откуда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

¶P(x,y)

 

¶ æ

¶u(x,y) ö

 

2u(x,y)

 

¶Q(x,y)

 

 

¶ æ ¶u(x,y) ö

2u(x,y)

 

 

 

 

=

 

ç

 

÷

=

 

 

 

 

 

;

 

 

 

=

 

ç

 

 

÷ =

 

.

 

¶y

 

¶x

 

 

¶y¶x

 

¶x

 

 

 

¶y

¶x¶y

 

 

¶y è

ø

 

 

 

 

 

 

 

 

 

¶x è

ø

 

Но по теореме 12.9

2u(x,y)

=

 

2u(x, y)

(две эти смешанные произ-

 

 

 

 

¶x¶y

 

 

 

 

 

 

 

 

 

¶y¶x

 

 

 

 

 

 

 

 

 

 

 

 

 

водные непрерывны в области D по условию теоремы), значит,

 

 

 

 

 

 

 

 

 

 

 

P(x,y)

 

=

Q(x,y).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

Достаточность. Для простоты проведем доказательство для слу- чая открытого прямоугольника (т.е. прямоугольника без гра ницы).

14. Дифференциальные уравнения первого порядка...

Пусть (x0,y0)ОD – произвольная фиксированная точка; (x,y) D – произвольная точка. Предположим, что нужная нам функция u(x,y) (т.е. такая, что du = Pdx +Qdy ) существует. Тогда

u(x,y) = P(x,y).

x

Проинтегрируем это равенство от x0 до x (для удобства переменную интегрирования тоже будем обозначать через x):

 

 

 

x

¶u(x,y)

 

x

 

 

 

 

x

x

 

 

ò

dx =

ò P(x,y)dx;

u(x,y)

 

= ò P (x,y)dx;

 

 

¶x

 

 

 

 

x0

 

 

 

 

 

 

 

 

 

 

x0

 

 

 

 

 

x0

 

 

 

 

 

x0

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

u(x,y) = ò P(x,y)dx +u(x0,y).

 

 

 

 

 

 

 

 

 

 

 

 

x0

 

 

 

 

Положим, что u(x0,y) = ϕ(y) , и используем второе равенство

u(x,y)

=

 

 

 

 

¶u(x,y)

 

æ x

ö¢

 

 

 

 

 

 

 

 

ç

÷

 

 

 

y

 

Q(x,y) :

 

 

 

 

 

 

=

ç ò P(x,y)dx ÷

 

+ j¢(y) =Q(x,y). Òàê êàê

 

 

 

 

¶y

 

 

 

 

 

 

P(x, y)

 

 

 

 

 

 

 

è x0

ø y

 

 

 

P(x, y)

è

непрерывны, то в этой формуле символы интегри-

 

 

 

 

 

y

 

 

 

 

 

 

 

 

 

 

 

 

рования и дифференцирования можно поменять местами (см. теорему

13.3.), тогда

òx ¶P(x,y)dx + j¢(y) = Q (x,y).

x0 ¶y

Последнююформулувсилуусловия(14.21)можнопереписатьвви де

x

¶Q(x,y)

 

 

 

 

 

x

 

 

ò

dx+ j¢(y)=Q (x,y);Q(x,y)

+j¢(y)=Q(x,y).

¶x

x0

 

 

 

 

 

 

 

x0

 

 

 

 

 

 

 

 

 

Отсюда Q(x,y)

Q(x0,y)

+ ϕ′

=

Q(x,y);

ϕ′

= Q (x,y).

 

(y)

 

(y)

 

Интегрируем это неравенство по y от y0 до y (для удобства переменную интегрирования тоже будем обозначать через y):

y

y

y

ò j¢(y)dy = òQ(x0,y)dy; j(y)

 

y

= òQ(x0,y)dy;

 

 

y0

 

y0

y0

y0

y

j(y)= òQ(x0,y) dy+j(y0).

y0

274

275

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

y

Таккакздесь ϕ(y0) = u(x0,y0) –постоянная,то j(y) = òQ(x0,y)dy +C è

 

x

 

y

 

y0

 

u(x,y)=

ò

P (x,y) dx+

ò

Q(x

0,y) dy+C.

(14.22)

 

 

 

 

 

x 0

y 0

 

 

Таким образом, если нужная нам функция u(x,y) существует, то она задается формулой (14.22). Теперь проверим, что эта формула действительно дает функцию u(x,y) , такую, что du = Pdx +Qdy , èëè

u

= P, u

= Q.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

В соответствии с теоремой 10.6

 

 

 

 

 

 

 

 

 

 

 

¶u(x,y)

 

 

æ x

 

ö¢

 

 

 

 

 

 

 

 

 

 

=

ç

ò

 

P (x,y) dx÷

 

= P(x,y).

 

 

 

 

 

 

 

¶x

 

 

 

 

 

 

 

 

 

ç

 

÷

 

 

 

 

 

 

 

 

 

 

 

è x0

ø x

 

 

 

 

 

 

По той же теореме и по теореме 13.3

 

 

 

 

 

 

 

¶u (x,y)

æ x

 

ö'

 

 

æ

y

 

 

ö'

 

x

¶P (x,y)

 

 

 

ç

ò

 

÷

 

 

ç

ò

Q(x0,y) dy

÷

 

ò

 

0,y).

 

¶y

=ç

P (x,y) dx ÷

+ç

÷

=

¶y

dx+Q(x

 

ç

 

÷

 

 

ç

 

 

÷

 

 

 

 

 

è x0

øy

 

è y0

 

øy

 

x0

 

 

 

Отсюда, в силу условия (14.21)

¶u(x,y)

x

¶Q(x,y)

 

 

x

 

 

= ò

¶x

dx +Q(x0

,y) =Q(x,y)

 

x0

+Q(x0,y) =

¶y

 

 

 

 

 

x0

 

 

 

 

 

 

=Q(x,y)Q(x0,y)+Q(x0,y) =Q(x,y).x

Формула (14.22) соответствует следующему принадлежащему пря - моугольнику D пути от точки (x0,y0) до точки (x,y) (рис. 79).

Аналогично,начинаясравенства u(x,y) = Q(x,y),можемполучить

другую формулу для u(x,y):

y

x

y

 

u(x,y) = ò P(x,y0)dx + òQ(x,y )dy +C ,

(14.23)

x0

y0

 

соответствующую пути, указанному на рис. 80.

Для областей D более сложного вида (чем прямоугольник со сторонами, параллельными осям координат) формулы (14.22) и (14.23) дают

14. Дифференциальные уравнения первого порядка...

 

(x,y)

(x,y)

 

 

(x0,y0)

 

(x0,y0)

 

 

Ðèñ. 79 Ðèñ. 80

u(x,y) во всех точках (x,y) , которые можно соединить с точкой (x0,y0) ломаной такой формы.

Теорема 14.2 и формулы (14.22) и (14.23) являются ответами на задачи, сформулированные в начале этого пункта.

Пример. Решить дифференциальное уравнение

(x3 y2x)dx + (y3 x2 y)dy = 0.

Ðå ø å í è å

Âэтом примере P(x,y) = x3 y2x, Q(x,y) = y3 x2y непрерывны вместе со своими частными производными на всей плоскости 0xy.

P(x,y)

= −2yx,

Q(x,y)

= −2xy , ò.å.

P(x,y)

=

Q(x,y)

. Значит, это

y

x

y

x

 

 

 

 

уравнение является уравнением в полных дифференциалах (к стати, отметим, что оно также является однородным).

Âçÿâ x0 = y0 = 0 , из формулы (14.22) получаем (достаточно знать одну функцию u(x,y) , поэтому берем C = 0 ):

x

y

 

x4

 

x

 

 

x2

 

x

 

y4

 

y

 

 

x4

 

x2 y2

 

y4

 

 

 

 

 

 

 

 

 

 

 

 

 

u(x, y) = ò

(x3 y2 x) dx + ò

( y3 0) dy =

 

y2

 

 

+

 

 

=

 

+

.

0

0

4

 

0

 

 

 

2

 

0

 

4

 

0

 

 

4

2

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x4

 

 

 

x2y2

 

 

y4

 

 

 

 

 

 

Общий интеграл уравнения имеет вид

 

 

 

 

 

 

 

 

+

 

 

= C , èëè

 

 

 

4

 

 

 

2

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x4 2x2 y2 + y4 = C.

14.4. Дифференциальные уравнения высших порядков

Рассмотрим следующую задачу: пусть точка движется вдоль о сиx0 с постоянным ускорением a. Найти x(t) – положение точки в произвольный момент времени t.

По условию примера x′′(t) = a , что является дифференциальным уравнением второго порядка. Далее имеем x¢(t) = òadt + c1 = at + c1;

276

277

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

x(t) = ò(at + c1)dt + c2 = a t22 + c1t + c2 – любое решение нашего уравнения.

Это решение зависит от двух произвольных постоянных c1 è c2. Для того чтобы решение было единственным, нужно задать положе ние точ-

ки x(t) и ее скорость x¢(t) в некоторый момент времени t0: x(t0) = x0 ,

x¢(t

0

) = x¢

, ãäå x

0

и x¢ – некоторые числа. Тогда

 

 

 

 

 

 

 

 

0

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

¢(

)

= at0 + c1

¢

¢

- at0 è x(t0) = a

t2

 

 

 

 

 

 

 

 

0

+ c1t0 + c2 = x0 Þ

 

 

t0

 

= x0 Þ c1 = x0

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t2

 

 

 

t2

 

 

 

 

 

 

 

 

 

 

t2

 

 

 

c = x - a

0

 

- c t

= x - a

0

- (x¢

- at

0

)t

0

= x

0

- x¢t

+ a

0

.

 

 

 

 

 

 

 

 

 

 

2

0

 

 

2

1 0

0

2

0

 

 

 

 

0 0

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Отталкиваясь от этого примера, по аналогии с разд. 14.2 привед ем следующие определения и формулировки.

Дифференциальное уравнение n-го порядка имеет вид

F (x,y,y¢,y¢¢,...,y(n) ) = 0.

Если из этого уравнения можно выразить старшую производн уюy(n) , то мы получим так называемое уравнение, разрешенное относ ительно старшей производной

y(n) = f (x,y,y¢,...,y(n1)),

(14.24)

где f – некоторая функция (n+1)-й переменной.

Определение 14.6. Задачей Коши для уравнения 14.24 называется задача

 

 

 

 

 

ìy(n) = f (x,y,y¢,...,y(n1) );

 

 

 

 

 

 

ï

 

= y0;

 

 

 

 

 

 

 

ïy(x0)

 

 

 

 

 

 

 

ïy¢(x ) = y¢ ;

 

 

 

 

 

 

 

í

0

 

0

 

 

 

 

 

 

 

ï.....................................

(14.25)

 

 

 

 

 

ï

 

 

 

 

 

 

 

 

 

 

ïy(n1)(x

) = y(n1)

,

 

 

 

 

 

 

î

 

0

0

 

 

ãäå

x ,y ,y¢

,...,y(n1)

– некоторые числа.

 

 

 

0

0

0

0

 

 

 

 

 

 

Теорема 14.3 (существования и единственности решения задачи Коши). Пусть функция f и ее частные производные первого порядка по всем аргументам, кроме x, непрерывны в некоторой области D

14. Дифференциальные уравнения первого порядка...

(n+1)-мерного пространства и точка (x ,y ,y¢

,...,y(n1))ОD . Тогда

0

0

0

0

задача Коши (14.25) имеет решение, и притом единственное. Это решение определено в некоторой окрестности точки x0 .

Доказательство этой теоремы (как и теоремы 14.1) здесь не при водится.

Замечание. В следующей главе будут рассмотрены так называемые линейные уравнения высших порядков, т.е. уравнения вида

y(n) + P1(x)y(n1) + P2(x)y(n2) + ...+ Pn(x)y = f (x),

где все коэффициенты Pi (x), i = 1, 2,...,n , и правая часть f (x) определены и непрерывны на некотором интервале (a, b). Такие уравнения могут быть записаны в виде

y(n) = f (x) - P1(x)y(n1) - P2(x)y(n2) - ...- Pn (x)y

и, естественно, удовлетворяют всем условиям теоремы 14.3 в об ластиD, в которой x (a, b) , а остальные переменные – любые. Особенностью данного случая является то, что (как можно доказать) единстве нное решение задачи Коши (14.25) будет при этом определено на всем интервал еa(, b) (а не только в окрестности точки x0).

Определение 14.7. Пусть выполняются условия теоремы 14.3. Функция

y = j(x,c1,c2,...,cn ),

(14.26)

ãäå c1, c2,..., cn – постоянные, называется общим решением уравнения (14.24) в некоторой окрестности точки x0 U (x0) , åñëè:

1. Ïðè "x ÎU (x0) и " наборе (c1,c2,...,cn)ÎC , где C – некоторое множество (в простых случаях c1, c2,..., cn будут любыми числами),

функция (14.26) является решением уравнения (14.24).

2. Любое решение уравнения (14.24), такое, что x ОU (x0) è (x,y,y¢,...,y(n1))ОD , получается из формулы (14.26) при некотором

наборе (c1,c2,...,cn)ÎC .

Точно так же, как в разд. 14.2, показывается, что в этом опреде-

лении условие 2 можно заменить на условие 2¢. Какие бы начальные условия y(x0) = y0 , y¢(x0) = y0¢ , ... , y(n1)(x0) = y0(n1), где точка

278

279

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

(x ,y ,y¢

,...,y(n1))ОD, ни задались, существует набор (c ,c ,...,c )ОC,

0

0

0

0

1

2

n

при котором функция (14.26) удовлетворяет этим начальным усло виям. Ниже будем проверять именно это условие.

Определение 14.8. Равенство вида Φ(x,y,c1,c2,...,cn) = 0 , неявно задающее общее решение, называется общим интегралом дифференциального уравнения (14.24).

14.5. Уравнения, допускающие понижение порядка

Одним из методов решения дифференциальных уравнений выс ших порядков является сведение их к дифференциальным уравне ниям меньшего (лучше всего – первого) порядка.

Рассмотрим два основных типа уравнений, допускающих пони жение порядка.

1. Уравнение не содержит явным образом искомую функцию y и, может быть, несколько ее первых производных, т.е. имеет вид

F (x,y(k),y(k+1),...,y(n) ) = 0.

(14.27)

Сделаем в этом уравнении замену y(k) = z , где z = z(x) – новая неизвестная функция (т.е. за новую неизвестную функцию берется производная наименьшего порядка, входящая в это уравнение). Т огда y(k+1) = z¢ , y(k+2) = z¢¢, ..., y(n) = z(nk) и (14.27) принимает вид F(x, z, z,..., z(nk)) = 0. Таким образом, мы сумели уменьшить (или понизить) порядок уравнения.

Пример. Решить задачу Коши для дифференциального уравнения

ìx2y¢¢ = y¢2;

ï

ïy(1) = 1;

í

ïïy¢(1) = 1. î 2

 

 

 

 

 

 

 

Ð å ø å í è å

 

 

 

 

 

 

 

Обозначая y¢ = z,

z = z(x) , имеем x2z¢ = z2.

Это уравнение первого

порядка является уравнением с разделяющимися переменны ми:

dz

 

 

dz

dx

. Интегрируя, получаем

 

 

 

 

 

 

 

x2 dx = z2; z2

= x2

 

 

 

 

 

 

 

 

 

 

dz

dx

 

1

1

1 1

 

 

1

+ c1x

 

x

 

 

ò z2 = ò x2

+ c1;

- z

= - x + c1;

z = x

+ c1

=

 

 

; z =

 

,

 

 

x

1+ c x

 

 

 

x

 

 

 

 

 

 

 

 

 

1

 

ò.å. y¢ =

 

 

.

 

 

 

 

 

 

 

 

 

 

 

1

+ c x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

280

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

14. Дифференциальные уравнения первого порядка...

Прежде чем интегрировать еще раз, найдем c1 из второго начального ус-

ловия. При x = 1 имеем

 

 

 

1

 

=

1

; 1+ c1 = 2; c1 = 1.

Значит, y¢ =

 

x

. Отсюда

1+ c1

2

 

 

1+ x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y =

x

dx + c =

x +1-1

dx + c = dx -

 

 

 

dx

+ c = x - ln |1

+ x | +c .

 

 

 

 

 

 

 

 

 

 

 

 

 

ò1+ x

 

 

 

2

ò

 

 

1+ x

 

 

 

2

ò

 

 

ò1+ x

2

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Подставляя сюда x = 1 , из первого начального условия находим по-

стоянную c2: 1- ln2 + c2 = 1;

c2 = ln2. Таким образом,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y = x - ln |1+ x | + ln2 = x - ln

 

1+ x

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

Если бы нужно было найти общее решение исходного уравнени я, то

y =

 

 

 

x

 

dx + c2

=

1 1+ c1x -1

dx + c2

=

1

 

 

 

dx -

1

 

dx

 

+ c2 =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ò1

+ c1x

c1

ò 1 + c1x

 

c1 ò

c1

ò1 + c1x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

1

 

 

 

d(1+ c1x)

 

 

1

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

x -

 

ò

 

 

+ c2

=

 

x -

 

ln |1+ c1x | +c2,

 

 

 

 

 

 

 

 

c

c2

 

1+ c x

c

c2

 

 

 

 

 

 

 

 

1

 

1

 

 

 

 

 

 

1

 

 

 

1

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

÷òî ïðè c1 ¹ 0 дает общее решение; если же c1 = 0 , òî y¢ = x è y = òxdx + c = x22 + c.

2. Уравнение не содержит явным образом независимую переме н- ную х, т.е. имеет вид

F(y, y, y′′,..., y(n)) = 0.

(14.28)

Сделаем в этом уравнении замену y¢ = z , где z = z(y), т.е. за новую независимую переменную возьмем y, а за новую независимую функцию – y¢ = z(y). Покажем, что при такой замене порядок уравнения понижается на единицу:

y¢¢ = (y¢)¢x = zx¢ = zy¢ × yx¢ = z¢z , ò.å. y¢¢ = z¢z ,

y¢¢¢ = (y¢¢)¢x = (z¢z)¢x = (z¢z)¢y × yx¢ = (z ¢¢z + z ¢2)z = z ¢¢z 2 + z ¢2z,

yIV = (y¢¢¢)¢x = (z¢¢z2 + z¢2z)¢y yx¢ = (z ¢¢¢z 2 + z¢¢2zz ¢ + 2z ¢z ¢¢z + z ¢3 )z =

=z¢¢¢z3 + 4z¢¢z¢z2 + z¢3z

èтак далее, т.е. порядок каждой производной становится на е диницу меньше.

281