Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Учебник 1 курс 2 семестр.pdf
Скачиваний:
91
Добавлен:
16.12.2013
Размер:
1.39 Mб
Скачать

в) характеристическое уравнение λ2 +3λ −4 = 0

имеет

 

корни λ1 =1 , λ2 = −4 , так что

y = C ex +С

2

у4 x – общее решение соответствующего однородного уравнения. Ищем частное реше-

1

 

 

 

 

 

 

 

 

 

 

 

ние исходного уравнения в виде

 

 

 

 

 

 

 

 

 

 

 

 

ϕ(x) = ( Asinx + Bcosx)ex .

 

 

 

 

Имеем:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+ B) cos x)e

x

;

 

 

 

 

 

 

 

ϕ (x) = (( A B) sinx +( A

 

 

 

 

 

 

 

 

′′

 

x

.

 

 

 

 

 

 

 

 

 

ϕ (x) = (2B sin x + 2 Acos x)e

 

 

 

Подставляя в уравнение, приводя подобные и деля на ex , получаем:

 

 

 

 

 

 

 

 

 

(5B A) sin x +(B +5A) cos x) = cos x ,

откуда A =

5

,

B = −

1

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

 

 

 

 

 

 

 

 

 

 

 

 

26

 

26

 

 

 

 

 

 

y = C1ex +C2 e4 x + 261 (5sinx cos x)ex

является общим решением исходного уравнения.

Вопросы и упражнения для самостоятельной работы

1.Найти общие решения уравнений:

а) x2 y′′− xy′ = 3x2 ;

б) (2x +1) y′′+(2x 1) y′−2 y = x2 + x .

2.Найти решения уравнений, удовлетворяющие указанным начальным условиям:

а)

y

′′

4 y

+3y = 0 ,

y(0) = 6 ,

 

=10

;

 

 

y (0)

 

б)

y

′′

2 y

+ y = 0 , y(2) =1 ,

= 3 ;

 

 

 

y (0)

 

в)

y

′′

2 y

+2 y = 0 ,

y(0) = 0 ,

 

=1 ;

 

 

 

y (0)

 

г)

y

′′

2 y

+3y = 0 ,

y(0) =1 ,

 

= 3 .

 

 

y (0)

3.Найти общие решения уравнений:

а)

y′′−4 y′+ 4 y = x2 ;

б) y′′+8 y′=8x ;

в) y′′−2 y′−3y = e4 x ;

г)

y′′+ 4 y′+ 4 y = 8e2 x ;

д)

y′′+ y = 4 sin x ;

е)

y′′+ 2 y′+5y = ex sin 2x .

4.Найти решения уравнений, удовлетворяющие указанным начальным условиям:

а)

y

′′

4 y

+ 4 y = e

2 x

, y(0)

= 2 ,

= 8 ;

 

 

 

 

y (0)

 

б)

y

′′

y

= −5e

x

(sin x +cos x) ,

 

= 5 .

 

 

 

y(0) = −4 , y (0)

5.Изменение количества товара на рынке пропорционально разнице между предложением и спросом: q′= 0,48(s d ) , а изменение цены пропорционально отклонению количества товара от ба-

зового уровня q = 50 : q′ = −0,1(q 50) . Функции спроса и предложения имеют вид: s =10 +20 p , d = 70 10 p . Требуется:

а) найти в общем виде p = p(t) и q = q(t). С каким периодом будет колебаться цена? При каких начальных значениях цены останутся неизменными?

б) найти p = p(t) и q = q(t), если в начальный момент p(0) =1,5 , q(0) = 44 . Найти значения p и q в момент времени t1 = 0,1. Найти моменты времени, когда p и q принимают экстремальные значения; вычислить значения p и q для этих моментов времени.

§5. Уравнения, допускающие понижение порядка

I.Автономные уравнения

Уравнения вида

F( y, y, y′′,..., y(n) ) = 0 ,

не содержащие переменную x , называются автономными и допускают понижение порядка, если принять за новую переменную неизвестную функцию y , а за новую функцию – p( y) = y.

II. Уравнения, не содержащие y(x)

F(x, y, y′′,..., y(n) ) = 0 .

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

p( y) = y.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

′′

 

(k )

, то следует выбрать новую

Если уравнение не содержит также и производные y , y ,...., y

 

функцию в виде p( y) = y(k +1) , при этом порядок уравнения понижается до n k .

 

III.

 

 

 

′′

 

(n)

 

 

 

Уравнения, однородные относительно y, y , y ,...., y

 

 

 

 

Уравнение вида

 

 

 

 

 

 

 

 

 

 

(n)

)

= 0

 

 

 

 

 

 

 

F(x, y, y ,..., y

 

 

 

 

 

 

называется

′ ′′

(n)

,

если

 

 

 

′ ′′

на

однородным относительно y, y , y ,...., y

 

одновременная замена y, y , y ...

′′

не меняет уравнения.

 

 

 

 

 

 

 

 

ky, ky , ky ...

 

 

 

 

 

 

 

 

Порядок такого уравнения понижается с помощью введения новой неизвестной функции

z :

y′ = yz .

 

 

 

 

 

 

 

 

 

 

IV. Уравнения, приводимые к дивергентному типу

Уравнение вида

dxd (F(x, y, y,..., y(n) )= 0

называетсяуравнениемдивергентноготипа, интегрированиекоторого

F(x, y, y,..., y(n) ) = C

приводит к понижению порядка.

П ри ме р ы

1.Решить уравнение yy′′+ (y)2 yy′ = 0.

 

 

 

Ре ш е н и е . Данное уравнение не содержит x , следовательно, является автономным. Соглас-

но

п. I,

сделаем

замену

переменной на y , а

неизвестной функции на p = y,

при этом

y

′′

=

d (y

)

dp

dp dy

 

по переменной y , а не x ).

 

 

dx

 

= dx = dy dx = p p (здесь указана производная p

 

 

 

 

 

Получим уравнение пониженной степени

 

 

 

 

 

 

 

 

 

 

 

ypp′+ p2 yp = 0 .

 

 

 

 

 

Разделим уравнение на yp :

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

p′+

 

p =1,

 

 

 

 

 

 

 

 

 

 

y

 

 

 

при этом необходимо проверить удовлетворяет ли

y = 0 и

p = y′ = 0 исходному уравнению (иначе

мы можем "потерять" решения). Подстановка y = 0

в первоначальное уравнение приводит к верному

тождеству, следовательно,

y(x) = 0 – частное решение; аналогично устанавливается, что

y′ = 0 так-

же удовлетворяет уравнению, т.е. y = C , где C произвольная постоянная, является еще одним частным решением уравнения (заметим, что при C = 0 y(x) = 0 ).

Полученное уравнение p'+ 1y p =1 было уже рассмотрено нами в § 3 пример 2 (замечание 2); его

решение имеет вид p = y ln y +Cy , где C произвольная постоянная. Возвращаясь к исходной переменной x и функции y , имеем

y′ = y ln y +Cy .

Действуядалее методом разделения переменных, получаем ответ y(x) = C1 exp(C2ex ) ,

где C1 и C2 произвольные постоянные.

В итоге общее решение уравнения записывается в виде

y(x) = C1 exp(C2ex ) ; y(x) = C3 .

2.Решить уравнение y(4) + y′′′ = x.

Ре ш е н и е . Данное уравнение не содержит y и производных y, y′′, следовательно, производим замену функции p = y′′′ (согласно п. II). Получаем линейное неоднородное уравнение первого порядка с постоянными коэффициентами

p′+ p = x ,

решение которого имеет вид p = Cex 1 + x , где C произвольная константа (убедитесь в этом).

Возвращаясь к первоначальной функции, получим

y′′′ = Cex 1 + x ,

трехкратное интегрирование которого приводит к окончательному ответу

y = C1ex + x4 x3 +C2 x2 +C3x +C4 , 4! 3!

где Ck произвольные постоянные ( k =1,4 ).

3.Решить уравнение yy′′+ 1 x2x2 (y)2 = 0.

Решение. Проверим однородность уравнения по y, yи y′′:

 

 

 

 

 

1 x2

 

2

 

 

 

 

2

 

 

 

 

1 x2

 

 

2

 

 

 

 

 

 

 

 

 

 

ky(ky )+

 

 

 

)

= k

 

 

 

 

′′

+

 

 

 

 

 

 

)

 

= 0.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

2 (ky

 

 

yy

 

 

 

 

x

2 (y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Следовательно, замена y ky, y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

и

 

ky ... не меняет уравнение, т.е. уравнение однородно по y, y

 

y

′′

. Применим согласно п. III замену

y

= yz ,

при этом

 

y

′′

 

 

 

 

= yz

2

+ yz

. Подставляя новые

 

 

 

 

= y z + yz

 

 

 

выражения для yи y′′ в уравнение, получим:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y2 z

+

1

 

z2

= 0.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Данное уравнение выполняется либо при y(x) = 0 , либо при

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z′+

 

1

 

z2 = 0 ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

решениекоторого легко получить методом разделения переменных:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z

=

 

 

 

x

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Cx 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где C произвольная константа.

Вернемся к функции

y :

 

y′ = yz

и после подстановки найденной

функции z получим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

yy= Cxx1 .

Данное уравнение также легко интегрируется. В итоге записываем окончательный ответ

 

 

 

x

 

1

 

 

 

 

 

 

 

y(x) = C

exp

 

+

 

ln

C

 

x 1

 

,

 

 

 

 

 

1

 

 

 

2

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

C2

 

C2

 

 

 

 

 

 

где C1

и C2 произвольные постоянные. Заметим,

что найденное ранее частное решение y(x) = 0

включено в ответ при C1 = C2 = 0 .

 

 

 

 

 

 

 

 

 

 

4.

Решить уравнение y′′ = (y)2 .

 

 

 

 

 

 

 

 

 

 

Ре ш е н и е . К данному уравнению можно применить методики пп. I, II, т.к. уравнение является автономным и не содержит функцию y . Мы же решим данное уравнение IV способом, т.е. при-

ведем его к дивергентному виду. Для этого разделим уравнение на y′ ≠ 0 , предварительно убедившись, что y′ = 0 y(x) = C является частным решением. В итоге получим

yy= y,

которое оказывается дивергентного типа, т.к. допускает следующую запись:

(ln y)= (y).

Проинтегрируем полученное уравнение:

ln y′ = y +C,

и после тождественных преобразований и переобозначений произвольной постоянной C получим y′ = Ce y , интегрирование которого методом разделения переменных приводит к ответу

y(x) = −ln C1x +C2 ,

где C1 и C2 произвольные постоянные. В окончательный ответ следует включить все найденные серии частных решений:

y(x) = −ln C1x +C2 , y(x) = C3 .

Вопросы и упражнения для самостоятельной работы

1.Верно ли, что при параллельном переносе вдоль оси Ox интегральная кривая автономного уравнения переходит в интегральную кривую того же уравнения? Ответ поясните.

2.Всегда ли однородное по y, y,..., y(n) уравнение имеет частное решение y(x) = 0 ? Поче-

му?

3.

Каким образом

и

насколько можно понизить порядок уравнения вида

 

dF( y(m) , y(m+1) ,..., y(n) )

= y(k ) ?

Рассмотрите случаи

 

dx

 

 

 

 

а) k < m +1 < n +1 ;

 

б) m +1 < k < n +1 .

