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

Лекции по обыкновенным дифференциальным уравнениям Ржавинская Белякова Жаркова

.pdf
Скачиваний:
176
Добавлен:
27.03.2016
Размер:
906.51 Кб
Скачать

æç y1(t)

Пусть Y (t) = ç ...

çè yn (t)

ö

÷÷ - решение системы (4.6), функции fi (t, y1(t),..., yn (t)), i =1,..., n имеют

÷

ø

непрерывные частные производные по всем аргументам до порядка n – 1 включительно. Подставляем Y(t) в (4.6), получим n тождеств:

ì dy1

= f (t, y (t),..., y

n

(t))

ï

dt

1

1

 

 

ï

 

 

...

 

 

 

 

í

dyn

 

 

 

 

 

 

ï

= f

n

(t, y (t),..., y

n

(t))

 

ï

dt

 

1

 

 

î

 

 

 

 

 

 

 

в частности,

 

 

 

 

 

 

 

 

dy1 = f

(t, y (t),..., y

n

(t)) .

(4.12)

 

 

 

 

 

 

 

 

dt

 

1

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Продифференцируем (4.12) по t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d

2 y

 

f

 

 

f dy

 

 

 

 

 

 

 

 

f

 

 

dy

n

 

 

 

 

 

 

 

1

=

1

+

1

 

1

+ ...+

 

 

1

 

 

 

 

=

 

 

 

 

dt2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

y1 dt

 

 

 

 

 

 

 

 

yn dt

 

 

 

 

 

f1

 

n

f1

 

dyk

 

 

 

f1

 

 

 

 

 

n

f1

 

 

 

 

=

+

å

 

 

=

+ å

fk

 

 

yk dt

 

 

 

 

 

 

 

 

t

 

k=1

 

 

 

 

t

 

 

 

 

 

k =1

yk

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

или

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d 2 y

 

f

 

n

 

f

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

=

 

1 + å

 

1

 

 

 

fk .

 

 

 

 

(4.13)

 

 

 

 

 

 

 

 

dt

2

 

yk

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

k =1

 

 

 

 

 

 

 

 

 

 

 

 

Обозначим правую часть в (4.13) через F2(t, y1, …, yn):

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f

 

 

 

 

n

 

 

f

 

 

 

 

 

 

 

 

 

 

 

F2 (t, y1,K, yn ) =

 

 

 

1

+

å

 

1

fk .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

k =1

y

k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Дифференцируем (4.13) по t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d3 y

 

F

 

 

n

 

F dy

 

 

F

 

 

 

n

F

 

 

1

 

=

 

 

2

+ å

2

 

 

 

k

 

=

 

 

 

 

2

+

å

 

2 fk .

 

dt

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

 

k =1

yk dt

 

 

 

 

t

 

k =1

yk

83

PDF created with pdfFactory Pro trial version www.pdffactory.com

Обозначим

F3(t, y1,K, yn ) =

F

n

F

fk

2

+ å

y

2

 

t

k =1

k

 

 

 

 

 

и т.д. n – 1 раз.

Получили n – 1 тождество

ì

 

 

dy1 = f (t, y (t),..., y

n

(t))

ï

 

 

dt

1

1

 

 

ï

 

 

 

 

 

 

 

 

d 2 y

 

 

 

 

 

ï

 

 

1

= F (t, y (t),..., y

n

(t))

 

 

 

í

 

dt2

2

1

 

 

 

 

...

 

 

 

ï

 

 

 

 

 

 

 

ïd n−1y

 

 

 

 

 

ï

 

 

1

= Fn−1

(t, y1(t),..., yn (t))

 

 

n−1

ï

dt

î

 

 

 

 

 

 

и еще одно уравнение

d n y1 = Fn (t, y1(t),..., yn (t)) . dtn

(4.14)

(4.15)

Систему (4.14) можно рассматривать как систему, неявно задающую функции y2(t), y3(t), …, yn(t). Предположив, что в рассматриваемой области якобиан

D( f1, F2, F3,..., Fn−1) ¹ 0 , систему (4.14) можем разрешить относительно y2, y3, …, yn, выра-

D( y2, y3,..., yn )

 

dy

 

d n−1y

 

зив их через переменные t, y ,

1

,...,

1

. Подставив эти выражения в уравнение (4.15),

dtn−1

1

dt

 

 

получим дифференциальное уравнение порядка n относительно y1(t):

d n y

 

dy

 

d n−1y

1

