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

Математика.-6

.pdf
Скачиваний:
2
Добавлен:
05.02.2023
Размер:
2.78 Mб
Скачать

 

d

 

(x, ) (x, )

 

(x, ) (x, )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

2

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dx

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Заметим,

что, в

силу леммы

Адамара,

lim

1

 

 

f

,

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

lim

2

f

. Поэтому, переходя в полученном

выше

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

уравнении к пределу при 0 , имеем

d

 

 

 

f

 

 

 

 

 

 

 

 

 

 

 

 

dx

 

 

 

 

 

f

.

(4.6)

 

 

 

 

 

 

Уравнение (4.6) называется уравнением в вариациях, а в теории автоматического управления уравнением чувствительности. Обоснование справедливости предельного перехода, а также непрерывности и дифференцируемости соответствующих функций можно посмотреть в [14].

4.2. Определение устойчивости по Ляпунову

В предыдущем пункте мы занимались зависимостью решения от начальных данных в том случае, когда решение определено на ограниченном отрезке. В данном пункте мы будем исследовать этот вопрос в случае, когда решение определено на бесконечном полуинтервале [a, ) . Большой

вклад в развитие данного вопроса внёс А.М.Ляпунов. Более подробное изложение можно найти в [10, 13, 14, 15, 18]

и других пособиях.

 

y 0 (x) y 0 (x), y

 

 

(x) T

 

Определение. Пусть

0

(x),..., y 0

 

 

 

 

 

1

 

 

2

n

 

 

решение системы уравнений (3.1)

 

 

 

 

 

 

 

f1 (x, y1 , y2 ,..., yn ),

 

 

 

 

y1

 

 

 

 

y

f

 

(x, y , y

 

,..., y

 

),

 

 

 

2

 

2

1

2

 

n

 

 

 

 

 

.......... .......... .......... .......

 

 

 

 

y

f

n

(x, y , y

2

,..., y

n

).

 

 

 

n

 

1

 

 

 

 

 

 

91

с начальными данными

y(x

0

) y 0

y 0

, y 0

,..., y 0

T

,

 

 

 

1

2

n

 

 

определённое на полуинтервале [a, ) . Это решение

будем называть устойчивым,

если для всякого 0

существует 0 ,

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

x1 x0

 

,

 

y1 y 0

 

 

yi1 yi0 2

 

( x1 [a, ) )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i 1

 

 

 

 

 

 

 

 

для любого решения y1 (x) y1

(x) y1

(x),...,

y1 (x) T

сис-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

2

 

 

n

 

темы

 

 

уравнений

(3.1)

с

начальными

данными

 

y(x ) y

y1 , y1 ,..., y1

T ,

определённого

также

на

1

1

1

2

 

 

n

 

 

 

 

 

 

 

 

 

полуинтервале

 

[a, ) ,

выполнено

неравенство

 

y1 (x) y0 (x)

 

, i 1,2,..., n , для всех

x из полуинтер-

 

 

 

i

 

 

 

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

вала

 

 

 

 

[a, ) .

 

 

 

 

Если

 

 

при

 

этом

 

lim

y1

(x) y 0 (x) 0, i 1,2,..., n ,

 

 

то

 

решение

 

x

i

 

 

 

 

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y 0 (x) y 0 (x), y 0

(x),..., y 0 (x) T

называется

 

асимптоти-

 

 

 

1

2

 

 

 

 

n

 

 

 

 

 

 

 

 

чески устойчивым.

 

 

 

 

 

 

 

 

 

Исследование устойчивости решений можно свести к исследованию устойчивости нулевого решения с нулевыми начальными данными некоторой другой системы диффе-

ренциальных уравнений. Покажем, как это сделать. Пусть

y 0 (x) y 0 (x), y 0

(x),..., y 0 (x) T

решение системы уравне-

1

2

n

 

 

 

 

ний (3.1)

с

начальными

данными

y 0 (x

0

) y0 ,

 

 

 

 

 

 

y(x) y1 (x), y2 (x),..., yn (x) T – любое решение этой же сис-

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

92

вектор-функции

y(x) y1 (x), y2 (x),..., yn (x) T

 

скалярного

аргумента вычисляется по формуле

 

 

 

y (x), y

2

(x),..., y

n

(x) T y