4.Решить автономные уравнения:

а) yy

′′

2

3

б) y

′′

= e

y

 

= (y )

(y ) ;

 

 

y .

5.Решить уравнения, не содержащие y :

а) xy

′′′

′′

б) y

(4)

+ y

′′

= x.

 

= y ;

 

 

6.Решить однородные по y, y,..., y(n) уравнения:

а) xyy

′′

2

б) xyy

′′

2

= 2 yy

;

 

x(y )

= yy ;

 

+ x(y )

 

в) x2 yy′′ = (y xy)2 .

7.Решить уравнения путем приведения к дивергентному виду:

а) yy

′′

);

б) yy

2

 

= y (1

+ y

 

+ (y ) =1.

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

Система дифференциальных уравнений относительно неизвестных функций y1(x),..., yn (x) вида

y

= a

y

+ a

y

2

+... + a

y

n

+ q (x),

1

11

1

12

 

1n

 

1

y2= a21 y1 + a22 y2 +... + a2n yn + q2 (x),

M

 

 

 

 

 

 

 

 

 

 

 

 

+ an2 y2

+... + ann yn + qn (x),

yn′ = an1 y1

где aki заданные числа ( k,i =1, n ), а q1(x),..., qn (x) заданные функции, называется системой ли-

нейных дифференциальных уравнений с постоянными коэффициентами. При qk (x) = 0 для всех

k =1, n говорят об однородной системе, в противном случае данную систему называют неоднород-

ной.

Метод подстановки

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

П ри ме р ы

1. Найти общие решения следующих систем линейных дифференциальных уравнений (однородной (а) и неоднородной (б):

а)

y′= −y + z,

б)

y′ = 5y 2z,

 

 

 

sin x.

 

z′= 4 y z,

 

z′ =13y 5z

Ре ш е н и е .

а) Дифференцируя первое уравнение системы, получаем: y'' = −y'+z' .

Исключаем из системы переменную z . Имеем: z = y′+ y , z′ = 4 y y′− y , так что y′′ = −y′+ 4 y y′− y . В результате получаем однородное линейное уравнение второго порядка

y′′+ 2 y′−3y = 0 .

Его общее решение имеет вид y = C1ex + C2 e3x . Теперь находим z :

z = y′+ y = C1ex 3C2 e3x + C1ex + C2 e3x = 2C1ex 2C2 e3x .

б) Дифференцируя первое уравнение, с учетом обоих уравнений находим y′′ = 5y′− 2z. Исключая из системы переменную z , получаем линейное неоднородное уравнение второго порядка

y′′+ y = 2sin x .

Его общее решение имеет вид y = C1sin x +C2 cos x x cos x . Теперь находим z :

z = −0,5y′+ 2,5y ;

z= (2,5C1 + 0,5C2 0,5x) sin x (0,5C1 2,5C2 + 2,5x 0,5) cos x .

2.Найти общее решение неоднородной системы линейных дифференциальных уравнений с постоянными коэффициентами

x(t) = 4x y z,

 

 

 

 

 

 

 

 

 

= x + 2 y z + 2,

 

 

 

 

 

 

 

 

y (t)

 

 

 

 

 

 

 

 

= x y + 2z.

 

 

 

 

 

 

 

 

 

z (t)

 

 

 

 

 

 

 

 

 

Р е ш е н и е . Выразим неизвестную

функцию

z(t)

из

первого

уравнения

системы:

z = 4x y xи подставим ее в оставшиеся уравнения:

x′′

= 6x′− y′−9x +3y,

Далее,

подставим

 

= x′−3x +3y +

2.

 

y

 

 

 

 

 

′′

+5x

6x = 2 . Итак, мы полу-

y (t) из последнего уравнения полученной системы в первое: x

 

 

чили неоднородное линейное дифференциальное уравнение 2-го порядка с постоянными коэффициентами относительно одной неизвестной функции x(t) . Решая его, получаем

x(t) = C1e2t +C2e3t 13 .

Подставим найденное решение x(t) в уравнение для y(t) : y′−3y = x′−3x + 2 ,

получим y′−3y = −C1e2t + 3 , общим решением которого является

y(t) = C3e3t 1 +C1e2t .

Для нахождения последней неизвестной функции z(t) подставим найденные x(t) и y(t) в выражение для z(t) : z = 4x y x. В результате получим

z(t) = C1e2t + e3t (C2 C3 ) 13 .

З а м е ч а н и я . Несмотря на универсальность метода подстановки, применение его к системам высших порядков (пример 2) может привести к неоправданно большим вычислениям. Данный метод удобнее всего применять к системам для двух неизвестных функций (см. Пример 1).