= F(t, y ,

1

,...,

1

) .

dtn

 

1

dt

 

dtn−1

Справедливо следующее утверждение, которое мы приводим без доказательства.

Теорема 4.2. Пусть дана система дифференциальных уравнений

dyi = f

(t, y ,..., y

n

), i =1,...,n .

(4.16)

dt

i

1

 

 

 

 

 

 

 

Функции fi (t, y1,..., yn ) имеют непрерывные частные производные по всем аргументам в некоторой области

G :t0 - a £ t £ t0 + a; yi0 - bi £ y £ yi0 + bi , i =1,..., n ;

J = D( f1, F2 , F3,..., Fn−1) ¹ 0 в области G.

D(y2 , y3,..., yn )

Тогда следствием системы (4.16) является уравнение

84

PDF created with pdfFactory Pro trial version www.pdffactory.com

d n y

 

dy

 

d n−1y

 

1

= F(t, y ,

1

,...,

1

), i =1,..., n .

(4.17)

dtn

 

1

dt

 

dtn−1

 

Обратное утверждение. По любому решению y1(t) уравнения (4.17) можно подобрать

æç y1(t) n – 1 функцию y2(t), …, yn(t), такую, что вектор-функция Y (t) = ç ...

çè yn (t)

системы (4.16).

(Доказательство см., например, [1, гл. 3, § 2].)

ö

÷÷ является решением

÷

ø

ì dx

= x - y .

В примере 4.2 íï dt

ï

dy

= y - 4x

î dt

 

В обозначениях теоремы 4.2: f1 = x y, f2 = y – 4x и J = f2 = -4 ¹ 0

Þ

систему (4.7)

 

 

 

x

Теор. 4.2

 

можно

свести к одному уравнению второго порядка относительно y(t)

(аналогично

~

=

f

= -1 ¹ 0 и систему (4.7) можно свести к одному уравнению относительно x(t)).

J

1

 

 

y

 

 

 

4.3. Системы линейных дифференциальных уравнений

Определение 4.2. Система дифференциальных уравнений называется линейной, если уравнения системы линейны относительно всех компонент неизвестной вектор- функции и их производных, т.е. имеют вид

dyi

 

n

 

 

+ f (t), i =1,...,n .

(4.18)

=

å

a (t)y

 

dt

 

i j

j

i

 

 

 

 

 

j=1

Обозначим

æ y (t) ö

 

ç

1

÷

,

Y (t) = ç

L

÷

ç

 

÷

 

è yn (t)ø

 

dY dt

çdy1(t)

çdt

=ç L

çdyn (t)

çè dtæ

ö

æ

f (t) ö

 

÷

 

÷÷ ,

ç

1

÷

; A =||ai j (t)||, i, j =1,K,n .

F = ç

L

÷

÷

ç

 

÷

 

÷

è

fn (t)ø

 

ø

 

 

 

 

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

dY

= AY + F .

(4.19)

dt

 

 

85

PDF created with pdfFactory Pro trial version www.pdffactory.com

 

æ

0 ö

 

Определение 4.3. Если в системе (4.19)

ç

 

÷

, система (4.19) называется сис-

F = çL÷

 

ç

0

÷

 

 

è

ø

 

темой линейных однородных дифференциальных уравнений.

Рассмотрим сначала систему линейных однородных дифференциальных уравнений.

 

 

 

 

 

dY = AY .

(4.20)

 

 

 

 

 

dt

 

Пусть C(1)

- совокупность всех вектор-функций, компоненты которых непрерывны

[a, b]

 

 

 

 

 

 

вместе со своими производными первого порядка на отрезке [a, b], C[a,b] -

совокупность

всех вектор-функций, компоненты которых непрерывны на отрезке [a, b].

 

Обозначим L[Y ] = dY - AY , L :C(1)

]

® C

[a,b]

- дифференциальный оператор.

 

 

