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

Линейные дифференциальные уравнения и системы линейных уравнений. Учебное пособие

.pdf
Скачиваний:
0
Добавлен:
12.08.2026
Размер:
1 Мб
Скачать

3. Линейныесистемы с постоянными коэффициентами

Решение. Будем отыскивать решение системы в виде

x 1ekt . y 2ekt , z 3ekt .

Составляем и решаем характеристическое уравнение

1 k

1

2

 

1 k 2 1 k 8 4

1 k 4

1 k 0,

 

4

1 k

0

0,

k1

1, k2

k3 1.

 

 

2

1

1 k

 

 

 

1) Строим собственный вектор, соответствующий простому корню k1 1. Решим систему (3.7), подставляя k 1 :

2 1 2 2 3 0,4 1 0,2 1 2 2 3 0.

Из последней системы находим 1

0 , 2 2 3 . Тогда собствен-

ный вектор

 

1 и первое частное решение X , соответствующие

A

 

 

 

 

 

 

 

 

1

корню k1 1, запишутся в виде

 

 

 

 

 

 

 

0

 

 

0

 

 

 

 

1

2

, X

 

2 et .

 

 

A

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

1

2) Построим два линейно независимых частных решения, соот-

ветствующих характеристическому числу

k 1, кратность кото-

рого 2 . При k 1 получаем матрицу системы (3.7):

0

1

2

 

4 2

0

 

A kE

.

 

2

1

0

 

 

 

Её размерность n 3 ,

ранг

r 2 . Число линейно независимых

собственных векторов

равно

m n r 1 .

Поскольку m , то

решение следует искать в виде

 

 

 

 

x a

a

t e t , y b

b t e t , z c

c t e t .

(3.11)

1

2

 

 

1

2

1

2

 

71

Нежельская Л. А. Линейные дифференциальныеуравнения

Для нахождения коэффициентов ai , bi , ci , i 1, 2 , подставим (3.11) в исходную систему и получим

a e t a

a t e t

b b t e t 2

c

c t e t

a

a

t e t ,

 

2

1

2

 

 

 

1

2

 

 

 

1

2

 

1

2

 

b e t b

b t e t 4

a

a t

e t

 

b

b t e t ,

 

 

 

 

 

2

1

2

 

 

 

1

2

 

 

 

1

2

 

 

 

 

c e t c

c t e t 2

a

a t e t

 

b b t e t c c t e t ,

 

 

1

2

 

 

 

1

2

 

 

 

1

2

 

1

2

 

2

 

 

 

 

 

 

 

 

или, в результате сокращения на множитель e t 0

и приведения

подобных членов, запишем

 

 

 

 

 

 

 

 

 

 

 

 

 

a

2

b b t

2c

 

2c t 0,

 

 

 

 

 

 

 

 

 

 

1

2

 

1

 

2

 

 

 

 

 

 

 

 

4a1 4a2t 2b1 2b2t b2 0,

 

 

 

 

 

 

 

2a

 

 

2a

t b

b t c

0.

 

 

 

 

 

 

 

 

1

 

2

 

1

2

 

2

 

 

 

 

 

Приравнивая коэффициенты при t0 и t1 в левой и правой частях каждого уравнения, приходим к системе

 

 

t0 : a

2

b

 

2c

0, t1 :

b

2c 0,

 

 

 

 

 

 

 

1

 

1

 

 

2

2

 

 

 

 

 

 

4a1 2b1 b2 0,

4a2 2b2 0,

 

 

 

 

 

2a1 b1 c2 0.

 

2a2 b2 0.

 

 

 

 

Поскольку кратность корня

k 1

равна 2 ,

то произволь-

ными будем считать два коэффициента a1

и a2 в решении (3.11).

Коэффициенты b1 , b2 ,

c1 ,

c2

выразим через a1 и a2 . Из последней

системы находим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

b1 2a1 a2 , b2 2a2 , c1 a1 a2 , c2 a2 .

 

Тогда решение (3.11) принимает вид

 

 

 

 

 

 

x a a

t e t , y

2a a

2

2a t e t ,

z a a

2

a

t e t .

1

2

 

 

 

 

1

 

2

 

1

 

2

 

В качестве линейно независимых частных решений, соответствующих характеристическим числам k2 k3 1, можно взять

72

3. Линейныесистемы с постоянными коэффициентами

 

1

 

 

 

 

 

te t

 

 

 

 

 

 

0

 

 

1

 

 

 

 

 

t

 

 

 

 

 

 

 

 

t

 

 

 

 

 

 

 

 