(x), y

(x),...,

y

(x) T .

1

 

 

1

2

 

n

 

Подробнее о дифференцировании вектор-функций можно

посмотреть

в

[3].

Рассмотрим

вектор-функцию

z(x) y(x) y0 (x) .

 

Тогда

y(x) z(x) y0 (x) ,

 

 

 

 

0 (x) – решения сис-

y (x) z (x) y0 (x) . Так как y(x) и y

темы уравнений (3.1), то подставляя их в эту систему, по-

лучаем

y (x) z (x) y0

 

f (x, y z) .

Поэтому

(x)

z f (x, y0 ) f (x, y z)

или,

что то

же самое,

z f (x, y z) f (x, y0 ) . Рассмотрим теперь решение полученного уравнения с начальными данными z(x0 ) 0 . Име-

ем z(x

0

) y(x

0

) y0 (x

0

) y(x

0

) y

0

0 . Поэтому

y(x

0

) y

0

 

 

 

 

 

 

 

 

 

 

и, если для системы дифференциальных

уравнений

(3.1) y f (x, y)

выполнены условия теоремы существова-

ния и единственности, то y(x)

совпадает с y0 (x) . Следова-

тельно z(x) 0 .

 

 

 

 

 

 

 

 

 

 

 

 

Определение. Решение y(x) 0 системы дифференциальных уравнений (3.1) с начальными данными y(x0 ) 0 будем называть точкой покоя этой системы.

Определение устойчивости точки покоя формулируется следующим образом.

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

выполнении

условий

 

x

 

,

 

y0

, i 1,2,..., n

 

 

 

 

 

 

1

 

 

 

i

 

( x1 [a, ) )

для

любого

 

решения

y(x) y1 (x), y2 (x),..., yn (x) T

системы уравнений (3.1)

93

с начальными данными

y(x ) y 0

y 0

, y 0

,..., y 0

T

,

 

1

1

2

n

 

 

определённого также на полуинтервале [a, ) , вы-

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

 

yi (x)

 

, i 1,2,..., n

для всех x

 

 

из

полуинтервала

 

[a, ) . Если

при этом

lim

yi (x) 0, i 1,2,..., n ,

 

то точка покоя

называется

x

 

 

 

 

 

 

асимптотически устойчивой.

Имеются различные методы исследования устойчивости невозмущённого движения, в том числе и точки покоя.

4.3. Метод функций Ляпунова

Определение. Функцию V ( y) V ( y1 , y2 ,..., yn ) назо-

вём не отрицательно определённой в содержащей на-

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

области

D Rn ,

если для

всякого

y D

V ( y) 0 .

Если, кроме того,

V ( y) 0

тогда и

только

тогда,

когда

y 0 ,

то

функция

V ( y) V ( y1 , y2 ,..., yn ) называется положительно опре-

делённой.

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

Теорема 4.3 (Ляпунова об устойчивости). Если для системы дифференциальных уравнений (3.1)

 

 

 

 

 

 

 

f1

(x, y1 , y2 ,..., yn ),

 

 

 

 

 

 

 

y1

 

 

 

 

 

 

 

y

f

 

(x, y , y

 

,..., y

 

),

 

 

 

 

 

 

 

2

 

 

2

 

 

1

2

 

n

 

 

 

 

 

 

 

 

.......... .......... .......... .......

 

 

 

 

 

 

 

y

f

n

(x, y , y

2

,..., y

n

),

 

 

 

 

 

 

 

n

 

 

 

 

1

 

 

 

существует

положительно

определённая

функция

V ( y) V ( y1 , y2 ,..., yn )

 

 

 

 

такая,

 

что

функция

 

n

V

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

W

 

y

 

f

1

(x, y , y

2

,...,

y

n

)

не положительно опреде-

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

i 1

 

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

94

лена, то точка покоя устойчива. Если функция W отрицательно определена, то точка покоя асимптотически устойчива.

Общего метода построения функций Ляпунова нет. Часто их берут в виде суммы чётных степеней координат.

Пример 1. Исследовать на устойчивость точку покоя системы

dx y,dt

dy x.

dt

В качестве функции Ляпунова возьмем положительно определенную функцию V (x, y) x2 y2 . Чтобы воспользоваться теоремой Ляпунова, составляем функцию