[a,b

 

 

 

dt

Утверждение. L - линейный оператор.

В самом деле, пусть

 

 

 

æ

(1)

(1)

 

ç

(t) =

ç

Y1(t)ÎC[a,b], Y2

(t)ÎC[a,b] , Y1

ç

 

 

 

ç

 

 

 

è

y1(1) (t)ö÷

L ÷ , Y2

÷

yn(1) (t)÷ø

æç y1(2) (t

(t) = ç L

ç

çè yn(2) (t

)ö÷ ÷ .

÷

)÷ø

Имеем

éêæç y1(1) (t)

L[Y1 +Y2 ] = Lêç êç

êç y(1) (t)

ëè n

æ d(y(1) + y(2) )ö

 

ç

1

1

÷

æ

 

dt

ç

 

÷

ç

= ç

 

...

÷

- Aç

ç d(yn(1) + yn(2) )÷

çç

ç

 

 

÷

è

 

 

ç

 

dt

÷

è

 

ø

 

+ y1(2) (t)ö÷ùú

... ÷ú = ÷ú

+ y(2) (t)÷ú n øû

(1)

 

(2)

ö

(y1

+ y1

)÷

 

 

 

÷

(1)

...

(2)

÷ .

(yn

+ yn

÷

)ø

Привлекая свойства производных, а также свойство дистрибутивности умножения матриц, приходим к равенству:

 

 

æ dy(1)

ö

æ dy(2)

ö

 

 

 

 

 

ç

1

÷

ç

1

÷

 

 

 

L[Y1

+Y2

ç

dt

÷

ç

dt

÷

- AY1

- AY2

 

] = ç ...

÷

+ ç ...

÷

=

 

 

ç

(1)

÷

ç

(2)

÷

 

 

 

 

 

ç dyn

÷

ç dyn

÷

 

 

 

 

 

ç

dt

÷

ç

dt

÷

 

 

 

 

 

è

ø

è

ø

 

 

 

= L[Y1] + L[Y2 ] .

Пусть λ - произвольное действительное число, Y (t)ÎC[(a1), b] .

86

PDF created with pdfFactory Pro trial version www.pdffactory.com

Рассмотрим LY ] .

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

 

 

 

æ d(ly ) ö

 

 

 

 

æ ly (t) ö

ç

 

1

÷

æ ly

ö

 

dt

 

ç

1

÷

ç

 

÷

ç

1

÷

=

L[lY ] = Lç ...

÷

= ç

 

...

÷

- Aç ...

÷

ç

 

÷

ç d(lyn ) ÷

ç

 

÷

 

èlyn (t)ø

ç

 

 

 

÷

èlyn ø

 

 

dt

 

 

 

 

 

è

 

ø

 

 

 

 

æ dy (t) ö

 

 

 

ç

 

1

÷

æ y

ö

 

dt

 

ç

 

÷

ç 1

÷

= l L[Y ].

= lç ...

÷

- l Aç ...

÷

ç dyn (t) ÷

ç

÷

 

ç

 

 

 

÷

è yn ø

 

 

dt

 

 

è

 

ø

 

 

 

Откуда заключаем, что L - линейный оператор.

Используя введенные обозначения, систему (4.20) запишем в следующем виде:

L[Y ] = 0.

(4.21)

Теорема 4.3. Пусть Y1(t), Y2(t), …, Yk (t) - решения системы линейных однородных дифференциальных уравнений L[Y ] = 0. Тогда при любых λ1, λ2, …, λk вектор-функция Y(t) = λ1Y1(t) + λ2Y2(t) + … + λkYk(t) является решением системы L[Y ] = 0.

Доказательство. Рассмотрим L[Y ] = L1Y1(t) + λ2Y2(t) +…+ λkYk(t)].

L[Y]

=

λ1L[Y1(t)] + λ2L[Y2(t)] +…+ λk L[Yk(t)].

 

L − линейный

 

Далее, так как Y1(t), Y2(t), …, Yk(t) - решения системы (4.21), L[Y1] = 0, L[Y2] = 0,…, L[Yk] = 0, и, таким образом, L[Y ] = 0 Þ Y(t) = λ1Y1(t) + + λ2Y2(t) +…+ λkYk(t) - решение системы

(4.21).

Итак, любая линейная комбинация решений системы линейных однородных дифференциальных уравнений является решением этой системы.

Теорема 4.4. Если система линейных однородных дифференциальных уравнений (4.21) с действительными коэффициентами aij(t) имеет решением комплекснозначную векторфункцию Y(t) = U(t) + i V(t), то U(t) и V(t) - решения системы (4.21).

Доказательство слово в слово повторяет доказательство аналогичного утверждения для одного дифференциального уравнения порядка n и мы его опустим.

Определение 4.4. Система вектор-функций Y1(t), Y2(t), …, Yk(t) называется линейно зависимой на отрезке [a, b], если найдутся такие числа α1,…, αk, не все равные нулю одновременно, что "t Î[a, b]

a1Y1(t) + ... + akYk (t) = 0 .

87