t

.

X2

2 e

 

, X3 1 2t e

 

 

 

1

t

2

e

 

 

 

 

 

 

 

1 t e

t

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

1

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3) Общее решение системы , согласно (3.5), имеет вид

 

 

x

0

 

 

 

1

 

 

 

 

 

0

1

 

 

 

 

 

t

 

 

 

t

C3

 

 

 

 

 

 

t

.

 

y C1 2 e

 

C2

2 e

 

 

1 t

2 e

 

 

z

1

 

 

 

1

 

 

 

 

 

 

1

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

 

 

 

 

 

 

 

 

 

 

 

x

2x y z,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x 2 y z,

 

 

 

 

 

 

 

 

 

 

 

 

 

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x y 2z.

 

 

 

 

 

 

 

 

 

 

 

 

 

z

 

 

 

 

 

 

 

 

Решение. Характеристическое уравнение имеет вид

 

 

 

 

 

 

2 k

 

1

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

2 k

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0 ,

 

 

 

 

 

 

 

 

 

 

1

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 k

 

 

 

 

 

 

 

 

2 k 3 1 1 1 2 k 1 2 k 1 2 k 0 , k3 6k 2 9k 4 0 .

Подбором находим корень k1 1; тогда характеристическое уравнение принимает вид k 1 k2 5k 4 0 или k 1 2 k 4 0 .

Таким образом,

характеристические

числа равны k1 k2 1,

k3 4 .

 

 

1) Для корня

k 1 кратности 2

определим число линейно

независимых собственных векторов. При k 1 система (3.7) принимает вид

73

Нежельская Л. А. Линейные дифференциальныеуравнения

 

 

2

 

3

0,

1

 

 

 

1 2 3 0,

 

2

3 0.

1

Размерность матрицы последней системы n 3 , ранг r 1. Число линейно независимых собственных векторов m n r 2 и совпадает с кратностью корня. Так как 1 2 3 , то в качестве

собственных векторов A 1 и A 2 можно взять

 

 

1

 

 

1

 

 

1 1

,

 

2

0 .

A

A

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

1

Два линейно независимых частных решения, соответствующих корню k 1 кратности 2 , запишутся в виде

 

 

 

 

 

 

 

1

 

 

 

 

 

1

 

 

 

 

 

X

 

 

1 et , X

2

 

 

0 et .

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

1

 

 

 

2) Для простого корня k3 4

 

находим собственный вектор

 

 

3 , решая систему

 

 

 

 

 

 

 

 

 

 

 

 

A

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

2

 

3

0,

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

1

2 2 3 0,

 

 

 

 

 

 

 

 

 

2 2 3 0.

 

 

 

 

 

 

 

1

 

 

Из последней системы находим

 

1 3

и

2 3 . Значит, соб-

ственный вектор

 

3

и частное решение

X3 , соответствующие

A

простому корню k3 4 , принимают вид

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

1

 

 

 

 

 

 

 

 

 

3

 

 

X3

 

 

4t

 

 

 

 

 

A

 

.

 

 

 

 

 

1 ,

1 e

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

1

 

 

74

3.Линейныесистемы с постоянными коэффициентами

3)Общее решение исходной системы, согласно (3.5), запишется

ввиде

x

 

1

 

 

1

 

1

 

 

 

 

 

 

 

t

 

 

t

 

4t

.

y

C1

1 e

 

C2

0 e

 

C3 1 e

 

 

 

 

0

 

 

 

 

 

 

 

 

z

 

 

 

 

1

 

1

 

 

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

x 3x y 2z,y 5y 3z,

z x 3y.

Решение. Составим и решим характеристическое уравнение

3 k

1

2

 

 

k 3

k 5 k 3 2 5 k 9 3 k 0,

 

0

5 k

3

 

0 ,

 

 

k

3

 

8k

2

22k 20

0.

 

1

3

k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Подбором находим корень

k1 2 ; при этом характеристическое

уравнение принимает

вид

k 2 k 2 6k 10 0 и уравнение

k2 6k 10 0 имеет

пару простых комплексных сопряжённых

корней k2,3 3 i .

Таким образом, все корни характеристического уравнения различные. Построим фундаментальную систему решений.

1) Строим собственный вектор, соответствующий корню k1 2 . Решая систему

 

2

 

2

3

0,

1

 

 

 

3 2 3 3

0,

 

 

 

 

2 3 0,

1 3 2

находим 1 3 , 2 3 . Тогда собственный вектор A 1 и частное решение X1 , соответствующие простому корню k1 2 , имеют вид

