Конспект лекций по обыкновенным дифференциальным уравнениям. Ч. 1. Учебное пособие
.pdf
|
|
|
|
. |
|
|
|
|
|
|
. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
x |
xs ( y ys ), |
y ys (x xs ) |
|
|
|
|
(d) |
||||||||||
Положим |
|
|
x a2 , |
y |
b2 , |
(x x ) U , |
y y |
s |
V |
|
|
||||||||||
|
|
|
|
s |
|
|
s |
|
|
|
|
|
|
|
s |
|
|
|
|
|
|
Тогда система (d) примет вид: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
. |
|
|
|
. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
U a2V , |
V |
b2U |
|
|
|
|
|
(e) |
|||||||
Продифференцируем по t, например, первое уравнение и подставим в него второе |
|||||||||||||||||||||
|
|
|
|
|
|
|
.. |
a2b2U 0 . |
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
U |
|
|
|
|
|
|
|
|||||||
Легко проверить, что общим решением этого уравнения является |
|
|
|||||||||||||||||||
|
|
|
|
|
U c1 cos(abt) c2 sin(abt) , |
|
|
|
|
(f) |
|||||||||||
где с1 и с2 - постоянные, определяемые из начальных условий (b). Из (e) |
|||||||||||||||||||||
получаем: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
V |
1 |
. |
b |
c sin(abt) c cos(abt) . |
|
|
|
|
|||||||||
|
|
|
|
U |
|
|
|
(g) |
|||||||||||||
|
|
|
|
|
a2 |
|
a |
|
1 |
|
|
|
|
|
2 |
|
|
|
|
|
|
Исключая из (f) и (g) время, имеем: |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
U 2 |
|
V 2 |
|
c2 |
c2 |
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
a2 |
b2 |
|
1 |
a2 |
2 |
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
Возвращаясь к исходным обозначениям, получим |
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
(x x )2 |
|
( y y |
|
)2 |
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
s |
|
|
|
|
s |
|
1, |
|
|
|
|
(h) |
||
|
|
|
|
|
|
c2 x |
|
c2 y |
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
s |
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
s |
|
|
|
|
|
|
|
|
|
|
|
||
где |
с2 с2 |
. |
Соотношение |
|
(h) |
– |
это |
|
|
|
|
|
|
||||||||
c2 1 |
2 |
|
y |
|
|
|
|
|
|||||||||||||
|
a2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
уравнение эллипсов (в зависимости от с2, т.е. |
|
|
|
|
|
|
|||||||||||||||
от |
начальных |
условий) |
с |
центром |
|
в точке |
y1 |
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(xs , ys ) ). Стрелки указывают направление, |
|
|
|
|
|
|
|||||||||||||||
соответствующее возрастанию времени. Из (d) |
|
|
|
|
|
|
|||||||||||||||
видно, что при |
x xs , растет y |
(т.к. >0) |
и |
|
|
|
x |
s |
x |
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
наоборот. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
11 |
|
|
|
|
|
|
|
|
|
|
Видно, что изменение популяций носит циклический характер: через
определенное время популяция возвращается к первоначальному уровню.
1.2.1 Уравнения 1-ого порядка, интегрируемые в квадратурах
Рассмотрим теперь некоторые уравнения 1-ого порядка, которые
интегрируются в квадратурах.
1) Уравнения с разделяющимися переменными.
Уравнением с разделяющимися переменными называется уравнение вида:
|
dy |
f (x)g( y) |
(13) |
|
dx |
||
|
|
|
|
будем считать, что функции f (x) , g( y) непрерывны при a x b |
,c y d , и |
||
g( y) 0 , т.к. если, например, в какой-то точке y1 , функция g( y1 ) 0 , то из (13)
следует, что y y1 есть решение. Запишем уравнение (13) в виде:
dy f (x)dx g( y)
и проинтегрируем
|
dy |
f (x)dx C . |
(14) |
|
g( y) |
||||
|
|
|
Это и есть общее решение уравнения (13).
Пример 7. Найдем решение следующего д.у. с разделяющимися переменными:
|
y y2 . |
|||||
Перепишем уравнение в виде: |
|
dy |
|
dx . |
||
|
y2 |
|||||
|
|
|
||||
После интегрирования получим |
|
1 |
x C . |
|||
y |
||||||
|
|
|
|
|||
Тогда решение уравнения имеет вид: |
y 1 (x C) , |
|||||
где С - произвольная постоянная. Кроме того, имеется еще решение y 0.
12
Найдем теперь решение следующей задачи (такую задачу называют задачей
Коши): y y2 , y(x0 ) y0 |
при y0 0 . |
Определив значение C из начального |
||||
условия, получим: |
|
|
|
|
|
|
|
y(x) |
|
|
|
y0 |
|
|
|
|
|
. |
||
|
1 y |
0 |
(x x ) |
|||
|
|
|
|
0 |
|
|
2)Уравнения, приводящиеся к уравнениям с разделяющимися переменными.
Втаких уравнениях переменные можно разделить, выполнив некоторую
подходящую замену.
а) Рассмотрим уравнение вида:
dydx f (ax by) ,
где a,b const .
Перейдем к новым переменным z ax by , в результате будем иметь
|
dz |
a b |
dy |
, т.е. |
dz |
a bf (z) |
||
|
|
|
dx |
|||||
|
dx |
|
|
dx |
|
|||
или |
|
|
|
|
|
|
|
|
|
|
|
|
dz |
dx , |
|||
|
|
|
|
|
||||
|
|
|
a bf (z) |
|||||
т.е. переменные разделились. Интегрируя, получим:
dz
x a bf (x) C .
б) Рассмотрим теперь следующее уравнение
dy |
|
y |
|
|
|
f |
|
|
(x, y) , |
|
|
|||
dx |
x |
|
||
где (x, y) однородная функция нулевого порядка. Такое уравнение называют
однородным д.у. 1-ого порядка.
Функция (x, y) называется однородной функцией нулевого порядка, если
(tx,ty) t 0 (x, y) (x, y)
13
Выполним замену y
x z , т.е. y xz , тогда
dy |
x |
dz |
z, x |
dz |
z f (z), |
dz |
|
dx |
. |
|
|
|
f (z) z |
|
|||||
dx |
|
dx |
|
dx |
|
|
x |
||
Таким образом, мы снова получили д.у. с разделенными переменными. После интегрирования получим:
|
dz |
|
|
|
|
|
|
|
|
|
dz |
|
||
|
|
|
|
|
|
|
|
|||||||
|
|
ln |
x |
ln C, x C exp |
|
|
. |
|||||||
|
|
|
||||||||||||
f (z) z |
|
|
|
|
|
|
|
|
|
f (z) z |
||||
|
|
|
|
|
|
|
|
|||||||
в) Уравнение вида |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
M (x, y)dx N(x, y)dy 0 |
|
|
||||||||||
будет являться однородным, если |
|
|
|
|
|
|
|
|||||||
|
|
|
dy |
|
M (x, y) |
|
y |
|
|
|||||
|
|
|
|
|
f |
|
, |
|
|
|||||
|
|
|
|
|
|
|
|
|||||||
|
|
|
dx |
|
|
|
N (x, y) |
|
x |
|
|
|||
т.е. |
M (x, y) , N (x, y) |
- |
однородные |
функции |
одного |
порядка, |
||
( M (tx,ty) t n M (x, y), |
N(tx,ty) t n N(x, y) ), |
и |
будет |
уравнением |
с |
|||
разделяющимися переменными, если |
|
|
|
|
|
|||
MM1(x)M2 ( y), N N1(x)N2 ( y) .
2.Линейные дифференциальные уравнения 1-го порядка.
Линейным дифференциальным уравнением 1-го порядка называется уравнение 1-го порядка, линейное относительно неизвестной функции и еѐ
производной.
|
|
|
dy |
p(x) y f (x) |
(1) |
|
|
|
|
|
|
||
|
|
|
dx |
|
|
|
где p(x), f (x) – непрерывны в области, где ищется решение уравнения (1). |
|
|||||
Если |
f (x) 0 , то |
уравнение (1) называется линейным однородным, |
в |
|||
противном |
случае это |
уравнение называется |
линейным неоднородным. |
В |
||
14
линейном однородном уравнении переменные разделяются, и его решение не составляет труда:
|
dy |
p(x) y 0, т.е. |
dy |
|
p(x)dx |
|
||||
|
dx |
y |
|
|||||||
|
|
|
|
|
|
|
|
|||
Таким образом, |
|
|
|
|
|
|
|
|
|
|
|
|
ln |
|
y |
|
p(x)dx ln C, |
C 0, |
|
||
|
|
|
|
|
||||||
|
|
|
|
|
|
|||||
и общее решение однородного уравнения имеет вид: |
|
|||||||||
|
|
y C exp{ p(x)dx}, C 0 |
(2) |
|||||||
При делении на y мы потеряли решение y=0, однако оно может быть включено в (2), если считать, что C может быть равно нулю.
Для интегрирования же неоднородных уравнений применяются метод вариации произвольной постоянной и метод подстановки, которые и будут рассмотрены ниже.
2.1 Метод вариации произвольной постоянной
Для интегрирования неоднородного уравнения (1) можно применять метод вариации произвольной постоянной.
В соответствии с этим методом вначале находится решение соответствующего однородного уравнения, которое задается формулой (2). Затем в полученном решении «варьируют постоянную», т.е. полагают C C(x) , и
подставляют это решение в исходное неоднородное уравнение (1). Тогда,
учитывая что:
dydx dCdx e p( x)dx C(x) p(x)e p(c)dx ,
будем иметь:
dCdx e p( x)dx C(x) p(x)e p( x)dx p(x)C(x)e p( x)dx f (x) ,
или
15
dCdx f (x)e p( x)dx .
Таким образом, мы получили линейное однородное д.у. относительно C(x),
интегрируя которое, находим
C(x) f (x)e p( x)dxdx C1 ,
и, следовательно,
y C(x)e |
p( x)dx |
C1e |
p( x)dx |
e |
p( x)dx |
f (x)e |
p( x)dx |
dx . |
(3) |
|
|
|
|
Итак из формулы (3) видно, что общее решение неоднородного линейного дифференциального уравнения 1-го порядка равно сумме общего решения соответствующего однородного уравнения
С1e p( x)dx
и частного решения неоднородного уравнения
e p( x)dx f (x)e p( x)dxdx ,
получающегося из (3) при С1 0 .
Пример 1. Методом вариации произвольной постоянной найдем решение неоднородного уравнения
|
|
|
|
|
|
|
|
|
|
|
|
dy |
|
y |
x2 . |
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
dx |
|
x |
|
|
|
|
|
|
|
|
|
|
|
|
|||
Вначале интегрируем соответствующее однородное уравнение: |
|||||||||||||||||||||||||||||
|
dy |
|
y |
0, т.е. |
dy |
|
|
dx |
, следовательно ln |
|
y |
|
ln |
|
x |
|
ln C; y Cx . |
||||||||||||
|
|
|
|
|
|
||||||||||||||||||||||||
|
|
|
|
|
|||||||||||||||||||||||||
|
dx |
|
x |
|
|
y |
|
x |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
dy |
|
dC |
|
|
|
||||||||||||||
Затем полагаем |
C C(x) , |
тогда y C(x) x, |
и |
|
x C(x) . Подставляя эти |
||||||||||||||||||||||||
|
|
||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
dx |
|
dx |
||||||||
выражения в исходное уравнение, получим: |
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||
|
|
|
|
|
dC |
x C(x) C(x) x2 , т.е. |
dC |
x, |
dC xdx , |
||||||||||||||||||||
|
|
|
|
|
|
|
|||||||||||||||||||||||
|
|
|
|
|
dx |
|
|
|
|
|
|
|
dx |
|
|
|
|
|
|
|
|
|
|
|
|
||||
и после интегрирования будем иметь:
16
C(x) x2 C1. . 2
Таким образом, общее решение имеет вид:
y C1x x3 . 2
Пример 2. Методом вариации произвольной постоянной найдем решение неоднородного уравнения:
dydx yctgx 2x sin x .
Интегрируя однородное уравнение, находим
dyy cossin xx dx , т.е. dydx yctgx 0, ln y ln sin x ln C, y C sin x .
Варьируем постоянную:
y C(x)sin x, y C (x)sin x C(x)cos x .
Подставим эти выражения в исходное уравнение:
C (x) sin x C(x) cos x C(x) cos x 2x sin x ,
получим: |
|
C(x) x |
2 |
C1 |
, |
C (x) 2x, |
|
||||
Откуда |
y x2 sin x C sin x . |
|
|||
|
|
1 |
|
|
|
2.2 Метод подстановки (Даламбера)
Для решения линейного неоднородного уравнения 1-ого порядка (1) можно применять также метод подстановки (Даламбера). В соответствии с этим методом в уравнении (1) применяют подстановку y(x) U (x)V (x) , тогда это уравнение примет вид:
V ( dUdx pU ) (U dVdx f ) 0.
Выберем функцию U так, чтобы первая скобка обратилась в нуль:
dUdx p(x)U 0, или dUU p(x)dx .
17
Обозначим через U0 (x) одно из частных решений этого уравнения, тогда для нахождения V имеем уравнение:
U0 (x) dVdx f (x) 0 .
Находим общее решение этого уравнения V V (x,C) . Теперь общее решение уравнения (1) будет иметь вид:
y(x) U0 (x)V (x, C) .
Пример 3. Методом подстановки найдем решение следующего
неоднородного уравнения:
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
y2 |
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
y |
2xy 3 . |
|
|
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
Перепишем это уравнение в виде: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
dx |
|
|
2x |
|
|
3 |
|
. |
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
dy |
|
|
|
y |
|
|
|
y2 |
|
|
|
|
|||||||
Оно линейно относительно x и |
|
dx |
|
. Положим x U ( y)V ( y) , тогда имеем: |
||||||||||||||||||||||||||||||||
|
dy |
|||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
V ( |
dU |
|
|
2 |
U ) (U |
dV |
|
3 |
) 0 . |
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
dy |
|
y |
|
|
|
|
|
|
|
|
dy |
y2 |
|
|
|||||||||||
Теперь |
dU |
|
2 |
U 0 , или |
dU |
|
|
|
2dy |
, |
т.е. U y2 C . |
|
|
|||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||
|
dy |
|
|
|
y |
|
|
|
|
U |
|
|
|
|
|
|
|
y |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
Положим U |
0 |
y2 |
, тогда |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
y2 |
dV |
|
|
3 |
0, или |
|
dV |
3 |
|
dy , т.е. V ( y, C) |
1 |
C , |
||||||||||||||||||||||
|
|
|
|
|
|
|
y4 |
|
y3 |
|||||||||||||||||||||||||||
|
|
|
|
|
dy |
|
y2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
откуда окончательно имеем
x Cy 2 1y .
2.3 Уравнения, сводящиеся к линейным
18
Многие нелинейные дифференциальные уравнения после замены переменных могут быть сведены к линейным.
1) Уравнение Бернулли
Уравнением Бернулли называется нелинейное относительно неизвестной
функции уравнение вида:
|
dy |
p(x) y f (x) yn , (n 1) , |
|
|||||||
|
dx |
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
или |
|
|
|
|
|
|
|
|
|
|
|
y n |
dy |
p(x) y1 n |
f (x) . |
(4) |
|||||
|
dx |
|||||||||
|
|
|
|
|
|
|
|
|
|
|
Выполним замену переменной |
y1 n z, |
|
|
|
|
|||||
тогда |
(1 n) y n |
dy |
|
dz |
. |
|
||||
|
|
|
||||||||
|
|
|
|
|
dx |
|
dx |
|
||
Подставив эти выражения в (4), получим линейное дифференциальное уравнение:
|
|
|
1 dz |
|
p(x)z f (x) . |
||||||||||||||
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
1 n dx |
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
Пример 4. Найдем решение следующего уравнения Бернулли: |
|||||||||||||||||||
|
|
|
|
dy |
|
|
y |
|
|
x2 |
. |
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
dx |
|
|
|
|
2x 2 y |
||||||||||
Перепишем его в виде: |
2 y |
dy |
|
y2 |
|
|
x2 , (n 1) . |
||||||||||||
|
|
|
|||||||||||||||||
|
|
dx |
x |
|
|
|
|
|
|
|
|
|
|
||||||
Выполним замену: y2 z. |
Тогда 2 y |
dy |
|
|
dz |
, и мы получим линейное уравнение |
|||||||||||||
|
|
|
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
dx |
|
|
dx |
|||||
из примера 1: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
dz |
|
z |
|
x2 , |
||||||||||
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
dx |
|
|
x |
|
|
|
|
|
||||||
которое может быть решено как методом вариации произвольной постоянной, так и методом подстановки.
19
Отметим, что уравнение Бернулли можно решать и непосредственно методом подстановки.
2) Уравнение Риккати
Уравнением Риккати называется нелинейное относительно неизвестной
функции уравнение вида:
dydx p(x) y q(x) y2 f (x) .
В общем виде это уравнение не интегрируется в квадратурах, однако,
заменой переменной может быть преобразовано в уравнение Бернулли, если известно одно частное решение y1 (x) этого уравнения.
Полагая y y1 |
z , получим: |
|
|
|||
|
|
|
|
|
|
|
|
|
y1 z p(x)(y1 z) q(x)(y1 z |
2 |
) f (x) , |
||
|
|
|
||||
но так как |
|
p(x) y1 q(x) y1 f (x) , то для z(x) получим уравнение Бернулли: |
||||
y1 |
||||||
|
|
|
|
2 |
|
|
z p(x) 2y1(x)q(x) z q(x)z2 0 .
Пример 5. Найдем решение следующего уравнения Риккати:
dydx y2 x22 .
Здесь легко подобрать частное решение: y1 1
x . Полагая y1 1
x z , получим:
|
|
|
|
|
1 |
|
|
|
|
|
|
1 |
|
1 2 |
2 |
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
y |
z |
x2 |
, т.е. |
z |
x2 |
|
(z x) |
x2 , |
||||||||||||
|
|
|
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
z |
|
|
|
|
||
или: |
|
|
|
|
|
z |
z |
|
2 x |
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|||||||||||
Последнее уравнение является уравнением Бернулли. Перепишем его в виде
z 2 1, z2 xz
Выполним замену:U 1z , тогда dUdx zz2 , а уравнение примет вид:
dUdx 2Ux 1; .
20
