Смекни!
smekni.com

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

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

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

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

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

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

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

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

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

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

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

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

t =y(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,

2=y1+y3, (2)

3=y1+y2.

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

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

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

y˝11−2y1=0; (4)

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

y3= 1y2. (5)

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

2= 1y2+y1 (6)

или

2 +y2= 1 +y1 (7)

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

y1=C1ex+C2e2x (8)

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

2+y2=−C1ex+2C2e2x+C1ex+C2e2x, (9)

или

2+y2=3C2e2x (10)

откуда

y2=C3ex+c2e2x (11)

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

y3= 1y2=−C1ex+3C2e2xC3exC2e2x . (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)