75

Нежельская Л. А. Линейные дифференциальныеуравнения

 

 

 

 

 

 

 

 

1

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 1

, X1 1 e2t .

 

 

 

 

 

 

 

 

 

 

 

A

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

1

 

 

 

 

 

 

 

 

2) Подставляя корень k2 3 i

в систему (3.7), находим

 

 

 

 

 

 

 

 

 

 

2 3 0,

 

 

 

 

 

 

 

 

 

 

 

 

i 1 2

 

 

 

 

 

 

 

 

 

 

 

 

2 i 2 3 3 0,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3 i 3 0.

 

 

 

 

 

 

 

 

 

 

1 3 2

 

 

 

 

 

 

Приведём матрицу системы к ступенчатому виду

 

 

 

 

1 3

3 i

 

 

1

3

3 i

1

3

3 i

 

 

 

 

 

 

 

3

 

 

0

2 i

3

 

 

0 2

i

3

 

 

0 1

 

 

.

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

 

i

1

2

 

 

 

 

0 1 3i

3

 

 

 

 

 

 

 

 

 

 

 

 

 

3i

0

0

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Выпишем связь между элементами собственного вектора:

2 23 i 3 , 1 3 2 3 i 3 .

Тогда собственный вектор можно задать в виде

 

 

 

 

 

 

2 i

 

 

 

2

 

 

 

 

A

3

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 i

Комплексное решение, соответствующее корню k2 3 i , есть

 

 

2 i

 

 

2 i

 

 

 

 

 

2 i

 

 

 

 

 

 

 

 

3

 

 

3 i t

 

 

 

3t

e

it

 

3

 

3t

cost i sin t

 

 

 

 

X

 

e

 

 

 

3 e

 

 

 

e

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

2 i

 

 

 

 

 

 

 

 

i

 

 

2

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2cost sin t

 

 

 

cost 2sin t

 

 

 

 

 

 

 

 

 

3cost

 

 

 

 

 

3sin t

 

 

3t

.

 

 

 

 

 

 

 

i

 

 

e

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2cost sin t

 

 

cost 2sin t

 

 

76

3. Линейныесистемы с постоянными коэффициентами

Взяв действительную и мнимую части от комплексного реше-

ния X , получим два линейно независимых действительных решения

 

 

 

2cost sin t

 

 

 

 

 

cost 2sin t

 

 

 

 

 

 

3cost

 

3t

,

 

 

 

3sin t

 

3t

.

 

 

 

X2 Re X

e

 

X3 Im X

e

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2cost sin t

 

 

 

 

 

cost 2sin t

 

 

3) Общее решение системы – линейная комбинация векторфункций Xi , i 1,3 :

x

1

 

 

2cost sin t

 

 

cos t 2sin t

 

 

 

 

 

 

2t

 

3cost

 

3t

 

3sin t

 

3t

.

y

C1 1 e

 

C2

e

 

C3

e

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z

1

 

 

2cost sin t

 

 

cost 2sin t

 

 

3.2.2. Матричный метод

Линейную однородную систему всегда можно проинтегрировать матричным методом. Это значит, что всегда можно найти сразу всю фундаментальную матрицу системы.

Пусть дана система

dxi

aji x j t ,

i 1, n ,

(3.12)

dt

n

 

 

 

 

j 1

 

 

 

 

где все a ji – постоянные действительные числа. Введём в рассмотрение две матрицы

 

x11 t

x12 t ... x1n t

 

a

 

a

 

...a

 

 

 

x

 

 

 

 

t

11

12

 

1n

 

 

X

t

x

t ... x

 

a

 

a

 

 

...a

 

,

21

 

22

 

2n

 

 

, A

21

 

22

 

2n

 

.................................

 

..................

 

 

 

 

t

 

 

 

 

 

 

 

a

 

 

...a

 

 

 

 

xn2

 

 

 

 

a

n1

n2

 

 

 

xn1

t ... xnn t

 

 

nn

 

где X – фундаментальная матрица системы (3.12), A – транспонированная матрица коэффициентов системы (3.12).

77

Нежельская Л. А. Линейные дифференциальныеуравнения

Тогда систему (3.12) можно записать в виде одного равносильного матричного уравнения

 

 

 

dX

XA .

(3.13)

 

 

 

dt

 

 

 

 

 

В (3.13) матрица

dX

называется производной от матрицы

X и

dt

 

 

 

 

 

находится дифференцированием всех элементов матрицы X . Всякая невырожденная матрица (с определителем, отличным

от нуля), обращающая уравнение (3.13) в тождество, называется