WV f1 (x, y1 , y2 ,..., yn ) . Имеем:

i 1 yin

W V f

1

V f

2

2x y 2y x 0 .

x

y

 

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

Пример 2. Исследовать на устойчивость точку покоя системы

dx y x 3 ,dt

dy x 3y 3 .

dt

В качестве функции Ляпунова возьмем положительно определенную функцию V (x, y) x2 y2 . Чтобы восполь-

зоваться

 

 

теоремой

Ляпунова, составляем функцию

 

n

V

 

 

 

 

 

 

 

 

 

W

 

y

 

f

1

(x, y , y

2

,..., y

n

) . Имеем:

 

 

 

 

 

 

 

 

1

 

 

 

 

 

i 1

 

i

 

 

 

 

2x y x3 2y x 3y3 2x4 6y 4 .

W V f

1

V f

2

 

 

x

 

 

y

 

 

 

 

 

95

Причем, W 2x4 6y4 0 тогда и только тогда, когда обе координаты нулевые, то есть x y 0 .

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

Пример 3. Исследовать на устойчивость точку покоя системы

dx

 

 

 

 

 

 

2

 

 

2

 

 

 

 

x

2 y x

 

 

y

 

,

dt

 

 

 

 

 

dy

 

 

y

 

x

3

y

 

 

 

 

 

x

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

 

2

 

 

2

 

 

 

 

 

Будем

искать функцию Ляпунова

в виде

V (x, y) ax2

by2 , где a и b произвольные положительные

константы.

Составляем

функцию

WV f1 (x, y1 , y2 ,..., yn ) . Имеем:

i 1 yin

 

 

 

 

 

 

 

y

 

 

x3 y

 

 

W V f

1

V f

2

2ax x 2 y x 2 y 2 2by x

 

 

 

 

 

 

 

 

x

y

 

 

 

2

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2ax2 by2 b 2a 2xy x3 y2

.

 

 

 

 

 

 

 

Первое слагаемое полученного выражения 2ax2

by2

не

положительно при любых значениях координат, причем2ax2 by2 0 тогда и только тогда, когда обе координаты нулевые, то есть x y 0 . А знак второго слагаемого не определен, то есть зависит от значений координат. Следо-

вательно, надо подобрать константы a

и b так, чтобы оно

обнулилось. Пусть b 2a , тогда

W 2ax2 by2

и

V (x, y) ax2 2ay2 . При любом положительном значении константы a функция V x, y положительно определена, а

функция W отрицательно определена. Следовательно, точка покоя данной системы асимптотически устойчива.

96

4.4. Устойчивость линейных систем

Рассмотрим линейную систему

дифференциальных

уравнений

 

y A(x) y b(x) .

(4.7)

Соответствующая однородная система будет иметь вид

y A(x) y .

(4.8)

Очевидно, что система (4.8) имеет тривиальное решение и, в случае непрерывно дифференцируемых коэффициентов, это решение единственно.

Пусть матрица A постоянна, то есть системы (4.7) и (4.8) есть системы с постоянными коэффициентами. Тогда общее решение системы (4.8) может быть записано в виде

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y Cl l e l x ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l 1

 

 

 

 

 

 

 

 

где l , l 1,2,..., n

– собственные числа, а l , l 1,2,..., n – со-

ответствующие им собственные векторы матрицы

A . Со-

ответственно,

y

i

 

e l x

, i 1,2,..., n ,

фундаментальная

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

система решений системы уравнений (4.8).

 

 

 

Исследуем устойчивость точки покоя системы (4.8).

Имеем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

n

 

 

 

e Re l iJm l x

 

 

 

 

y

 

 

 

Cl

 

l

 

e l x

 

 

Cl

 

l

 

 

 

 

 

 

 

 

 

 

 

l 1

 

 

 

 

 

 

 

 

 

l 1

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

max Re l x

 

Cl

 

l

eiJm l x

eRe l x

 

Cl

 

l

eRe l x Me 1 l n

,

l 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l 1

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где M

 

Cl

 

 

 

 

 

l

 

 

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

если действительные

 

 

 

 

 

 

 

 

l 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

части Re l , l 1,2,..., n, собственных

чисел

l , l 1,2,..., n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

max Re l x

0 . Поэтому

