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

Дифференциальные уравнения. Учебное пособие-3

.pdf
Скачиваний:
0
Добавлен:
15.08.2026
Размер:
1 Мб
Скачать

ТЕОРИЯ УСТОЙЧИВОСТИ

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

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

 

dyi

 

(t, y , y

 

, ...,

y

 

),

i

 

.

(1)

 

2

n

1, n

 

 

 

dt

 

i

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Решения i (t),

i

 

, системы(1) устойчивыпоЛяпунову. Если

1, n

 > 0 можно найти ( ) > 0 такое, что для всякого yi (t) – решения (1), начальные условия которого удовлетворяют соотношениям:

yi (t0 ) i (t0 )

( )

yi (t) i (t)

,

i

1, n

t t0 .

(2)

Определение. Если число из предыдущего определения можно выбиратьнезависящимотt, тотакаяустойчивостьназываетсяравно-

мерной в области T.

Отметим, что если система (1) удовлетворяет условиям теоремы о непрерывной зависимости решения от начальных условий, то в определении вместо t t0 может быть записано t0 t tn , т. е. при близких

начальных условиях решения остаются близкими на конечном отрезке от t0 до tn .

Определение. Если при сколь угодно малом  > 0 хоть для одного решения yi (t) второе неравенство из (2) не выполняется, то решение

i (t), i 1, n, называется неустойчивым.

61

Определение. Если i (t) не только устойчиво, но и удовлетво-

ряет соотношению

 

yi (t) i (t)

 

0

при t при условии, что

 

 

yi (t0 ) i (t0 )

 

( ) ,

то решение i (t) называется асимптотически

 

устойчивым.

Замечание. Из асимптотической устойчивости не следует устойчивость решения.

ОБЩИЕ ТЕОРЕМЫ ОБ УСТОЙЧИВОСТИ СИСТЕМ ЛИНЕЙНЫХ ДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИЙ (СЛДУ)

Пусть дана СЛДУ

 

 

 

 

 

dy

A(t) y(t) f (t) ,

(3)

 

dt

 

 

 

 

 

где A(t) и f (t) – непрерывные функции аргумента t ( , ) :

 

 

 

 

dy

A(t) y(t)

(4)

 

 

 

 

 

 

 

dt

 

– соответствующая однородная система.

Определение. СЛДУ (3) называется устойчивой (или вполне не-

устойчивой), если все решения этой системы y(t) соответственно

устойчивы (неустойчивы) по Ляпунову при t .

Замечание. Необходимо отметить, что решения СЛДУ либо все устойчивы, либо все неустойчивы.

Теорема. ДляустойчивостиСЛДУ(3) прилюбойправойчасти f (t)

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

 

Доказательство

 

 

Пусть (t) – некоторое устойчивое

решение системы (1),

t0

t . Значит, > 0 ( ) > 0: y(t) – решения системы (1) при

t0

t будет справедливо неравенство

 

 

 

y(t) (t)

 

 

 

,

(5)

 

 

 

 

62

если выполняется

 

 

 

 

y(t0 ) (t0 )

 

 

 

.

(6)

 

 

 

Но если y и – решения (3), то

 

 

 

 

x(t) y(t) (t)

(7)

– решение однородной системы (4).

Любое решение однородной системы всегда можно представить

в виде (7). В cсоответствии с (5) и (6) получаем, что

 

 

 

x(t)

 

 

 

, если

 

 

 

 

 

x(t0 )

 

 

 

. Это выполняется x(t), значит, выполнится

и для х0 0,

 

 

 

т. е. тривиальное решение однородной системы (4) устойчиво по Ляпунову при t .

Замечание. Из доказательства следует, что устойчивость тривиального решения системы (4) вытекает из устойчивости хотя бы одного решения линейной неоднородной системы (3) при любой правой части f (t), которая может быть даже тождественно равна нулю.

Пусть х0 0 – тривиальное решение однородной системы (4). И пусть х0 устойчиво по Ляпунову при t , тогда, если x(t) – про-

извольное решение однородной системы (4),

 

 

 

x(t0 )

 

 

 

( , t0 ) выполня-

 

 

 

 

ется при условии, что

 

 

 

x(t)

 

 

 

. Пусть (t) ,

 

y(t)

 

– некоторые произ-

 

 

 

 

 

 

вольные решения системы (3), тогда из y(t0 ) (t0 ) будет следовать, что y(t) (t) t0 t . Последние два неравенства означают, что (t) будут устойчивыми по Ляпунову при t .

Следствие. СЛДУ устойчива, когда устойчиво хотя бы одно решениеэтойсистемы, ивполненеустойчива, еслинеустойчивоеенекоторое решение.

Определение. СЛДУ (3) называется равномерно устойчивой, если все решения этой системы равномерно устойчивы при t .

Теорема. СЛДУ (3) равномерно устойчива тогда и только тогда, когда тривиальное решение соответствующей ей однородной системы (4) равномерно устойчиво при t .

Определение. СЛДУ (3) называется асимптотически устойчивой,