Матричное представление

Запишем общий вид линейной системы дифференциальных уравнений в матричном виде. Пусть

y

 

(x)

 

a

 

a

L a

 

 

 

q

(x)

 

 

1

 

 

 

 

11

 

12

 

1n

 

1

 

 

 

y2

(x)

,

a21

a22

L a2n

,

q2

(x)

,

Y =

 

M

 

A =

M

 

M

O

M

 

Q =

M

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

 

L a

 

 

 

 

y

n

(x)

 

a

n1

n2

 

 

 

q

(x)

 

 

 

 

 

 

 

 

 

 

nn

 

n

 

 

 

где Y – столбец неизвестных функций,

A – матрица заданных числовых коэффициентов и Q

столбец правых частей (неоднородность). Тогда матричная форма записи системы приобретает вид

Y ′ = AY +Q .

Общим решением данной системы, согласно теории линейных уравнений, является сумма частного решения Y1 неоднородной системы и общего решения Y0 однородной системы:

Y =Y0 +Y1 .

Рассмотрим методы построения Y0 и Y1 .

Однородная система

Для решения однородной системы Y ' = AY необходимо произвести следующее действие.

1. Найти собственные числа λk матрицы A как решения характеристического уравнения

 

 

a

−λ

a

L

a

 

 

 

 

 

11

 

12

 

1n

 

 

 

 

 

a21

a22 −λ L a2n

 

= 0 .

 

det

 

M

M

O

M

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

an1

an2

 

 

 

 

 

 

 

L ann −λ

 

2.

В зависимости от кратности каждого корня λ , найти все линейно независимые собствен-

ные векторы V1,K,Vm как решение линейной алгебраической системы

 

a

−λ

a

L

a

 

 

 

 

11

 

 

12

 

1n

 

 

 

 

a21

 

a22 −λ L a2n

 

 

 

 

 

M

 

M

O

M

V

= 0

 

 

 

 

 

 

 

 

 

 

 

 

an1

 

an2

L ann −λ

 

 

 

 

 

 

 

и записать соответствующие частные решения:

а) случай корня λ кратности 1. Ему соответствует лишь один собственный вектор V1 , следова-

тельно, частное решение, отвечающее данному корню λ , имеет вид C1V1eλx , где C1 – произ-

вольная постоянная (пример 3);

б) случай корня λ кратности k . Ему соответствует ровно k линейно независимых собственных векторов V1,K,Vk . Тогда соответствующее частное решение, отвечающее данному кор-

ню λ , имеет вид

C1V1eλx + C2V2eλx +... + CkVk eλx ,

где C1,C2 ,K,Ck – произвольные постоянные (пример 4).

в)* случай корня λ кратности k , но ему соответствует m (m < k) линейно независимых собственных векторов V1,K,Vm . В этом случае решение, отвечающее данному λ , необходимо искать в специальном виде:

y = (a + b x + c x2

+... + d xk m )eλx ,

 

 

 

 

1

 

 

 

 

 

M

 

 

 

 

 

 

 

 

 

 

 

yn = ( p + r x + s x2 +... + h xk m )eλx ,

 

 

 

 

 

 

y

 

 

где коэффициенты a, b, c, d,..., p, r, s, h требуют определения путем подстановки

 

1

 

в од-

 

M

 

 

 

 

 

 

 

 

 

yn

 

нородную систему (пример 5).

3. В окончательный ответ для общего решения однородной системы Y0 необходимо записать сумму найденных частных решений, соответствующих каждому λk (согласно пункту 2).

П ри ме р ы

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

а)

y(x) = 4 y + 3z,

б)

y(x) = 4 y 5z,

 

z (x) = 3y 4z,

 

z (x) =5y 4z.

Ре ш е н и е :

а) матрица данной системы имеет вид

4

3

 

A =

 

 

.

 

3

4

 

 

 

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

 

4 −λ

3

 

= 0

 

λ2 25 = 0 ,

 

 

 

3

4 −λ

 

решением которого являются собственные числа λ = ±5 .

2.Найдем собственные векторы, отвечающие найденным собственным значениям. Для λ = 5 получаем следующую систему:

 

 

 

1

3 v

 

=

 

0

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

.

 

 

 

 

 

 

 

 

3

 

9

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

v2

 

 

 

0

 

 

 

 

Выбирая в качестве свободной переменной

v2

и

 

полагая

v2 =1 , получаем v1 = 3 , т.е.

V

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

– собственный вектор, отвечающий λ = 5 . Соответствующее частное решение, со-

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

гласно пункту 2а, записываем в виде

C

 

3

 

 

 

 

 

 

 

 

 

 

 

e5x .

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

Аналогично для λ = −5

получим: V

 

 

1

и решение C

 

 

1

 

=

 

 

2

 

 

e5x .

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

3

3.Общим решением однородной системы а) является

y(x)

= C

 

3

 

 

1

 

 

 

e5x + C

 

e5x .

 

 

 

1

 

 

 

2

 

z(x)

 

 

 

1

 

 

3

 

 

4

 

5

 

 

б) Матрица данной системы имеет вид A =

 

 

 

.

 

 

 

 

 

5

 

4

 

 

 

 

 

 

 

 

 

 

1. Вычисление собственных значений приводит к λ2 +9 = 0 , λ = ±3i .

