Матика 2 курс / Дифференциальные уравнения 2б
.pdf121
ее характеристического уравнения.
ПРИМЕР П. 1.1. Найти собственные значения и собственные векторы матрицы
|
|
|
|
|
A |
|
|
0 |
1 |
|
|
|
|
|
|
|
|
|
|
. |
|
|
|
|
|
|
|
|
|
1 |
0 |
|
|
|
|
|
|
|
|
|
|
||
Характеристическое уравнение матрицы имеет вид |
||||||||||
|
|
|
|
|
1 |
|
2 1 0. |
|||
|
|
|
|
|
||||||
|
|
|
|
1 |
|
|
||||
|
|
|
|
|
|
|
|
|
||
Решая |
его, |
находим |
собственные значения матрицы |
|||||||
1 1; |
2 1. Для того, чтобы найти собственные векторы, под- |
|||||||||
ставим значения 1 |
и последовательно в систему |
|||||||||
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
(a11 )x1 a12 x2 0 |
|||||||
|
|
|
|
(a22 )x2 0. |
||||||
|
|
|
a21x1 |
|||||||
Подставляя 1 1, имеем систему уравнений |
||||||||||
x1 x2 |
0 |
|
|
|
|
|
|
|||
|
x1 x2 0 |
|
или |
x1 x2 , |
|
|||||
|
|
|
|
|
|
|
|
|
которая имеет множество решений. Обозначая x1 C , где C - про-
извольная |
постоянная, |
получаем |
x2 C . |
Тогда в |
качестве соб- |
||||||
ственного |
вектора |
можно |
взять |
простейший |
вектор |
||||||
|
1 |
|
|
t11 |
|
|
|
|
|
|
|
e |
|
|
или e |
|
|
При подстановке 2 |
1 получаем |
||||
1 |
|
|
|
1 |
|
. |
|||||
|
1 |
|
|
t21 |
|
|
|
|
|
|
|
систему двух одинаковых уравнений |
|
|
|
||||||||
|
|
|
|
|
x1 x2 |
0 |
или x1 x2. |
|
|
Поэтому в качестве второго собственного вектора можно взять век-
|
|
|
|
|
|
|
|
|
122 |
|
|
|
1 |
или |
|
t12 |
|
||
|
e |
|
|
e |
|
|
. |
||
тор |
2 |
|
1 |
|
2 |
t |
22 |
|
|
|
|
|
|
|
|
|
|
Справедлива теорема: матрицы, связанные преобразованием
подобия (эквивалентные матрицы), имеют одинаковые собственные значения. Другими словами, собственные значения матрицы явля-
ются основной характеристикой матрицы, которая не изменяется при переходе от одного базиса к другому. При этом простейший вид J , к которому может быть приведена матрица A при выборе соответствующего базиса, определяется характером собственных значений матрицы A и называется жордановой формой J матрицы
A. Матрицы A и J связаны преобразованием подобия
J T 1AT; |
A TJT 1 . |
|
(П. 1.5) |
|||
Для матрицы второго порядка характеристическое уравнение |
||||||
имеет вид |
|
|
|
|
|
|
|
a11 |
a12 |
|
0 |
, |
(П. 1.6) |
|
|
|||||
|
a21 |
a22 |
|
|||
|
|
|
|
|
|
2 (a a ) a a |
22 |
a a |
0 |
. |
||||||
|
11 |
22 |
11 |
21 |
12 |
|
|||||
Обозначив a11 a22 |
SpA- след матрицы A, |
||||||||||
a11a22 |
a21a12 |
|
a11 |
a12 |
|
DetA - определитель матрицы A, |
|||||
|
|
||||||||||
|
|
|
a21 |
a22 |
|
|
|
|
|
|
|
получим |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
SpA DetA 0 . |
|
|
(П.1.7) |
||||
Корни уравнения (П.1.7) записывают в виде: |
|
|
|
123
1,2 |
|
SpA |
(SpA)2 4DetA |
. |
(П. 1.8) |
|
2 |
||||
|
|
|
|
|
При этом возможны случаи:
0. Корни 1 2 действительные и раз-
личные; |
|
|
|
|
|
2) |
Sp2А - 4 DetA = 0. Корни |
|
|
действительные и |
|
|
|
1 |
2 |
0 |
|
одинаковые; |
|
|
|
|
|
3) |
Sp2А - 4 DetA < 0. Корни |
i |
являются ком- |
||
|
|
1,2 |
|
|
|
плексными. |
|
|
|
|
|
Для любой действительной матрицы |
A второго порядка су- |
ществует действительная матрица Т (DetT ≠ 0 ), переводящая мат-
рицу A к жордановой форме J . При этом вид матрицы J определя-
ется характером собственных значений матрицы A, согласно (П.
1.8) по таблице П.1.1
Таблица П.1.1. Действительные жордановы формы матрицы
n 2
Характер соб- |
Вид матрицы J |
||||||
ственных значений |
|||||||
|
|
|
|
||||
|
|
|
|
|
|
||
|
|
1 |
0 |
|
|||
|
|
|
|
|
|||
|
1 |
2 |
|
0 |
|
|
|
|
|
|
|||||
|
|
|
|
|
2 |
|
|
|
|
0 |
1 |
|
|||
|
|
|
|
||||
1 |
2 |
0 |
|
0 |
|
|
|
|
|
|
|||||
|
|
|
|
|
0 |
|
124
1,2 i |
|
|
|
0 |
||
|
|
|
|
, |
||
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
Способы нахождения матрицы перехода Т.
Способ 1. Найдите собственные значения и собственные век-
торы матрицы A. При этом координаты собственных векторов
|
t11 |
|
e |
|
|
1 |
|
|
|
t21 |
|
|
|
t12 |
|
и |
e |
|
|
1 |
|
|
|
|
|
t22 |
|
образуют столбцы матрицы T |
t |
|
t |
|
A |
||
11 |
12 , переводящей матрицу |
||||||
|
|
|
|
|
|
|
|
|
|
|
t21 |
t22 |
|
||
|
|
|
|
|
|
|
|
из базиса E1 i , |
j к базису |
E2 e1, e2 . Продолжая пример (П. 1.1), |
|||||
можно написать |
|
|
|
|
|
|
|
|
|
T |
1 |
|
1 |
|
|
|
|
|
|
|
. |
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
1 |
|
|
|
|
|
1 |
|
1 |
|
|
С учетом того, |
|
|
|
|
|
|
|
что T 1 1/ 2 |
|
|
|
, в базисе из собственных |
|||
|
|
|
1 |
|
1 |
|
|
|
|
|
|
|
|
векторов исходная матрица имеет диагональный вид
J T |
1 |
1 0 |
1 1 |
1 |
1 |
0 |
|||||
1AT 1/ 2 |
|
|
|
|
|
|
|
|
|
. |
|
|
|
1 |
1 |
|
1 |
0 |
|
|
|
0 |
|
|
|
|
1 |
1 |
|
1 |
Способ 2. Найти собственные значения матрицы 1, |
2 . |
По |
||
характеру собственных значений установить вид матрицы |
J |
со- |
||
гласно таблице П.1.1. Далее, решая матричное уравнение |
|
|
|
|
J T 1AT, |
TJ AT , |
(П. 1.9) |
||
найти матрицу T . |
|
|
|
|
|
|
|
125 |
|
|
|
Продолжение примера (П.1.1). Так как 1 1; |
2 1, то |
|||
|
1 |
0 |
|
||
J |
|
. |
|
||
|
0 |
|
|
||
|
1 |
|
|
a |
|
b |
|
|
|
|
Берем за матрицу T |
|
, тогда уравнение (П.1.9) имеет вид |
|||||
|
|
|
|
|
|
|
|
|
c |
|
d |
|
|
|
|
a |
b 1 |
0 |
0 |
1 a |
b |
|
|
|
|
|
|
|
|
|
; |
|
|
0 |
|
|
|
|
|
c |
d |
1 |
1 |
0 c |
d |
|
|
|
a |
b |
c |
d |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
. |
|
|
|
c |
d |
a |
b |
|
|
Приравнивая элементы матрицы, стоящие на одинаковых местах,
получаем систему решений
a c, |
|
b d. |
|
Выбирая простейшие решения a 1, |
b 1, получаем |
||
T |
1 |
1 |
|
|
|
. |
|
|
|
|
|
|
1 |
1 |
Правильность нахождения Т можно проверить подстановкой в уравнение J= T-1AT.
126
ПРИЛОЖЕНИЕ 2.
ДИФФЕРЕНЦИАЛЬНЫЕ УРАВНЕНИЯ И СИСТЕМЫ УРАВНЕНИЙ ПОРЯДКА n 2
Дифференциальным уравнением порядка |
n называют урав- |
нение вида |
|
F (t, y, y(1) , y(2) ,..., y(n) ) 0 . |
(П.2.1) |
Общее решение уравнения содержит n произвольных посто- |
|
янных и имеет вид |
|
y (t,C1,C2,...,Cn ) . |
(П.2.2) |
Для определения произвольных постоянных C1,C2 ,. , Cn и
нахождения частного решения задают n начальных условий
y t |
0 |
y ; y(1) |
(t |
0 |
) y (1) |
; |
y(2) t |
0 |
y 2 |
, , |
y(n 1) t |
0 |
y |
n 1 |
. (П.2.3) |
|
0 |
|
0 |
|
|
0 |
|
|
|
0 |
|
В некоторых случаях порядок уравнения может быть понижен
путем введения новой переменной. Так если уравнение явно не со-
держит функцию y , то порядок легко понижается введением z y 1 t . Если уравнение явно не содержит независимую пере-
менную t , то для понижения порядка делаем подстановку y 1 z y dydt ;
y 2 dydz dydt z dydz .
ПРИМЕР П. 2.1. Дифференциальное уравнение порядка n 3
127
y 3 cos t
явно не содержит искомую функцию. Решение легко находится пу-
тем трехкратного интегрирования
|
y 2 |
sin t C ; |
|
|
|||||||
|
|
|
|
|
|
1 |
|
|
|
||
|
y 1 cos t C t C |
2 |
; |
||||||||
|
|
|
|
|
|
1 |
|
|
|||
|
|
|
C t 2 |
|
|
|
|
|
|||
|
y sin t |
1 |
|
C2t C3. |
|||||||
|
2 |
|
|||||||||
|
|
|
|
|
|
|
|
|
|||
При нулевых начальных условиях |
y 0 y 1 (0) y 2 0 частное ре- |
||||||||||
шение имеет вид y t sin t. Покажите это самостоятельно. |
|||||||||||
ПРИМЕР П. 2.2. Уравнение y 2 |
y 1 |
|
t |
||||||||
t |
|||||||||||
|
|
|
|
|
|
|
|
|
|||
явно не |
содержит переменную |
|
у. |
Поэтому подстановкой |
|||||||
z y 1 t , |
z(1) (t) y(2) (t) оно приводится к линейному уравнению |
||||||||||
первого порядка относительно переменной |
|
z |
|||||||||
|
z 1 |
z |
|
t. |
|
|
|||||
|
t |
|
|
||||||||
|
|
|
|
|
|
|
|
|
Покажите самостоятельно, что общее решение уравнения имеет вид
|
t 3 |
C t 2 |
|
||
y |
|
|
1 |
C2. |
|
3 |
2 |
||||
|
|
|
|||
ПРИМЕР П. 2.3. Уравнение |
y(2) y 3 явно не содержит не- |
зависимую переменную t. Поэтому подстановка
|
|
|
128 |
|
|
|
|
|
y |
1 z y |
dy |
; |
|
y |
2 z |
dz |
|
dt |
|
dy |
||||||
|
|
|
|
|
|
приводит исходное уравнение к уравнению с разделяющимися пе-
ременными zdz y 3dy ,
решением которого является функция:
z |
|
(2C1 y 2 1) |
|
dy |
|
|
|
|
. |
||
|
y |
dt |
Окончательно, общее решение исходного уравнения имеет вид:
(2C1 y 2 1)1 (t C2 ). 2C1
ПРИМЕР П. 2.4. В уравнении третьего порядка
( y(2) )2 y(1) y(3) |
|
1 |
|
( y(1) )2 |
t2 |
||
|
правую и левую часть уравнений можно представить как производ-
ные некоторых функций
|
|
y |
2 1 |
|
|
1 |
|
1 |
|||
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
C1 ; |
|||
|
|
|
|
|
|||||||
|
|
y |
(1) |
|
|
|
t |
|
|
||
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
y 2 |
|
1 |
C |
||
|
|
|
|
||
y(1) |
|
|
t |
1 . |
Полученное уравнение второго порядка не содержит явно перемен-
ную у и приводится к уравнению первого порядка подстановкой y 1 z t ; y 2 z 1 t .
129
Покажите самостоятельно, что общее решение уравнения имеет вид
y(t) C2e |
|
|
|
|
t |
|
|
|
e |
C1t |
|
|
|
C1t |
tdt C3 |
C2 |
|
e |
C1t |
|
|
) |
|
C3. |
|||
|
|
C |
|
|
2 |
|
|||||||
|
|
|
|
|
1 |
|
|
|
C1 |
|
|
|
Общее решение однородного линейного уравнения порядка n с
постоянными коэффициентами:
|
a |
n |
y(n) |
a |
n 1 |
y(n 1) |
... a y(1) a |
0 |
y 0 |
(П.2.4) |
|||||
|
|
|
|
|
|
|
1 |
|
|
|
|
||||
записывается |
|
|
как |
линейная |
комбинация |
базисных |
функций |
||||||||
y1 t , y2 t , , yn t , |
|
образующих |
фундаментальную |
систему |
|||||||||||
решений, и имеет вид |
|
|
|
|
|
|
|
|
|
|
|||||
yoo C1y1(t) C2 y2 (t) ... Cn yn (t). |
|
(П.2.5) |
|||||||||||||
При этом каждому простому действительному корню характе- |
|||||||||||||||
ристического уравнения |
i |
в общем решении соответствует |
|||||||||||||
слагаемое Ci e it . |
|
|
|
|
|
|
|
|
|
|
|
||||
Каждому действительному корню кратности |
r в общем ре- |
||||||||||||||
шении соответствует слагаемое |
|
|
|
|
|
|
|
||||||||
P |
|
t e 1t |
C C t C |
|
t r 1 |
e i t . |
|
||||||||
r 1 |
|
|
|
|
|
|
0 |
1 |
|
r 1 |
|
|
|
||
Каждой паре комплексно-сопряженных |
корней 1,2 i |
||||||||||||||
соответствует слагаемое |
|
|
|
|
|
|
|
|
|
e t (C1 cos t C2 sin t) .
Если комплексно-сопряженные корни имеют кратность r , то в общем решении присутствует слагаемое
130
e t Pr 1 t cos t Qr 1 t sin t .
ПРИМЕР П. 2.5. Найти решение уравнения
y 3 y 1 0
при условии y 0 3; y 1 0 1; y 2 0 1. Все корни характери-
стического уравнения действительные и различные
3 1 1 0;
1 0; 2 1; 3 1.
Поэтому общее решение имеет вид
yoo C1 C2et C3e. t
Покажите, что с учетом начальных условий частное решение имеет вид
yoo 2 e t .
ПРИМЕР П. 2.6. Для уравнения
y 3 3y 2 3y 1 y 0
характеристическое уравнение
3 32 3 1 1 3 0
имеет один действительный корень 1 |
кратности r 3. По- |
этому |
|
yoo (C0 C1t C2t 2 )e t .
ПРИМЕР П. 2.7. Для уравнения y 4 y 0 характеристическое
уравнение