если все ее решения асимптотически устойчивы.

63

Теорема. Система (3) асимптотически устойчива тогда и только тогда, когдатривиальноерешениесоответствующейсистемы(4) асимптотически устойчиво при t .

УСТОЙЧИВОСТЬ ОДНОРОДНЫХ СИСТЕМ ЛИНЕЙНЫХ ДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИЙ

Пусть дана однородная СЛДУ (4). Если матрица A(t) непрерывна,

то устойчивость системы (4) эквивалентна ограниченности всех решений СЛДУ (4).

Теорема. Однородная СЛДУ (4) устойчива по Ляпунову тогда и только тогда, когда каждое решение y(t) ограничено на полуоси

t t0 , .

Теорема. Однородная СЛДУ (4) асимптотически устойчива тогда и только тогда, когда все ее решения y(t) 0 при t .

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

Пусть

dx

Ax(t) ,

(8)

dt

 

 

где А – квадратная матрица (n × n) с постоянными коэффициентами. Будем искать решение СЛДУ (8) в виде

x(t) eAtС,

(9)

где С – вектор-константа.

 

x(t0 ) x0 , тогда с учетом началь-

Пусть заданы начальные условия

ных условий решение (9) примет вид

 

 

x(t

0

) x

eAt0 С x(t) eA(t t0 ) x .

(10)

 

0

 

0

 

64

Пусть 1, 2 , ..., m – собственные значения матрицы А, m n и e1, e2 , ..., em – собственные векторы, соответствующие собственным значениям. Тогда матрицу А можно представить в виде

A S diag J1( 1), J2 ( 2 ), ..., Jm ( m ) S 1 ,

где S – ортогональная матрица; Ji – соответствующие клетки Жор-

дана [3].

Тогда решение системы (8) имеет вид

x(t) S diag eJ1( 1)(t t0 ) , eJ

2 ( 2 )(t t0 ) ,

..., eJm ( m )(t t0 ) S 1x(t

0

) ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

exp(Jk ( k )(t t0 ))

 

 

 

 

 

 

 

 

 

 

 

 

(t t0 )

 

 

(t t0 )

2

 

 

 

(t t0 )

k 1

 

 

 

 

exp(

k

(t t

0

)) E

I

 

 

I

2

...

 

I

k 1

,

1!

 

 

(k 1)!

 

 

 

1

2!

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где I j – единичная матрица со сдвинутой на j единиц вверх главной

диагональю.

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

щественные части: Re j 0, j 1, n , причем собственным значениям,

имеющимнулевуювещественнуючасть, соответствуюттолькопростые элементарные делители, т. е. единичная жорданова клетка.

Доказательство

 

 

Пусть собственные значения матрицы A:

j j i j ,

j 1, p

(р – кратность собственных значений), причем j 0 и k i k ,

k 1, q (q – кратность собственных значений). То есть максимальное

число клеток Жордана матрицы А: m = p + q. Решение (8) можно представить в виде

p

q

y(t) e jt (cos jt isin jt)Pj (t)

(cos kt isin kt)Ck , (11)

j 1

k 1

где Pj (t) – полиномиальные векторы; Ck

– векторы постоянных.

65

 

Покажем,

 

что решение (11)

 

ограничено.

Так

как

j 0 ,

 

cos t isin t

 

1, а Pj (t) имеют ограниченную степень, которая ниже

 

 

кратности собственного значения

 

j

, то e jt P (t) 0

при

t .

 

 

 

 

 

 

j

 

 

Ck – ограниченныеконстанты. Такимобразом, | y(t) | M , гдеМ– кон-

станта, т. е. решение ограниченное и, следовательно, устойчивое, а значит, и СЛДУ (8) устойчива.

ПустьСЛДУ (8) устойчива. Покажем, что всесобственные значения матрицы А имеют неположительные вещественные части.

Докажем от противного. Пусть существует такое собственное зна-

чение s

i , у

которого 0. Система (8) однородная, и, следо-

вательно,

частное

решение,

соответствующее s , примет вид

e t (cos t i sin t)С,

 

 

 

C

 

 

 

0.

Решение (t) возрастает неограни-

 

 

 

 

ченно, что противоречит условию теоремы об устойчивости. Значит, наше предположение неверно и Re s 0.

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

т. е.

A S diag J

( ),

J

2

(

2

), ..., J

m

(

m

) S 1,

det(s) 0. И

некото-

 

 

1

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

рому

s с нулевой действительной частью соответствует жорданова

клетка размерности ls 1,

 

тогда частное решение имеет вид

 

 

 

 

 

 

 

 

 

 

 

(t) S diag 0,

0,

..., exp (Js ( s )t), ..., 0, 0 S 1 .

(12)

 

ВсоответствиисопределениемрешенияСДУ S 1 (t)S – тожечаст-

ное решение СДУ (8), его норма:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

S 1 (t)S

 

 

 

 

 

diag 0, 0, ...,

 

exp (Js ( s )t), ...,

