А. Р. Лакерник. Высшая математика.Краткий курс
.pdf5. Дифференциальные уравнения
шение получается из (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|= e−2x−2c; 3 -2u = ±e−2xe−2c; |
||||
2 ln | 3 |
||||||
|
|
|
|
u = ± 1 e−2ce−2x + |
3 . |
|
|
|
|
|
2 |
2 |
|
Заменяя здесь ± |
1 e−2c |
на c, имеем: u = ce−2x |
+ 3 |
, откуда |
||
|
|
|
2 |
|
2 |
|
y= 12 æç3x - ce−2x - 23 + 5ö÷ = - 2c e−2x + 23 x + 47 ,
èø
или, заменяя - 2c íà c,
y = ce−2x + 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(n−1)), |
(14.24) |
где f – некоторая функция (n+1)-й переменной.
Определение 14.6. Задачей Коши для уравнения 14.24 называется задача
|
|
|
|
|
ìy(n) = f (x,y,y¢,...,y(n−1) ); |
|
||||
|
|
|
|
|
ï |
|
= y0; |
|
|
|
|
|
|
|
|
ïy(x0) |
|
|
|||
|
|
|
|
|
ïy¢(x ) = y¢ ; |
|
|
|||
|
|
|
|
|
í |
0 |
|
0 |
|
|
|
|
|
|
|
ï..................................... |
(14.25) |
||||
|
|
|
|
|
ï |
|
|
|
|
|
|
|
|
|
|
ïy(n−1)(x |
) = y(n−1) |
, |
|
||
|
|
|
|
|
î |
|
0 |
0 |
|
|
ãäå |
x ,y ,y¢ |
,...,y(n−1) |
– некоторые числа. |
|
|
|||||
|
0 |
0 |
0 |
0 |
|
|
|
|
|
|
Теорема 14.3 (существования и единственности решения задачи Коши). Пусть функция f и ее частные производные первого порядка по всем аргументам, кроме x, непрерывны в некоторой области D
14. Дифференциальные уравнения первого порядка...
(n+1)-мерного пространства и точка (x ,y ,y¢ |
,...,y(n−1))ОD . Тогда |
||
0 |
0 |
0 |
0 |
задача Коши (14.25) имеет решение, и притом единственное. Это решение определено в некоторой окрестности точки x0 .
Доказательство этой теоремы (как и теоремы 14.1) здесь не при водится.
Замечание. В следующей главе будут рассмотрены так называемые линейные уравнения высших порядков, т.е. уравнения вида
y(n) + P1(x)y(n−1) + P2(x)y(n−2) + ...+ Pn(x)y = f (x),
где все коэффициенты Pi (x), i = 1, 2,...,n , и правая часть f (x) определены и непрерывны на некотором интервале (a, b). Такие уравнения могут быть записаны в виде
y(n) = f (x) - P1(x)y(n−1) - P2(x)y(n−2) - ...- 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(n−1))ОD , получается из формулы (14.26) при некотором
наборе (c1,c2,...,cn)ÎC .
Точно так же, как в разд. 14.2, показывается, что в этом опреде-
лении условие 2 можно заменить на условие 2¢. Какие бы начальные условия y(x0) = y0 , y¢(x0) = y0¢ , ... , y(n−1)(x0) = y0(n−1), где точка
278 |
279 |
5. Дифференциальные уравнения
(x ,y ,y¢ |
,...,y(n−1))О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(n−k) и (14.27) принимает вид F(x, z, z′,..., z(n−k)) = 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