2. Возникновение пары комплексно-сопряженных корней всегда упрощает дальнейшее нахождение соответствующего им решения. В этом случае достаточно найти собственный вектор, отвечающий какому-либо одному из найденных λ (в нашем случае, например, для λ = 3i ), и в полученном решении выделить реальную и мнимую части, которые и будут соответство-

вать сразу обоим корням:

4

3i

5 v

 

 

0

. В качестве свободной пере-

λ = 3i

 

 

1

 

=

 

 

 

 

5

4

 

 

 

0

 

 

 

 

3i v2

 

 

 

 

менной выберем v

и для простоты положим v

= 5 , тогда

v

 

= 4

3i V

 

 

5

 

2

=

 

 

. Част-

 

 

 

1

 

 

 

 

 

1

 

 

 

 

 

1

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3i

ное решение, соответствующее

V , задается выражением

 

5

 

 

 

 

 

 

 

 

 

 

 

ei3x , от которого вычислим

 

 

 

 

 

 

1

 

 

 

 

 

4

3i

 

 

 

 

 

 

реальную и мнимую части:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5

 

 

5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Re

 

ei3x = Re

 

(cos3x +i sin 3x) =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

3i

 

3i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5cos3x

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4cos3x +3sin 3x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5

 

 

 

5sin 3x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Im

 

ei3x =

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3i

 

3cos 3x + 4 sin 3x

 

 

 

 

 

 

3.Общее решение системы б) приобретает следующий вид:

y(x)

= C

 

5cos 3x

 

+C

 

5sin 3x

 

 

 

 

 

 

 

 

.

 

 

1

 

 

 

 

2

 

 

z(x)

 

 

4 cos 3x +3sin 3x

 

 

3cos 3x + 4 sin 3x

4. Решить однородную систему линейных уравнений из примера 2:

x(t) = 4x y z,

y(t) = x + 2 y z,z(t) = x y + 2z.

Ре ш е н и е . Найдем собственные значения матрицы системы как решение характеристического уравнения

 

4 −λ

1

1

 

 

 

 

1

2 −λ

1

= 0 .

 

1

1

2 −λ

 

Разложим определитель левой части уравнения по первой строке и после несложных упрощений получим (λ −3)(−λ2 +5λ −6) = 0 . Первый сомножитель дает корень λ = 3 , второй сомножитель – λ = 2

и λ = 3 . Таким образом, мы получили два различных собственных значения, причем корень λ = 3 встречается два раза, следовательно, его кратность k = 2 . Перейдем к определению собственных векторов.

Рассмотрим случай λ = 2 кратности 1. Соответствующая система для собственного вектора V

приобретает вид

 

 

 

 

 

 

 

2 1

1 v1

 

 

0

 

1

0

 

 

 

0

 

 

1 v2

 

=

.

 

1

 

 

 

 

0

 

 

1

0 v3

 

 

 

Здесь мы выделили один из возможных базисных миноров (определитель максимального порядка, не равный нулю). Таким образом, переменная v1 становится свободной. Положим v1 =1 , тогда остав-

1

шиеся компоненты собственного вектора легко находятся: v2 = v3 =1 или V1 = 1 , что соответствует

1

1 2t

частному решению C1 1 e .

1

Перейдем к случаю λ = 3 кратности k = 2 . Действуя аналогично предыдущему случаю, получим следующую систему:

[1]

1

1 v

 

 

0

 

1

1

1

 

 

0

 

 

1 v2

 

=

.

 

1

1

 

 

 

0

 

 

1 v3

 

 

 

Заметим, что двойная кратность корня λ = 3 в данном случае приводит к появлению двух свободных переменных v2 и v3 , не вошедших в выделенный базисный минор. Выберем сначала v2 =1, v3 = 0 ,

 

 

 

1

 

 

 

тогда v1

=1 и в совокупности V1

 

1

 

. Далее выберем

v2 = 0, v3 =1 , для которых v1 =1 , и соответ-

=

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

1

 

 

ствующий второй линейно независимый собственный вектор V2

 

0

 

. Количество линейно незави-

=

 

 

 

1

 

 

 

 

 

 

симых собственных векторов совпадает с кратностью корня k = m (случай 2б), следовательно, част-

 

 

 

1

 

 

1

 

ное решение, соответствующее данному корню λ = 3

, записывается в виде C2

 

1

 

 

0

 

 

e3t +C3

 

e3t .

 

 

 

0

 

 

1

 

 

 

 

 

 

 

Общим решением системы будет являться сумма полученных частных решений

x(t)

 

1

 

1

 

1

y(t)

 

= C

1 e2t +C

2

 

1

e3t +C

3

 

0

e3t .

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

1

 

z(t)

 

 

1

 

 

 

 

 

 

З а м е ч а н и е . Если в полученном ответе произвести переобозначение констант C1 C1,C2 C3,C3 C2 C3 , то новая форма записи ответа в точности совпадет с полу-

ченной ранее в примере 2 частью решения, отвечающего однородной системе (покажите это).

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

x(t) = x + z,

y(t) = y + 2z,z(t) = 2x y +3z.

Ре ш е н и е . Составим характеристическое уравнение для определения собственных значе-

ний:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 −λ

0

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

1 −λ

2

 

= 0 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

1

3 −λ

 

 

 

 

 

 

 

 

 

 

