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

Сандаков Решение ОДУ

.pdf
Скачиваний:
26
Добавлен:
17.05.2024
Размер:
758 Кб
Скачать
☆

плексных) решений системы (6.5), подчиняющейся условиям:γ2 = γ3 =... = γrm −1 = 0 , γ1 ≠ 0 , т.е. находят максимальное число ли-

нейно независимых вещественных (или комплексных) решений системы

 

 

 

 

 

= 0;

 

 

(A −λm E)

γ1

 

(6.7)

 

 

E)

γ

 

=

γ

.

(A −λ

m

0

 

 

 

 

 

1

 

 

При этом очевидно, что γ0 также отлично от нуля. Если число

вновь найденных линейно независимых решений системы (6.7) вместе с числом линейно независимых решений системы (6.6) окажется меньше rm, то разыскивают максимальное число линейно независимых вещественных (или комплексных) решений системы

(6.5), подчиняющейся условиям γ3 = γ4 =... = γrm −1 = 0 , γ2 ≠ 0 , т.е.

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

(A

(A(A

−λm

−λm

−λm

E)γ

2

= 0;

 

 

 

 

 

 

 

;

(6.8)

E)γ1

= 2γ

2

 

 

 

 

 

 

E)γ0

= γ1.

 

 

При этом, как легко видеть, векторы γ0 и γ1 также отличны от нуля.

Если число вновь найденных линейно независимых решений системы (6.8) вместе с числом линейно независимых решений систем (6.7) и (6.6) окажется меньше rm, то указанный процесс построения линейно независимых решений вида (6.4) продолжают.

Утверждается, что на каком-то k-м шаге (0 ≤ k ≤ rm – 1) мы получим rm линейно независимых решений системы (6.1) вида (6.4), при этом все rm линейно независимых решений (6.4) системы (6.1)

окажутся такими, что γk +1 = γk +2 =... = γrm −1 = 0 . Пример 6.1. Найти общее решение системы

y1′(x) =3y1 − y2 + y3 ;y2′(x) = y1 + y2 + y3 ;

y3′(x) = 4y1 − y2 + 4y3.

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

31

det(A −λE) =

 

3 −λ

−1

1

 

 

 

 

1

1−λ

1

= 0

 

 

4

−1

4 −λ

 

или

det(A −λE) = (3 −λ)(1−λ)(4 −λ) −2(1−λ) =

= (1−λ)(λ2 −7λ+10) = (1−λ)(λ −2)(λ −5) = 0,

т.е. это уравнение имеет три различных вещественных корня λ1 = 1,

λ2 = 2, λ3 = 5.

Координаты собственных векторов h1, h2, h3 , отвечающих соб-

ственным значениям λ1, λ2, λ3, находятся из систем, соответственно,

2h

−h

+ h = 0,

h

−h

+ h

 

= 0,

−2h

−h

 

+ h

 

= 0,

 

11

21

31

 

12

22

32

 

 

 

13

23

33

 

h11

 

+ h31 = 0,

h12

−h22

+ h32

= 0,

h13

−4h23

+ h33

= 0,

4h −h +3h = 0;

4h −h +

2h = 0;

4h −h

−h = 0.

 

11

21

31

 

12 22

 

32

 

13

23

 

33

 

 

Из этих систем видим, что в качестве собственных векторов можно взять

 

h

 

1

 

 

 

h

 

−1

 

 

h

 

1

 

11

 

 

 

 

;

12

 

 

2

 

;

13

 

 

 

h1

= h21

= 1

 

h2

= h22

=

 

h3

= h23

= 1

.

 

h

 

 

−1

 

 

h

 

3

 

 

 

h

 

3

 

31

 

 

 

 

 

 

32

 

 

 

 

 

 

33

 

 

 

Следовательно, общее решение заданной системы дифференциальных уравнений имеет вид

 

y

 

1

 

 

−1

 

1

 

 

 

1

 

 

 

 

x

 

2

 

2 x

 

 

5x

.

y

(x) = y2

= c1 1

e

 

+c2

e

 

+c3 1

e

 

 

y

 

 

−1

 

3

 

 

3

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

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

y1′(x) = 2y1 − y2 ;y2′(x) = y1 + 2y2 .

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

вид

32

det(A −λE) =

 

2 −λ

−1

 

= 0

 

 

 

 

1

2 −λ

 

 

или

det(A −λE) = (2 −λ)2 +1 = 0.

Это уравнение имеет два комплексно-сопряженных корня λ1 = 2 + i,

λ2 = 2 – i.

Найдем собственный вектор h1,

отвечающий собственному зна-

чению λ1 = 2 + i. Вектор h1

является нетривиальным решением си-

−i

−1

h11

 

0

 

 

Очевидно, что в качестве собственного

стемы 1

−i

h

 

= 0 .

 

 

21

 

 

 

 

 

 

 

 

 

 

 

вектора h1

можно взять вектор

h1 = i .

 

 

 

 

 

 

 

 

 

 

 

1

 

Собственному значению λ2 = 2 – i будет соответствовать соб-

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

= h1

= −i . Следовательно, векторфункции

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

i

 

(2+i) x

i

 

2 x

(cosx +isinx),

 

y1(x) = e

 

 

 

= e

 

 

 

 

 

1

 

 

 

1

 

 

 

 

 

 

(x)

 

−i

 

(2−i) x

−i

2 x

(cosx −isinx)

 

y2

=

e

 

 

=

e

 

 

 

 

 

1

 

 

 

 

1

 

 

 

образуют фундаментальную (комплексную) систему решений, а вектор-функции

 

 

−sinx

2 x

;

 

 

cosx

2 x

u

(x) = Re y1

(x) =

e

 

v

(x) = Im y1

(x) =

e

 

 

 

cosx

 

 

 

 

 

sinx

 

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

 

−sinx

2 x

cosx

2 x

.

y(x) = c1

 

e

 

+c2

e

 

 

cosx

 

 

sinx

 

 

33

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

 

 

 

y′(x) = y

− y

;

(6.9)

 

 

1

1

2

 

 

 

y2′(x) = y1 +3y2 .

 

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

 

 

1−λ

−1

 

= 0

или λ2

– 4λ + 4 = 0

 

 

 

 

 

1

3 −λ

 

 

 

 

 

 

имеет кратный корень λ1 = λ2 = 2. Будем искать линейно независимые решения системы в виде

y(x) = (γ0 + γ1x)e2 x ,

(6.10)

где γ0 и γ1 – подлежащие определению постоянные вектор-столбцы. Подставляя решение вида (6.10) в заданную систему, имеем

γ1e2 x +2(γ0 + γ1x)e2 x = A(γ0 + γ1x)e2 x .

Приравняв коэффициенты при одинаковых степенях х, получим для определения вектор-столбцов γ0 ,γ1 систему

 

 

 

= 0,

 

 

(A −2E)γ

 

 

(6.11)

 

1

 

.

(A −2E)γ

0

= γ

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

А. Положим γ1 = 0 , тогда система (6.11) примет вид

 

 

 

 

 

 

= 0

 

 

 

 

 

(6.12)

или

(A −2E)γ0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

−1

−1

γ

10

 

=

0

 

 

 

 

 

 

 

 

 

 

 

 

.

 

 

 

 

1 1

 

γ20

 

 

 

0

 

 

 

 

 

Единственным линейно независимым решением этой системы

 

 

−1

. Следовательно, согласно (6.10)

является вектор-столбец γ0

=

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

−1

 

 

2 x

системы (6.9).

получим ненулевое решение y1

=

 

1

e

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= 0,

 

≠ 0;

 

Б. Рассмотрим систему

(A −2E)

γ

 

 

γ

1

Первое уравне-

 

 

 

 

 

 

1

 

 

 

.

 

 

(A −2E)

γ

0

= γ

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ние системы имеет такой же вид, что и (6.12), следовательно,

34

γ1 = γγ11 = −11 .

21

Подставив γ1 во второе уравнение системы, получим

−1

−1

γ

10

 

−1

 

1

1

 

 

 

=

 

 

 

γ20

1

 

или γ10 + γ20 = 1, т.е. за γ0 можно взять вектор-столбец

γ0 = 0 .1

Таким образом, согласно (6.10) получим следующее ненулевое решение системы (6.9):

 

0

 

−1

 

2 x

−1

 

2 x

.

y2

(x) =

 

+

1

x e

 

=

e

 

 

1

 

 

 

 

1+ x

 

 

 

 

 

 

 

 

 

 

 

Очевидно, что y1(x) и y2 (x) – линейно независимые решения си-

стемы (6.9) и поэтому образуют фундаментальную систему решений. Значит, общее решение системы (6.9) имеет вид

 

−1

2 x

−1

 

2 x

.

 

y(x) = c1

 

 

e

 

+c2

e

 

 

 

 

 

1

 

1+ x

 

 

 

 

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

 

 

 

 

 

 

 

4

−1

−1

 

′

 

 

 

 

 

2

 

 

 

(6.13)

y (x) = Ay(x) , где

A = 1

 

−1 .

 

 

 

 

 

 

−1

 

 

 

 

 

 

 

 

 

1

 

2

 

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

 

det(A −λE) =

 

4 −λ

−1

−1

 

 

 

 

 

 

 

 

1

 

2 −λ

−1

 

= 0

 

 

 

 

1

 

−1

2 −λ

 

 

или

det(A −λE) = (4 −λ)(2 −λ)2 +(2 −λ) = (2 −λ)(λ2 −6λ+9) = = (2 −λ)(λ−3)2 = 0

имеет корни λ1 = 2 и λ2 = λ3 = 3.

35

Рассмотрим простой корень λ1 = 2. Найдем собственный вектор, отвечающий этому собственному значению λ1 = 2. Он является нетривиальным решением системы

 

 

2

 

 

−1

−1 h

 

 

0

 

h11 −h31 = 0;

 

 

 

1

 

 

0

 

 

 

11

 

=

 

0

 

или

 

 

 

 

 

−1

h21

 

 

 

 

 

 

 

 

 

−1 0

 

 

 

 

 

h11 −h21 = 0.

 

 

1

 

 

h

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

31

 

 

 

 

 

 

 

Следовательно, в качестве собственного вектора h1 можно взять

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

. Тогда первое ненулевое решение системы (6.13)

вектор h1

1

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

2 x

.

 

 

 

 

 

 

 

 

имеет вид y1(x)

1 e

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

Рассмотрим корень λ2 = λ3 = 3 кратности два. Найдем собственный вектор, отвечающий этому собственному значению. Он является нетривиальным решением системы

1

−1

−1 h

 

0

 

 

−1

1

 

 

 

или x1 − x2 − x3 = 0 .

1

−1 h2

 

= 0

 

−1

 

 

 

 

 

1

−1 h3

 

0

 

Эта система имеет два линейно независимых решения:

 

 

1

 

 

 

1

 

=

 

 

=

 

 

h2

1

 

и h3

0 .

 

 

0

 

 

1

 

Следовательно, для собственного значения λ = 3 его алгебраическая и геометрическая кратности совпадают. Поэтому согласно вышеизложенному функции

 

1

 

 

 

1

 

 

 

 

1

 

 

 

2 x

,

 

 

3x

,

 

3x

y1

(x) = 1 e

 

y2

(x) = 1

e

 

y3

(x) = 0 e

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

0

 

 

 

1

 

 

образуют ФСР системы (6.13). Тогда общее решение системы (6.13) имеет вид

36

y

 

1

 

 

1

 

1

 

(x) = c

1 e2 x + c

1

 

+c

0

e3x .

o.o

1

 

 

2

 

 

3

 

 

 

 

 

1

 

 

0

 

1

 

 

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

 

 

 

 

4

 

−1

 

0

 

 

 

′

 

 

 

 

1

 

 

 

 

(6.14)

y (x) = Ay(x) , где A = 3

 

 

−1 .

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

1

 

 

1

 

 

 

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

4 −λ

−1

0

 

 

3

1−λ

−1

= 0

1

0

1−λ

 

или (4 −λ)(1−λ)2 +1+3(1−λ) = −λ3 +6λ2 −12λ+8 = (2 −λ)3 = 0 име-

ет корень λ1 = λ2

= λ3

= 2 кратности три.

 

Будем искать линейно независимые решения системы в виде

 

 

y(x) = (γ0 + γ1x + γ2 x2 )e2 x ,

(6.15)

где γ0 , γ1, γ2 – подлежащие определению постоянные вектор-

столбцы. Подставляя решение вида (6.15) в уравнение (6.14) и приравнивая коэффициенты при одинаковых степенях х, получим для

определения вектор-столбцов γ0 , γ1, γ2 систему

 

 

(A −2E)γ

2 = 0;

 

 

 

 

 

 

 

 

 

 

(6.16)

 

 

(A −2E)

γ1 = 2γ2 ;

 

 

 

 

 

 

 

 

 

 

(A −2E)

γ

0 = γ1.

 

 

= 0 , тогда система (6.16) примет вид

А. Положим γ1

= γ2

 

 

 

 

 

= 0

(6.17)

 

 

(A −2E)γ0

 

 

2 −1

 

0

γ10

 

0

или

 

 

 

 

 

 

 

 

 

3 −1

−1

γ20

= 0 .

 

 

1 0

−1

γ

30

 

0

 

 

 

 

 

 

 

 

37

Единственным линейно независимым решением этой системы

является вектор-столбец γ0 = 1 . Следовательно, согласно (6.15)

21

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 x

 

системы (6.14).

получим ненулевое решение y1

= 2 e

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

=

0, тогда система (6.16) примет вид

Б. Положим γ1 ≠ 0,

γ2

 

 

 

 

 

 

 

 

=

0;

 

 

 

 

 

 

(A

−2E)γ

 

 

 

 

 

(6.18)

 

 

 

 

 

1

 

 

 

 

.

 

 

 

 

(A

−2E)γ

0

= γ

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

Согласно (6.17) имеем

γ1 =

 

 

. Подставив

γ1 во второе уравнение

2

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

−1 0

 

γ10

1

 

 

системы (6.18), получим

 

 

 

−1

 

 

 

 

 

 

 

 

 

 

. Единственным

 

3

 

−1

 

γ20

= 2

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

−1

 

γ30

1

 

 

 

 

1

 

 

ненулевым решением этой системы являетсявектор

 

 

.

γ0

= 1

 

 

 

 

 

 

 

0

 

Таким образом, за второе ненулевое решение системы (6.14) вида (6.15) можно взять вектор-функцию:

 

 

1

 

1

 

 

 

 

 

 

 

 

 

 

2 x

.

y2

(x) =

1

 

+ 2 x e

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

1

 

 

 

В. Рассмотрим, наконец, систему

(A −2E)γ

2

= 0;

 

 

 

 

 

 

 

;

(A −2E)γ1

= 2γ

2

 

 

 

 

 

 

(A −2E)γ0

= γ1.

 

Очевидно, что в качестве γ2 можно взять вектор

(6.19)

 

1

γ2

=

2 .

 

 

 

 

1

38

Подставив γ2 во второе уравнение системы (6.19), получим

2 −1

 

0

γ11

 

 

 

2

 

 

 

 

 

 

 

 

 

 

=

 

 

3 −1

−1

γ21

 

4 .

1 0

 

−1

γ

31

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

Отсюда без труда находим

 

 

γ11

 

 

 

 

 

 

 

 

 

 

2

 

 

 

=

γ

 

 

 

= 2 .

 

γ

1

21

 

 

 

 

γ

 

 

 

 

 

 

 

 

 

 

31

 

0

 

 

И, наконец, подставляя вектор γ1

в последнее уравнение системы

(6.19), получим для определения γ0

систему

2 −1

 

0

 

γ10

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

=

 

 

3 −1

−1

 

γ20

 

2 .

1 0

 

−1

 

γ

30

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

γ10

 

 

1

 

Откуда без труда находим

 

=

γ

 

=

0 .

γ

 

 

 

 

0

 

 

γ

20

 

 

 

 

 

 

 

 

 

30

 

1

 

 

Таким образом, согласно (6.15) получим третье ненулевое решение

y3

 

1

2

1

 

(x) =

0

+ 2 x + 2 x2

e2 x .

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

0

 

1

 

Решения y1(x) , y2 (x) , y3 (x) образуют фундаментальную си-

стему решений системы. Поэтому согласно утверждению 6.2 общее решение системы (6.14) имеет вид

 

1

 

 

 

 

 

 

1

 

 

1

 

 

 

 

 

 

1

 

2

1

 

 

 

 

 

 

2 x

 

 

 

 

 

+

 

 

 

2 x

+c3

 

 

 

 

2

 

2 x

.

y

(x) = c1 2 e

 

 

+c2

1

 

2 x e

 

 

0 +

2 x

+ 2 x

 

e

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

0

 

1

 

 

 

 

1

 

0

1

 

 

 

 

или

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y(x) =

c1

+c2

+c3

 

 

c2

+2c3

 

 

c3

 

 

 

 

 

 

 

 

 

2c

+c

 

+

2c

+

2c

x + 2c

x2

e2 x .

 

 

 

 

 

 

 

 

 

 

c 1+c

2

 

 

 

2 c

3

 

 

c 3

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

3

 

 

 

 

2

 

 

 

 

3

 

 

 

 

 

 

 

39

Задачи для самостоятельного решения

Вариант 1:

x′ = 2x −2 а) y′ = −2x +z′ = −2y;

