Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Моделирование процессов средствами дифференциальных уравнений. Учебное пособие.pdf
Скачиваний:
0
Добавлен:
12.08.2026
Размер:
808 Кб
Скачать

Решая систему, получаем:

 

 

 

6

 

 

 

 

C1

=

 

 

 

 

 

 

37

 

 

 

 

 

 

 

 

5

 

 

 

 

(6cos15 sin15)e

6cos15

 

 

 

 

 

 

C2

=

 

 

 

 

 

.

 

 

37sin15

 

 

 

 

 

 

 

Вычисляя с точностью до двух знаков после запятой, получаем

C1 0,16C2 ≈ −31,94.

Окончательно,

у(t) et (0,16cos3t 31,94sin 3t) 0,16cos3t +0,03sin 3t

искомое отклонение в любой момент времени t.

Ответ: у(t) et (0,16cos3t 31,94sin 3t) 0,16cos3t + 0,03sin 3t .

2.СИСТЕМЫ ДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИЙ

2.1.Пусть заданы функции Φ и Ψ переменных t, x и y . Нор-

мальной системой двух дифференциальных уравнений первого порядка называется система вида

dx = Φ(t, x, y)dt

dy = Ψ(t, x, y).dt

Пара функций x = x(t), y = y(t) называется решением системы диф-

ференциальных уравнений вида, если при подстановке в каждое из уравнений эти функции x(t), y(t) обращают его в тождество.

Класс вида

x = x(t,C1,C2 )y = y(t,C1,C2 )

называется общим решением системы, если

1) при всех значениях произвольных постоянных C1, C2 соответствующая пара функций {x, y} является решением системы;

92

2) для заданных начальных условий

x(t0 ) = x0 , y(t0 ) = y0

постоянные C1, C2 могут быть однозначным образом найдены так, что

соответствующее решение удовлетворит начальным условиям.

С точки зрения механики решить систему – значит восстановить за-

 

 

dx

dy

кон движения точки по известному вектору скорости v =

 

i +

 

j .

 

 

 

 

dt

dt

2.2. Ограничимся алгоритмом нахождения общего решения так называемой линейной однородной системы с постоянными коэффициентами, т.е. системы вида

dx

 

= ax +by

 

 

 

 

 

 

 

 

 

a, b, p, q = const .

(22)

dt

 

dy

 

= px + qy,

 

 

 

 

 

dt

 

 

 

Используем обозначения

x′ = dxdt , y′ = dydt ,

и будем считать для определенности p 0 . Выразим х из второго уравнения системы (22) в виде

x =

1

( y′−qy) .

(23)

p

Теперь второе уравнение в (22) продифференцируем по переменной t и подставим в него x′ = ax +by и х (из (23)). Последовательно получим

y′′ = px′+ qy;

y = p(ax +by) + qy;

y′ = a( y′−qy) + pby +qy.

Тогда система принимает вид

y′′−(a +q)y′+(aq pb)y = 0

 

 

1

 

(24)

 

( y′−qy).

x =

 

 

p

 

 

 

 

93

Первое соотношение в (24) – это однородное дифференциальное уравнение второго порядка с постоянными коэффициентами. Легко проверить, что его характеристическое уравнение может быть записано следующим образом:

a −λ

b

λ

 

(25)

= 0 .

p

q

 

 

В соответствии с характером корней λ1, λ2 найдём фундаментальную систему решений y1 и y2 , а значит и общее решение

y = C1y1(t) +C2 y2 (t) .

Затем из равенства (23) находим х. В результате будет найдено общее решение системы (22).

2.3. При решении конкретных систем вида (22) можно сразу поступать следующим образом:

а) составить характеристическое уравнение вида (25);

б) в соответствии с его корнями λ1, λ2 построить ФСР {y1(t), y2(t)} и общее решение

y = C1 y1(t)

x = 1p ( y′−

+C2 y2 (t)

qy).

Пример. Найти общее решение системы

 

x

= −6x +4y

 

 

 

 

 

 

 

 

 

 

 

 

y′ = 9x 6y.

 

 

Решение. Запишем характеристическое уравнение (25) с

a = −6, b = 4, p = 9, q = −6

 

 

 

 

 

 

 

 

 

6 −λ

 

 

4

 

 

= 0

 

 

 

 

 

 

 

 

9

6 −λ

 

 

 

т.е. (6 −λ)2 36 = 0 , откуда λ+6 = ±6 и λ1 = 0,λ2 = −12 .

Следовательно,

 

 

 

 

 

 

 

 

y = e0

=1; y

2

= e12t

;

1

 

 

 

 

 

 

 

y = C +C

2

e12t .

 

 

1

 

 

 

 

 

 

 

94

Далее, из второго уравнения системы

 

 

 

1

( y′+6y) .

 

 

 

 

 

x =

 

 

 

 

 

 

 

9

 

 

Поскольку

 

 

 

 

 

 

 

 

 

 

 

y′ = (C1 +C2e12t )= −12C2e12t ,

 

то

 

 

 

 

 

 

 

 

 

 

 

x =

1

(12C2e12t +6C1

+6C2e12t ) или x =

2

(C1 C2e12t ).

9

