- •Вопросы и упражнения для самостоятельной работы
- •Вопросы и упражнения для самостоятельной работы
- •Вопросы и упражнения для самостоятельной работы
- •Вопросы и упражнения для самостоятельной работы
- •Вопросы и упражнения для самостоятельной работы
- •Вопросы и упражнения для самостоятельной работы
- •Вопросы и упражнения для самостоятельной работы
- •Вопросы и упражнения для самостоятельной работы
- •Ответы и указания к упражнениям
в) характеристическое уравнение λ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 e−4 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 = 8e−2 x ; |
|
д) |
y′′+ y = 4 sin x ; |
е) |
y′′+ 2 y′+5y = e−x 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 = Ce−x −1 + x , где C − произвольная константа (убедитесь в этом).
Возвращаясь к первоначальной функции, получим
y′′′ = Ce−x −1 + x ,
трехкратное интегрирование которого приводит к окончательному ответу
y = C1e−x + 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′ = Cxx−1 .
Данное уравнение также легко интегрируется. В итоге записываем окончательный ответ
|
|
|
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 e−3x . Теперь находим z :
z = y′+ y = C1ex − 3C2 e−3x + C1ex + C2 e−3x = 2C1ex − 2C2 e−3x .
б) Дифференцируя первое уравнение, с учетом обоих уравнений находим 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 |
|
|
e−5x . |
||||||||||
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
−3 |
|
|
|
|
|
|
|
|
−3 |
3.Общим решением однородной системы а) является
y(x) |
= C |
|
3 |
|
|
1 |
||
|
|
|
e5x + C |
|
e−5x . |
|||
|
|
|
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 |
|
|
e−5x + |
|
|
|
|
|||||||||
|
|
|
|
|
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 |
Действуя далее аналогично рассмотренным случаям, получаем ответ и для этого частного решения: