- •Глава 1. МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ. ПОНЯТИЙНЫЙ АППАРАТ
- •1. МОДЕЛЬ КАК НЕОПРЕДЕЛЯЕМОЕ ПОНЯТИЕ
- •2. МАТЕМАТИЧЕСКАЯ МОДЕЛЬ: ХАРАКТЕРИСТИКИ И ЭТАПЫ МОДЕЛИРОВАНИЯ
- •3. МОДЕЛИРОВАНИЕ СИСТЕМ
- •4. КЛАССИФИКАЦИЯ СИСТЕМ. СИСТЕМНЫЙ ПОДХОД
- •5. ТРЕБОВАНИЯ К МОДЕЛЯМ
- •Глава 2. УРАВНЕНИЯ ПЕРВОГО ПОРЯДКА
- •1. ОСНОВНЫЕ ПОНЯТИЯ
- •2. УРАВНЕНИЯ С РАЗДЕЛЯЮЩИМИСЯ ПЕРЕМЕННЫМИ
- •3. ОДНОРОДНОЕ УРАВНЕНИЕ ПЕРВОГО ПОРЯДКА
- •4. ЛИНЕЙНЫЕ ДИФФЕРЕНЦИАЛЬНЫЕ УРАВНЕНИЯ ПЕРВОГО ПОРЯДКА
- •5. ТИПОВЫЕ ЗАДАНИЯ ПО ПОСТРОЕНИЮ МАТЕМАТИЧЕСКИХ МОДЕЛЕЙ СРЕДСТВАМИ УРАВНЕНИЙ ПЕРВОГО ПОРЯДКА
- •Глава 3. УРАВНЕНИЯ ВТОРОГО И ВЫСШИХ ПОРЯДКОВ
- •1. УРАВНЕНИЯ, ДОПУСКАЮЩИЕ ПОНИЖЕНИЕ ПОРЯДКА
- •3. ЛИНЕЙНЫЕ ОДНОРОДНЫЕ УРАВНЕНИЯ ВТОРОГО ПОРЯДКА
- •6. ТИПОВЫЕ ЗАДАНИЯ ПО РЕШЕНИЮ ЛИНЕЙНЫХ ДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИЙ ВТОРОГО ПОРЯДКА
- •Глава 4. КРАЕВЫЕ ЗАДАЧИ И СИСТЕМЫ ДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИЙ
- •1. КРАЕВЫЕ ЗАДАЧИ
- •2. СИСТЕМЫ ДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИЙ
- •3. ТИПОВЫЕ ЗАДАНИЯ ПО РЕШЕНИЮ СИСТЕМ ДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИЙ
- •Глава 5. ПРИБЛИЖЁННОЕ РЕШЕНИЕ ДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИЙ
- •1. РАЗЛОЖЕНИЕ РЕШЕНИЯ В СТЕПЕННОЙ РЯД
- •2. МЕТОД ЭЙЛЕРА ПРИБЛИЖЁННОГО РЕШЕНИЯ ЗАДАЧИ КОШИ ДЛЯ ОБЫКНОВЕННОГО ДИФФЕРЕНЦИАЛЬНОГО УРАВНЕНИЯ ПЕРВОГО ПОРЯДКА
- •3. ТИПОВЫЕ ЗАДАНИЯ ПО ПРИБЛИЖЁННОМУ РЕШЕНИЮ ДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИЙ
- •ЗАКЛЮЧЕНИЕ
Решая систему, получаем:
|
|
|
6 |
|
|
|
|
C1 |
= |
|
|
|
|
|
|
37 |
|
|
|
|
|||
|
|
|
|
5 |
|
|
|
|
|
(6cos15 −sin15)e |
−6cos15 |
|
|||
|
|
|
|
|
|||
C2 |
= |
|
|
|
|
|
. |
|
|
37sin15 |
|
||||
|
|
|
|
|
|
||
Вычисляя с точностью до двух знаков после запятой, получаем
C1 ≈ 0,16C2 ≈ −31,94.
Окончательно,
у(t) ≈e−t (0,16cos3t −31,94sin 3t) −0,16cos3t +0,03sin 3t
искомое отклонение в любой момент времени t.
Ответ: у(t) ≈e−t (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 |
= e−12t |
; |
|||||
1 |
|
|
|
|
|
|
|
||
y = C +C |
2 |
e−12t . |
|
|
|||||
1 |
|
|
|
|
|
|
|
||
94
Далее, из второго уравнения системы
|
|
|
1 |
( y′+6y) . |
|
|
||||||
|
|
|
x = |
|
|
|
|
|||||
|
|
|
9 |
|
|
|||||||
Поскольку |
|
|
|
|
|
|
|
|
|
|||
|
|
y′ = (C1 +C2e−12t )′ = −12C2e−12t , |
|
|||||||||
то |
|
|
|
|
|
|
|
|
|
|
|
|
x = |
1 |
(−12C2e−12t +6C1 |
+6C2e−12t ) или x = |
2 |
(C1 −C2e−12t ). |
|||||||
9 |
3 |
|||||||||||
|
|
|
|
|
|
|
|
|
|
|||
Итак, |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
x = |
2 |
(C1 −C2e−12t ) |
|
|
||||||
|
|
|
|
|
||||||||
|
|
|
3 |
|
|
|
|
|
|
|||
|
|
y = C |
+C |
2 |
e−12t . |
|
|
|||||
|
|
|
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′ = (1− N )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
