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

1









ЛИНЕЙНЫЕ СИСТЕМЫ ДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИЙ  С  ПОСТОЯННЫМИ   КОЭФФИЦИЕНТАМИ

§  I. МЕТОД ЭЙЛЕРА



211.   Предварительные  замечания.  В  этой  главе  мы  будем изучать линейные системы уравнений:



   







Или                                                                                                             (1′)                                                                                                                                                                                                                                   









Где коэффициенты akl=(k, l=1,2,…,n) – постоянные вещественные числа, а fk(x) (k=1,2,…,n) – функции от х, непрерывные в интервале (a,b).

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

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









В силу теоремы о построении общего решения, для построения общего решения системы  (2)  достаточно построить хоть одну фундаментальную систему решений.

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

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

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



212. Построение фундаментальной системы решений и общего решения однородной линейной системы в случае различных корней характеристического уравнения. По аналогии с однородным линейным уравнением с постоянными коэффициентами будем искать частное решение системы  (2) в виде

          	    (3)

	nи  - некоторые постоянные числа, причем числа n не равны нулю одновременно, ибо в противном случае мы получили бы очевидное нулевое решение, которое не может входить в состав фунда- ментальной системы и, следовательно, не может быть использовано для построения общего решения.

Обратим особое внимание на то, что число  мы берем одно и то же для всех функций, составляющих решение.

Подставляя функции (3) в систему (2), сокращая на еx  и пере- нося все члены направо, получим для определения чисел  следующую систему:  

                  





(4)







Нас интересует ненулевое решение этой системы. Такое решение существует лишь при условии, что определитель системы равен нулю, т. е. при условии

   



                                                (5)                                                                                      

                                                                                                             



 Уравнение (5) называется характеристическим уравнением системы

(2), его корни— характеристическими числами, а определитель   ) -характеристическим определителем.

      Рассмотрим сначала случай, когда все характеристические числа     

n    различны. В этом случае имеем:  i)=0, но   Вследствие этого ранг матрицы 













Составленной из коэффициентов системы











которая получается из системы (4) после замены в ней на iравен n-1.

    Действительно, вычисляя ’(), имеем:   