решением этого уравнения или интегральной матрицей.

Напомним, что показательной функцией eA

матрицы A назы-

вается сумма ряда

 

 

 

 

 

 

 

eA E

A

 

A2

 

A3

... ,

(3.14)

 

1!

 

2!

 

3!

 

 

где E – единичная матрица той же размерности, что и матрица A . Ряд сходится для любой матрицы A .

Интегральной матрицей уравнения (3.13) является матрица ви-

да

X eAt , X 0 E .

(3.15)

Каждый элемент интегральной матрицы (3.15) представляет собой степенной ряд по степеням t , сходящийся при всех t . Все эти ряды можно просуммировать, так что интегральная матрица (3.15) будет состоять из элементарных функций.

Рассмотрим возможные случаи матрицы системы

1) Матрица A каноническая, при этом чисто диагональная

A k1, k2 ,...,kn ,

так что всем характеристическим числам k1,k2 ,...,kn матрицы A соответствуют простые элементарные делители

k k1, k k2 ,...,k kn ,

78

3. Линейныесистемы с постоянными коэффициентами

тогда как среди характеристических чисел k1,k2 ,...,kn могут быть и равные между собой. Пользуясь формулой (3.14), находим

eAt e k1,k2 ,...,kn t e k1t,k2t,...,knt ek1t ,ek2t ,...,eknt .

 

 

 

В этом случае интегральная матрица (3.15) принимает вид

X ek1t ,ek2t

,...,eknt .

(3.16)

 

 

 

 

 

Матрица X чисто диагональная,

с функциями

ekit , i

 

, на

1, n

главной диагонали.

Замечание. Матрица A чисто диагональная, если:

а) все собственные значения ki простые, т.е. имеют кратность

i 1;

б) среди собственных значений ki имеются кратные, но каждо-

му кратному корню, кратность которого

i 1,

соответствует

i

одномерных жордановых клеток.

 

 

 

 

 

 

 

 

 

 

 

 

2) Матрица A каноническая, но квазидиагональная

 

 

A

I

1

k

,I

2

k

2

,...,I

s

k

s

,

 

2

...

s

n .

 

 

 

 

1

 

 

 

 

 

 

 

1

 

 

 

 

 

 

Здесь I i

ki

 

– жордановая клетка, есть матрица размерности

i ,

имеющая вид

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k

 

 

0

 

0 ...

0

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

ki

 

0 ...

0

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

I i

ki

0

 

1

ki ...

0

0

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

..........................

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

0

0 ...

1 k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

 

 

 

 

 

 

 

Заметим,

что

 

I1 ki ki .

 

Матрица

 

I i

ki

 

соответствует элемен-

тарному делителю

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k ki i , i 1, s .

79

Нежельская Л. А. Линейные дифференциальныеуравнения

Среди характеристических чисел k1,k2 ,...,kn могут быть и равные

между собой.

В этом случае, пользуясь формулой (3.14), находим интегральную матрицу

X eAt

I

1

k

,I

2

k

 

,...,I

s

k

 

t

 

 

I

1

k

t,I

2

k

 

t,...,I

s

k

 

t

 

e

1

 

 

 

2

 

 

 

 

s

 

e

1

 

 

 

 

2

 

 

 

 

s

 

 

 

 

 

 

 

 

I 1 k1 t

,e

I

2 k2 t

,...,e

I s ks t

,

 

 

 

 

 

 

 

 

(3.17)

 

 

 

e

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ekit

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

...

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

tekit

 

 

 

 

 

 

 

 

 

ekit

 

 

 

 

 

 

...

0

 

 

 

 

 

 

 

 

eI i ki t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

..................................................

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t i 1

 

 

 

k t

 

 

t i 2

 

 

 

 

k t

 

 

 

 

 

k t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

 

 

 

 

 

 

 

 

 

 

i

 

 

 

 

 

 

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 !e

 

 

 

 

2

!e

 

 

...

e

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3)Матрица A произвольного вида.

Вэтом случае её всегда можно привести к каноническому виду при помощи преобразования подобия, т.е. можно построить такую

невырожденную матрицу S det S 0 , чтобы матрица

 

B SAS 1

(3.18)

была канонической.

В матричном уравнении (3.13) выполним подстановку

X ZS , (3.19)

где Z – новая неизвестная интегральная матрица (решение нового дифференциального матричного уравнения). Получаем

ddtΖS ZSA , ddtZ ZSAS 1

или

dZ

ZB .

(3.20)

dt

 

 

80

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