0,

0

 

 

 

 

 

exp (Js ( s )t)

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Пусть

 

 

 

A

 

 

 

max

aij

 

, тогда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 i n

j 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

S 1 (t)S

 

 

 

 

 

 

 

S 1

 

 

 

 

 

(t)

 

 

 

 

 

 

 

S

 

 

 

,

 

 

 

 

 

 

(13)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

66

 

exp(Js ( s )t)

 

 

 

 

1

 

t

 

t2

...

 

tls 1

 

 

tls 1

.

(14)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1!

2!

 

(ls 1)!

(ls 1)!

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Воспользуемся (14) и получим оценку:

 

 

 

 

 

 

 

 

 

(t)

 

 

 

 

 

tls 1

 

 

 

.

 

 

 

 

 

(15)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(ls 1)!

 

 

 

S

 

 

 

 

 

S 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Знаменатель соотношения

(15)

ограничен,

 

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

при

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

Замечание. Устойчивая линейная однородная система ДУ с постоянной матрицей А равномерно устойчива относительно начального мо-

мента t0 .

Теорема. ЛинейнаяоднороднаясистемаДУ(8) спостояннойматрицей А асимптотически устойчива тогда и только тогда, когда все собственные значения матрицы А имеют отрицательные вещественные части.

ТОЧКИ ПОКОЯ

Определение. Пусть дана система ДУ

dxi

f

(t,

x , ...,

x

), i

 

.

1, n

 

 

 

 

 

dt

 

i

 

 

1

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Если эта система имеет решение xi xi0 const,

i

 

, то траектория

1, n

этого

решения сходится в n-мерном

пространстве

 

 

к точке

M (x10

, x20 , ..., xn0 ), и эта точка называется точкой покоя для данной

системы ДУ.

 

 

 

 

 

 

 

 

 

 

 

Если интерпретировать решение ДУ

dxi

f

(t, x , ..., x

 

)

как неко-

 

 

 

 

dt

i

 

1

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

торую траекторию в n-мерном пространстве, зависящую от параметра t (времени), то решение системы, состоящее из констант, будет точкой покоя этой системы. Естественно, что тривиальное решение (x0 0) –

точка покоя данной системы.

67

ТИПЫ ТОЧЕК ПОКОЯ

Пусть дана однородная система ДУ с постоянной матрицей, и соответственно вид решений ДУ будет определяться видом собственных значений матрицы. Рассмотрим типы точек покоя на примере двумерного пространства:

 

 

 

 

 

dY

 

 

x

 

 

 

 

 

 

 

AY, Y .

 

 

 

 

 

dt

 

 

 

 

 

 

 

y

I.  Собственные значения 1

, 2

матрицы А – действительны и раз-

личны

( 1 2 ).

Тогда решение

можно записать в виде Y (t)

C e 1te C

2

e 2te .

 

 

 

1

1

 

2

 

 

 

 

1.   1 0; 2 0 – тогда точка покоя будет устойчивым узлом, поскольку все точки траектории, находящиеся в начальный момент t0

в любой -окрестности начала координат, при достаточно большом t стремятся к точкам, принадлежащим сколь угодно малой -окрестности начала координат.

2.   1 0; 2 0 – неустойчивый узел.

68

3.   1 0; 2 0 – тогда точка покоя называется седлом.

II.  Собственные значения – комплексно-сопряженные:

1,2 p iq.

1.   p 0; q 0 – тогда точка покоя называется устойчивым фоку-

сом.

2.   p 0; q 0 – неустойчивый фокус.

69

3.   p 0; 1,2 iq – точка устойчива – центр.

III.   1 2 – кратные корни.

 

 

 

 

 

 

 

Y (t) C ee t C

(et f )e t ,

е – собственный вектор, f – присоеди-

1

2

 

 

 

 

 

 

 

 

 

 

 

ненный вектор.

 

 

 

 

 

 

 

 

 

 

 

1.   0 , тогда при t Y(t) 0 , т. е. точка покоя будет асимп-

тотически устойчивой, это будет устойчивый узел.

 

 

 

 

2.   0 , тогда точка покоя – неустойчивый узел.

 

 

 

 

МЕТОД ЛЯПУНОВА ИССЛЕДОВАНИЯ СДУ

 

 

 

НА УСТОЙЧИВОСТЬ

 

 

 

 

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

dxi

f

(t, x , x

, ..., x ),

i

 

 

. Будем ис-

1,

n

 

 

 

 

 

dt

i

1

2

n

 

 

 

 

следовать на

 

 

 

 

решение

или

точку покоя

устойчивость тривиальное

xi 0, i

 

. Если с ростом t

точки всех траекторий будут прибли-

1, n

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

 

Пусть (x1

, x2 , ..., xn ) – расстояние от некоторой точки траектории

 

 

 

n

 

до

начала координат,

xi2 (t) .

Знак производной функции

 

 

 

i 1

 

 

n

 

 

 

 

xi2 (t)

определяется

знаком

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

i 1

70

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]