PDF created with pdfFactory Pro trial version www.pdffactory.com

Система вектор-функций Y1(t), Y2(t), …, Yk(t) называется линейно независимой на отрезке [a, b], если из равенства

a1Y1(t) + ...+ akYk (t) = 0

следует, что α1 = … = αk = 0.

Приведем без доказательства следующее утверждение.

Теорема 4.5 (о структуре общего решения системы линейных однородных диффе-

ренциальных уравнений). Пусть дана система линейных однородных дифференциальных уравнений

dY

= A(t)Y ,

(4.22)

dt

 

 

где A = ||ai j (t)||, i, j =1,K,n , функции ai j (t) непрерывны на отрезке [a, b], и пусть Y1(t), Y2(t),

…, Yn(t) - линейно независимые на отрезке [a, b] решения системы (4.22).

Тогда общее решение системы линейных однородных дифференциальных уравнений (4.22) имеет вид

Yобщ = С1Y1(t) +K+ СnYn (t) ,

где С1, …, Сn - произвольные действительные числа.

Пример 4.3.

ì dx

= y

 

ï

dt

.

(4.23)

í

ï

dy

= -x

 

î dt

 

Частные решения:

x1(t) = cos t ; y1(t) = -sin t ; x2 (t) = sin t ; y2 (t) = cost (убеждаемся непо-

средственно подстановкой). В векторной форме

 

æ

cost ö

 

æ sin t ö

Y (t) = ç

÷

; Y (t) = ç

÷ .

1

ç

÷

2

ç

÷

 

è

- sin t ø

 

ècost ø

Вектор-функции Y1(t) и Y2(t) - линейно независимы на любом отрезке [a, b]. В самом деле, пусть λ1 Y1(t) + λ2 Y2(t) = θ, или

ì

l cost + l

2

sin t = 0

(4.24)

í

1

cost = 0.

î

- l sin t + l

2

 

1

 

 

 

88

PDF created with pdfFactory Pro trial version www.pdffactory.com

Определитель системы линейных однородных алгебраических уравнений (4.24)

D =

cos t

sin t

= cos2 t + sin2 t = 1 ¹ 0 система (4.24) имеет единственное решение λ1 = λ2 =

 

- sin t

cost

 

0 и вектор-функции Y1(t) и Y2(t) линейно независимы на любом отрезке [a, b]. Тогда по теореме 4.5 общее решение системы (4.23)

Y

= С

æ

cost ö

+ С

æ sin t ö

ç

÷

ç

÷ .

общ

1

ç

÷

 

2 ç

÷

 

 

è

- sint ø

 

ècost ø

Определение 4.5. Любые n линейно независимых частных решений системы линей-

ных однородных дифференциальных уравнений называются фундаментальной системой решений.

В примере 4.3 вектор-функции Y1

æ

cost ö

и Y2

æ sint ö

- фундаментальная система

ç

÷

ç

÷

(t) = ç

÷

(t) = ç

÷

 

è

- sint ø

 

ècost ø

 

решений системы дифференциальных уравнений (4.23).

4.4. Системы линейных однородных дифференциальных уравнений с постоянными коэффициентами

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

 

 

 

dY

= A×Y ,

(4.25)

 

dt

 

 

 

где A = ||ai j ||, i, j =1,K,n , aij - действительные числа.

Так как постоянные aij непрерывные на любом отрезке [a, b] функции, то применима сформулированная выше теорема 4.5.

Частные решения системы (4.25), составляющие фундаментальную систему решений, ищем в виде

æ a1 ö

= ç ÷ kt . (4.26)

Y (t) ç L÷e

çèan ÷ø

Подставляем (4.26) в систему уравнений (4.25)

ì ka1ekt = a11a1ekt + a12a2ekt +K+ a1nanekt

ï

 

ekt = a a ekt + a

 

a

ekt +K+ a

 

a

ekt

ïka

22

2n

í

2

21 1

2

 

n

 

ï

 

 

 

...

 

 

 

 

 

ïka

ekt = a

a ekt + a

n2

a

ekt +K+ a

nn

a

ekt

î

n

 

n1 1

2

 

n

 

89

PDF created with pdfFactory Pro trial version www.pdffactory.com

или

ì(a - k)a + a a

2

+K+ a

a

n

= 0

 

ï

11

1

 

 

12

 

 

1n

 

 

 

ïa21a1 + (a22 - k2 )a2 +K+ a2nan = 0

(4.27)

í

 

 

 

 

...

 

 

 

 

 

 

 

ï

 

 

 

 

 

 

 

 

 

 

 

 

ïa

a + a

n2

a

2

+K+ (a

nn

- k)a

n

= 0.

 

î

n1 1

 

 

 

 

 

 

 

 

 

Система (4.27) - система n линейных однородных алгебраических уравнений относи-

тельно α1, …, αn. Она имеет нетривиальные решения, если

 

a11 - k

a12

...

a1n

 

 

 

 

 

D =

a21

a22 - k2

...

a2n

= 0 .

(4.28)

...

...

...

...

 

 

 

 

an1

an2

...

ann - kn

 

 

Уравнение (4.28) называется характеристическим уравнением для системы дифференциальных уравнений (4.25), а его корни - характеристическими корнями системы

(4.25).

Случай 1. Характеристические корни действительны и различны, тогда их ровно n. Пусть это числа k1, …, kn. Подставляя ki, i = 1, …, n в (4.27), находим соответствующие

a(i) ,...,a(i) , i =1,...,n .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

æa(i) ö

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ç 1

÷

kit

 

линейно независимы на любом отрезке [a, b].

 

 

Векторы Yi =

ç

÷

, i =1,...,n

 

 

ç ...

÷e

 

 

 

 

ça(i) ÷

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

è n

ø

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Покажем это. Допустим, существуют β1, …, βn, такие, что b1Y1(t) + ... + bnYn (t) = 0 . Тогда

выполняются следующие равенства:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ìb a(1)ek1t + b a(2)ek2t

+...+ b a(n)eknt = 0

 

 

 

 

 

 

ï

1 1

 

 

2

1

 

 

n

1

 

 

 

 

 

 

 

 

 

ïb a(1)ek1t + b

a(2)ek2t

+...+ b

a(n)eknt = 0

(4.29)

 

 

 

 

 

í

1 2

 

 

2

2

 

...

n

2

 

 

 

 

 

 

 

 

ï

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ï

(1)

e

k1t

 

(2)

e

k2t

 

(n)

e

knt

= 0.

 

 

 

 

 

 

îb1an

 

+ b2an

 

+ ...+ bnan

 

 

Так как система функций ek1t ,ek2t ,...,eknt линейно независима на любом отрезке [a, b], из равенств (4.29) следует n2 равенств

ìb a(1)

= 0, b

a(2)

= 0,...,b

a(n)

= 0,

 

ï

1

1

 

2

1

 

n

1

 

 

ïb a(1)

= 0, b

2

a(2)

= 0,...,b

n

a(n)

= 0,

(4.30)

í

1

2

 

2

 

2

 

ï

 

 

 

 

...

 

 

 

 

ï

 

(1)

 

 

(2)

 

 

(n)

= 0.

 

îb1an

= 0, b2an

= 0,...,bnan

 

90

PDF created with pdfFactory Pro trial version www.pdffactory.com

 

Так как Y1

θ, первый столбец равенств в (4.30) означает, что

β1

= 0 (среди a(i) ,...,a(i)

хотя бы одно число ≠ 0), аналогично заключаем, что β2 = … = βn =

 

1

n

 

0. Тогда система вектор-функций Y1, ..., Yn линейно независима на любом отрезке [a, b].

 

 

æa(1)

ö

 

 

æa(1)

ö

 

 

 

 

ç 1

÷

k1t

 

ç 1

÷

knt

 

Вектор-функции

Y1 =

ç

÷

, ..., Yn =

ç

÷

составляют фундаментальную систему

ç ...

÷e

 

ç ...

÷e

 

 

 

ça(1)

÷

 

 

ça(1)

÷

 

 

 

 

è n

ø

 

 

è n

ø

 

 

решений системы дифференциальных уравнений (4.25). Тогда по приведенной теореме 4.5 общее решение системы уравнений (4.25) имеет вид

æa(1)

ö

 

ç 1

÷

k1t

ç

÷

Yобщ = C1ç ...

÷e

 

ça(1)

÷

 

è n

ø

 

 

æa(1)

ö

 

ç 1

÷

+...+ Cn

ç ...

÷eknt .

 

ç

÷

 

ça(1)

÷

 

è n

ø

Пример 4.4. Проинтегрировать систему дифференциальных уравнений

ì

x& = 5x + y

.

(4.31)

í

= 10x + 2y

îy&

 

 

Решение. Записываем соответствующее характеристическое уравнение

 

5 - k

1

 

= 0 .

 

 

 

10

2 - k

 

 

Или

k2 – 7k = 0.

Его корни: k1 = 0, k2 = 7. Частные решения ищем в виде

 

 

 

 

 

 

æ a

ö

0t

æ a

 

ö

æ b

ö

 

7t

 

 

 

 

 

 

 

Y1(t) =

ç 1

÷

ç 1

 

÷

ç 1

÷

 

.

 

 

 

 

 

 

ç

÷e

 

= ç

÷ ,

Y2 (t) = ç

÷e

 

 

 

 

 

 

 

 

 

èa2

ø

 

èa2

ø

èb2

ø

 

 

 

 

 

Подставляем Y1 в систему (4.31)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ì 5a1 + a2 = 0

 

 

 

 

 

 

 

 

 

 

 

 

 

í

 

 

 

 

= 0,

 

 

 

 

 

 

 

 

 

 

 

 

 

î10a1 + 2a2

 

 

 

 

 

 

 

 

 

æ

a

ö

 

 

 

 

 

 

 

æ

 

1

ö

 

 

откуда α2

= –5α1

и Y1

ç

1

÷

 

 

 

 

 

 

 

ç

 

 

÷

 

 

(t) = ç

- 5a

÷

. Пусть α1 = 1, тогда Y1 = ç

- 5÷ .

 

 

 

 

è

1

ø

 

 

 

 

 

 

 

è

 

 

ø

 

 

Аналогично находим Y2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ì- 2b + b

2

= 0

Þ

 

 

 

 

æ b

ö

7t

 

 

 

 

 

1

 

b2 = 2b1 и Y2 (t) =

ç

 

1

÷

.

 

 

 

í10b - 5b

2

= 0

 

 

ç2b

÷e

 

 

 

 

î

1

 

 

 

 

 

 

 

è

 

1

ø

 

 

91

PDF created with pdfFactory Pro trial version www.pdffactory.com

Полагаем β1 = 1, тогда Y2

æ

1

ö

7t

.

ç

 

÷

= ç

2

÷e

 

 

è

ø

 

 

Итак, общее решение системы (4.31)

æ

1 ö

æ

1

ö

7t

.

ç

 

÷

ç

 

÷

Yобщ = C1ç

- 5

÷

+ C2ç

2

÷e

 

è

ø

è

ø

 

 

Случай 2. Среди характеристических корней есть комплексные. Так как коэффициенты в системе уравнений (4.25), а следовательно, и в уравнении (4.28) действительные числа, комплексные корни входят сопряженными парами.

Пусть k1 = α + iβ и k2 = α – iβ являются корнями характеристического уравнения (4.28).

æ a1 ö

В этом случае вместо комплекснозначных вектор-функций Y = ç... ÷e(α+iβ)t и

1 ç ÷ çèan ÷ø

 

 

æ b

ö

 

 

 

~

 

Y

=

ç 1

÷

e(α−iβ)t

 

 

= ReY

ç ...

÷

в фундаментальную систему решений включаем вектор-функции Y

2

 

 

 

 

1

1

 

 

ç

÷

 

 

 

 

 

 

 

èbn

ø

 

 

 

 

 

и

~

 

 

 

~

~

линейно независимы.

Y2

= ImY1 (теорема 4.4). Самостоятельно докажите, что Y1

и Y2

Тогда по теореме 4.5 общее решение системы уравнений (4.25) имеет вид

 

 

 

æa(1)

ö

 

 

æa(n−2)

~

~

+ C3

ç 1

÷

k

t

ç 1

Yобщ = C1Y1

+ C2Y2

ç ...

÷e 3

 

+ ... + Cn ç ...

 

 

 

ç

÷

 

 

ç

 

 

 

èa(n1)

ø

 

 

èa(nn−2)

ö

÷

÷eknt .

÷

ø

Пример 4.5. Проинтегрировать систему дифференциальных уравнений

ìdx

= x + y,

ï

 

 

 

í dt

(4.32)

ï

 

dy

= -2x + 3y.

î dt

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

 

1- k

1

 

= 0

 

 

 

- 2

3 - k

 

 

или

k2 - 4k - 5 = 0 Þ k1,2 = 2 ± i .

Частные решения ищем в виде

æ a

ö

(2+i)t

 

ç 1

÷

.

Y(t) = ç

÷e

 

èa2

ø

 

 

92

PDF created with pdfFactory Pro trial version www.pdffactory.com