- •ГЛАВА 2. ВТОРОЙ МЕТОД ЛЯПУНОВА
- •Под вторым методом Ляпунова понимают совокупность приемов и средств исследования устойчивости решений систем дифференциальных уравнений при помощи специальных функций Ляпунова.
- •§ 1. Знакоопределенные функции
- •Рассмотрим вещественную непрерывную числовую функцию
- •Следует заметить, что в качестве W(x) из приведённого определения иногда удается использовать функцию вида
- •О п р е д е л е н и е 2. Положительно или отрицательно определенную функцию именуют знакоопределенной.
- •Пример 1. Рассмотрим функцию
- •Если |a| < 1, то она является положительно определенной, так как
- •причем
- •При a = 1 функция V1 имеет вид
- •З а м е ч а н и е 1. Если функция V = V(x) непрерывна в начале координат и V(0) = 0, то, очевидно, V допускает б.м.в.п.
- •З а м е ч а н и е 2. Для того чтобы установить, что функция V(t, x) имеет б.м.в.п., достаточно убедиться в существовании непрерывной в точке х = 0 функции W(x), удовлетворяющей следующим двум условиям:
- •Пример 2. Функция
- •на основании замечания 2 допускает б.м.в.п. в силу
- •тогда как функция
- •Упражнения
- •1) Выяснить, являются ли данные функции положительно определенными:
- •2) Определить, допускают ли указанные функции б.м.в.п.:
- •§ 2. Теорема Ляпунова об устойчивости
- •Рассмотрим систему дифференциальных уравнений
- •Будем считать, что система (1) имеет нулевое решение ξ = 0, т.е.
- •есть градиент функции V.
- •Пример 1. Пусть
- •Тогда производная функции V в силу данной системы имеет вид
- •Введем сферу
- •где произвольно выбранное ε удовлетворяет неравенствам 0 < ε < H. Очевидно, Sε – непустой компакт и функция W непрерывна на нем. Введём
- •Рассмотрим функцию
- •По условию теоремы
- •что невозможно.
- •если только ||x(t0)|| < δ. Полученное означает, что нулевое решение ξ = 0 системы (1) устойчиво по Ляпунову ■
- •Для линейной однородной системы устойчивость нулевого решения влечет ограниченность каждого решения этой системы. Поэтому имеет место следующее
- •существует положительно определенная функция V(t, x), для которой производная в силу системы неположительна, то все решения этой системы определены и ограничены при t ≥0.
- •Пример 1. Рассмотрим систему
- •и исследуем на устойчивость её нулевое решение х = 0, y = 0. Примем
- •Имеем V > 0 для всех x,y, одновременно не обращающихся в нуль, и
- •Это говорит о том, что функций, удовлетворяющих условия теоремы 1, может быть бесконечно много.
- •Проблема: построить функцию Ляпунова V(t, x), т.е. такую функцию, которая удовлетворяет условиям теоремы 1.
- •Упражнения
- •3) С помощью теоремы 1 установить устойчивость нулевого решения уравнения
- •§ 3. Теорема Ляпунова об асимптотической устойчивости
- •Вновь обратимся к системе
- •По условию теоремы
- •С учетом обозначения
- •Поэтому полученная выше оценка сверху для v(t) принимает вид
- •Отсюда при достаточно большом t следует
- •В итоге установлено равенство
- •Следовательно, предельное соотношение (#) истинно ■
- •Следствие 1. Если для линейной однородной системы
- •Заметим, что равенства нулю всех частных производных функции V вытекают из условия положительной определённости и предположения о непрерывной дифференцируемости этой функции в начале координат (в этом можно легко убедиться, предположив противное).
- •Рассмотрим систему
- •Пример 2. Исследуем на устойчивость нулевое решение уравнения
- •Если a = 0, то, очевидно, нулевое решение устойчиво, но не является асимптотически устойчивым.
- •где m+1 – четное. Следовательно, если a < 0, то функция V удовлетворяет всем требованиям последней теоремы, а значит нулевое решение асимптотически устойчиво.
- •Упражнения
- •1) С помощью надлежащей функции Ляпунова установить асимптотическую устойчивость нулевого решения уравнения
- •2) Подобрать функцию Ляпунова и с её помощью установить асимптотическую устойчивость нулевого решения системы
- •§ 4. Теоремы Ляпунова о неустойчивости
- •10. Первая теорема Ляпунова о неустойчивости.
- •Теорема 1. Пусть для системы
- •1) V(t,x) допускает б.м.в.п.;
- •Тогда нулевое решение системы (1) неустойчиво.
- •□ Согласно условию 3) имеем
- •Пусть
- •Это означает, что нулевое решение системы (1) неустойчиво по Ляпунову ■
- •20. Вторая теорема Ляпунова о неустойчивости.
- •Упражнения
- •2) Проверить, что указанные функции являются функциями Ляпунова для систем
- •имеющих неустойчивое нулевое решение.
- •§ 5. Экспоненциальная устойчивость
- •выполнено неравенство
- •Как показывает следующий результат, экспоненциальная и асимптотическая устойчивость в классе линейных однородных систем с постоянной матрицей совпадают.
- •Лемма 2. Линейная однородная система
- •с постоянной матрицей А асимптотически устойчива тогда и только тогда, когда эта система экспоненциально устойчива (т. е. экспоненциально устойчивым является всякое её решение).
- •Всякое решение линейной системы (4) может быть представлено в виде
- •Введём числа
- •Очевидно,
- •Из курса линейной алгебры известно, что
- •Далее
- •А значит,
- •с непрерывной на промежутке [0, +∞) матрицей A(t). Для таких систем нередко вводят несколько более «жёсткое» понятие экспоненциальной устойчивости, чем это было сделано выше.
- •Нетрудно понять, что из экспоненциальной устойчивости системы (5) вытекает её асимптотическая устойчивость (но не обратно). В случае, когда матрица A(t) = A постоянна, оба указанных понятия эквивалентны.
- •В следующей теореме представлен критерий экспоненциальной устойчивости системы (5).
- •Поскольку
- •Рассмотрим подынтегральное выражение из (9):
- •Используя условие 1), для него можно записать
- •В таком случае из равенства (9) вытекают неравенства
- •Аналогичным образом поступим и с правым неравенством:
- •Из последних двух неравенств легко вывести (6) ■
- •которое именуют матричным уравнением Ляпунова.
- •имела единственное решение относительно P и чтобы это решение представляло собой симметричную положительно определённую матрицу.
- •Упражнения
- •1) Проверить, являются ли экспоненциально устойчивыми нулевые решения следующих уравнений
- •§ 6. Устойчивость квазилинейных систем
- •Линейная система
- •называется системой линейного приближения системы (1).
- •то нулевое решение системы (1) асимптотически устойчиво по Ляпунову.
- •Рассмотрим функцию
- •С учётом (*) для нее получаем представление
- •Получим
- •Напомним, что x – вектор-столбец. Благодаря
- •имеем
- •если только
- •Из доказанной теоремы вместе с теоремой 1 предыдущего параграфа получаем следующее утверждение.
- •Следствие 1. Рассмотрим нелинейную систему
- •Следствие 2. Рассмотрим нелинейную автономную систему
- •и можно применить доказанную выше теорему 1 ■
- •20. Теорема Ляпунова о неустойчивости квазилинейной системы.
- •Теорема 2. Пусть дана нелинейная система
- •□ Пусть
- •Выполним замену
- •где вектор y, вообще говоря, комплексный. Имеем
- •Поэтому
- •Положим
- •и перепишем систему (*) следующим образом
- •Переходим к комплексно-сопряженной системе:
- •Очевидно,
- •Поэтому для функции
- •получаем представление
- •Полагая
- •Если теперь вернуться к исходной переменной, то получим
- •откуда
- •Пример 1. Исследовать на устойчивость нулевое решение системы второго порядка
- •где a – const. Воспользовавшись формулой Маклорена, найдем разложения правых частей системы и перепишем её в виде квазилинейной системы
- •Упражнения
- •Исследовать на асимптотическую устойчивость нулевые решения следующих систем:
- •§ 7. Теорема Зубова
- •Рассмотрим автономную дифференциальную систему
- •в которой
- •Относительно вектор-функции f будем предполагать следующее:
- •□ Ограничимся доказательством лишь части «достаточность» (доказательство второй части можно найти в [7]). По условию функция V(x) – отрицательно определена, а её производная по времени в силу системы (1) положительно определена, так как
- •С этой целью выполним замену, осуществляющую изменение переменной времени по правилу
- •где правые части уравнений непрерывны и ограничены относительно x, так как
- •вытекающим из условия 3) доказываемой теоремы ■
- •З а м е ч а н и е 1. Уравнение в частных производных
- •из последнего условия теоремы 1 называют уравнением Зубова.
- •З а м е ч а н и е 2. На основе теоремы 1 разработан метод Зубова построения оценки области асимптотической устойчивости, общая идея которого состоит в следующем. Пусть система (1) представима в виде
- •где A − постоянная матрица,
- •Существует вычислительный алгоритм построения оценок области притяжения для случая, когда в качестве V(x) выбрана положительно определенная квадратичная форма.
- •2) отрицательно определена;
- •4) удовлетворяет равенству
- •Как видим, все условия теоремы Зубова выполнены, поэтому множество A является областью притяжения нулевого решения.
- •Упражнения
- •1) Построить область асимптотической устойчивости для уравнения
- •2) Построить область асимптотической устойчивости для системы
ГЛАВА 2. ВТОРОЙ МЕТОД ЛЯПУНОВА_______________________________________________________90
Отсюда в результате интегрирования в пределах от t до t получаем
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
V(t, x(t)) ≥ V(t, |
x) + γ(t − t) , |
|
|
|
|||||||||||||
что противоречит установленной ранее ограниченности функции V(t, x) в об- |
||||||||||||||||||||||||||
ласти |
{(t,x) |t ≥ 0, |
|
|
|
x |
|
|
|
≤ 0 }. |
|
|
|
|
|
|
|
|
|
и произвольно выбранного |
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
|
Таким образом, для фиксированного |
0 > 0 |
|
|||||||||||||||||||||||
δ (0, 0 ) указано ненулевое решение |
x(t) , для которого при некотором t1 > |
|
|
|||||||||||||||||||||||
t |
||||||||||||||||||||||||||
имеют место неравенства |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
0 < |
|
|
|
< δ, |
|
|
x(t1) |
|
> |
0 . |
|
|
|
||||
|
|
|
|
|
|
|
|
|
x(t) |
|
|
|
|
|
|
|||||||||||
Это означает, что нулевое решение системы (1) неустойчиво по Ляпунову ■ |
||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
& |
|
= a x |
m |
. Теперь рассмотрим случай, |
||||||||
|
Пример 1. Вернемся к уравнению x |
|
|
|||||||||||||||||||||||
когда |
a > 0 и m – нечетное число. Функция Ляпунова V (x) = x2 |
при a > 0 |
||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
удовлетворяет требованиям теоремы 1. Следовательно, нулевое решение данного уравнения неустойчиво. К аналогичному результату приходим и в случае, когда m – четное (в этом случае следует взять V2 (x) = a x ).
20. Вторая теорема Ляпунова о неустойчивости. |
|
|
|
|
|
|
|
|
|
||||||||||||
Теорема 2. Пусть существует функция |
V(t,x) C(1,1)tx (Ζ), обладающая |
||||||||||||||||||||
следующими свойствами: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
1) |
функция V(t, x) |
ограничена сверху в некоторой |
0 -окрестности |
||||||||||||||||||
|
точки x = 0, т.е. существует M > 0, такое, что V(t, x) ≤ M |
|
|
|
для всех |
||||||||||||||||
|
t ≥ 0 и всех x, удовлетворяющих ; |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
2) |
существует число |
|
|
≥ 0 такое, что для любого δ > 0 |
найдётся точ- |
||||||||||||||||
|
t |
||||||||||||||||||||
|
ка x , удовлетворяющая неравенствам |
|
x |
|
|
|
x) > 0 ; |
|
|
|
|
|
|
|
|
||||||
|
|
|
< δ, V(t, |
|
|
|
|
|
|
|
|
||||||||||
3) |
производная функции V(t, x) в силу системы (1) представима в виде |
||||||||||||||||||||
|
dV |
|
|
= λV(t, x) + W(t,x) |
|
для всех (t, x) Z , |
(2) |
||||||||||||||
|
|
|
|||||||||||||||||||
|
dt |
|
(1) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
где λ – некоторая положительная константа, а W(t, x) является не- |
||||||||||||||||||||
|
которой непрерывной неотрицательной функцией при t ≥ 0 , |
|
|
|
x |
|
|
|
≤ H . |
||||||||||||
|
|
|
|
|
Тогда нулевое решение ξ = 0 системы (1) неустойчиво по Ляпунову.
□ Предположим напротив, что нулевое решение системы (1) устойчиво.
Тогда для любого ε (0, |
0 ) и t = |
|
можно подобрать число δ > 0 такое, что |
||||||||||||
t |
|||||||||||||||
из неравенства |
|
|
|
< δ |
всякий раз будет следовать |
||||||||||
x(t) |
|||||||||||||||
|
|
|
|
|
x(t; |
|
|
|
|
|
t ≥ |
|
. |
||
|
|
|
|
|
t, x(t)) |
< ε < 0 |
t |
ГЛАВА 2. ВТОРОЙ МЕТОД ЛЯПУНОВА_______________________________________________________91
В таком случае в соответствии с условием 1) доказываемой теоремы вдоль решения x(t; t, x(t)) функция V(t, x) ограничена сверху при всех t ≥ t .
Благодаря условию 2) существует точка x , x < δ, для которой выполнено V(t, x) > 0 . Функция V = V(t, x(t; t, x)) в силу условия 3) удовлетворяет линейному дифференциальному уравнению
dVdt = λV + W
с начальным условием V(t, x(t; t, x)) = V(t, x) > 0 . Интегрируя это уравнение в пределах от t до t, с учётом неотрицательности функции W(t, x) получим неравенство
V(t, x(t; t, x) ≥ V(t, x) eλ (t−t) t ≥ t .
Это неравенство вместе с условием λ > 0 противоречит установленной в начале доказательства ограниченности сверху функции V(t, x) при всех t ≥ t ■
Упражнения
1) С помощью надлежащей функции Ляпунова установить неустойчивость ну-
левого решения уравнения |
& |
−cos |
4 |
x (t ≥ 0). |
x =1 |
|
2)Проверить, что указанные функции являются функциями Ляпунова для систем
i)x& = x5 + y3 ,
y& = x3 + y5 ,
|
|
& |
4 |
, |
|
|
|
|
|
|
|
|
ii) |
x = x |
|
|
|
|
|
|
|
|
|
||
|
|
|
2 |
|
2 |
|
|
2 |
|
|
|
|
|
|
& |
|
y |
− x |
, |
|
|
||||
|
y = 2x |
|
|
|
|
|
||||||
|
x = 2y −xy, |
|
|
|
|
|
||||||
iii) |
|
& |
|
|
|
|
|
|
|
|
|
|
& |
|
+ x |
2 |
+2y |
2 |
, |
||||||
|
y = 2x |
|
|
имеющих неустойчивое нулевое решение.
( V = x4 − y4 )
(V = x y)
(V = x + y)
§ 5. Экспоненциальная устойчивость
10. Экспоненциальная устойчивость решений систем общего вида.
Рассмотрим систему дифференциальных уравнений
ГЛАВА 2. ВТОРОЙ МЕТОД ЛЯПУНОВА_______________________________________________________92
|
dx |
= f(t,x) |
(t ≥ 0), |
(1) |
|
dt |
|||
|
|
|
|
|
где f(t, 0) = 0 для всех t ≥ 0, причем f(t, x) |
непрерывна в Z = {(t, x) | t ≥ 0 , |
||x|| ≤ H} при некотором H > 0.
О п р е д е л е н и е 1. Нулевое решение ξ = 0 системы (1) называют экспоненциально устойчивым, если существует такое h > 0, что для любого решения x = x(t) этой системы с начальными данными в области
Z0 = {(t0 , x0 ) | t0 ≥ 0, x0 < h < H }
выполнено неравенство
|
|
|
|
|
|
|
|
|
|
|
|
x(t) |
|
≤ Ν |
|
x(t0 ) |
|
e |
−α(t−t0 ) |
|
t ≥ t0 , |
|
|
(2) |
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||
в котором некоторые положительные числа |
N и |
α |
не зависят от выбора x(t). |
|||||||||||||||||||||||||||||||||||||||
Лемма 1. Из экспоненциальной устойчивости следует асимптотическая |
||||||||||||||||||||||||||||||||||||||||||
устойчивость нулевого решения системы (1). |
|
|
|
|
|
ε > 0 . |
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||
□ В самом деле, выберем произвольное |
Для любых начальных |
|||||||||||||||||||||||||||||||||||||||||
данных |
(t |
0 |
, x |
0 |
) Z |
0 |
, |
удовлетворяющих неравенству |
x |
0 |
< |
ε |
N |
, |
рассмотрим |
|||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
решение |
x = x(t) = x(t; t0 , x0 ) , x(t0 ) = x0 , |
|
и положим |
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
δ:= |
ε |
. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Ν |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
Тогда для введённого δ из (2) вытекает |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
x(t) |
|
|
|
≤ Ν |
ε |
e−α(t−t0 ) < ε |
|
|
|
|
|
t ≥ t0 . |
|
|
|
|
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||
|
|
|
|
|
|
|
|
Ν |
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ξ = 0 |
|
|
устойчиво по Ляпунову. Более того, |
||||||||||||||||||||
Следовательно, нулевое решение |
|
|
||||||||||||||||||||||||||||||||||||||||
имеет место равенство |
|
|
lim x(t) = 0 ■ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
t→+∞ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
Аналогичным образом определяется экспоненциальная устойчивость лю- |
||||||||||||||||||||||||||||||||||||||||||
бого (в том числе ненулевого) решения. А именно, решение |
ξ(t) |
системы (1) |
||||||||||||||||||||||||||||||||||||||||
экспоненциально устойчиво, если при некотором |
h |
> |
0 |
|
любое решение |
|||||||||||||||||||||||||||||||||||||
x(t) = x(t; t0 , x0 ) , такое, что |
t0 ≥ 0, |
|
x(t0 ) −ξ(t0 ) |
|
|
|
|
< h < H , |
удовлетворяет нера- |
|||||||||||||||||||||||||||||||||
|
|
|
||||||||||||||||||||||||||||||||||||||||
венству |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|| x(t) −ξ(t)|| ≤ Ν ||x(t0 ) |
|
− ξ(t0 ) || e |
−α(t − t0 ) |
|
t ≥ t0 |
|
(3) |
||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||
при некоторых положительных постоянных |
N и |
α. |
|
|
|
|
|
|
|
|
Как показывает следующий результат, экспоненциальная и асимптотическая устойчивость в классе линейных однородных систем с постоянной матрицей совпадают.
ГЛАВА 2. ВТОРОЙ МЕТОД ЛЯПУНОВА_______________________________________________________93
Лемма 2. Линейная однородная система
dx |
= A x |
(4) |
|
dt |
|||
|
|
с постоянной матрицей А асимптотически устойчива тогда и только тогда, когда эта система экспоненциально устойчива (т. е. экспоненциально устойчивым является всякое её решение).
□ Достаточность. Из экспоненциальной устойчивости системы (4) следует экспоненциальная устойчивость её нулевого решения. В соответствии с леммой 1 нулевое решение системы (1) асимптотически устойчиво. Тогда в силу результатов § 2 главы 1 асимптотически устойчивой является и система (4).
Необходимость. Ранее было установлено (см. § 4 главы 1), что нулевое решение системы (4) асимптотически устойчиво тогда и только тогда, когда все
собственные значения λj матрицы |
|
А |
имеют отрицательные вещественные |
|||||||
части, т.е. Re λj < 0, j = 1,2,...,n . Положим (т.е. введём |
α) |
|
||||||||
max Re λj |
< −α < 0 . |
|
|
|||||||
j |
|
|
|
|
|
|
|
|||
Тогда (см. гл. 1, § 12 [2]) для матричной экспоненты при некотором N > 0 |
||||||||||
справедлива оценка |
|
|
|
|
|
|
|
|||
|
|
|
etA |
|
|
|
≤ N e−αt |
t ≥ t0 . |
(*) |
|
|
|
|
|
Всякое решение линейной системы (4) может быть представлено в виде
x(t) = e(t−t0 ) A x(t0 ) ,
где начальный момент времени t0 произволен. Следовательно, с использованием (*) получаем
|| x(t) ||=|| e(t−t0 ) A x(t0 ) ||≤|| e(t−t0 ) A || || x(t0 ) ||≤
≤ N || x(t0 )|| e |
−α(t − t0 ) |
t ≥ t0 . |
|
||
Если в это неравенство вместо x(t) |
подставить разность x(t) − ξ(t) , где ξ(t) |
– произвольное решение системы (4) (такая подстановка правомерна, так как данная разность также является решением системы (4)), то получим (3). Полученное будет означать экспоненциальную устойчивость системы (4) ■
Пример 1. Если A = A(t) , то лемма 2 может быть неверна. Рассмотрим, например, скалярное линейное дифференциальное уравнение