книги / Математика. Дифференциальные уравнения
.pdfб) Найдем общее решение соответствующего однородного уравнения. Составляем характеристическое уравнение:
k4 + 2k2 +1 = 0 ,
|
|
|
|
( |
) |
|
|
|
|
|
|
|
|
k2 +1 2 = 0 . |
|
|
|
Тогда, k1,2 = ±i, k3,4 |
= ±i . |
|
|
|
||||
Согласно п. 6 §2 |
|
|
|
|
|
|
||
y = e0 x cos x = cos x, y |
2 |
= e0 x sin x =sin x, y = x cos x, |
y |
4 |
= xsin x . |
|||
1 |
|
|
|
3 |
|
|
||
|
y =C1 cos x +C2 sin x +C3 x cos x +C4 xsin x |
|
|
|
||||
|
|
|
|
|
или |
|
|
|
|
|
=C1 cos x +C2 sin x + x(C3 cos x +C4 sin x) . |
|
|
|
|||
|
y |
|
|
|
||||
f (x) =5sin x − xcos x |
или f (x) = e0 x (5sin x − x cos x) . |
Так как α ±iβ = 0 ±i |
являются корнем характеристического |
уравнения, кратности λ = 2 (см. замечание 1, гл. II, §3, п. 3), то |
|
y* = ( Ax + B)sin x +(Cx + D)cos x x2 . |
|
|
|
71
Глава IV
СИСТЕМЫ ДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИЙ
§1. Общие понятия
С системами дифференциальных уравнений встречаются при изучении процессов, для описания которых одной функции недостаточно. Так, решение задач динамики криволинейного движения приводит к системе трех дифференциальных уравнений, в которых неизвестными функциями являются проекции движущейся точки на оси координат, а независимой переменной – время.
Решение задач электротехники для двух электрических цепей, находящихся в электромагнитной связи, потребует решения системы двух дифференциальных уравнений. Количество подобных примеров легко можно увеличить.
Обычно приходится иметь дело не с произвольными системами дифференциальных уравнений, а с так называемыми нормальными системами.
|
|
Основные формулы |
|
|
Определения и замечания |
|||||||||||||
|
dy1 |
|
|
= f1 (x, y1 , y2 ,..., yn ) |
(4.1) |
|
Система дифференциальных уравне- |
|||||||||||
|
|
|
|
|
|
|
|
ний первого порядка, относительно |
||||||||||
|
dx |
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
искомых функций |
y1 (x), y2 (x), ..., |
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
dy2 |
|
|
= f |
2 (x, y1 , y2 ,..., yn ) |
|
|
|
|
|
||||||||
1. |
|
|
|
|
|
|
|
yn (x) , причем правые части уравне- |
||||||||||
|
|
|
|
|
||||||||||||||
dx |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
... |
|
|
|
|
|
|
|
|
|
|
|
|
|
ний не содержат производных иско- |
|||
|
|
dyn |
|
|
|
|
|
|
|
|
|
|
мых функций, называется системой |
|||||
|
|
= f |
|
(x, y |
, y |
|
,..., y |
|
) |
|
|
|||||||
|
n |
2 |
n |
|
|
в нормальной |
форме или нор- |
|||||||||||
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
1 |
|
|
|
|
|
||||||
|
dx |
|
|
|
|
|
|
|
|
|
|
|
|
|
мальной системой. |
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
y1 = ϕ1 (x,C1,C2 ,...,Cn ) |
|
(4.2) |
|
(4.2) – общее решение нормальной |
|||||||||||||
|
|
|
||||||||||||||||
|
|
|
= ϕ2 |
(x,C1,C2 ,...,Cn ) |
|
|
системы (4.1), |
где |
C1 ,C2 ,...,Cn – |
|||||||||
2. |
y2 |
|
|
произвольные постоянные. |
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
... |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
= ϕn |
(x,C1 ,C2 ,...,Cn ) |
|
|
|
|
|
|||||||||
|
yn |
|
|
|
|
|
||||||||||||
72 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Основные формулы |
|
Определения и замечания |
y1 (x0 ) = y10 |
(4.3) |
(4.3) – начальные условия, при по- |
|
|
(x0 ) = y20 |
|
мощи которых из общего решения |
y2 |
|
выделяется частное. |
|
|
|
|
|
... |
|
|
|
yn (x0 ) = yn0 |
|
|
|
ϕ1 ( x0 ,C1,C2 ,...,Cn ) = y10 |
(4.4) |
Подставляя начальные условия в об- |
|
|
( x0 ,C1,C2 ,...,Cn ) = y20 |
|
щее решение (4.2), получаем систему |
ϕ2 |
|
алгебраических уравнений для опре- |
|
|
|
|
|
... |
|
|
деления произвольных постоянных. |
ϕn (x0 ,C1 ,C2 ,...,Cn ) = yn0 |
|
|
§2. Метод исключения неизвестных
Для интегрирования системы (4.1) можно применить метод, с помощью которого данная система, содержащая n уравнений относительно n искомых функций, сводится к одному уравнению n-го порядка относительно одной неизвестной функции. Этот метод называется методом исключения неизвестных.
Покажем его применение на примере. Пусть дана система уравнений:
dx |
=3x |
+ y, |
|
dt |
|
||
|
|
|
|
dy |
= −x |
+ y. |
|
|
|
||
dt |
|
|
|
Из всех решений выделить то решение, которое удовлетворяет |
|||
условию: x(0) =1, y (0) = 0. |
|
|
|
Решение |
|
|
|
Дифференцируя первое из уравнений системы по t, |
находим: |
||
d 2 x =3 dx |
+ dy |
(4.5) |
|
dt2 |
dt |
dt |
|
|
|
|
73 |
Подставляя в это равенство выражение dydt из второго уравне-
ния системы, получим:
d 2 x |
=3 dx |
− x + y . |
(4.6) |
|
dt2 |
||||
dt |
|
|
Заменяя, наконец, функцию y ее выражением из первого уравнения системы
y = dx |
−3x , |
(4.7) |
dt |
|
|
приходим к линейному однородному уравнению второго порядка относительно одной неизвестной функции:
d 2 x |
=3 dx |
− x + dx −3x |
||
dt2 |
|
dt |
dt |
|
или |
|
|
|
|
|
d 2 x |
−4 dx + 4x = 0. |
||
|
|
|||
|
dt2 |
dt |
||
Составляем характеристическое уравнение |
||||
|
|
k2 −4k + 4 = 0, |
||
k = 2, x = e2t , x =te2t . |
||||
1,2 |
|
1 |
2 |
Общее решение:
x =C1e2t +C2te2t
или
x = e2t (C1 +C2t ) .
Дифференцируя равенство (4.9), находим:
dxdt = 2e2t (C1 +C2t) +C2e2t
(4.8)
(4.9)
74
|
или |
|
dx |
= e2t (2C1 + 2C2t +C2 ) . |
(4.10) |
dt |
|
|
Подставляя выражения (4.9) и (4.10) в равенство (4.7), получим:
y = e2t (2C1 + 2C2t +C2 ) −3e2t (C1 +C2t )
или
y = −C1e2t +C2e2t (1−t) .
Общее решение системы:
|
|
2t |
|
|
2t |
|
|
x =C1e |
+C2te |
, |
|
||||
|
|
|
|||||
y = −C e2t +C |
e2t |
(1−t ). |
|||||
|
1 |
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
Решим задачу Коши. Пусть x(0) =1 . Тогда
C1e0 +C2 0 e0 =1,
следовательно,
C1 =1 .
Пусть y (0) = 0. Тогда:
−C1e0 +C2e0 1 = 0,
−C1 +C2 = 0,
C2 =1 .
Искомое решение системы
|
2t |
(1 |
+t ) |
|
x = e |
|
|||
|
. |
|||
|
|
|
|
|
y = −te2t |
|
|
||
|
|
|
|
|
(4.11)
75
ПРИЛОЖЕНИЯ
Задача 1. Найти общий интеграл дифференциального уравнения.
1.01. ( x +1)3 dy −( y −2)2 dx = 0
1.02. sec2 x sec ydx +ctg x sin ydy = 0 1.03. x 1+ y2 + yy′ 2 + x2 = 0
1.04. 3x2 ydy −2xy2dx = 6xdx −6 ydy
1.05. ( xy + x ) y′− y = 0 1.06. ye2 xdx +(e2 x +5)dy = 0
1.07. |
2x+y −3x−2 y y′= 0 |
1.08. |
x 4 + y2 dx + y 3 + x2 dy = 0 |
1.09.2(1+ex ) yy′= ex
1.10.4xdx −3ydy =3x2 ydy −2xy2 dx
1.11. yy′ |
1− x2 |
+1 = 0 |
|
1− y2 |
|
1.12.3ex sin ydx =(ex −1)sec ydy
1.13.(1+ x2 ) y3dx −( y2 −1)x3dy = 0
1.14.5 + y2 dx − ydy = x2 ydy
1.15. 2x + 2xy2 + 2 − x2 y′= 0
1.16.y ln y + xy′= 0
1.17.y′sin2 x ln y + y = 0
1.18.( xy − x)2 dy + y(1− x)dx = 0
76
1.19.y = 2(1+ x2 y′)+ xy′
1.20.( y −1) xdy − y2 ln xdx = 0
1.21.(4 +ex ) y′= 2 yex
1.22.x + xy2 + 1− x2 y′= 0
1.23.2 ydy −6xdx =3xy2dx −2 yx2 dy
1.24.y (1+ln y) + xy′= 0
1.25.y3 ln2 xdx −( y −1) xdy = 0
1.26.2x2 yy′+ y2 = 2
1.27.5 + y2 dx + 4(x2 y + y)dy = 0
1.28.y′cos x −( y +1)sin x = 0
1.29. 2(x2 y + y)dy + 3 + y2 dx = 0
1.30.4 + y2 dx − ydy = x2 ydy
Задача 2. Найти общий интеграл дифференциального уравнения.
2.01. (x2 + y2 )dx −2xydy = 0 2.02. y − xy′= y ln xy
2.03. xdy − ydx = x tg xy dx
2.04.3y′= y2 +8 y + 4
x2 x
2.05. xy′ =3 2x2 + y2 + y
2.06. y − xy′= x + yy′
77
2.07. y = x(y′− x ey )
2.08. ( y2 −3x2 )dy + 2xydx = 0
2.09.xy′+ y ln xy −1 = 0
2.10.ydx +(2 xy − x)dy = 0
2.11.y′(2x2 −2xy) = x2 + 2xy − y2
2.12.4 y′= y2 +10 y +5
x2 x
2.13.x3 y′= xy2 − x2 y
2.14.xy′= y + xsin xy
2.15.xy′− y(ln y −ln x) = 0
2.16. y′cos |
y |
− |
y |
|
cos |
y |
|
+1 = 0 |
|||||
x |
|
x |
|||||||||||
|
|
|
|
x |
|
|
|
||||||
|
y = xy′− xe |
y |
|
|
|
|
|
|
|||||
2.17. |
x |
|
|
|
|
|
|
||||||
2.18. |
y (2 |
x − |
y )dx + xdy = 0 |
||||||||||
2.19. y − xy′= 2(x+ yy′) |
|||||||||||||
2.20. |
y′= |
x2 + xy −3y2 |
|
||||||||||
|
|
|
|
||||||||||
|
|
|
x2 −4xy |
|
|||||||||
2.21. |
y − xy′= xsec |
y |
|
|
|||||||||
x |
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
||||
2.22. |
xy′= 3x3 +10 yx2 |
|
|||||||||||
|
|
|
2 y2 +5x2 |
|
2.23. ( x − y) y′= x + y
78
2.24.xy′− y =(x + y)ln x +x y
2.25.dydx − xy (1+ln y −ln x) = 0
2.26. xy′− y =3 x2 + y2
2.27.x(2 y2 + 4x2 ) y′=3y3 +8yx2
2.28.y′(2x − y) = x + 2 y
2.29.2 y′ = y2 +8 y +8
x2 x
2.30.xy′− y = x tg xy
Задача 3. Найти решение задачи Коши.
3.01. |
y′+ |
2xy |
|
= |
2x2 |
|
|||||
1+ x2 |
|
|
|
|
|||||||
|
|
|
|
1+ x2 |
|||||||
3.02. |
y′sin x − y cos x =1 |
||||||||||
3.03. |
xy′+ y + xe−x2 |
= 0 |
|
||||||||
3.04. |
xy′− y − x3 = 0 |
|
|
|
|||||||
3.05. |
y′− |
y |
|
= tg |
|
x |
|
||||
sin x |
|
|
|
|
|||||||
|
|
2 |
|
|
|||||||
3.06. |
x( y′− y) = ex |
|
|
|
|||||||
3.07. |
y′+ |
2xy |
|
= 2 |
|
|
|
|
|||
1 − x2 |
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
||
3.08. |
xy′+( x +1) y =3x2e−x |
||||||||||
3.09. |
y′+ 2xy = xe−x2 sin x |
y(0) = 23
|
π |
= |
1 |
||
y |
4 |
|
|
||
2 |
|||||
|
|
|
y(1) = 12 е−1 y(2) = 4
y π = 02
y(1) = 0
y(0) = 2
y(1) = 0 y(0) =1
79
3.10.y′− x2+y1 = ex (x +1)2
3.11.y′− y ctgx = 2xsin x
3.12.y′+ y tgx = cos2 x
3.13.xy′+ y = ln x +1
3.14.y′+ 2 yx−1 = x−2
3.15.(x2 −1) y′− xy = x3 − x
3.16.(xy′−1)ln x = 2 y
3.17.y = x( y′− xcos x)
3.18.( x + 2) y′− y =(x2 + 2x)(x + 2)
3.19.y′− 2xx−2 5 y =5
3.20.xy′+ y = xsin x
3.21.( x +1) y′+ y = x3 + x2
3.22.y′+ 2yx = x2
3.23.y′− y cos x = sin 2x
3.24.y′− x2+y1 =(x +1)3
3.25.y′− 2xy = −lnxx
3.26.y′− y ctgx = sin1 x
80
y(0) =1
|
π |
= 0 |
||
y |
2 |
|
||
|
|
|
|
|
|
π |
= |
1 |
|
y |
4 |
|
2 |
|
|
|
|
||
y(1) = 0 |
|
|||
y(2) =1 |
|
|||
y( |
2 ) =1 |
y(e) = 0
y π = 02
y(−1) = 32 y(2) = 4
y(π) = π1
y(0) = 0 y(1) =1
y(0) = −1 y(0) = 12 y(1) =1
y π =12