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

Учебное пособие 800446

.pdf
Скачиваний:
7
Добавлен:
01.05.2022
Размер:
2.77 Mб
Скачать

~

~

~

 

 

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

 

 

0 ,0

Если в точке 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