3

 

 

 

 

 

 

 

 

 

 

Итак,

 

 

 

 

 

 

 

 

 

 

 

 

 

x =

2

(C1 C2e12t )

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

y = C

+C

2

e12t .

 

 

 

 

 

1

 

 

 

 

3. ТИПОВЫЕ ЗАДАНИЯ ПО РЕШЕНИЮ СИСТЕМ ДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИЙ

3.1. Типовое задание с решением. Пусть движение материальной точки на плоскости описывается системой дифференциальных уравнений

x′ = 4x yy′ = 5x +2y.

Здесь t – время, (x(t); y(t)) – координаты точки в момент t, (x(t); y(t)) –

скорость точки в момент t. Найти закон движения {x(t) и y(t)}.

Решение. Имеем систему линейных однородных уравнений первого порядка. Построим характеристическое уравнение

4 −λ

1

 

= 0

5

2 −λ

 

или (4 −λ)(2 −λ) +5 = 0 , так что λ2 6λ+13 = 0 , откуда λ1,2 = 3 ±2i.. Следовательно,

y1 = e3t cos 2t; y2 = e3t sin 2t; y = e3t (C1 cos2t +C2 sin 2t) .

95

Второе уравнение системы запишем в виде

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x =

 

1

( y′−2y) .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5

 

 

 

 

Теперь

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y′ = 3e3t (C1 cos 2t +C2 sin 2t) +e3t (2C1 sin 2t +2C2 cos 2t) ,

или

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y′ = e3t ((3C1 +2C2 ) cos 2t +(3C2 2C1)sin 2t) .

 

 

Тогда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x =

1

(e3t ((3C + 2C

2

)cos2t +(3C

2

 

2C )sin 2t) 2e3t (C cos2t +C

2

sin 2t)) ,

 

 

5

 

1

 

 

 

 

 

 

 

 

 

1

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

следовательно,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x =

1

e3t ((C + 2C

2

)cos2t +(C

2

2C )sin 2t) .

 

 

 

 

 

 

 

 

 

 

 

 

5

 

 

1

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ответ:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x =

 

1

e3t ((C

+2C

2

) cos 2t +(C

2

2C )sin 2t)

 

 

 

 

 

 

 

 

 

 

 

 

 

5

 

1

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y = e3t (C

cos 2t +C

2

sin 2t).

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

3.2. В следующих задачах для самостоятельного решения число N от 1 до 20 задаётся преподавателем.

Задание. Пусть движение материальной точки на плоскости описывается системой дифференциальных уравнений

x′ = a11x +a12 y

y′ = a21x +a22 y.

Здесь t время, (x(t); y(t)) – координаты точки в {x(t) и y(t)} момент t, ((x(t); y(t)) – скорость точки в момент t. Найти закон движения.

1.

x′ =

2x + y

 

 

 

y′ = (N 2)x +(N 1) y

x′ = 5x +(N 5) y

2.y′ = 3x +(N 3) y

3.x′ = − x +(N 2 +4) yy′ = − x +3y

4.x′ = (1N )x + yy′ = (N +1)x + y

x′ = (2N +1)x + y

5.y′ = − x +(2N 1) y

x′ = (2 + N)x + Ny

6.y′ = − Nx +(2 N) y

96

x′ =

7.y′ =

x′ =

8.y′ =

x′ =

9.y′ =

10.x′ =y′ =

11.x′ =y′ =

x′ =

12.y′ =

x′ =

13.y′ =

x′ =

14.y′ =

15.x′ =y′ =

16.x′ =y′ =

x′ =

17.y′ =

x′ =

18.y′ =

(N 2)x y 5x +(N +2) y

Nx 2y

N 2 x +3Ny

(N +2)x + Ny

(2 N)x +(4 N) y

(N 1)x +(6 N) y Nx +(5 N) y

x +(4 + N 2 ) y

x +5y

3Nx 10y N 2 x + Ny

Nx y

(N 2 +4)x Ny

(N + 2)x (N + 4) y 2x + 4y

(N +5)x + Ny

(N +10)x (N +5) y

(N 2)x Ny

Nx (2 + N ) y

(N +3)x +3y 3x +(N 3) y

4x y

(N 2 +16)x +4y

x′ =

19.y′ =

20.x′ =y′ =

21.x′ =y′ =

22.x′ =y′ =

23.x′ =y′ =

24.x′ =y′ =

x′ =

25.y′ =

26.x′ =y′ =

27.x& =y& =

28.x′ =y′ =

29.x′ =y′ =

30.x′ =y′ =

5x (N +5) y

2x (N +2) y

(N +4)x + Ny

(N +8)x (N +4) y

(N +2)x 2y

4x +(N +2) y

6x +2y

(N 6)x +(N 2) y

(N +1)x + y 2x +(N 1) y

(N 2)x + Ny

(6 N)x +(4 N) y

(3 N )x y

(N 2 +1)x + (3 + N) y

Nx + Ny

(5 N)x +(5 N) y

(N +2)x (N +2) y 2x 2y

(N 3)x 13y x +(N +3) y

(N + 5)x + 5y 5x + (N 5) y

Nx Ny

(6 + N)x +(6 + N) y

97

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]