матрицы A отрицательны, то lim e 1 l n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

97

точка покоя системы (4.8) асимптотически устойчива, и, следовательно, устойчива на полуинтервале [0, ) .

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

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

dx

1

1

 

 

 

a1 x a2 y,

 

 

dt

 

,

 

dy

 

a12 x a22 y.

 

dt

 

 

с невырожденной матрицей системы, то есть определитель

det A

a11

a12

матрицы системы отличен от нуля.

 

a2

a2

 

 

1

2

 

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

венных чисел матрицы системы

 

A E

 

 

a1

a1

0

 

 

1

2

 

 

 

 

 

 

a2

a2

 

 

 

 

 

 

 

 

1

2

 

и найти его корни 1

и 2 , то есть собственные числа мат-

рицы A .

 

 

 

 

 

 

 

 

Возможны следующие случаи:

 

 

1) Собственные числа 1 и 2 вещественны и различ-

ны, тогда:

 

 

 

 

 

 

 

 

а) если 1 и 2

отрицательны, то точка покоя асим-

птотически устойчива (устойчивый узел);

 

 

б) если 1 и 2

положительны, то точка покоя неус-

тойчива (неустойчивый узел);

 

 

в) если 1 и 2

имеют разные знаки, то точка покоя

неустойчива (седло);

 

 

98

2) Собственные числа 1 и 2 вещественны и кратны, тогда:

а) если 1 2 отрицательны, то точка покоя асимптотически устойчива (устойчивый узел);

б) если 1 2

положительны, то точка покоя неус-

тойчива (неустойчивый узел);

 

3) Собственные

числа 1 и

2 комплексные,

1,2 p qi , тогда:

а)

если

p 0

и

q 0 , то точка покоя асимптотиче-

ски устойчива (устойчивый фокус);

б)

если

p 0

и

q 0 , то точка покоя неустойчива

(неустойчивый фокус);

в)

если

p 0

и

q 0 , то точка покоя устойчива

(центр).

Пример 1. Исследовать на устойчивость точку покоя системы

dx 5x y,dt

dy 2x y.

dt

Составляем характеристическое уравнение:

 

A E

 

 

 

5

1

 

2 6 7 0 .

 

 

 

 

 

 

 

 

 

 

 

2

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

И находим его корни,

1 3 2 ,

2 3 2 . Так как

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

99

4.6. Устойчивость по первому приближению

Рассмотрим систему дифференциальных уравнений (3.1) в нормальной форме

 

 

 

 

 

 

f1 (x, y1 , y2 ,..., yn ),

 

 

 

 

 

y1

 

 

 

 

 

y

 

f

 

(x, y , y

 

,...,

y

 

),

 

 

 

 

 

 

2

 

2

1

2

 

 

n

 

 

 

 

 

 

.......... .......... .......... .......

 

 

 

 

 

y

 

f

n

(x, y , y

2

,...,

y

n

),

 

 

 

 

 

 

n

 

1

 

 

 

 

или, что то же самое, в векторной форме y f (x, y) .

Пусть y(x) 0

– точка покоя системы (3.1). Предполагая,

что функция

f (x, y)

дифференцируема, линеаризуем её

(см. [3]) в окрестности точки y 0,0,...,0 T ,

то есть пред-

ставим в виде

f (x, y) A(x) y b(x, y) , где A(x) производная

матрица

функции

 

 

 

f (x, y)

 

 

 

с

элементами

ai

(x)

fi

, i 1,2,..., n, j 1,2,...n , b(x, y)

– бесконечно малая

 

j

 

y j

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

вектор-функция более

высокого

порядка

малости, чем

y y1, y2 ,..., yn T . Тогда система дифференциальных уравнений (3.1) запишется в виде

y A(x) y b(x, y) .

(4.9)

Отбрасывая в системе (4.9) члены более высокого порядка малости, чем y y1, y2 ,..., yn T , можем переписать её в форме

y A(x) y .

(4.10)

Система дифференциальных уравнений (4.10) называется системой первого приближения для системы дифференциальных уравнений (3.1). Иногда удаётся судить об устойчивости точки покоя системы (3.1) по устойчивости точки покоя для системы (4.10).

Пусть система дифференциальных уравнений (3.1) является автономной, то есть функции fi , i 1,2,..., n не зави-

100