
- •СОДЕРЖАНИЕ
- •От автора
- •§1. Основные определения
- •§2. Интегрирование в квадратурах дифференциальных уравнений
- •§3. Теоремы существования и единственности решения задачи Коши
- •§1. Существование и единственность решения
- •§2. Интегрирование в квадратурах
- •§1. Дифференциальные уравнения высших порядков
- •§2. Интегрирование уравнений с помощью степенных рядов. Асимптотика
- •§3. Системы дифференциальных уравнений в нормальной форме
- •§1. Непрерывная зависимость решений от параметров
- •§2. Дифференцирование по параметрам и начальным значениям
- •§3. Уравнения в вариациях. Линеаризация
- •§1. Линейные однородные уравнения n-го порядка
- •§2. Линейные неоднородные уравнения
- •§1. Линейные однородные системы
- •§2. Линейные неоднородные системы
- •§3. Функции от матриц. Матричная экспонента
- •ГЛАВА 7. ЛИНЕЙНЫЕ УРАВНЕНИЯ ВТОРОГО ПОРЯДКА
- •§1. Линейное уравнение второго порядка
- •§2. Краевая задача. Функция Грина
- •§3. Собственные значения и интегральные уравнения
- •§4. Интегрирование уравнений с помощью рядов Фурье
- •§5. Уравнение Бесселя
- •ГЛАВА 8. ЭЛЕМЕНТЫ ТЕОРИИ УСТОЙЧИВОСТИ
- •§1. Устойчивость по Ляпунову
- •§2. Устойчивость линейных систем
- •§3. Устойчивость по первому приближению
- •§4. Нелинейные системы. Функция Ляпунова. Теоремы Ляпунова
- •§5. Теоремы Барбашина и Красовского
ГЛАВА 6. СИСТЕМЫ ЛИНЕЙНЫХ ДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИЙ
§1. Линейные однородные системы
Системой линейных однородных дифференциальных уравнений называется система следующего вида:
n |
|
|
|
|
|
|
|
y'i aik (t) yk , |
i |
|
|
|
(1.1) |
||
1,n, |
|
||||||
k 1 |
|
|
|
|
|
|
|
где функции aik (t) непрерывны на некотором интервале T . |
|
|
|||||
Линейные системы допускают более простую форму записи. |
|
Введем |
|||||
|
T |
|
|
a11 (t) a1n (t) |
|
||
следующие обозначения: y(t) (y1(t),..., yn(t)) |
A(t) |
|
|
. Тогда |
|||
, |
.................. |
|
|||||
|
|
|
|
|
|
|
|
систему (1.1) можно записать в виде матричной формы |
an1 (t) ann (t) |
|
|||||
|
|
|
|
||||
y' A(t)y. |
|
|
|
|
|
|
(1.2) |
Решением системы (1.2) называется непрерывно дифференцируемая векторфункция y(t), которая при подстановке в систему обращает все уравнения в
тождества. Если векторы y1(t),...,ym(t) |
являются решениями системы (1.2), то |
m |
(t) , где C1,..., Cm – действительные |
любая линейная комбинация Ci yi |
|
i 1 |
|
произвольные постоянные, вновь является решением системы (1.2).
Если требуется найти решение системы (1.2), удовлетворяющее условию y(t0 ) y0 , то говорят, что для системы (1.2) поставлена задача Коши и
записывают ее в виде |
|
|
y' A(t)y , |
y(t0 ) y0 . |
(1.3) |
Замечание. Из общей теоремы Коши о существовании и единственности решения нормальной системы дифференциальных уравнений сразу следует, что, если элементы матрицы A(t) – функции aik (t) непрерывны на отрезке
T [t0,t*], то на нем существует единственное решение задачи Коши (1.3). Теорема Коши не только дает условия существования и единственности
решения задачи |
Коши. |
При выполнении этих |
условий решение |
|||||
y(t) (y (t), , y |
n |
(t)) |
определено |
на отрезке t |
,t* , т.е. в концевых точках t t |
0 |
и |
|
1 |
|
|
0 |
|
|
|
||
t t* функции y (x), i 1,2, ,n, имеют односторонние производные. |
|
|
||||||
|
|
i |
|
|
|
|
|
|
Теорема |
|
Коши дает |
достаточные условия |
существования |
|
и |
единственности, т.е. задача Коши может иметь единственное решение и при невыполнении какого-либо из требований к элементам матрицы
73
Напомним, |
что векторы y1(t),..., yn(t) |
называются линейно-независимыми на |
||||||
интервале |
T , если тождество |
C1y1(t) Cnyn(t) 0 выполняется тогда и |
||||||
только тогда, когда все коэффициенты Ci равны нулю одновременно. |
|
|||||||
Пусть |
векторы 1(t),..., n(t) |
непрерывны на интервале |
T . |
|||||
|
|
1, , n |
|
11(t) 1n(t) |
|
|
|
|
|
|
|
|
|
||||
Функциональный определитель |
|
|
|
, где ik (t) – |
||||
|
|
|
|
|
n1(t) nn(t) |
|
|
|
координаты вектора k (t) , называется определителем Вронского, либо
вронскианом векторов 1(t),..., n(t) .
Теорема 1. Если определитель Вронского решений y1(x), , yn(x) системы
(1.2) тождественно равен нулю на интервале T , тогда эти решения линейно зависимы на T . Если определитель Вронского решений y1(x), , yn(x)
уравнения (1.2) не равен нулю ни в одной точке интервала T , тогда эти решения линейно независимы на T .
Замечание. Если 1(t),..., n(t) – произвольные векторы, то из равенства нулю их определителя Вронского, вообще говоря, не следует их линейная
зависимость. Действительно, рассмотрим два вектора |
1(t) ( 0 1)T , |
||||||||
2(t) (t 0 )T . Векторы линейно-независимы, например, на |
отрезке |
1,1 , |
|||||||
так как условие |
C (t) C |
2 |
(t) 0 выполняется только при |
C 0 и |
C |
2 |
0. |
||
|
1 |
1 |
2 |
|
1 |
|
|
Несмотря на это, определитель Вронского рассматриваемых векторов имеет нулевую строку и поэтому тождественно равен нулю. Предположив, что эти векторы являются решениями некоторой системы второго порядка. Получаем противоречие с результатом теоремы.
Часто бывает полезен другой критерий линейной независимости
произвольных векторов. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
Теорема 2. Для того чтобы произвольные |
векторы i |
(t), i 1, |
2,…,n были |
|||||||||||||||||
линейно независимы на a,b , необходимо и достаточно, |
чтобы определитель |
|||||||||||||||||||
|
|
|
|
1, , n |
|
( 1, 1)( 1, 2) ( 1, n) |
|
|
|
|||||||||||
|
|
|
|
|
|
|
||||||||||||||
Грама этих векторов |
|
|
..................................... |
|
|
был отличен |
||||||||||||||
|
|
|
|
|
|
|
( n, 1)( n, 2) ( n, n) |
|
|
|
||||||||||
|
, , |
|
0, |
|
|
|
|
|
|
|
|
b |
|
(t), |
|
(t))dt |
|
|||
от нуля: |
n |
где |
( |
i |
, |
j |
) ( |
i |
j |
– скалярное |
||||||||||
1 |
|
|
|
|
|
|
a |
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
произведение функций i(t) , j(t) i, |
j 1,n |
на отрезке. |
|
|
|
|
|
74
1.1. Фундаментальная матрица |
|
|
|
|
Определение. Фундаментальной |
системой решений y1(t),..., yn(t) системы |
|||
(1.2) называются любые n линейно независимых решений системы (1.2). |
||||
Если коэффициенты aij(t), |
i, j 1,2, ,n, в |
системе |
уравнений (1.2) |
|
непрерывны на t0,t* , |
то она |
имеет систему |
решений |
y1(t), , yn(t) , |
фундаментальную на t0,t* |
. |
|
|
|
Определение. Матрица W(t), столбцами которой являются координаты векторов, образующих фундаментальную систему решений, называется фундаментальной матрицей системы (1.2).
Определитель матрицы W(t) – это определитель Вронского системы n линейно независимых решений системы (1.2). Определитель не равен нулю,
поэтому матрица W(t) имеет обратную W 1 t . |
|
|
||
Пусть каждый из линейно независимых векторов |
y1(t),...,yn(t) является |
|||
решением |
матричного уравнения (1.2). Эту совокупность n |
векторных |
||
уравнений можно кратко записать в следующем виде: |
|
|
||
|
W (t) |
A(t) W (t) . |
|
(1.4) |
Столбцами матрицы W являются координаты векторов |
y1(t),..., yn(t). Каждый |
|||
столбец |
матрицы W (t) равняется |
произведению |
матрицы |
A t на |
одноименный столбец матрицы W(t) . Уравнение (1.4) называется матричным уравнением, сопоставленным векторному уравнению (1.2).
Очевидно, что фундаментальная матрица W(t) есть решение матричного уравнения (1.4).
Теорема 3. Линейная однородная система всегда имеет фундаментальную систему решений, а значит, и фундаментальную матрицу. По заданной системе n линейно независимых векторов y1(t),...,yn (t) можно найти единственную систему (1.2), для которой эти векторы образуют фундаментальную систему решений.
Пусть матрица W t , столбцами |
которой |
являются |
координаты |
этих |
|
векторов, удовлетворяет системе (1.4). Требуется найти матрицу |
A t . |
||||
Умножая тождество (1.4) справа на |
обратную |
матрицу |
W 1 t , |
получаем |
|
W'(t)W 1(t) A(t)W(t)W 1(t) A(t), где |
A t – |
матрица |
искомой |
системы |
|
уравнений. |
|
|
|
|
|
Замечание. Множество решений линейной однородной дифференциальной системы образует линейное пространство функций. Любая фундаментальная система решений является базисом этого пространства. Существует бесконечно
75
много фундаментальных систем решений однородной системы, переходящих одна в другую с помощью невырожденного линейного преобразования.
Общее решение линейной однородной системы (1.2) будет иметь вид
y W(t)C , |
(1.5) |
где W (t) – фундаментальная матрица, а C |
– произвольный постоянный |
вектор. Правую часть (1.5) удобно представлять как скалярное произведение строки (y1(t),...,yn(t)) векторов, столбцы координат которых образуют матрицу
W t , на столбец C : y(t) C1y1(t) Cn yn (t). Получается линейная комбинация векторов фундаментальной системы, коэффициентами которой служат координаты вектора C . По аналогии строилось общее решение линейного однородного уравнения n-го порядка.
1.2. Матрицант. Матрица Коши. Пусть фундаментальная матрица W(t)
удовлетворяет условию W (t0 ) E , т.е. векторы фундаментальной системы имеют единичную матрицу начальных значений. Тогда фундаментальная матрица W (t) называется матрицантом системы (1.2). Матрицант – аналог нормальной фундаментальной системы решений линейного однородного уравнения n-го порядка. Матрицант, очевидно, является решением следующей однородной дифференциальной системы:
|
d |
|
|
|
|
|
W (t) A(t) W (t), W (t0 ) |
E . |
(1.6) |
|
|
|||
|
dt |
|
|
|
Пусть W (t) – произвольная фундаментальная матрица. Тогда матрица |
|
|||
|
|
K(t,t0 ) W(t)W 1(t0 ), |
|
(1.7) |
зависящая от двух аргументов t и t0 , также является матрицантом системы
(1.2). Действительно, матрица K(t,t0 ) удовлетворяет следующей системе:
d |
|
|
|
K(t,t0 ) A(t) K(t,t0 ) , K(t0 ,t0 ) E. |
(1.8) |
|
||
dt |
|
Поскольку K'(t,t0) (W(t) W 1(t0))' W'(t)W 1(t0) A(t)W(t)W 1(t0) A(t)K(t,t0) и
при t t0 , выполнено начальное условие K (t0 ,t0 ) W(t0 )W 1(t0 ) E.
Принято также другое название матрицы, определяемой формулой (1.7), – матрица Коши системы (1.2).
Замечание. Матрица Коши K(t,t0 ) определяется по формуле (1.7)
единственным образом, несмотря на то, что матрица W (t) – произвольная фундаментальная матрица. Это следует из того, что матричное уравнение (1.2) с единичной матрицей начальных значений имеет единственное решение.
76
Если известна матрица Коши, то решение однородной задачи Коши
y' A(t)y , y(t0 ) y0 имеет следующий вид: y K(t,t0 )y0 . Справедливо и обратное. Если для произвольного вектора начальных значений y0 решение однородной задачи Коши y' A(t)y , y(t0 ) y0 может быть записано в форме
y K(t,t0 )y0 , то входящая в нее матрица является матрицей Коши. Замечание. Полученное представление является аналогом формулы, выражающей решение задачи Коши для линейного однородного уравнения n-го порядка через нормальную фундаментальную систему решений.
Поясним физический смысл матрицы Коши. Для этого рассмотрим следующую однородную дифференциальную систему:
|
|
|
|
d |
y A(t) y,t [t0,t |
* |
], |
|
||||||
|
|
|
|
|
|
|
|
(1.9) |
||||||
|
|
|
|
|
|
|||||||||
|
|
|
|
dt |
|
|
|
|
|
|
|
|||
|
|
|
|
y(t |
0 |
) y . |
|
|
|
|||||
|
|
|
|
|
|
|
0 |
|
|
|
||||
Пусть |
вектор |
Ki (t, ) |
является |
решением |
|
следующего |
векторного |
|||||||
дифференциального |
уравнения: |
|
|
d |
Ki(t, ) A(t) Ki(t, ) с начальным условием |
|||||||||
|
|
|
||||||||||||
|
|
|
ei (0,..., |
|
dt |
|
|
|
|
|||||
K i(t0,t0) ei. Здесь |
0,1.., 0,..., 0 )T – нулевой вектор с единицей |
|||||||||||||
на i-м |
|
|
|
|
|
|
|
|
|
i |
K(t, ) – |
|
|
|
месте. |
Обозначим |
через |
n n-матричную |
функцию, |
составленную из векторов Ki(t, ) как из столбцов, где i принимает значения от 1 до n.
Функция K(t, ) однозначно определена для любого значения аргумента t. Эта функция абсолютно непрерывна и удовлетворяет следующей системе дифференциальных уравнений:
d |
|
||
|
|
K(t, ) A(t) K(t, ), |
(1.10) |
|
|||
dt |
|||
|
|
|
|
K(t0,t0 ) E. |
|
||
Отсюда следует физический смысл матрицы Коши K(t, ) |
(функции Коши). |
||
Каждый i-й столбец функции |
Ki(t, ) представляет собой отклик в момент |
времени t однородной системы (1.9), которую возмутили в момент времени t0
единичным вектором y0 ei . Поэтому функции Коши часто называют еще функциями точечного источника (функциями Грина).
1.3. Линейные однородные системы с постоянными коэффициентами.
Построение ФСР линейной однородной системы с постоянными коэффициентами сводится к решению характеристического уравнения для матрицы системы и построению соответствующих собственных векторов.
77

|
|
d y1 |
|
7y 3y |
2 |
, |
|
|
|||||
|
|
dt |
1 |
|
||
|
|
|
|
|
|
|
Пример 1. Найти общее решение и ФСР системы уравнений |
|
|
|
|
|
|
d y2 |
|
6y1 4y2 , |
||||
|
|
|
|
|||
|
dt |
|
|
|
|
|
d y1(t) |
7 3 |
y1(t) |
||||||
т. е. |
|
|
y |
|
|
6 4 |
|
|
. Выпишем общее решение системы и |
|
|||||||||
|
dt |
2 (t) |
|
|
y2 |
(t) |
найдем векторы ФСР. Общее решение считается через собственные числа и соответствующие им векторы матрицы A. Получаем: A x x ; Ax Ex 0;
(A E)x 0; |
det(A E) 0. Отсюда |
7 3 |
0, |
или 2 |
11 10 0, |
|||
6 |
4 |
|||||||
|
|
|
|
|
|
|||
1 1, 2 10. |
Собственные |
числа являются |
действительными |
и простыми. |
||||
Найдем собственный вектор |
x1 , соответствующий собственному значению 1. |
Для этого решаем систему линейных алгебраических уравнений |
|
(A 1E) |
x1 0. |
|||||||
Получаем x |
|
1 |
|
|
1. Для второго вектора аналогично x |
2 |
|
1 |
2 10. |
|
1 |
|
|
|
1 |
|
|
|
|
|
|
|
|
2 |
|
|
|
|
1 |
|
|
Выписываем общее решение однородной дифференциальной системы в виде
y(t) C1 |
|
1 |
|
e |
t |
C2 |
|
1 |
10t |
. В матричной форме относительно |
C |
и C |
|
|
|
|
|
e |
2 |
||||||||
|
|
|
|
|
|
|
|
1 |
|
||||
|
|
2 |
|
|
|
|
1 |
|
|
|
|
|
|
et |
|
получаем y(t) |
|
t |
2e |
|
e10t |
C1 |
|
10t |
|
. Векторы ФСР |
e |
C2 |
|
et |
|
e10t |
|
|
|
t |
и 10t |
. |
2e |
|
e |
|
dy |
1 |
|
4y1 3y2 , |
|
|
|
|
||
|
|
|||
dt |
|
|
||
Пример 2. Найти общее решение и ФСР системы уравнений |
|
|
|
|
dy2 |
|
3y1 4y2 , |
||
|
|
|
|
|
dt |
|
|
d y1(t)
т.е.
dt y2 (t)
|
4 |
3 |
y (t) |
. Решаем характеристическое уравнение |
||
|
|
4 |
|
1 |
|
|
|
3 |
y2 (t) |
|
относительно матрицы A: det(A E) 0. Отсюда следует |
4 |
3 |
0; |
|||
3 |
4 |
|||||
|
|
|
|
|||
(4 )2 |
9; 4 3i; |
4 3i . Найдем собственный вектор x , |
||||
|
1,2 |
|
|
1 |
|
78
соответствующий собственному значению 1. Для этого решаем систему |
|
|
|
|||||||||||||||||||||||||||||||
линейных алгебраических уравнений (A 1E) |
|
x1 0. Получаем |
|
|
|
|
||||||||||||||||||||||||||||
1 |
|
|
|
4 3i. Для второго вектора аналогично x2 |
1 |
2 4 3i. |
|
|
|
|||||||||||||||||||||||||
x1 1 |
|
|
|
|
|
|||||||||||||||||||||||||||||
i |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
i |
|
|
|
|
|
|
|
|
|
Общее |
решение |
нашей |
системы |
имеет |
вид |
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||
|
1 |
(4 3i)t |
1 |
(4 3i)t |
. Переписывая в матричной форме |
|
|
|
||||||||||||||||||||||||||
y(t) c1 |
e |
|
|
c2 i e |
|
|
|
|
||||||||||||||||||||||||||
|
i |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
относительно c1 |
|
|
|
|
|
|
|
|
|
|
|
e(4 3i)t |
e(4 3i)t |
|
C1 |
|
|
|
|
|||||||||||||||
и c2 , получаем y(t) |
|
(4 3i)t |
|
|
(4 3i)t |
|
. Векторы |
|
|
|
||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ie |
|
|
|
ie |
|
|
|
C2 |
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
e(4 3i)t |
|
e(4 3i)t |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
ФСР соответственно будут |
(4 3i)t |
|
и |
|
(4 3i)t . Выделим в общем |
|
|
|
||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
ie |
|
|
|
|
ie |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
решении действительную и мнимую части. Имеем |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||
|
|
e(4 3i)t |
|
e4t |
cos3t ie4t sin3t, ie(4 3i)t |
e4t sin3t ie4t cos3t. |
|
|
|
|||||||||||||||||||||||||
Отсюда |
|
получаем |
y |
(t) |
e4t cos3t |
e4t sin 3t |
c1 |
|
Выписываем |
векторы |
||||||||||||||||||||||||
|
|
4t |
|
e |
4t |
|
|
|
|
|
. |
|||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
e |
|
sin 3t |
|
cos3t |
|
c |
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
||
|
4t |
|
|
|
|
|
e4t sin3t |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
ФСР e |
|
cos3t |
|
и |
|
|
|
. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
e4t sin3t |
|
|
e4t cos3t |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
dy1 |
|
|
5y |
y |
2 |
, |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
Пример 3. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|||||||
Найти общее решение и ФСР системы уравнений dt |
|
|
|
|
||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
dy2 |
y1 3y2. |
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
dt |
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
d y1(t) |
5 |
1 |
y1(t) |
|
|
|
|
|
|
|
|
|
|||||||||
Перепишем |
систему |
в виде |
|
|
|
|
|
|
|
|
|
|
|
. |
Решаем уравнение |
|||||||||||||||||||
|
|
|
|
|
|
|
||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
dt |
y2(t) |
1 |
3 y2(t) |
|
|
|
|
|
|
|
|
|
|||||||||
det(A E) 0: |
|
5 |
1 |
|
|
0; |
2 8 16 0; |
|
1 2 |
4. |
Если 1 |
– корень |
||||||||||||||||||||||
|
|
|
||||||||||||||||||||||||||||||||
|
|
|
|
|
|
1 |
3 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
характеристического уравнения кратностью m, то этому корню соответствует
решение |
x |
p |
(t)e 1t , x |
2 |
p |
2 |
(t)e 1t , |
…, |
x |
n |
p |
n |
(t)e 1t , |
где p (t), |
p |
2 |
(t), …, |
||||||
|
|
1 |
1 |
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|||||
pn (t) – многочисленные |
степени |
не |
|
выше |
|
m 1. Таким |
образом, |
||||||||||||||||
действительному |
|
корню |
|
4 кратностью |
два |
соответствуют |
решения |
||||||||||||||||
y |
e4t (a t a ), |
y |
2 |
e4t (b t b ). Дифференцируя y |
и y |
2 |
, |
получаем тождества |
|||||||||||||||
1 |
|
1 |
2 |
|
|
1 |
|
2 |
|
|
|
|
|
|
1 |
|
|
|
|
|
|
79
dy1 |
ae4t 4(at a |
)e4t , |
dy2 |
be4t 4(bt b )e4t . Значения |
x , |
x |
, |
dx1 |
, |
dx2 |
|||||
dt |
dt |
||||||||||||||
1 |
1 |
2 |
|
1 |
1 |
2 |
1 |
2 |
|
dt |
dt |
подставляем в нашу систему дифференциальных уравнений. После сокращения на e4t получаем два тождества
a1 4(a1t a2 ) 5(a1t a2 ) (b1t b2 ); b1 4(b1t b2 ) a1t a2 3(b1t b2 ).
Приравнивая коэффициенты при t |
и свободные члены, получаем следующие |
|||||||||||||||
системы алгебраических уравнений: |
4a 5a b, |
a 4a 5a b , |
||||||||||||||
|
1 |
|
1 |
1 |
|
1 |
2 |
2 |
2 . Тогда |
|||||||
|
|
|
|
|
4b1 a1 3b1, |
b1 4b2 a2 3b2. |
||||||||||
следует, что |
a1 b1, |
a2 b2 |
a1 |
b1. |
Полагая |
|
a1 C1, |
a2 |
C2 |
(C1,C2 – |
||||||
произвольные постоянные), находим b1 C1, |
b2 C2 |
C1 . Следовательно, |
||||||||||||||
|
|
x e4t (Сt С |
), |
x |
2 |
e4t (Сt С |
2 |
С ). |
|
|
||||||
|
|
1 |
1 |
2 |
|
|
|
1 |
|
|
1 |
|
|
|
||
|
|
|
|
|
|
4 t |
|
4 t |
|
C |
|
|
|
|
||
В матричном |
виде имеем y(t) |
e t |
|
e |
|
1 |
. |
Соответствующими |
||||||||
4 t |
(t 1) |
4 t |
|
|||||||||||||
|
|
|
|
e |
e |
C2 |
|
|
|
|
e4tt |
|
e |
|
векторами ФСР будут |
4t |
|
|
и |
e |
|
(t 1) |
e |
4t
4t
.
Замечание. Систему можно решить методом исключения. Действительно,
выразив из первого уравнения |
x2 и продифференцировав, подставим затем |
||||||||||||||||||||||||
значения |
x |
|
и |
|
x |
во второе |
|
уравнение. |
В |
|
результате получим |
линейное |
|||||||||||||
|
|
|
2 |
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
однородное уравнение второго порядка относительно x1. |
|
||||||||||||||||||||||||
Выражаем |
x |
|
из первого уравнения: |
x |
2 |
5x |
x . |
Дифференцируя, |
получаем |
||||||||||||||||
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
1 |
|
|
|
|
|
||
x |
5x x . Подставляем полученные значения во второе уравнение системы |
||||||||||||||||||||||||
|
2 |
1 |
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
и находим общее решение для x1: |
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
5x |
x |
x 15x 3x , |
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
1 |
1 |
|
|
1 |
|
|
1 |
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
x1 8x1 16x1 0. |
|
|
|
|
|
|||||||
Решая характеристическое уравнение, получаем 1 2 |
|
4. Отсюда общее |
|||||||||||||||||||||||
решение для |
|
y: |
y e4t(C C t) |
. Теперь подставляем |
x |
|
в выражение для x : |
||||||||||||||||||
|
|
|
|
|
1 |
|
1 |
|
1 |
2 |
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
2 |
x |
2 |
5x x |
, |
x |
2 |
5e4t (C C |
t) 4e4t (C C |
t) e |
4tC |
2 |
. В итоге имеем общее |
||||||||||||||
|
1 |
1 |
|
|
|
|
1 |
2 |
|
|
|
|
|
|
1 |
2 |
|
|
|
|
|||||
решение для |
x : |
x |
2 |
e4t (C |
С |
2 |
C |
2 |
t). Окончательно |
|
|
|
|
||||||||||||
|
|
|
|
|
2 |
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
e4tt |
e4t |
C1 |
|
e4tt |
e |
|||||||
y(t) |
e |
4t |
(t 1) |
e |
4t |
|
. Векторы ФСР |
|
4t |
(t |
|
и |
|
|
|
C |
|
e |
|
1) |
e |
||||
|
|
|
|
2 |
|
|
|
|
|
4t
4t
.
80