Вариант 2:

y,

y −2z,

x′ = б) y′ =z′ =

2y +2z,

2x +3y − z,

2x − y +3z;

x′ = 4x −5y +7z, в) y′ = x −4y +9z,z′ = −4x +5z.

x′= x −2y,

x′ = 7x + y + z,

 

 

2y −2z,

 

7y + z,

а) y′= −2x +

б) y′ = x +

 

 

3z;

 

 

 

z′= −2y +

z′ = x + y +7z;

Вариант 3:

 

 

 

 

 

x′=3x +

2y,

x′ = x + 2y + 2z,

 

4y −2z,

 

 

 

а) y′ = 2x +

б) y′ = 2x + y + 2z,

 

 

5z;

 

 

 

z′ = −2y +

z′ = 2x + 2y + z;

Вариант 4:

 

 

 

 

 

x′=5x −2y −2z,

x′=17x −2y −2z,

 

 

 

 

 

4z,

а) y′= −2x +

6y, б) y′= −2x +14y −

 

 

 

 

 

 

z′= −2x + 4z;

z′= −2x −4y +14z;

Вариант 5:

 

 

 

 

 

x′= 7x −4y,

x′=8x + 4y − z,

 

 

5y + 4z,

 

7y +

4z,

а) y′= −4x +

б) y′ = 4x −

 

3z;

 

4y +

8z;

z′= 4y +

z′ = −x +

Вариант 6:

 

 

 

 

 

x′= −2y,

 

 

x′= 2x − y + 2z,

 

 

2y −2z,

 

 

3z,

а) y′= −2x +

б) y′=5x −3y +

 

 

5z;

 

 

 

z′= −2y +

z′ = −x −2z;

 

x′= 4x −5y +7z, в) y′= x −4y +9z,z′ = −4x +5z.

x′ = −x −2y + 2z, в) y′ = −2x − y + 2z,z′ = −3x −2y +3z.

x′= −3x + 2y + 2z, в) y′= −3x − y + z,

z′= −x + 2y.

x′=3x −3y + z, в) y′ =3x −2y + 2z,z′ = −x + 2y.

x′ = 2x + y, в) y′ = x +3y − z,

z′ = −x + 2y +3z.

40