Учебное пособие 800446
.pdf~ |
~ |
~ |
|
|
1 |
|
|
|
1 |
|
x |
|
|
|
|
|
|
|
|
|
||||
y |
y1 |
y2 |
|
|
|
cosx |
|
|
e |
|
|
, |
|
а |
общее |
решение |
этого |
уравнения |
||||||
|
2 |
4 |
|
|
|
|||||||||||||||||||
|
~ |
|
1 |
|
|
|
|
1 |
|
|
|
x |
|
|
x |
|
|
x |
|
|
|
|
||
y |
y Y |
|
|
cosx |
|
|
e |
|
|
|
C1e |
|
C2xe |
|
. |
|
|
|
||||||
2 |
|
4 |
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Задачи для самостоятельного решения
Решить уравнения:
1.y 5y 4y 0. Ответ: y C1ex C2e4x.
2.y 6y 9y 0. Ответ: y e3x(C1 C2x).
3.y 8y 25y 0. Ответ: y e 4x (C1 cos3x C2 sin x).
4.y 3y 2y 0. Ответ: y C1ex C2e2x.
5.y 4y 4y 0. Ответ: y e2x C1 C2x .
6.y 2y 2y 0. Ответ: y ex(C1 cos x C2 sin x).
7.y 4y 3y 0. Ответ: y C1ex C2e3x.
8.y 4y 0. Ответ: y C1e2x C2e 2x.
9.y 4y 0. Ответ: y C1 C2e 4x.
10.y 3y 2y 0. Ответ: y C1e 2x C2e x.
11.y 2ay a2y 0. Ответ: y (C1x C2)e ax.
12.y 2y 5y 0. Ответ: y e x(C1 cos2x C 2 sin2x)
13.y y 0. Ответ: y C1ex C2e x.
14.y y 0. Ответ: y C1cosx C2 sin x.
15. |
y 2y y e |
x |
. |
Ответ: |
y (C1x C2)e |
x |
|
1 |
e |
x |
. |
||||
|
|
|
4 |
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
16. |
y 3y 2y ex. |
Ответ: |
y C1e2x (C2 |
x)ex. |
|
|
|||||||||
17. |
y y 2y 6x2. |
Ответ: |
y C ex |
C e 2x |
3(x2 |
x 1.5). |
|||||||||
|
|
|
|
|
1 |
2 |
|
|
|
3 |
|
|
|
|
|
18. |
y 3y 9x. |
|
|
Ответ: |
y C1 C2e |
3x |
|
|
x |
2 |
x. |
||||
|
|
|
|
2 |
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
19.y 2y xe x. Ответ: y C1ex2 C2e x2 (x 2)e x.
20.y 4y 8x3. Ответ: y C1e2x C2e 2x 2x3 3x.
21. |
y 5y 6y 13sin3x. Ответ: |
y C1e2x C2e3x |
1 |
5cos3x sin3x . |
|||||||||
|
|||||||||||||
|
|
|
|
|
|
3x |
6 |
|
3x |
|
|||
22. |
y y 2,5y 25cos2x. Ответ: |
y e x |
/2 |
|
|
|
|
|
|
||||
C |
cos |
|
|
C sin |
|
|
6cos2x 8sin2x. |
||||||
|
|
||||||||||||
|
|
|
|
1 |
2 |
|
2 |
|
2 |
|
|||
|
|
110 |
|
|
|
|
|
|
|
|
|
|
|
23. |
y 4y 3sin 2x. |
|
Ответ: y C1 cos2x C2 sin2x |
3 |
xcos2x. |
|
|
|
||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
4 |
|
|
|
|
|
|
|
|
|
|
|
|||||
24. |
y 4y sin 2x. |
|
|
Ответ: y C1cos2x C2 sin2x |
1 |
xcosx2x. |
|
|
||||||||||||||||||||||||
|
|
|
|
|
|
|||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
4 |
|
|
|
1 |
|
|
|
|
|
||||
25. |
y y xcosx. |
|
|
Ответ: y C1cosx C2 sinx |
|
xcosx |
x |
2 |
sinx. |
|||||||||||||||||||||||
|
|
4 |
4 |
|
||||||||||||||||||||||||||||
26. y 2y 3y e x cos x. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
Ответ: y ex(C 1cos |
|
|
x C2 sin |
|
|
x) |
e x |
|
(5cosx 4sin x). |
|||||||||||||||||||
|
|
|
|
|
2 |
|
2 |
|||||||||||||||||||||||||
|
|
|
|
|
|
|
||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
41 |
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
5x |
|
|
|
2x |
|
|
|
|
|
2 |
|
|
|
|
|
|
9 |
|
|||||||
27. |
|
|
|
Ответ: y C1e |
C2e |
|
|
|
|
|
|
sin x |
|
|
cosx. |
|||||||||||||||||
y 3y 10y sinx 3cosx. |
|
|
|
|
|
|
|
13 |
|
|
||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
39 |
|
|||||||||
|
3x |
|
2 |
|
|
x |
|
|
|
2x |
|
|
|
3x |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
28. y 3y 2y e |
(x |
|
x). |
Ответ: y C1e |
|
C2e |
|
|
|
|
e |
|
|
|
|
|
(C1 cosx C2 sin x). |
|||||||||||||||
|
|
|
2 |
|
|
|
|
|
||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
29.y 2y 2y ex(2cosx 4xsinx). Ответ: y ex(C1cosx C2 sinx) x2ex cosx.
30.y 4y ex ( 4x 4)cosx (2x 6)sinx .
31.y y x 2ex.
32.y 2y y 3ex x 1.
33.y 6y 8y ex e2x.
34.y 4y x e 4x.
Ответ: y C1e2x C2e 2x ex(xcosx sin x).
Ответ: y C |
|
cosx C |
2 |
|
sinx x ex. |
|
|
|
||||||||||||||||||
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
3 |
|
|
x |
|
2 |
|
|
|
|
|
x |
|
|
|
|
|
|
|
|
|
|
|
|||
Ответ: |
y |
|
|
e |
|
x |
|
x 3 e |
|
(C1x C2). |
|
|
||||||||||||||
|
|
|
|
|
|
|
||||||||||||||||||||
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
1 |
|
x |
|
|
1 |
2x |
|
|
|
|
2x |
|
|
|
|
4x |
|
|
|
||||||
Ответ: |
y |
|
|
e |
|
|
|
|
xe |
|
C1e |
|
|
C2e |
|
. |
|
|
||||||||
|
|
|
2 |
|
|
|
|
|
|
|||||||||||||||||
|
3 |
|
|
|
|
|
|
|
|
|
x2 |
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
x |
|
|
x |
|
|
|
|
|
|||||||
Ответ: |
y e 4x C |
|
|
|
|
|
|
|
|
|
|
C |
2 |
. |
||||||||||||
4 |
|
8 |
|
16 |
||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
y 4y e2x 3e 2x. |
|
|
|
x |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
3 |
|
|
|
|
|
|||||||||||||
35. |
Ответ: |
y e2x C |
|
|
|
e 2x C |
|
|
|
|
|
|
x . |
|
||||||||||||||||||||||||
4 |
|
4 |
|
|
||||||||||||||||||||||||||||||||||
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|||||||||||||||
|
y 9y 4sin3x x. |
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
x |
|
|||||||||||||
36. |
Ответ: y C sin3x C |
|
|
|
|
|
|
x cos3x |
|
. |
|
|
||||||||||||||||||||||||||
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||
|
|
|
|
|
1 |
|
|
|
2 |
|
3 |
|
|
|
|
|
|
|
|
|
|
|
|
|
9 |
|
|
|
||||||||||
|
y 3y x |
3 |
2. |
Ответ: y C1 C2e |
3x |
|
|
x4 |
|
|
|
|
|
x |
3 |
|
|
|
|
x2 |
20 |
x. |
||||||||||||||||
37. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
12 |
9 |
|
|
|
9 |
|
|
|
|
|||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
x |
|
|
|
|
|
|
27 |
|
|||||||||||||||
38. |
y 3y 1. |
|
Ответ: |
y C1 C2e |
3x |
|
|
|
. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
3 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
x2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
x |
|
|
|||||||||||||
39. |
y y xcosx. |
Ответ: |
y C |
|
|
|
|
|
sinx C |
|
|
|
|
|
|
|
cosx. |
|||||||||||||||||||||
4 |
|
|
|
|
|
4 |
|
|||||||||||||||||||||||||||||||
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
40. |
y y xe |
x |
. |
Ответ: |
y C1 C2e |
x |
|
e |
x x |
|
|
|
3 |
|
|
|
|
|
|
|
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
. |
|
||||||||||||||||||||||
|
|
|
|
|
|
2 |
|
4 |
|
|
||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
111
5. СИСТЕМЫ ДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИЙ
5.1.Нормальные системы дифференциальных уравнений
1.Понятие о нормальной системе. Линейные системы. Автономные системы
Нормальная система обыкновенных дифференциальных уравнений содержит n неизвестных функций одного переменного и состоит из n уравнений первого порядка, каждое из которых содержит производную только одной из функций и разрешено относительно этой производной:
|
dy1 |
f (x,y ,y |
2 |
, ,y |
n |
), |
||
|
||||||||
dx |
1 |
1 |
|
|
||||
|
|
|
|
|
|
|
||
dy2 |
f2(x,y1,y2, ,yn), |
|||||||
|
|
|||||||
|
||||||||
dx |
|
|
|
|
|
(5.1) |
dyn fn(x,y1,y2, ,yn);dx
где y1, y2,..., yn - искомые функции от независимой переменной x, f1,..., fn -
заданные функции, зависящие от (n 1) переменных. Число уравнений, входящих в систему (5.1), называется порядком этой системы. То есть, согласно этому определению, система (5.1) есть система n- го порядка.
К нормальному виду может быть приведена система уравнений высших порядков, каждое из которых разрешено относительно старшей производной. Для этого должны быть введены дополнительные функции, как в рассмотренном ниже примере.
|
Уравнение движения материальной точки массы m под действием силы |
|
|
dr |
|
Fимеет вид:
dt,rt,
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
d2 |
r |
|
|
|
|
|
|
dr |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
m |
|
|
|
F t,r, |
|
, |
|
(5.2) |
||||||||
|
|
|
|
|
|
|
|
|
|
|
dt2 |
|
|
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
dt |
|
dr |
|
|
|
|
|||||
где |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
- время, |
|
|
v |
- скорость |
|
r |
x(t)i |
y(t)j z(t)k радиус - вектор точки, t |
|
|
||||||||||||||||||||||||||
|
|
|||||||||||||||||||||||||||||
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
dt |
|||||
точки, |
d |
r |
|
|
a |
|
- ускорение. Записанное уравнение выражает второй закон |
|||||||||||||||||||||||
|
|
|
|
|||||||||||||||||||||||||||
|
|
|
|
|
|
|||||||||||||||||||||||||
|
|
|
|
|
||||||||||||||||||||||||||
|
|
|
dt2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
Ньютона. |
|
|
|
|
Спроектируем |
|
|
действующую |
на точку |
|
|
силу на оси декартово |
||||||||||||||||||
системы координат: |
|
|
112 |
|
|
|
|
|
|
|
|
|
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
dr |
|
|
|
|
|
|
|
dr |
|
|
|
|
|
|
|
|
dr |
|
|
|
|
|
|
|
|
dr |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
F t,r, |
|
|
|
|
f1 |
t,r, |
|
|
i |
f2 |
t,r, |
|
|
j |
f3 |
t,r, |
|
|
k. |
||||||||||||||||
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||
|
|
|
dt |
|
|
|
dt |
|
|
dt |
|
|
dt |
Подставив разложения r и F в уравнение (5.2), приравняв проекции на координатные оси. Получим систему трех уравнений:
|
d2x |
|
|
|
|
|
m |
|
|
|
|
f1(t,x, y,z,x,y,z), |
|
dt2 |
|
|||||
|
|
|
|
|
||
|
d2y |
|
|
|
|
|
|
|
|
|
|
||
m |
|
|
|
|
f2(t,x, y,z,x, y,z), |
(5.3) |
dt |
2 |
|
||||
|
|
|
|
|
|
|
|
d2z |
|
|
|
|
|
|
|
|
f3(t,x,y,z,x, y,z), |
|
||
m |
|
2 |
|
|
||
|
dt |
|
|
|
|
|
|
|
|
|
|
|
где точка над буквой обозначает производную по t. Каждое из уравнений системы (5.3) является дифференциальным уравнением второго порядка. Для
приведения системы к нормальному виду введем функции u x, |
v y , |
|||||||||||
w z - проекции |
скорости |
dr |
на координатные |
оси. Получим |
систему |
|||||||
dt |
||||||||||||
уравнений вида |
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|||||
|
|
|
|
x u, |
y v, |
z w, |
|
|
||||
|
|
|
|
|
f1(t,x, y,z,u,v,w), |
|
|
|||||
|
|
|
|
mu |
|
|
||||||
|
|
|
|
|
f2(t,x, y,z,u,v,w), |
|
|
|||||
|
|
|
|
mv |
|
|
||||||
|
|
|
|
|
f3(t,x, y,z,u,v,w). |
|
|
|||||
|
|
|
|
mw |
|
|
||||||
|
|
При этом |
обычно задаются |
начальное положение |
точки |
|||||||
|
r |
(t0) |
r0 x0, y0,z0 и начальная скорость |
dr |
(t0) |
r0 |
u0,v0,w0 . |
|
||||
|
|
|
|
|||||||||
|
|
|
|
|
|
|
dt |
|
|
Подобным образом и одному дифференциальному уравнению высшего порядка можно сопоставить систему дифференциальных уравнений первого порядка.
Система уравнений (5.1) называется линейной системой дифференциальных уравнений, если правые части зависят линейно от некоторых функций y1, y2, , yn :
113
|
dy1 |
|
a |
|
(x)y a (x)y |
2 |
a |
|
(x)y |
n |
f (x), |
||||||||||||||||
|
|
|
|
||||||||||||||||||||||||
dx |
11 |
|
|
1 |
12 |
|
|
|
|
1n |
|
|
|
1 |
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
dy2 |
|
a21(x)y1 |
a22 |
(x)y2 a2n(x)yn f2(x), |
|||||||||||||||||||||||
|
|
|
|||||||||||||||||||||||||
|
|||||||||||||||||||||||||||
dx |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(5.4) |
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
dyn |
|
a |
n1 |
(x)y a |
n2 |
(x)y |
2 |
a |
nn |
(x)y |
n |
f |
(x), |
|||||||||||||
|
|
||||||||||||||||||||||||||
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
n |
|
|||||||||||
dx |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
где aij (x) и fi (x) |
|
(i, j 1,2, ,n) |
- заданные функции от x. |
||||||||||||||||||||||||
Если правые части системы (5.1) не зависят явно от независимой |
|||||||||||||||||||||||||||
переменной x, т.е. если система имеет вид |
|
|
|
|
|
|
|
|
|||||||||||||||||||
|
|
|
|
|
|
|
|
dy1 |
|
|
f |
(y , y |
, , y |
), |
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
|
|
|
|
|
|
|
dx |
|
|
|
1 |
|
1 |
|
|
2 |
|
n |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
dy2 |
|
|
f2 |
(y1, y2, , yn), |
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
|
|
|
|
|
|
|
dx |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(5.5) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
|
|
|
|
|
|
|
|
dyn |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
f |
|
(y , y |
, , y |
|
), |
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
dx |
|
|
|
n |
1 |
|
|
2 |
n |
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
то она называется автономной или стационарной системой. Если в линейной системе (5.4) f1(x) f2(x) fn(x) 0, то имеем линейную автономную систему.
2. Решение системы. Геометрическое и механическое истолкования нормальной системы.
Совокупность n функций y 1(x), , y n (x), определенных и непрерывно дифференцируемых в интервале (a,b), называется решением системы (5.1) в этом интервале, если она обращает все уравнения системы в тождества:
1(x) f1(x, 1(x), 2(x), , n(x)),
|
(x) f |
|
(x, (x), |
|
(x), , |
|
(x)), |
|||
|
2 |
2 |
n |
|||||||
2 |
|
|
|
|
1 |
|
|
|||
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
(x) f |
|
(x, (x), |
|
(x), , |
|
(x)), |
|||
|
n |
2 |
n |
|||||||
n |
|
|
|
|
1 |
|
|
|||
справедливые при всех значениях x |
из интервала (a,b). |
|||||||||
Пример. Система уравнений |
|
|
|
|
|
|||||
|
|
dy |
7y 6z, |
|
|
|||||
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|||||
|
|
dx |
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|||
|
|
dz |
6y 7z, |
|
|
|||||
|
|
|
|
|
|
|
||||
|
|
dx |
|
|
|
|
|
114
имеет решение y1 ex, z1 ex. В этом можно убедиться, подставив эти функции в уравнения системы:
ex 7ex 6ex, ex 6ex 7ex.
Мы получили тождества, справедливые при x ( ; ) . Эта же система имеет и другие решения. Например, решением, является другая совокупность
функций: |
y2 e13x, |
z2 e13x при |
x ( ; ) . Решением этой системы будет |
||
являться |
и двухпараметрическое |
семейство |
функций: |
y Aex Be13x , |
|
z Aex Be13x , A |
и B – произвольные |
постоянные. |
Действительно, |
||
подставив эти функции в уравнения системы, получим: |
|
||||
|
Aex 13Be13x 7(Aex Be13x) 6( Aex Be13x), |
|
|||
|
Aex 13Be13x 6(Aex Be13x) 7( Aex Be13x). |
|
|||
Эти |
уравнения |
удовлетворяются тождественно относительно x при |
любых числовых значениях параметров A и B.
Рассмотренный пример системы двух уравнений наводит на мысль, что всякая система двух уравнений имеет решение, содержащее два произвольных параметра.
Перейдем к геометрическому истолкованию нормальной системы. |
|
|||||||
Изучая |
дифференциальные |
уравнения |
1-го |
порядка |
dy |
f (x, y) |
и |
|
|
||||||||
|
|
|
|
|
|
dx |
|
|
рассматривая |
x, y |
как прямоугольные координаты, |
мы установили, что такие |
|||||
уравнения задают |
на плоскости |
XOY некоторое |
поле направлений и |
что |
||||
направление |
касательной в каждой точке |
интегральной кривой y (x) |
совпадает с направлением поля в этой точке. Аналогичную геометрическую интерпретацию можно дать и нормальной системе n дифференциальных уравнений (5.1).
Будем рассматривать (x,y1, , yn) как координаты точки в (n 1)- мерном пространстве x, y1, , yn. Тогда решению y 1(x), y 2(x), , y n(x)
соответствует |
кривая в введенном (n 1)- мерном пространстве. Эта кривая |
|||
называется интегральной кривой системы (5.1). |
f2(x, y1, y2, , yn), |
|
||
Пусть |
функции |
f1(x,y1,y2, ,yn), |
, |
fn(x, y1, y2, , yn) в системе (5.1) определены и конечны в некоторой области
D изменения |
переменных x, y1, , yn. В |
каждой точке |
M(x, y1, , yn) D |
||
проведем |
отрезок |
прямой |
с направляющим вектором |
||
s |
1, f1(x, y1, y2, , yn), , fn(x, y1, y2, , yn) . Тогда получим |
поле направлений |
|||
в области D |
(n 1)- мерного пространства. Интегральная кривая в каждой |
||||
своей точке касается вектора поля |
|
|
|||
|
|
(dx:dy1 :dy2 : :dyn 1: f1 : f2 : fn). |
|
||
|
|
|
115 |
|
|
Если в точке M0(x0, y1(x0), , yn(x0)) все правые части системы (5.1)
или некоторые из них обращаются в неопределенные выражения вида
то в этой точке поле направлений не определено и через такую точку не проходит ни одна интегральная кривая.
Рассмотрим механическое истолкование нормальной системы.
Если дифференциальные уравнения описывают движение механической системы, то независимую переменную обозначают t, искомые функции через
x1, ,xn, |
а |
производные |
|
этих |
функций по времени |
t |
обозначают xk |
||||||||
(k 1,2, ,n). |
В этих |
обозначениях |
нормальная система |
уравнений |
(5.1) |
||||||||||
записывается так: |
x1 X1(t,x1,x2, ,xn), |
|
|
|
|||||||||||
|
|
|
|
|
|
||||||||||
|
|
|
|
|
X |
|
(t,x ,x |
|
, ,x |
|
), |
|
|
|
|
|
|
|
x |
2 |
2 |
2 |
n |
|
|
|
|||||
|
|
|
|
|
1 |
|
|
|
|
|
(5.6) |
||||
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
x |
n |
X |
n |
(t,x ,x |
2 |
, ,x |
n |
). |
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
||||
В |
евклидовом |
пространстве |
|
с |
прямоугольными |
координатами |
|||||||||
x1,x2, ,xn |
решение системы x1 |
1(t), ,xn n(t) |
определяет |
закон |
движения или движение по некоторой траектории с течением времени t. Если
ввести n |
- мерный |
вектор |
rn x1(t),x2(t), ,xn(t) , то его производная по |
|||||||
времени r |
(t) x (t), ,x |
n |
(t) интерпретируется как скорость движения точки в |
|||||||
n |
1 |
|
|
xk (t) (k 1,2, ,n)- |
|
|||||
n - мерном пространстве, |
проекции вектора скорости на |
|||||||||
координатные оси. |
|
|
|
|
|
|
|
|
|
|
При такой интерпретации, естественной |
в физических и механических |
|||||||||
задачах, система (5.6) или в векторной форме |
|
|||||||||
|
|
|
|
|
|
rn |
|
(t, |
|
(5.7) |
|
|
|
|
|
|
F |
rn ) , |
где |
F |
(t, |
rn) X1(t,x1,...,xn), ,Xn(t,x1,...,xn) , обычно называется динамической |
|||||||
системой, пространство с координатами x1,x2,...,xn называется |
фазовым |
|||||||||
пространством, кривая |
r |
|
r |
(t) называется фазовой траекторией. |
В случае |
|||||
n 2 фазовое пространство называется фазовой плоскостью. |
|
|
|
|||||||
|
|
Если функции X k (k 1,2,...,n) в системе (5.6) или вектор-функция |
|
в |
||||||
|
|
F |
векторном уравнении (5.7) не зависят от времени t, то система (5.6), или (5.7),
называется автономной динамической системой.
Взаимосвязь между траекторией и движением состоит в том, что
траектория |
|
есть проекция движения (расположенного |
в пространстве |
|||
t,x1,x2,...,xn ) |
в пространство x1,x2,...,xn . |
|
|
|
||
Таким |
|
образом, |
решение |
системы |
(5.6) |
- |
x1 1(t),x2 |
2(t), ,xn n(t), |
понимается |
как |
параметрические |
||
|
|
|
116 |
|
|
|
уравнения траектории движения (в пространстве x1,x2,...,xn ). Эти уравнению определяют траекторию как однопараметрическое геометрическое место точек; они показывают, как происходит движение точки по траектории с течением времени. В соответствии с этим решение системы x1 1(t),x2 2(t), ,xn n(t) обычно в механике называют движением.
Основной задачей интегрирования системы (5.6) или (5.7) является нахождение всех движений, определяемых этой системой и изучение свойств движений.
Пример. Пусть задана система уравнений x y, y x. Можно убедиться непосредственной подстановкой в уравнения, что семейство
функций x C1cos(t C2), |
y C1sin(t C2) |
является решением заданной |
системы уравнений. На фазовой плоскости |
XOY фазовыми траекториями |
|
является семейство окружностей x2 y2 C2 |
с центром в начале координат. |
|
|
1 |
|
Придавая параметру t числовые значения 0 t |
2 , будем иметь перемещение |
точки по определенной траектории в направлении, показанном на рис. 11.
y
s1 s2
x
Рис. 11. Фазовые траектории системы
Отметим, что уравнение траектории на фазовой плоскости не зависит от C2 . Поэтому все движения при фиксированном C1 происходят по одной и той
же траектории. При C1 0 (x2 y2 0) фазовая траектория состоит из одной точки (точки О), называемой в этом случае точкой покоя системы.
5.2.Задача Коши. Достаточные условия существования
иединственности решения задачи Коши. Общее, частное
иособое решения
1.Задача Коши для системы вида (5.1) ставится так: среди всех решений
системы найти такое решение y1 1(x), , yn n(x), в котором функцииk (x) (k 1,2, ,n) принимают заданные числовые значения y1(0), y2(0), , yn(0)
при заданном числовом значении x x0 независимой переменной:
117
1(x0) y1(0), , n(x0) yn(0).
Совокупность чисел x0, y1(0), , yn(0) называется начальными данными решения y1 1(x), , yn n(x); число x0 называют начальным значением
независимой переменной; числа y1(0), , yn(0) называются начальными
значениями искомых функций (начальные значения решения). Геометрический смысл задачи Коши заключается в том, что среди всех
интегральных кривых системы разыскивается та, которая проходит через
заданную точку (x0, y1(0), , yn(0)). |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
Остановимся |
на |
механическом |
истолковании |
задачи |
|
Коши |
для |
|||||||||||||||||||
нормальной системы (5.6): |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
xk |
Xk (t,x1, ,xn) |
(k 1,2, ,n). |
|
|
|
|
|
|
|
|
||||||||||
|
|
Задача Коши состоит в том, чтобы из всех движений, определяемых этой |
||||||||||||||||||||||||||
динамической |
системой, |
найти |
|
такое |
движение |
x1 1(t), |
|
x2 |
2(t),..., |
|||||||||||||||||||
x |
n |
|
n |
(t), |
в |
котором |
при |
t |
t |
0 |
x |
x(0),x |
2 |
x |
(0) |
, ,x |
n |
x(0), |
т.е. |
|||||||||
|
|
|
|
|
|
|
|
|
|
|
1 |
1 |
|
|
|
2 |
|
|
|
|
n |
|
||||||
движущаяся точка в заданный момент времени t |
t0 |
находится в заданной |
||||||||||||||||||||||||||
точке фазового пространства M0(x1(0),x2(0), ,xn(0)). |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
Совокупность чисел t0,x1(0), ,xn(0) |
называется |
начальными |
данными |
|||||||||||||||||||||||
движения; |
условия |
x |
x |
(0), ,x |
n |
x(0) |
называются начальными условиями |
|||||||||||||||||||||
|
|
|
|
|
|
|
|
1 |
1 |
|
n |
|
|
|
|
|
|
|
|
|
(x(0) |
|
(0)) - |
|||||
этого движения; t |
0 |
- начальным моментом времени; |
точка |
M |
0 |
, ,x |
||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
n |
||
называется начальной точкой. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
Пример. |
Найти решение автономной динамической системы |
|
|
|
||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
x y, |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
y x, |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
удовлетворяющее начальным условиям: t0 0, x0 2, y0 |
1. |
|
|
|
|
|
|
|
|
|||||||||||||||||||
|
|
Решение. |
|
До сих |
пор мы не излагали |
|
систематические методы |
получения решений систем дифференциальных уравнений. Представленная в примере система обладает особенностью, позволяющей найти ее решение без привлечения специальных методов. Сложив уравнения системы, а затем, вычитая из 1-го уравнения второе, получим следующую систему уравнений:
d(x y) |
||
|
|
(x y), |
|
||
|
dt |
|
|
d(x y) |
|
|
|
|
|
|
(x y). |
|
dt |
Видим, что новая система представляет собой систему независимых
118
уравнений относительно функций |
(x y) |
и |
(x y) . |
Интегрируя уравнения, |
||||||||||||||||||||||||
получим: |
|
|
d(x y) |
|
|
|
|
~ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
t |
C1, |
|
|
|
|
|
|
|
~ |
|
|
|
|
|
t |
|
|
||||
|
|
|
|
(x y) |
|
|
|
|
|
y) |
|
, |
|
|
, |
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ln(x |
|
t C1 |
|
x y C1e |
|
|
||||||||
|
|
|
|
|
d(x y) |
|
|
|
|
~ |
|
|
|
|
|
|
|
~ |
|
|
|
t |
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ln(x y) t C |
|
, |
|
|
|
|
|
|||||||
|
|
|
|
|
(x y) |
|
t C2, |
|
|
|
|
|
|
|
|
2 |
|
x y C2e . |
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
Из |
начальных |
|
|
условий, |
полагая |
t 0, |
|
x(0) 2, |
y(0) 1, |
получим |
||||||||||||||||
C1 3, C2 1. Подставив найденные постоянные в выражения для |
(x y) и |
|||||||||||||||||||||||||||
(x y), |
|
сложив, |
|
а |
|
затем |
|
вычтя |
|
их |
друг |
из |
друга, |
получим |
||||||||||||||
x |
1 |
(3e |
t |
e |
t |
), y |
1 |
(3e |
t |
e |
t |
). |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
2 |
|
|
2 |
|
|
|
|
|
|
|
|
x2 y2 3. Это гипербола на |
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
Уравнение фазовой траектории имеет вид: |
||||||||||||||||||||||||||
плоскости XOY . |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
2. |
|
Достаточные условия |
существования |
|
|
и единствен-ности решения |
задачи Коши. В этом разделе мы приведем без доказательства и в упрощенной формулировке основную теорему существования и единственности решений для системы (5.1).
Теорема (Пикара Э., 1856-1941, французский математик, почетный член АН СССР, 1924).
Пусть дана система уравнений (5.1)
dyk |
f |
k |
(x,y , , y |
n |
) (k 1,2, ,n) |
|
|||||
dx |
1 |
|
|||
|
|
|
|
и поставлены начальные |
условия |
y y(0) |
,y y(0), ,y |
y(0) |
при |
x x |
0 |
. |
|||||||||||||||
|
|
|
|
|
|
|
|
1 |
1 |
2 |
|
|
|
2 |
|
|
n |
n |
|
|
|
||
Пусть функции fk (x, y1, , yn) |
(k 1,2, ,n) |
|
определены |
в замкнутой |
|||||||||||||||||||
ограниченной области D (n 1)- мерного пространства: |
|
|
|
|
|
||||||||||||||||||
|
x x |
0 |
|
|
a, |
|
y |
y(0) |
|
b, , |
|
y |
|
y |
(0) |
|
b, |
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
1 |
1 |
|
|
|
|
|
n |
|
n |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
a и b - заданные положительные числа с точкой M0(x0, y1(0), ,yn(0)) внутри
D и удовлетворяют в D условиям:
функции |
fk (x, y1, , yn) |
(k 1,2, ,n) |
непрерывны по всем аргументам и |
|||||
потому ограничены |
|
|
|
|
|
|||
|
|
|
fk (x, y1, , yn) |
|
M |
(k 1,2, ,n), |
||
|
|
|
|
|||||
где M - положительное число, |
M (x, y1, , yn ) D ; |
|||||||
функции |
fk (x, y1, , yn) |
(k 1,2, ,n) имеют ограниченные частные |
||||||
производные |
|
|
|
|
|
|
|
|
|
|
fk (x, y1, , yn) |
K |
(k,l 1,2, ,n), |
||||
|
|
|
||||||
|
|
|
yl |
|
|
|
|
|
|
|
|
|
|
119 |
|