где ll (алгебраическое дополнение элемента all - определителя Так как   то из (8) видим, что хоть один из определи- телей (п—1)-го порядка, именно один из ll (i, отличен от нуля, так что ранг рассматриваемой матрицы равен п — 1.

Поэтому одно из уравнений системы (7) есть следствие осталь-ных и эта система имеет ненулевое решение, определенное с точ-ностью до произвольного множителя пропорциональности Ai:

i1 = Aimi1, i2 = Aimi2,…, in = Aimin (i=1.2,…, n).     (9)

Например, в качестве ik можно взять алгебраические дополнения элементов любой строки определителя l(i, если не все они равны нулю. В самом деле, так как сумма произведений элементов какой-либо строки определителя l(i на алгебраические дополнения эле-ментов другой строки равна нулю, а сумма произведений элементов строки на их алгебраические дополнения равна самому определителю (i т. е. снова равна нулю, то ясно, что, заменив в системе (7) неизвестные k взятыми алгебраическими дополнениями, мы  получим  тождества.

Фиксируя в формулах (9) множитель Ai, мы получим определен-ное решение системы (7).

Подставляя теперь в (3) вместо  последовательно характеристи-ческие числа  i ,а вместо   n — соответствующие им решения системы (7), определяемые формулами (9) при фиксированных множителях Ai, получим п решений:







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

   Если при этом все корни n  вещественны, то все решения (10) тоже будут вещественными.

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

Поэтому, в силу теоремы о построении общего решения, формулы







Дают общее решение системы (2) в области







Если характеристические числа различные, но среди них есть комплексные, то последние входят сопряженными парами. Пусть      a + ib и а – ib — простые корни характеристического уравнения. Корню a+ib соответствует согласно формуле (3) решение

      y1=e(a+ib)x,  y2=e(a+ib) , … , yn=ne(a+ib).                          (13)

здесь n – комплексные числа. Полагая 

1=11+i21,   2=12+i22, … ,n=1n+i2n ,

Получаем решение 

      y1=(11+i21) e(a+ib)x, y2=(12+i22) e(a+ib)x, … , yn=(1n+i2n) e(a+ib)x.

   Это   решение   комплексное. Отделяя в нем вещественные и мнимые части,  мы    получим,    согласно свойствам решений однородной системы, два  вещественных решения:





Эти решения, очевидно, линейно независимы в интервале . Нетрудно убедиться, что сопряженный корень а—ib не порождает новых вещественных линейно независимых частных решений.

Таким образом, если все характеристические числа — различные и вещественные, то мы получаем соответствующие им вещественные линейно независимые частные решения в виде (10). Если же все характеристические числа — различные, но среди них есть комплексные, то последние обязательно входят сопряженными парами и каждой паре таких характеристических чисел соответствуют два линейно независимых частных решения вида (15). Всего мы получим п вещественных частных решений.  Все  эти решения  линейно  независимы  в  интервале .

В самом деле, предположим обратное. Тогда, написав соответствующую систему соотношений между этими решениями и перейдя в ней от тригонометрических функций к 'показательным, мы получили бы, что системы функций вида (10), где 1,2, … ,n —различные числа, оказались бы линейно зависимыми.

Общее решение системы (2) в области (12) представляет собою линейные комбинации построенных п вещественных линейно независимых частных решений с произвольными постоянными коэффициентами.

 

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





Решая характеристическое уравнение 

	              или      2-10+9=0;

находим: 1=1, 2=9, так что характеристические числа различные и вещественные.

Составляем систему для определения чисел  и  соответствующих характеристическому числу 1 = 1. Матрица коэффициентов этой системы получается из матрицы



5-      4    

   4   5-заменой  на 1=1, так что искомая система будет иметь вид



                       

Здесь, как и следовало ожидать, второе уравнение является следствием первого (оно даже совпадает с первым уравнением) и его можно было и не выписывать. Полагая ,   находим , находим 2= −1

Таким  образом,  характеристическому  числу    Х1=1   соответствует  решение:

	y1= ех, z1 = −ex.	(20)

Аналогично,     решая     систему,     соответствующую    характеристическому числу  2=9:

                                         (21)

находим: 1=1, 2=1 так что этому характеристическому числу соот-ветствует решение:

	y2=e9x, z2=e9x	(22)

Мы получили фундаментальную систему решений:

	       (23)

	

	Беря линейную комбинацию,  получаем  общее решение:

	                           



Пример 2. Рассмотрим систему:

               



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









Имеет комплексные сопряженные корни λ1 = 2 + i, λ2 = 2 – i. Найдем решение, соответствующее λ1. Это решение имеет вид y = γ1e(2+i)x, 

z = γ2e(2+i)x. Числа γ1 и γ2 ищем из системы

         



Полагая γ1=1, находим γ2 = - i, так что искомым решением будет







Это решение комплексное. Отделяя в нем вещественные и мнимые части, получим два вещественных решения







Эти решения составляют фундаментальную систему решений, так что общим решением будет









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

     

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



Имеет различные и притом вещественные корни λ1 = 2,  λ2 = 3, λ3=6,  так что фундаментальная система решений имеет вид (10). Найдем сначала частное решение вида



Соответствующее характеристическому числу λ1 = 2. В качестве чисел γ11, γ22, … , γ1n можно взять алгебраические дополнения элементов первой строки определителя

                              

который получается из характеристического определителя Δ (λ) заменой  λ на λ1=2. Получаем



или (деля на 2)

 

Подставляя эти значения γ1k в (33), получим



Аналогично найдем, что в качестве чисел γ2k , γ3k , соответствующих характеристическим числам λ2 = 3, λ3=6, можно взять γ21 = 1, γ22 = 1, γ23 = 1, γ31 = 1, γ32 = -2, γ33 = 1. Фундаментальной системой решений будет



Так что общее решение имеет следующий вид





Случай наличия кратных корней характеристического уравнения.

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

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

Заметим, прежде всего, что если 1 есть простое характеристическое число, то независимо от того, будут среди остальных характеристических чисел встречаться кратные или нет, ему всегда соответствует одно частное решение вида:

	                       y1=1ex , y2=2ex , … , yn=nex     (38)

где 1, 2, … ,n — некоторые постоянные числа, определяемые с точностью до постоянного множителя.

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

При этом, так же как и для линейного однородного уравнения n-го порядка, оказывается, что одному характеристическому числу кратности k соответствует k линейно независимых частных решений.

Теорема. Если 1 есть характеристическое число кратности k, то ему соответствует решение вида

	y1=P1(x) ex , y2=P2(x) ex , … , yn=Pn(x) ex                                                 (39)

где P1(x), P2(x), … , Pn(x) суть полиномы от х степени не выше чем k−1, имеющие в совокупности k произвольных коэффициентов, так что среди всех коэффициентов всех этих полиномов k коэффициентов являются произвольными, а все остальные выражаются через них.

В частности может случиться, что все эти полиномы вырождаются в постоянные числа. В таком случае k-кратному характеристическому числу 1 будет соответствовать решение вида

	y1=1ex , y2=2ex , … , yn=nex	(40)

Однако здесь k из коэффициентов 1, 2, … , n являются произвольными, в то время как для простого характеристического числа произвольным является только один из них.

Практически при нахождении решения, соответствующего характеристическому числу 1 нужно искать решение в виде (39), считая P1(х), Р2(х), ..., Рп(х) полиномами (k−1)-й степени с неопределенными коэффициентами и, подставляя (39) в (2), выразить все коэффициенты через k из них, которые остаются произвольными.

Полагая поочередно один из этих произвольных коэффициентов равным единице, а остальные равными нулю, мы построим k линейно   независимых   решений,   соответствующих   характеристическому числу 1. Все эти частные решения будут составлены из произведений показательной функции ex на полиномы от х, степени которых не превышают k−1. Если же полиномы в формулах (39) вырождаются в постоянные числа, то мы получим  линейно независимых частных решений такого же вида, как и в случае простого корня характеристического уравнения.

Если 1 — вещественное характеристическое число, то построенные выше k линейно независимых решений будут вещественными.

Если же система (2) имеет комплексное характеристическое число a + ib кратности k, то оно имеет сопряженное характеристическое число а—ib той же кратности.

Построив k линейно .независимых комплексных решений, соответствующих характеристическому числу a + ib, и отделив в них вещественные и мнимые части, мы получим 2k, вещественных линейно независимых частных решений.

В общем случае каждому простому вещественному характеристическому числу соответствует одночастное решение, каждой паре простых сопряженных комплексных характеристических чисел соответствует два вещественных линейно независимых решения, вещественному характеристическому числу кратности k соответствует k вещественных линейно независимых частных решений, а каждой паре сопряженных комплексных характеристических чисел кратности k соответствует 2k вещественных линейно независимых частных решений. Всего получается п вещественных решений. Все эти решения линейно независимы в интервале(-∞,+∞), так что они образуют фундаментальную систему решений. Взяв линейные комбинации решений этой фундаментальной системы по столбцам с одними и теми же произвольными постоянными С1, С2, ..., Сп, мы получим общее решение системы (2) в области (12).



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

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

Указанный выше вид фундаментальной системы решений дает возможность сделать некоторые заключения об устойчивости нулевого решения однородной системы (2)* .

Теорема об асимптотической устойчивости (в смысле Ляпунова)  нулевого  решения однородной линейной системы  с постоянными коэффициентами. 

Рассмотрим систему: 



              



где akl — постоянные . вещественные числа, а х1, x2, ..., xn — неизвестные функции от времени t.

Теорема. Если все характеристические числа системы (41) отрицательные или имеют отрицательную вещественную часть, то нулевое решение

	x1 ≡ 0, х2 ≡ 0, ..., хп 5≡ 0	(42)

асимптотически устойчиво в смысле Ляпунова при t→+∞ причем начальные возмущения можно брать любыми.

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

 Теорема о неустойчивости нулевого решения однородной линейной системы с постоянными коэффициентами. Если хоть одно из характеристических чисел системы (41) положительно или имеет положительную вещественную часть, то нулевое решение (42)   неустойчиво в смысле Ляпунова при t→+∞.

 Приведение однородной линейной системы к системе с постоянными коэффициентами при помощи замены независимой переменной. Рассмотрим систему:





Введем вместо х новую независимую переменную t по формуле

	t =x).                                                            	(44)

Тогда получим систему:

	

	                   (k=1, 2, … , n).                                                                               (45)

Предположим, что коэффициенты этой системы постоянны, т. е.

         

Тогда pkl (х) имеют вид

       pkl(x)=aklψ΄(x),                                                        (47)

т. е. рк1 (х)  представляют собой произведения постоянных чисел на одну и ту же функцию от х.

Обратно,  если   коэффициенты pkl (x)   обладают этим  свойством, т. е. если

	                    pkl(x)=aklφ(x),	(48)

то, положив

               t=ψ(x)=∫φ(x)dx,                                             (49)

мы получим систему с постоянными коэффициентами akl. 



Пример 1.     Пусть дана система:

	

	Здесь условие (48) выполнено, причем φ(x)=1/x

	    Поэтому    подстановка

	   t=∫φ(x)dx=∫1/xdx=ln x   (x>0)	(51)

или

	x = et	(52)

приводит данную систему к системе с постоянными коэффициентами:

                         



Интегрируя эту систему, получаем:

                                         (54)

Поэтому общим решением системы (50) будет:

                                                       (55)

Отсюда видно, что решения системы (50) могут иметь особенность только в точке х=0, которая является единственной особой точкой этой системы. (В точке х=0 не выполнены условия теоремы существования). Наряду с такими решениями существует целое семейство решений  y1=Cx2, y2=Cx2, y3=Cx2 голоморфных в окрестности особой точки х = 0. Заметим, однако, что среди них (и вообще) нет решений, в которых функции у1, y2 и y3 стремились бы к пределам, не равным одновременно нулю, когда х стремится к особой точке х = 0.

	Интегрирование неоднородной линейной системы с постоянными   коэффициентами   методом   вариации   произвольных постоянных.   Рассмотрим  теперь   неоднородную   линейную систему с постоянными коэффициентами





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

Замечание. Если в системе (56) функции fk (x) представляют собою произведения показательной функции (с вещественным или комплексным показателем) на полином от х, то для построения общего решения этой системы можно вместо применения метода вариации произвольных постоянных найти частное решение методом неопределенных коэффициентов и прибавить его к общему решению соответствующей однородной системы. Тогда мы получим общее решение системы (56).





     Интегрирование  линейной  системы  с  постоянными  коэффициентами   при   помощи   приведения   ее   к  уравнению  п-го порядка (метод исключения). 

Применим к системе

      

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

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

y΄1=y2+y3,

y΄2=y1+y3,                                                         (2)

y΄3=y1+y2.

Дифференцируя  первое   уравнение  и,  пользуясь   вторым  и   третьим,  получаем:

	y˝1=2y1+y2+y3.	(3)

Но   y2 + у3 = у΄1. Поэтому

	                 y˝1− y΄1−2y1=0;	(4)

Исключим у3. Из первого уравнения системы (2) имеем:

	y3= y΄1−y2.	(5)

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

	y΄2= y΄1−y2+y1	(6)

или

	y΄2 +y2= y΄1 +y1	(7)

Таким   образом  система   (2)   приводится   к двум  линейным   уравнениям (4)  и   (7)  с неизвестными функциями у1 и у2 второго и первого порядка. Интегрируя уравнение (4), находим:

	y1=C1e−x+C2e2x	                        (8)

Подставляя это значение у1 в (7), получаем:

                        y΄2+y2=−C1e−x+2C2e2x+C1e−x+C2e2x,	                      (9)

или

	                                   y΄2+y2=3C2e2x	                         (10)

откуда

                                  y2=C3e−x+c2e2x                                                      (11)

Теперь находим у3:

                      y3= y΄1−y2=−C1e−x+3C2e2x−C3e−x−C2e2x .                      (12)



Общее решение системы (2) имеет вид: 

    

Метод Даламбера. Знание k (k<n) независимых первых интегралов нормальной системы п-го порядка дает возможность понизить порядок этой системы на k единиц. Если же мы знаем п независимых первых интегралов, то мы имеем общий интеграл.

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

Рассмотрим этот метод в случае линейной системы двух уравнений:

    

Умножим второе уравнение на некоторое число k и сложим почленно с первым. Получим:



Или





Выберем k так, чтобы 

                                                                               (17)

 Или

       a12+ka22=k(a11+ka21).                                                             (18)

тогда уравнение (16) можно переписать в виде

        

Это линейное уравнение первого порядка с искомой функцией y+kz. Интегрируя его, найдем:



          y+kz=e(a11+ka21)x{C+∫[f1(x)+kf2(x)]e−(a11+ka21)xdx}.                  (20)



если корни уравнения (18)различные и вещественные, то, обозначив их через  k1, k2, будем иметь два первых интеграла в неявной форме:





Разрешая систему (21) относительно С1 и С2, найдем общий интеграл системы (14), а разрешая относительно y и z, найдем общее решение этой системы.

Если корни уравнения (18) кратные: k1=k2, то формула (20) дает только один первый интеграл:

          y+k1z=e(a11+k1a21)x{C+∫[f1(x)+k1f2(x)]e−(a11+k1a21)xdx}.                      (22)

Подставляя значение г/, найденное отсюда, во второе уравнение системы (14), получим одно линейное уравнение первого порядка с неизвестной функцией z.

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



       

Составляем уравнение для k:

4+5k=k(5+4k).                                                 (24)

Отсюда k1,2=±1. поэтому первыми интегралами будут:        

                                                         

Решая систему (25) относительно y и z, получим общее решение системы (23).



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

       

Здесь мы имеем:

	4 — 2k = k{2 — k),  k2 — \k + 4 = 0.                     (27)

	

Это уравнение имеет двукратный корень k1,2 =2.   Поэтому   метод Даламбера дает возможность найти только один первый интеграл.

Умножая второе из уравнений (26) на 2 и складывая почленно с первым, получаем 



Отсюда находим первый интеграл системы (26)



Используя этот первый интеграл, мы можем переписать второе из уравнений (26) в виде



Отсюда

   

Поэтому







Общее решение системы (26) имеет следующий вид







   



ЛИНЕЙНЫЕ  СИСТЕМЫ  С   ПОСТОЯННЫМИ КОЭФФИЦИЕНТАМИ. СОДЕРЖАЩИЕ     ПРОИЗВОДНЫЕ                                                         ВЫШЕ ПЕРВОГО  ПОРЯДКА

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

Пример. Проинтегрировать систему :



Эта система приводится к одному уравнению четвертого порядка:

         y(4)−k4y=0                                                                                  (2)

        Отсюда:

        y=C1ekx+C2e−kx+C3coskx+C4sinkx                                              (3)

       Поэтому:                    ввв

	Метод Даламбера. 

	  Метод   Даламбера,   изложенный   ранее, распространяется и на линейные системы уравнений, содержащие производные выше первого порядка.

Рассмотрим систему двух уравнений:

   

Умножая второе уравнение на k, складывая почленно с первым уравнением и выбирая k из условия:

       a12+ka22=k(a11+ka21)                                             (6)

получаем:

   

Это линейное неоднородное уравнение с постоянными коэффициентами относительно y+kz. Интегрируя его, найдем:

    y+kz=φ(x, C1, C2, … , Cn).                                         (8) 

если корни уравнения (6) различные относительно y и z, получим общее решение системы (5).

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





1