Раскладывая определитель по третьему столбцу, получим (3 −λ)(1 −λ)2 = 0 , откуда λ = 3

крат-

ности 1 (простой корень) и λ =1

кратности k = 2 . Далее находим собственные векторы. Для λ = 3

 

 

2

0

1 v1

 

0

, переменная v – свободная, выберем v

 

 

 

= 1 ,

получаем

 

 

 

v

2

 

=

0

 

3

=1

и получим v

 

 

0

2

2

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

1

2

 

 

2

1

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0 v3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

v2 =1 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

В итоге соответствующий собственный вектор имеет вид V1 =

 

 

. Для дальнейшего удобства

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

умножим его на 2: V1 = 2 . Следовательно, частное решение, отвечающее корню λ = 3 , принимает

2

 

 

1

 

 

 

 

 

 

 

 

 

 

вид C

 

2

e3t .

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

0 1 v1

 

 

0

 

Переходим к следующему корню λ =1 кратности k = 2 :

 

0

0

 

 

 

0

 

, выбираем

 

2 v2

 

=

 

 

 

 

 

 

2

 

 

 

 

0

 

 

 

 

 

 

 

1

2 v3

 

 

 

 

1

свободную переменную v1 =1 и получаем V1 = 2 . Мы нашли все линейно независимые собствен-

0

ные векторы (всего один – V1 ), отвечающие λ =1 . Однако кратность корня k = 2 превышает количе-

ство собственных векторов m =1 , следовательно, согласно пункту 2в, мы должны искать соответствующее решение в виде

x(t)

a +bt

 

 

 

 

 

 

1 t

.

y(t)

= c + dt e

 

 

 

 

 

 

z(t)

p + qt

 

 

Подставим искомое решение в систему и получим:

a +b +bt t

c + d + dt e

p + q + qt

 

a + p +(b + q)t

 

 

c + 2 p +(d + 2q)t

 

=

et .

 

 

 

 

2a c +3p +(2b d +3q)t

Приравнивая подобные слагаемые в каждом уравнении, получим систему алгебраических уравнений для определения неизвестных констант a,b, c, d, p и q :

a +b = a + p

 

p = b

 

 

+ q

 

 

q = 0

 

b = b

 

 

 

= c + 2 p

Решая ее, получим

 

d = 2b

.

c + d

d = d + 2q

 

 

d = d

 

p + q = 2a c +3p

 

 

c = 2(a + b)

 

 

 

 

0 = 2b 2b

 

q = 2b d +3q

 

 

Заметим, что четвертое и шестое уравнения системы выполнены тождественно, а оставшиеся

уравнения выражают неизвестные

p, q, d

и c через произвольные константы a и b , которые естест-

венно переобозначить в качестве a = C2 ,

b = C3 . Выпишем полученный ответ для корня λ =1 крат-

ности 2:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

C2 +C3t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2(C2 +C3 ) + 2C3t et .

 

 

 

 

 

 

 

C3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Общее решение однородной системы принимает вид суммы

 

x(t)

 

1

 

C2 +C3t

 

 

y(t)

 

= C

 

2

e3t +

2(C

2

+C

3

) +

2C t

et .

 

 

1

 

 

 

 

 

 

3

 

 

 

 

 

2

 

 

 

C3

 

 

z(t)

 

 

 

 

 

 

 

 

Однородная система. Метод присоединенных векторов*

В последнем примере 5 мы разобрали решение однородной системы, в которой количество собственных векторов для некоторого корня ( λ =1 ) меньше его кратности (случай 2в), что привело к достаточно громоздким вычислениям. Возникающих в этом случае трудностей можно избежать, применяя метод присоединенных векторов.

Пусть λ – собственное значение кратности k >1, и ему соответствует один собственный вектор

V1 , т.е. V1

удовлетворяет уравнению AV1 = λV1 (случай 2в). Определим присоединенные векторы

V2 ,V3,K,Vk

как последовательные решения следующих алгебраических систем:

 

AV2

= λV2 +V1,

 

AV3

= λV3 +V2 ,

 

M

 

AVk = λVk +Vk 1.

Частное решение, соответствующее данному корню λ , принимает вид

C V eλt +C

 

(

t

V +V )eλt +C

(

t2

V +

t

V +V )eλt +K

 

 

 

 

 

 

1

1

 

2

 

1!

1

2

3

2!

1

1!

2

3

 

 

K+C

 

(

tk

V

+

tk 1

V +K+

t

V

+V )eλt .

 

 

 

k!

(k 1)!

 

 

 

 

 

k

 

1

 

2

 

 

1!

k 1

k

Несмотря на кажущуюся громоздкость, применять этот метод намного проще, чем решать задачу на отыскание неопределенных коэффициентов, как это было сделано в примере 5.

П р и м е р 6*. Решить однородную систему линейных дифференциальных уравнений из примера 5, используя метод присоединенных векторов.

Ре ш е н и е . В примере 5 мы получили следующие собственные значения и векторы: λ = 3

 

 

1

 

1

 

 

кратности 1, ему соответствует вектор V

=

 

2

 

и решение C

 

2

e3t ; и

λ =1 кратности 2, ему соот-

1

 

 

 

 

1

 

 

 

 

 

 

 

2

 

 

 

2

 

 

 

 

 

 

 

 

 

 

1

ветствует только один собственный вектор V1 = 2 . Для построения решения необходимо найти еще

0

один недостающий вектор V2 , но уже не собственный (линейно независимых с V1 больше нет), а

 

0

0 1 v1

 

 

1

 

 

присоединенный:

 

0

0

 

 

 

2

 

. Здесь выписано уравнение

AV2 = λV2 +V1 или

 

2 v2

 

=

 

 

 

2

 

 

 

 

0

 

 

 

 

 

1

2 v3

 

 

 

 

 

( A −λE)V2 =V1 . Данная система является неоднородной, для нахождения какого-либо частного решения положим свободную переменную v1 = 0 и найдем оставшиеся базисные переменные: v2 = 2 ,

 

 

0

 

v3 =1 . Следовательно, V2

 

2

 

=

.

 

 

1

 

 

 

 

 

 

 

1

 

 

1

 

0

 

 

Запишем соответствующее частное решение

C2

 

2

 

t

 

2

 

 

2

 

t

. Общее решение систе-

 

e

 

+C3

t +

e

 

 

 

 

0

 

 

 

0

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

мы представится в виде суммы

 

 

 

 

 

 

 

 

 

 

 

 

 

x(t)

 

1

 

 

C2 +C3t

 

 

y(t)

 

= C

 

2

e3t +

 

2(C

2

+C

3

) +

2C t

et ,

 

 

1

 

 

 

 

 

 

 

3

 

 

 

 

 

2

 

 

 

 

C3

 

 

z(t)

 

 

 

 

 

 

 

 

 

что полностью совпадает с ответом в примере 5, полученным стандартным способом.

Неоднородная система

Частное решение неоднородной системы линейных дифференциальных уравнений строится по тем же правилам, что и для одного уравнения, с разницей лишь в том, что в случае совпадения кор-

ней определяемое решение не умножается на xs ( s – кратность собственного значения), но при этом на s повышается степень многочлена, содержащего искомые коэффициенты. Рассмотрим специфику данного метода на конкретных примерах.

П ри ме р ы

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

 

x

,

 

 

 

 

 

 

 

 

 

 

 

 

 

y(x) = 4 y +3z + e

 

б)

y (x) = 4 y +3z +2,

 

 

 

 

 

 

 

 

 

а)

 

x

 

 

5x

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z (x) = 3y 4z 2e ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ре ш е н и е .

 

 

 

 

 

 

 

 

y

y

 

y

 

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

 

. Об-

 

 

=

 

 

+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z

z

0

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

щее решение однородной системы

y

было получено в примере 3а, таким образом, остается найти

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z 0

 

 

 

 

 

 

 

 

 

 

 

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

y

 

Q

 

e

x

 

 

 

. Столбец правых частей для данной системы

=

 

можно представить

 

 

 

 

 

 

 

 

2ex

 

 

 

 

z

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

в виде

 

1

1 x

. При экспоненте стоит столбец из многочленов нулевой степени, а множитель 1

Q =

e

 

 

 

 

 

 

 

 

 

 

 

 

при x

 

2

 

 

 

 

 

 

 

 

 

в показателе экспоненты не совпадает с собственными значениями: 1 ≠ ±5 , следовательно, ча-

стное решение необходимо искать в том же виде, что и вектор Q :

y

 

a

 

1 x

, где a и b – кон-

 

 

=

e

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z

 

b

 

 

 

 

 

 

 

 

 

1

 

 

 

 

станты, требующие определения. Подставим искомое частное решение в неоднородную систему а)

a

 

4a +3b +1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ex

=

ex Сокращая на экспоненту и приравнивая соответствующие компоненты век-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

b

 

3a 4b 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

торного уравнения, в итоге получим a =

 

1

24

и b = 3

8

. Запишем общее решение неоднородной сис-

темы:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y(x)

 

3

 

 

 

1

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

+C

 

 

24 ex .

 

 

 

 

 

 

 

 

= C

 

 

e5x

2

 

 

e5x +

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

z(x)

 

1

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

8

 

 

 

 

 

 

 

Система б). Общее решение однородной системы, как и в случае а), было получено ранее в при-

мере 3а. Найдем частное решение данной неоднородной системы. Вектор правых частей

 

2

 

Q =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5x

необходимо представить в следующем виде:

 

 

 

 

 

 

 

 

 

 

 

e

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

0

 

2

 

 

 

0

 

 

+Q

 

.

 

 

 

 

 

 

 

Q =

+

 

5x

= e0 x +

e5 x = Q

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

0

e

 

 

 

 

0

 

 

 

1

 

 

 

 

 

 

 

 

 

Мы видим, что правая часть состоит из двух простейших слагаемых, разнесенных по экспонен-

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

y

 

y

+

y

,

 

 

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z

 

z

I

z

II

 

 

 

 

 

 

1

 

 

 

 

 

где первый столбец отвечает за неоднородность Q

 

2

 

, второй соответствует Q

 

=

0

 

 

 

 

=

 

e0 x

2

 

e5 x .

 

 

1

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

Для

y

получаем аналогичную случаю а) ситуацию: при e0 x

стоит столбец констант – мно-

 

 

 

 

 

 

 

 

 

 

 

 

z I

 

 

 

 

 

 

гочленов нулевой степени

 

2

 

, а множитель 0 в показателе при x

не совпадает с собственными

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

числами 0 ≠ ±5 . Таким образом, частное решение

y

 

ищем в виде

y

a

 

 

a

 

 

 

 

 

 

 

=

e0 x

=

, а под-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z I

 

 

 

z I

b

 

b

 

 

становку для определения a и b осуществляем не во всю неоднородную систему

y

 

y

 

 

 

= A

 

+Q , а

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z

 

 

z

 

 

 

y

y

 

0

 

4a +3b + 2

 

 

 

y

 

 

8

 

 

 

 

 

 

 

=

 

25 .

только в соответствующую ее часть

 

 

= A

+Q

:

 

=

 

. В итоге

 

 

 

 

 

 

 

 

 

1

 

 

 

3a 4b

 

 

 

 

 

 

6

 

 

 

z

 

z

 

0

 

 

 

 

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

25

y

, отвечающего за неоднородность Q

 

0

 

 

 

 

 

 

 

 

 

 

Для

 

2

=

e5 x , получаем несколько иную ситуацию.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z II

 

 

 

 

 

 

1

 

 

 

 

при x

совпадает с

При экпоненте стоит столбец многочленов нулевой степени, однако множитель 5

одним из собственных значений, причем кратность этого корня s =1 . Следовательно,

мы должны

увеличить на s =1 степень нулевого многочлена, и частное решение

y

представится в виде

 

 

 

 

 

 

 

z II

 

y

a + bx

. Подстановку осуществляем в систему с правой частью Q

 

:

 

 

=

e5x

2

 

 

 

 

 

 

z II

c + dx

 

 

 

b +5a +

5bx

 

4a + 3c + (4b

+ 3d )x

 

 

 

 

e5x =

 

 

 

 

 

 

 

e5x .

 

 

 

 

 

 

 

 

 

 

 

 

 

d +5c +

5dx

 

3a 4c +1 + (3b 4d )x

Приравнивая соответствующие слагаемые, получаем систему

 

 

b +5a = 4a +3c

 

 

 

 

a = 3(c d )

 

 

5b = 4b +3d

 

 

 

 

 

 

b = 3

 

 

 

 

 

 

 

 

10 .

 

 

 

= 3a 4c +1

 

 

 

d +5c

 

 

 

c любое

 

5d

= 3b 4d

 

 

 

 

 

 

d = 1

 

 

 

 

 

 

 

 

10

Выберем какое-либо c , например, пусть c = d =

1

 

,

следовательно, a = 0 . Таким образом, мы на-

y

 

 

0.3x

 

 

10

 

 

 

 

 

 

 

 

 

 

 

 

 

 

шли последнее неизвестное:

 

=

 

 

e5x . Сумма общего решения однородной системы и

 

 

 

 

 

 

 

 

 

 

 

 

 

z II

 

0.1 +0.1x

 

 

 

 

 

 

двух найденных частных решений будет являться общим решением системы б).

8.Найти общее решение неоднородной системы линейных уравнений

 

 

x

,

y (x) = 4 y 5z + 2xe

 

 

x

 

 

 

+sin x.

z (x) = 5y 4z + 2e

 

Ре ш е н и е .

Представим

правую

часть

 

 

2xex

 

в

виде

Q =

 

x

 

 

 

 

 

 

 

2e

 

 

 

 

 

 

 

 

 

 

+sin x

 

 

 

2x

 

0

 

0

 

= Q

+Q

 

. Первое простейшее слагаемое отвечает за множитель

Q =

 

ex +

 

 

sin x +

 

cos x

 

2

 

2

 

 

1

 

 

0

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

e1 x , второе (выделенное в квадратных скобках) соответствует e±i x . Общее решение системы будем

искать в виде

y

y

y

y

. Общее решение однородной системы было получено в при-

 

 

=

 

+

 

+

 

 

 

 

 

 

 

 

 

 

 

 

z

z 0

z I

z II

 

мере 3б. Найдем частные решения. Для

y

, соответствующего правой части Q

 

2x

1 x

, получа-

 

 

=

 

e

 

 

 

1

 

2

 

 

 

 

z I

 

 

 

 

 

ем следующее: столбец при экспоненте состоит из многочленов первой и нулевой степени, следовательно, выбираем максимальную, т.е. первую степень, а множитель при x в экспоненте не совпадает

с собственными значениями 1 ≠ ±3i ( λ = ±3i

см. пример 3б), следовательно

y

a

+bx

 

 

 

 

=

 

 

ex . После

 

 

 

 

 

 

 

 

 

 

 

 

z I

c + dx

 

 

подстановки в систему с соответствующей правой частью Q , получаем решение:

y

 

1 + x

 

 

=

4

ex .

 

1

 

 

 

 

 

5

+ x

 

 

 

 

z

I

 

 

 

Для

y

, отвечающего за Q

 

 

0

 

 

0

 

cos x

и

 

 

2

=

 

sin x +

 

cos x , получаем: столбцы, стоящие при

 

 

 

 

 

1

 

 

0

 

 

 

 

z II

 

 

 

 

 

 

 

 

sin x , состоят из многочленов нулевой степени, а ±i ≠ ±3i , следовательно, второе частное решение ищем в виде

y

a

b

a cos x +bsin x

 

 

=

cos x +

sin x =

.

 

 

 

 

 

 

 

 

z 2

c

d

 

ccjsx + d sin x

Действуя далее аналогично рассмотренным случаям, получаем ответ и для этого частного решения: