Санкт-Петербургский государственный университет
Кафедра управления медико-биологическими системами
Воробьёва Анна Алексеевна
Выпускная квалификационная работа бакалавра
Исследование условий абсолютной устойчивости
некоторых классов нелинейных систем
010400
Прикладная математика и информатика
Научный руководитель,
доктор физ.-мат. наук,
профессор
Александров А. Ю.
Санкт-Петербург
2016
Содержание
1. Введение
3
2. Основные понятия и теоремы
2.1. Система с переключениями . . . . . . . . . . . . . . . . . . . . .
2.2. Система с запаздыванием . . . . . . . . . . . . . . . . . . . . . .
5
5
7
3. Обзор литературы
9
4. Система с переключениями
4.1. Постановка задачи . . . . . . . . . . . . . . . . . . . . .
4.2. Полученные результаты . . . . . . . . . . . . . . . . . .
4.2.1. Необходимые условия абсолютной устойчивости
4.2.2. Достаточные условия абсолютной устойчивости
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
11
11
11
12
15
5. Системы с запаздыванием
5.1. Постановка задачи . . . . . . . . . . . . . . . . . . . . . . . . . .
5.2. Условия существования диагональных функционалов Ляпунова — Красовского . . . . . . . . . . . . . . . . . . . . . . . . . . .
5.2.1. Системы с якобиевыми матрицами . . . . . . . . . . . . .
5.2.2. Условия диагональной устойчивости матриц третьего порядка со специальной структурой . . . . . . . . . . . . .
20
20
6. Заключение
36
2
21
21
28
1. Введение
Одной из важнейших задач современной теории управления является анализ устойчивости различных классов нелинейных систем [11]. В данной работе исследуются условия устойчивости для некоторых классов нелинейных
систем. Предполагается, что функции, входящие в правые части рассматриваемых уравнений, удовлетворяют ограничениям секторного типа [1]. Такие
функции называются допустимыми. Можно выделить две основные задачи.
Первая задача заключается в исследовании условий асимптотической устойчивости нулевого решения для одного класса нелинейных систем с переключениями. Системы с переключениями широко используются в задачах управления технологическими процессами, механическими системами, а так же во
многих других областях [1, 10]. Система такого типа представляет собой гибридную динамическую систему, включающую в себя семейство непрерывных
по времени подсистем, а также закона переключения между ними [1]. Переключения могут быть обусловлены внешними возмущениями или же задаваться управлением. В каждый момент времени активна только одна из подсистем семейства. Зачастую при проектировании управляемой системы необходимо обеспечить её устойчивость для любых допустимых переключений [1, 10].
Данная проблема является актуальной, если изначально закон переключений
неизвестен. Необходимым, но не достаточным условием устойчивости системы
является устойчивость каждой подсистемы семейства [1]. Для исследования
необходимых условий к соответствующим линейным подсистемам применяется критерий Гурвица [4]. Для получения достаточных условий устойчивости
эффективным средством является прямой метод Ляпунова. Существование
общей для всех подсистем функции Ляпунова гарантирует асимптотическую
устойчивость равномерную относительно закона переключения [1].
Во второй задаче рассматриваются условия асимптотической устойчивости
нулевого решения для некоторого класса нелинейных систем с постоянным запаздыванием. Далеко не всегда реальные динамические системы могут быть
описаны системами обыкновенных дифференциальных уравнений, поэтому зачастую при моделировании динамических процессов используются системы с
запаздыванием [2]. Запаздывание возникает тогда, когда динамика системы
зависит не только от ее состояния в текущий момент времени, но и от поло-
3
жения системы в некоторый предыдущий момент времени [2]. Такие системы
используются для моделирования биологических сообществ, применяются в
экологии, экономике, инженерии и многих других областях [2, 3, 9]. Устойчивость систем с запаздыванием также является актуальной проблемой теории
управления. В данной работе рассматриваются как системы с запаздыванием без переключений, так и системы с запаздыванием и с переключениями.
Анализ устойчивости основан на использовании прямого метода Ляпунова и
функционалов Ляпунова — Красовского [2]. Найдены условия существования
диагональных функционалов Ляпунова — Красовского.
Выполнение полученных условий гарантирует, что нулевые решения изучаемых систем будут являться асимптотически устойчивыми при любых допустимых нелинейностях, переключениях и при произвольном неотрицательном
запаздывании.
4
2. Основные понятия и теоремы
В данном разделе приведем основные понятия и теоремы, которыми будем
оперировать.
Определение 1 [1]. Постоянная матрица A называется метцелеровой, если
все её внедиагональные элементы неотрицательны.
Определение 2 [9]. Матрица A называется устойчивой, если она является
гурвицевой, то есть все ее собственные числа имеют отрицательные вещественные части.
Определение 3 [9]. Матрица A называется D-устойчивой, если матрица DA
устойчива при любой положительной диагональной матрице D.
2.1. Система с переключениями
В общем случае нелинейная система дифференциальных уравнений с переключениями представима в виде
ẋ(t) = Fσ (t, x).
(2.1)
Здесь x = (x1 , . . . , xn )T – вектор состояния системы; σ = σ(t) – кусочно-постоянная функция, которая задает закон переключений [1].
Определение 4 [1]. Будем говорить, что закон переключения является допустимым, если он задается кусочно-постоянной функцией σ(t) : [0; +∞) →
Q = 1, . . . , N , где σ(t) является непрерывной справа в точках разрыва, и к тому же любой конечный интервал времени содержит конечное число моментов
переключений.
В таком случае, в каждый момент времени система (2.1) характеризуется
одной из подсистем семейства
ẋ(t) = Fs (t, x),
s = 1, . . . , N.
(s)
(s)
Предполагается, что вектор-функции Fs (t, x) = (F1 , . . . , Fn )T , s = 1, . . . , N,
непрерывны в области
t > 0,
kxk < H
(0 < H 6 +∞).
5
(2.2)
В данной работе будут рассматриваться системы с переключениями вида
ẋ(t) = Pσ f (x(t)),
(2.3)
где f (x) = (f1 (x1 ), . . . , fn (xn ))T ; σ = σ(t)
допустимый закон переключения,
—
n
(s)
σ(t) : [0; +∞) → Q = 1, . . . , N ; Ps = pij
— постоянная матрица для
i,j=1
каждого фиксированного s = 1, . . . , N .
Определение 5 [1]. Будем говорить, что скалярные функции fi (xi ) являются
допустимыми, если они непрерывны при |xi | < H (0 < H 6 +∞) и удовлетворяют условиям секторного типа, то есть xi fi (xi ) > 0 при xi 6= 0, i = 1, . . . , n.
Из свойств функций f1 (x1 ), . . . , fn (xn ) следует, что система (2.3) имеет нулевое решение.
Определение 6 [1]. Система (2.3) абсолютно устойчива, если ее нулевое решение асимптотически устойчиво при любых допустимых переключениях, а
также при любых допустимых нелинейностях.
Для исследования устойчивости системы типа (2.3) в данной работе применяется прямой метод Ляпунова.
Теорема Ляпунова об асимптотической устойчивости [11]. Рассмотрим систему
ẋ = G(t, x),
где вектор-функция G(t, x) определена и непрерывна в области (2.2), причем G(t, 0) ≡ 0 при t > 0, то есть данная система имеет нулевое решение.
Пусть существует заданная и непрерывно дифференцируемая на множестве
(2.2) функция V (t, x), которая обладает следующими свойствами:
1) V (t, x) положительно определена;
2) V (t, x) допускает бесконечно малый высший предел;
3) производная V (t, x) в силу рассматриваемой системы отрицательно определена.
Тогда нулевое решение этой системы асимптотически устойчиво.
Замечание. Здесь и далее под k · k понимается евклидова норма вектора.
Определение 7 [9]. Система (2.3) называется диагонально устойчивой, если
существует положительно-определенная диагональная матрица Λ, такая, что
для любого s = 1, . . . , N выполняется матричное неравенство Ляпунова, то
6
есть
ΛPs + PsT Λ < 0,
s = 1, . . . , N.
2.2. Система с запаздыванием
В данной работе будут рассматриваться системы с постоянным запаздыванием. Система с постоянным запаздыванием представима в виде функционального дифференциального уравнения
ẋ(t) = F (t, x(t), x(t − τ )).
(2.4)
Здесь x = (x1 , . . . , xn )T – n-мерный вектор состояния системы; вектор-функция
F (t, x, y) определена в области t > 0, kxk < H, kyk < H, кусочно-непрерывна
по t и непрерывна по x и y; τ – постоянное неотрицательное запаздывание [2].
Каждое решение x(t, t0 , ϕ) системы (2.4) при t > t0 определяется начальным
моментом времени t0 и начальной функцией ϕ(θ) [2].
Будем считать, что начальные функции принадлежат пространству
C([−τ, 0], Rn ) непрерывных вектор-функций ϕ(θ) : [−τ, 0] → Rn с равномерной нормой kϕkτ = maxθ∈[−τ,0] kϕ(θ)k [2]. Через xt (t0 , ϕ) обозначим отрезок
решения: xt (t0 , ϕ) : θ → x(t + θ, t0 , ϕ), θ ∈ [−τ, 0]. Если начальный момент
времени и начальная функция очевидны из контекста или несущественны, то
аргументы t0 и ϕ в данных обозначениях будем опускать.
Для исследования устойчивости нулевого решения систем типа (2.4) применяется прямой метод Ляпунова и функционалы Ляпунова — Красовского [2,
6].
Определение 8 [13]. Функцией Хана называется непрерывная неубывающая
функция χ(s), такая что χ(s) > 0 при s > 0 и χ(0) = 0.
Теорема Ляпунова—Красовского об асимптотической устойчивости [2]: Если существует непрерывно дифференцируемый функционал
V (t, xt ) : R × C → R, такой, что можно указать функции Хана χ1 (s), χ2 (s),
χ3 (s), такие, что
χ1 (kϕ(0)k) 6 V (t, ϕ) 6 χ2 (kϕkτ ),
7
при t > 0, kxt kτ < H, то нулевое решение системы (2.4) асимптотически
устойчиво.
В данной работе будут рассматриваться системы с запаздыванием вида
ẋ(t) = Af (x(t)) + Bf (x(t − τ )).
(2.5)
Здесь A и B – постоянные матрицы, f (x) = (f1 (x1 ), . . . , fn (xn ))T – заданная и
непрерывая при kxk < H векторная функция.
Будем говорить, что f (x) является допустимой функцией, если она удовлетворяет условиям из определения 4.
Определение 9 [2]. Система (2.5) абсолютно устойчива, если ее нулевое решение асимптотически устойчиво при любом положительном запаздывании, а
также при любых допустимых нелинейностях.
Определение 10 [6].
Система (2.5) диагонально устойчива, если
существуют
диагональные
положительно-определённые
матрицы
P = diag{p1 , . . . , pn } и Q = diag{q1 , . . . , qn }, при которых матрица
S=
AT P + P A + Q P B
BT P
−Q
отрицательно определена.
8
!
3. Обзор литературы
Базовые понятия и методы теории устойчивости систем обыкновенных дифференциальных уравнений подробно рассматриваются в [11]. В [14] вводится понятие абсолютной устойчивости и приводится метод построения функций Ляпунова для исследования абсолютной устойчивости нелинейных систем.
Также в этой книге ставится задача абсолютной устойчивости нелинейных систем с переключениями и приводятся методы решения этой задачи. Задача
исследования диагональной устойчивости в экологических моделях типа Лотки — Вольтерры рассматривается в [9]. Некоторые результаты в исследовании
условий D-устойчивости для матриц размерности 4 × 4 приведены в работе
[12]. В [5] приводятся условия устойчивости для систем вида
ẋ(t) = Pσ f (x),
где σ(t) : [0; +∞) → Q = 1, . . . , N – допустимый закон переключений; f (x)
– допустимые нелинейности, а Ps – метцелерова матрица для любого s =
1, . . . , N , а также получены некоторые результаты для систем с неметцелеровыми матрицами.
В [10] рассматриваются основные проблемы устойчивости систем с переключениями, а именно:
1) найти условия гарантирующие, что переключающаяся система асимптотически устойчива для любого закона переключения;
2) найти такие классы законов переключения, для которых система с переключениями будет являться асимптотически устойчивой;
3) построить закон переключения, который обеспечивает асимптотическую
устойчивость нулевого решения системы.
В данной работе рассматривается первая из этих трех проблем для некоторых классов нелинейных систем.
Для систем с запаздыванием основные методы исследования устойчивости
подробно рассмотрены в [2]. В [6] вводится понятие диагональной устойчивости
для линейных систем с запаздыванием, а также предлагается метод анализа
устойчивости таких систем. В статье [7] получен критерий диагональной устойчивости для систем с запаздыванием, но проблема в том, что в общем случае
9
этот критерий оказывается трудно проверяемым и не дает конструктивных
условий устойчивости.
Основной целью настоящей работы является нахождение таких классов систем, для которых можно получить конструктивно проверяемые условия диагональной устойчивости, а также нахождение этих условий.
10
4. Система с переключениями
4.1. Постановка задачи
Пусть задана система
ẋ(t) = Pσ f (x(t)).
(4.1)
Здесь x(t) = (x1 (t), . . . , xn (t))T ; f (x) = (f1 (x1 ), . . . , fn (xn ))T , скалярные функции fi (xi ) являются допустимыми нелинейностями; σ = σ(t) — допустимый
n
(s)
закон переключения, σ(t) : [0; +∞) → Q = 1, . . . , N ; Ps = pij
— поi,j=1
стоянная матрица для каждого фиксированного s = 1, . . . , N . Тогда в любой
момент времени такая система характеризуется одной из подсистем
ẋ(t) = Ps f (x(t)),
s = 1, . . . , N.
(4.2)
Основной целью данной работы является исследование условий, при которых система (4.1) будет абсолютно устойчива. Для получения необходимых
условий абсолютной устойчивости используется критерий Гурвица, а для исследования достаточных условий применяется прямой метод Ляпунова. Функцию Ляпунова будем строить в виде
V (x) = 2
n
X
Z
λi
xi
fi (u)du,
(4.3)
0
i=1
где λi > 0 для всех i = 1, . . . , n.
4.2. Полученные результаты
В данной работе будем рассматривать системы с матрицами размерности
[3 × 3]. Будем считать, что семейство систем (4.2) представимо в следующем
виде
As 0 0
−d1 a12 a13
(4.4)
0 Bs 0 ẋ(t) = a21 −d2 a23 f (x(t)), s = 1, . . . , N.
0 0 Cs
a31 a32 −d3
11
Представление системы в таком виде обусловлено тем, что во многих задачах математического моделирования [3, 9] особый интерес представляет Dустойчивость подобных систем. Например, при рассмотрении модели "хищник
– жертва" в экологии, часто ставится вопрос, каким должно быть межвидовое
взаимодействие, чтобы положение равновесия сообщества оставалось устойчивым при любых коэффициентах естественного прироста, соответствующих
отдельным видам вне взаимодействия [9].
В соответствии с определением D-устойчивости будем рассматривать случай, когда элементы As , Bs и Cs являются положительными для любого s.
Для существования функции Ляпунова в виде (4.3) также необходимо, чтобы
выполнялось d1 , d2 , d3 > 0.
Для того чтобы нулевое решение этой системы было асимптотически устойчивым для любых допустимых функций необходимо, чтобы оно было асимптотически устойчиво и для соответствующей линейной системы. Поэтому в
данном случае для выяснения необходимых условий абсолютной устойчивости
будем проверять выполнение критерия Гурвица. Для применения критерия
Гурвица, семейство линейных систем, соответствующее семейству (4.4), может
быть приведено к виду
−d1 a12 a13
As As As
a21 −d2 a23
x(t),
ẋ(t) =
B
B
B
s
s
s
a
a32 −d3
31
Cs
Cs
Cs
s = 1, . . . , N.
(4.5)
Будем рассматривать частный случай, когда a13 = a31 = 0.
Замечание. Результаты для случаев a23 = a32 = 0 и a12 = a21 = 0 могут
быть получены аналогичным образом. При том все найденные условия устойчивости для случая a13 = a31 = 0 могут быть переформулированы для двух
других случаев с точностью до обозначений.
4.2.1. Необходимые условия абсолютной устойчивости
Для того чтобы система с переключениями, соответствующая семейству
(4.5), была абсолютно устойчива необходимо, чтобы была абсолютно устой-
12
чива каждая индивидуальная подсистема семейства (4.5).
Рассмотрим систему
−d1 a12
0
A
A
a21 −d2 a23
x(t).
ẋ(t) =
B
B
B
a32 −d3
0
C
C
(4.6)
Исследуем условия, необходимые для асимптотической устойчивости нулевого
решения такой системы.
Необходимые условия: метцелерова матрица.
Для начала рассмотрим случай, когда матрица системы (4.6) является метцелеровой. В этом случае условия выполнения критерия Гурвица эквивалентны условиям Севастьянова — Котелянского [5], то есть −∆1 > 0, ∆2 > 0, −∆3 >
0, где ∆i — главный диагональный минор порядка i матрицы системы (4.6).
Исходя из этих условий, можно сформулировать утверждение.
Утверждение 1: Если матрица системы (4.6) является метцелеровой, то
для того чтобы нулевое решение этой системы было асимптотически устойчиво, необходимо выполнение неравенства
a32 a23 d1 + a12 a21 d3 < d1 d2 d3 .
Необходимые условия: неметцелерова матрица.
Теперь рассмотрим случай, когда матрица системы (4.6) не является метцелеровой. В таком случае условия Севастьянова — Котелянского применить
нельзя, значит, будем применять критерий Гурвица напрямую. Полученные
условия сформулируем в виде утверждения.
Утверждение 2: Для того чтобы нулевое решение системы (4.6) являлось
асимптотически устойчивым, необходимо выполнение условий
13
C 2 (d1 B + d2 A)
C 2 (d1 B + d3 A)
a
a
+
a
a
<
d
d
+
+
32 23
2 3
2 (d B + d C) 12 21
2 (d B + d C)
A
A
3
2
3
2
B 2 (d1 C + d3 A) 2BC (d1 d2 d3 )
+
,
+
2 (d B + d C)
A
A
(d
B
+
d
C)
3
2
3
2
a a d + a a d < d d d .
32 23 1
12 21 3
Доказательство:
(4.7)
1 2 3
Рассмотрим матрицу Гурвица для системы (4.6)
α 1 0
Γ = γ β α ,
0 0 γ
где
d1 d2 d3
+
+ ,
A B
C
d1 d3 d1 d2 − a21 a12 d2 d3 − a23 a32
+
+
,
β=
AC
AB
BC
d1 d2 d3 −d3 a21 a12 −d1 a23 a32
γ=
.
ABC
Для того чтобы нулевое решение системы (4.6) было асимптотически устойчивым, необходимо, чтобы главные миноры матрицы Гурвица были положительны. Учитывая, что A, B и C, d1 , d2 и d3 положительны, очевидно, что ∆1 > 0.
Условие ∆2 > 0 эквивалентно неравенству
α=
C 2 (d1 B + d2 A)
C 2 (d1 B + d3 A)
a32 a23 + 2
a12 a21 < d2 d3 + 2
+
A (d3 B + d2 C)
A (d3 B + d2 C)
B 2 (d1 C + d3 A) 2BC (d1 d2 d3 )
+ 2
+
.
A (d3 B + d2 C)
A (d3 B + d2 C)
А ∆3 > 0 тогда и только тогда, когда
a32 a23 d1 + a12 a21 d3 < d1 d2 d3 .
Что и требовалось доказать.
Таким образом, для того чтобы система с переключениями, соответствую-
14
щая семейству (4.5), была абсолютно устойчивой, необходимо, чтобы условия
вида (4.7) выполнялись для каждой индивидуальной подсистемы этого семейства.
Отметим, что вышеуказанные условия являются лишь необходимыми условиями абсолютной устойчивости системы с переключениями. Для того чтобы
получить достаточные условия устойчивости следует воспользоваться прямым
методом Ляпунова. Чтобы доказать абсолютную устойчивость равномерную
относительно закона переключения, достаточно построить общую функцию
Ляпунова для всех индивидуальных подсистем, соответствующих рассматриваемой системе [1].
4.2.2. Достаточные условия абсолютной устойчивости
Учитывая, что a13 = a31 = 0, семейство систем (4.4) будет иметь вид
As 0 0
−d1 a12
0
0 Bs 0 ẋ(t) = a21 −d2 a23 f (x(t)),
0 0 Cs
0 a32 −d3
s = 1, . . . , N.
(4.8)
Функцию Ляпунова будем строить следующим образом:
x1
Z
V (x) = 2Ã
Z
f1 (u)du + 2B̃
0
x2
Z
f2 (u)du + 2C̃
0
x3
f3 (u)du.
0
Здесь Ã, B̃, C̃ — положительные постоянные.
Тогда производная функции V (x) будет иметь вид
Ã
B̃
Ã
a12 + a21
0
−2 A d1
As
Bs
s
Ã
B̃
B̃
B̃
C̃
V̇ (x) = f T (x(t))
a
+
a
−4
d
a
+
a
2
23
32 f (x(t)).
A 12 B 21
B
B
C
s
s
s
s
s
B̃
C̃
C̃
0
a23 + a32
−2 d3
Bs
Cs
Cs
(4.9)
Если существуют такие положительные Ã, B̃ и C̃, что функция V̇ (x) является
отрицательно определенной для любых s = 1, . . . , N , то система с переключениями, соответствующая семейству систем (4.8), абсолютно устойчива [1].
15
Для исследования отрицательной определенности матрицы каждой системы
семейства (4.9) будем использовать критерий Сильвестра [4].
2Ã
Условие −∆1 =
d1 > 0 всегда выполняется. Рассмотрим неравенства
As
∆2 = 4
−∆3 = 8
à B̃
d1 d2 −
As Bs
Ã
B̃
a12 + a21
As
Bs
à B̃ C̃
C̃
d1 d2 d3 − 2 d3
As Bs Cs
Cs
−2
Ã
d1
As
!2
> 0,
B̃
Ã
a12 + a21
As
Bs
!2
C̃
B̃
a23 + a32
Bs
Cs
!2
−
> 0.
(4.10)
Отметим, что выполнение условия ∆2 > 0 будет следовать из выполнения
условия −∆3 > 0. Таким образом можно утверждать, что если выполняется
неравенство (4.10), то система с переключениями, соответствующая семейству
систем (4.8), будет являться абсолютно устойчивой.
à B̃ C̃
, получим
Разделим левую и правую часть неравенства (4.10) на 2
As Bs Cs
d1
B̃ Cs 2
C̃ Bs 2
à Bs 2
B̃ As 2
a23 +2d1 a23 a32 +d1
a32 +d3
a12 +2d3 a12 a21 +d3
a21 < 4d1 d2 d3 .
C̃ Bs
B̃ Cs
B̃ As
à Bs
B̃
Ã
As
Bs
Сделаем замену переменных. Пусть
= α,
= β,
= xs ,
= ys . В
Bs
Cs
B̃
C̃
новых переменных условие примет вид
d3
α 2
xs
β
ys
a12 + 2d3 a12 a21 + d3 a221 + d1 a223 + 2d1 a23 a32 + d1 a232 < 4d1 d2 d3 ,
xs
α
ys
β
что эквивалентно
r
d3
α
a12 +
xs
r
xs
a21
α
s
2
+ d1
β
a23 +
ys
r
ys
a32
β
!2
< 4d1 d2 d3 .
(4.11)
Достаточным условием абсолютной устойчивости системы с переключения-
16
ми, соответствующей семейству систем (4.8), является выполнение неравенства
(4.11) для всех s = 1, . . . N . То есть если для любого s выполняется,
s
β
a23 +
ys
r
ys
a32
β
!2
< 4d2 d3 −
r
t(xs , α) =
α
a12 +
xs
r
d3
max t(xs , α),
d1 s=1,...,N
xs
a21
α
2
,
то система с переключениями, соответствующая семейству систем (4.8), является абсолютно устойчивой. Однако α и β могут быть произвольными положительными постоянными. Таким образом самая широкая область ограничений
для левой части неравенства будет при
r
t̃ = min max
α>0 s=1,...,N
α
a12 +
xs
r
xs
a21
α
2
.
Точками, подозрительными на максимум по xs , будут
x̃m = min xs
s=1,...,N
и
x̃M = max xs .
s=1,...,N
Обозначим t̃m (α) = t(x̃m , α), t̃M (α) = t(x̃M , α). Отметим, что при всех s =
1, . . . , N будет выполняться одно из неравенств
"
t(xs , α) 6 t̃m (α),
t(xs , α) 6 t̃M (α).
Пример функций t̃m (α) и t̃M (α) изображен на рис.1. Если изобразить графики
t(xs , α), то получим, что все они лежат в заштрихованной области.
17
Рис. 1. Графики функций t̃m (α) и t̃M (α).
Приходим к тому, что α – величина, которая доставляет минимум функции
r
max
s=1,...,N
α
a12 +
xs
r
xs
a21
α
2
– абсцисса точки пересечения графиков функций t̃m (α) и t̃M (α). Таким образом находим, что оптимальная
x̃m x̃M .
a12
Проводя аналогичные рассуждения относительно
s
q(ys , β) =
β
a23 +
ys
r
ys
a32
β
и решая задачу поиска
q̃ = min max q(ys , β),
β
s=1,...,N
18
!2
,
находим, что подозрительными на максимум по ys будут точки
ỹm = min ys
s=1,...,N
и
ỹM = max ys .
s=1,...,N
А абсциссой точки пересечения q̃m (β) = q(ỹm , β) и q̃M (β) = q(ỹM , β) будет
ỹm ỹM .
a23
Таким образом можно сформулировать теорему.
Теорема. Для того чтобы система с переключениями, соответствующая семейству систем (4.8), была абсолютно устойчива, достаточно, чтобы для каждого s выполнялось неравенство
1 + xs
d3 |a12 a21 | √
xs
где xs =
2
1 + ys
+ d1 |a23 a32 | √
ys
Bs
As
, ys =
.
Bs
Cs
19
2
< 4d1 d2 d3 ,
5. Системы с запаздыванием
5.1. Постановка задачи
Пусть задана система
ẋ(t) = Af (x(t)) + Bf (x(t − τ )).
(5.1)
Здесь x(t) = (x1 (t), . . . , xn (t))T — n-мерный вектор состояния системы; A и B —
постоянные матрицы; f (x) = (f1 (x1 ), . . . , fn (xn ))T , скалярные функции fi (xi )
являются допустимыми нелинейностями; τ — постоянное неотрицательное запаздывание. Из сделанных предположений следует, что у изучаемой системы
существует решение x(t) ≡ 0.
Уравнения такого вида широко используется для моделирования систем автоматического управления и нейронных сетей [15, 16].
Исследуем условия, при выполнении которых можно гарантировать абсолютную устойчивость системы (5.1).
Для решения поставленной задачи используется прямой метод Ляпунова.
Пусть
Z xi (t)
Z t
n
n
X
X
V (xt ) = 2
pi
fi (u)du +
qi
fi2 (xi (θ))dθ,
(5.2)
i=1
0
i=1
t−τ
где p1 , . . . , pn и q1 , . . . , qn — положительные коэффициенты.
Нетрудно проверить, что при любых положительных значениях p1 , . . . , pn и
q1 , . . . , qn существуют заданные при s ∈ [0, H) функции Хана χ1 (s) и χ2 (s)
такие, что
χ1 (kx(t)k) 6 V (xt ) 6 χ2 (kxt kτ )
при kxt kτ < H.
Продифференцируем функционал (5.2) в силу системы (5.1). Получим
V̇ (xt ) =
f (x(t))
f (x(t − τ ))
!T
20
S
f (x(t))
f (x(t − τ ))
!
.
Здесь матрица
T
S=
A P + PA + Q PB
BT P
−Q
!
,
а P = diag{p1 , . . . , pn }, Q = diag{q1 , . . . , qn }.
В соответствии с определением 10, для диагональной устойчивости системы
(5.1) необходимо и достаточно, чтобы cуществовали диагональные положительно определенные матрицы P и Q, при которых матрица S отрицательно
определена.
Известно (см. [6]), что условие отрицательной определенности матрицы S
эквивалентно условию отрицательной определенности матрицы
R = AT P + P A + Q + P BQ−1 B T P.
(5.3)
Критерий диагональной устойчивости систем такого вида установлен в работе [7]. Однако в общем случае проверка условий указанного критерия довольно
затруднительна. Поэтому большой интерес представляет задача нахождения
классов матриц A и B, для которых можно получить конструктивно проверяемые условия существования требуемых матриц P и Q.
5.2. Условия существования диагональных
функционалов Ляпунова — Красовского
5.2.1. Системы с якобиевыми матрицами
Будем рассматривать случай, когда A является якобиевой матрицей [8], т. е.
a1
c1
0
A=.
..
0
0
d1 0 · · ·
0
0
a2 d 2 · · ·
0
0
c 2 a3 · · ·
0
0
,
.. .. . . .
..
..
. .
.
.
0 0 · · · an−1 dn−1
0 0 · · · cn−1 an
21
а матрица
b1
0
B=
...
0
b2 · · · bn
0 ··· 0
.
.. . . . ..
.
.
0 ··· 0
(5.4)
Тогда в выражении (5.3)
P
n
b2i p21
qi
i=1
−1 T
P BQ B P = 0
..
.
0
0 0 ··· 0
0 0 · · · 0
.
.. .. . . . ..
. .
.
0 0 ··· 0
Кроме того, пусть для элементов матрицы A имеют место соотношения
ai < 0, dj cj < 0,
i = 1, . . . , n,
j = 1, . . . , n − 1.
(5.5)
Матрицы такого рода широко используются для моделирования трофических
цепей [3].
Для того чтобы получить конструктивно проверяемые условия диагональной устойчивости, будем рассматривать достаточные условия. Пусть
di
pi+1 = − pi ,
ci
i = 1, . . . , n − 1,
(5.6)
где p1 — произвольное положительное число. Тогда матрица R принимает диагональный вид, и для ее отрицательной определенности необходимо и достаточно, чтобы выполнялись неравенства
n
X
b2i p21
+ 2a1 p1 + q1 < 0,
q
i
i=1
q + 2a p < 0, i = 2, . . . , n.
i
i i
Чтобы получить наиболее широкую область для ai , bi , ci , di , следует выбирать qi достаточно близкими к −2ai pi при i = 1, . . . , n. Используя предельный
22
переход, приходим к неравенству
n
X
b2i p1
q1 b21 p1
> 2a1 + +
.
2a
p
p
q
i
i
1
1
i=2
Учитывая соотношение (5.6), имеем
n
i−1
2a
p
q
i
1
1
i=2
j=1
Для того чтобы получить наиболее широкую область допустимых значений
параметров необходимо минимизировать
правую
часть неравенства по p1 , q1 .
2
q1
q1 b1 p1
Очевидно, что если b1 6= 0, то min
+
= 2 |b1 | при
= |b1 |. В случае,
p1
q1
p1
когда b1 = 0, область допустимых значений параметров системы будет тем
шире, чем меньше значение величины q1 /p1 . Исходя из этого, можно сформулировать теорему.
Теорема. Пусть A — матрица Якоби, а B имеет вид (5.4). Пусть, кроме
того, элементы матрицы A удовлетворяют условиям (5.5). Если справедливо
соотношение
n
i−1
> 4(a1 + |b1 |),
a
i
i=2
j=1
то система (5.1) абсолютно устойчива.
Теперь рассмотрим случай, когда матрицы
a1 0 · · · 0
0
c 1 a2 · · · 0
0
A=
,
..
..
... ... . . .
.
.
0 0 · · · cn−1 an
0
0
B=
...
0
0 ···
0 ···
.. . . .
.
0
0
..
.
b
0
,
..
.
(5.7)
0 ··· 0 0
где ai < 0, i = 1, . . . , n; cj 6= 0, j = 1, . . . , n−1; b 6= 0. Таким образом, исследуем
систему с замкнутой петлёй обратной связи. При этом наличие запаздывания
может быть обусловлено задержками, возникающими при формировании и
передаче управляющего сигнала [17].
Теорема. Пусть матрицы A и B определяются по формулам (5.7). Тогда
23
для диагональной устойчивости системы (5.1) необходимо и достаточно, чтобы
выполнялось неравенство
|a1 . . . an | > |b c1 . . . cn−1 |.
(5.8)
Доказательство: В данном случае матрицу (5.3) можно представить в форe + Q,
e где
ме R = R
p21 b2
p 2 c1
0
2p1 a1 + q
n
p2 c1
2p2 a2 p3 c2
e
R=
0
p3 c2 2p3 a3
..
..
..
.
.
.
0
0
0
···
···
···
...
0
0
0
..
.
· · · 2pn an + qn
,
e = diag{q1 , . . . , qn−1 , 0}. Получаем, что система (5.1) диагонально устойа Q
чива тогда и только тогда, когда существуют диагональная положительноопределённая матрица P = diag{p1 , . . . , pn } и число qn > 0, при которых матe отрицательно определена.
рица R
e выполнение условий критерия Сильвестра, имеем
Проверяя для R
p21 b2
∆1 = 2a1 p1 +
< 0,
qn
∆2 = 2a2 p2 − p22 c21 > 0,
(−1)j ∆j = (−1)j (2aj pj ∆j−1 − p2j c2j−1 ∆j−2 ) > 0,
j = 3, . . . , n − 1,
(−1)n ∆n = (−1)n ((2an pn + qn )∆n−1 − p2n c2n−1 ∆n−2 ) > 0.
Из первых (n − 1)-го неравенства данной системы следует, что
p1 = −k1
a1
qn ,
b2
p2 = k2
a2
∆1 ,
c21
pj = kj
24
aj ∆j−1
, j = 3, . . . , n − 1,
c2j−1 ∆j−2
(5.9)
где 0 < ks < 2, s = 1, . . . , n − 1. Тогда
a21
∆1 = k1 (k1 − 2) 2 qn ,
b
∆2
a22
= k2 (2 − k2 ) 2 ∆1 ,
∆1
c1
a2j ∆j−1
∆j
= kj (2 − kj ) 2
,
∆j−1
cj−1 ∆j−2
j = 3, . . . , n − 1.
Подставляя полученные выражения в последнее неравенство системы (5.9),
приходим к условию
qn2
+ 2an pn qn +
p2n b2
n−1
Y
c2j
2 < 0.
k
(2
−
k
)a
j
j
j
j=1
(5.10)
Левую часть неравенства (5.10) можно рассматривать как квадратный трехчлен относительно qn или относительно pn . Нетрудно видеть, что если такой
трехчлен имеет вещественные корни, то они одного знака. Таким образом для
того чтобы существовали положительные числа pn и qn , удовлетворяющие соотношению (5.10), необходимо и достаточно, чтобы дискриминант квадратного
уравнения, соответствующего неравенству (5.10), был положителен, что эквивалентно выполнению неравенства
n
Y
i=1
a2i
n−1
Y
c2j
.
>b
k
(2
−
k
)
j
j
j=1
2
(5.11)
Неравенство (5.11) задает наиболее широкую область допустимых значений
параметров системы в случае, когда
kj = 1,
j = 1, . . . , n − 1.
Таким образом, приходим к условию (5.8).
Теорема доказана.
Система с переключениями
Теперь рассмотрим случай, когда помимо запаздывания в системе имеют
место допустимые переключения. В каждый момент времени система описы-
25
вается одной из подсистем семейства
ẋ(t) = As f (x(t)) + Bs f (x(t − τ )),
s = 1, . . . , N.
(5.12)
Пусть
(s)
a1
(s)
c1
As = .
..
0
0 ···
(s)
a2 · · ·
.. . . .
.
0
0
0
..
.
(s)
0
0
..
.
(s)
· · · cn−1 an
0 · · · b(s)
Bs = ... . . . ... ,
0 ··· 0
,
(5.13)
(s)
где ai < 0, i = 1, . . . , n.
Для того чтобы такая система была диагонально устойчива, достаточно,
чтобы выполнялись неравенства
c
As + Bs ξ < 0,
cs + Br
A
T
θ < 0,
s = 1, . . . N,
(5.14)
s, r = 1, . . . N,
(5.15)
cs – метцелегде ξ = (ξ1 , . . . , ξn ), θ = (θ1 , . . . , θn )T – положительные векторы; A
рова матрица, соответствующая матрице A; Bs – матрица, элементы которой
соответствуют элементам матрицы B, взятым по модулю [18].
Распишем системы неравенств (5.14) и (5.15) для семейства систем вида
(5.12) с матрицами (5.13). Имеем
(s)
a1 ξ1 + |b(n) |ξn < 0,
(s)
(s)
(s)
(s)
(s)
(s)
|c1 |θ2 + a1 θ1 < 0,
(s)
(s)
|c2 |θ3 + a2 θ2 < 0,
(s)
(s)
|c3 |θ4 + a3 θ3 < 0,
|c1 |ξ1 + a2 ξ2 < 0,
|c2 |ξ2 + a3 ξ3 < 0,
···
···
|cn−1 |(s) ξn−1 + an(s) ξn < 0,
26
|b(r) |θ1 + a(s)
n θn < 0,
откуда следует, что
(s) (s)
(s)
|a1 |ξ1
(s)
|b |ξn <
<
|a1 a2 |
(s)
|c1 |
(s)
i=1 |ai |
Qn−1 (s) ξn ,
j=1 |cj |
Qn
ξ2 < · · · <
(s) (s)
(r)
|b |θ1 <
|a(s)
n |θn
<
|an an−1 |
(s)
|cn−1 |
(s)
i=1 |ai |
Qn−1 (s) θ1 .
j=1 |cj |
Qn
θn−1 < · · · <
Исходя из этого, можно сформулировать теорему.
Теорема. Для того чтобы система с допустимыми переключениями и постоянным запаздыванием, соответствующая семейству систем вида (5.12), была
диагонально устойчива, достаточно, чтобы выполнялось неравенство
(r)
max |b | <
r=1,··· ,N
(s)
i=1 |ai |
min Qn−1 (s) .
s=1,··· ,N
j=1 |cj |
Qn
Отметим, что полученные условия для систем с запаздыванием и с переключением являются достаточными, но не необходимыми. Приведем пример.
Пример: Пусть n = 3, N = 2,
−1 0
0
A1 = 2 −2 0 ,
0 0.02 −1
−10 0
0
A2 = 10 −20 0 ,
0 0.25 −10
0 0 1
B1 = 0 0 0 ,
0 0 0
0 0 100
B2 = 0 0 0 .
0 0 0
В данном случае
(s) (s) (s)
(r)
max |b | = 100;
r=1,2
min
s=1,2
27
|a1 a2 a3 |
(s) (s)
|c1 c2 |
= 50,
то есть условия теоремы не выполняются. Но если взять матрицы
0.01 0
0
P = 0 0.01 0 ,
0
0 50
0.00001
0
0
Q= 0
0.00001 0 ,
0
0
10
то матрица R будет отрицательно определена. Значит, система с переключениями будет абсолютно устойчивой.
5.2.2. Условия диагональной устойчивости матриц третьего
порядка со специальной структурой
Будем рассматривать систему вида (5.1) с матрицами
a1 0 0
0 0 b1
A = c1 a2 0 , B = 0 0 b2 ,
0 c 2 a3
0 0 0
где a1 , a2 , a3 < 0.
e + Q,
e где
В данном случае матрицу R можно представить в форме R = R
b1 b2 p1 p2
p21 b21
p
c
+
2p
a
+
2
1
1
1
q3
q3
e=
p 2 b2
bbpp
R
p2 c1 + 1 2 1 2 2p2 a2 + 2 2
q3
q3
0
p3 c 2
0
p3 c 2
q1 0 0
e=
, Q
0 q 2 0 .
0 0 0
2p3 a3 + q3
Получаем, что система (5.1) диагонально устойчива тогда и только тогда,
когда существуют диагональная положительно-определённая матрица P =
e отрицательно опреdiag{p1 , p2 , p3 } и число q3 > 0, при которых матрица R
делена.
e ОтРассмотрим выполнение условий критерия Сильвестра для матрицы R.
метим, что от условия ∆3 < 0 можно перейти к условию
∆3
=
p23
a3 q 3
2 + 2
p3 p3
28
∆2 − c22 ∆1 < 0.
Так как ∆1 и ∆2 от p3 не зависят, для получения наиболее широкой области
допустимых значений параметров следует выбирать p3 таким образом, чтобы
∆3 /p23 было минимальным, откуда следует, что
p3 = −
q3
.
a3
Тогда для того чтобы существовали положительные числа p1 , p2 , p3 и q3 , при
e будет отрицательно определенной, необходимо и достаточкоторых матрица R
но, чтобы существовали положительные числа p1 , p2 и q3 при которых будет
отрицательно определена матрица
p21 b21
b1 b2 p1 p2
p2 c1 +
0
2p1 a1 + q
q3
3
2 2
b
b
p
p
p
b
1
2
1
2
2
2
p2 c1 +
2p
a
+
c
2 2
2 .
q
q
3
3
a23
0
c2
−
q3
Теперь рассмотрим условие
p21 b21
∆1 = 2p1 a1 +
< 0.
q3
(5.16)
Неравенство (5.16) выполняется, если
p1 = −k1
a1
q3 ,
b21
0 < k1 < 2.
Тогда
∆1 = −k1 (2 − k1 )
a21
q3 .
b21
Рассмотрим неравенство
∆2 =
p22 b22
2p2 a2 +
q3
b1 b2 p1 p2
∆1 − p2 c1 +
q3
29
2
> 0.
Подставляя выражения, полученные для ∆1 и p1 , приходим к условиям
−2a2 k1 (2 −
a21
k1 ) 2 q3
b1
− p2 k1 (2 −
a2 b2
a1 b2 c1
2k1 12 2 + c21 − 2k1
b1
b1
a21 b22
k1 ) 2
b1
2 !
a1 b 2
+ c1 − k1
> 0,
b1
p2 < −2k1 (2 − k1 )
a21 a2
q3 ,
b21
2k1 a21 b22 + c21 b21 − 2k1 a1 b1 b2 c1 p2 < −2k1 (2 − k1 )a21 a2 q3 ,
p2 = −
k1 k2 (2 − k1 )a21 a2
q3 ,
2k1 a21 b22 + c21 b21 − 2k1 a1 b1 b2 c1
0 < k2 < 2.
Подставляя выражение для p2 в ∆2 , получим
∆2 =
k12 k2 (2 − k1 )2 (2 − k2 )a41 a22
q32 .
2
2
2
2
2
b1 (2k1 a1 b2 + c1 b1 − 2k1 a1 b1 b2 c1 )
Теперь рассмотрим условие
∆3 = −
a23
∆2 − c22 ∆1 < 0,
q3
которое эквивалентно неравенству
a21 c22
k12 k2 (2 − k1 )2 (2 − k2 )a41 a22 a23
+ k1 (2 − k1 ) 2 < 0,
− 2
b1 (2k1 a21 b22 + c21 b21 − 2k1 a1 b1 b2 c1 )
b1
откуда получаем
c22
k1 k2 (2 − k1 )(2 − k2 )
<
.
a21 a22 a23
(2k1 a21 b22 + c21 b21 − 2k1 a1 b1 b2 c1 )
Чтобы область допустимых значений параметров системы была наиболее
широкой, надо найти, при каких k1 , k2 правая часть неравенства достигает
своего максимума. Нетрудно видеть, что оптимальное значение k2 = 1. Тогда
a21 a22 a23
2k1 a21 b22 + c21 b21 − 2k1 a1 b1 b2 c1
<
,
k1 (2 − k1 )
c22
30
a21 b22
2
c1 b 1
k1 −
a21 a22 a23
a1 b2
<
− a21 b22 ,
2
k1 (2 − k1 )
c2
2
c1 b1
k1 −
a22 a23
a1 b 2
< 2 2 − 1.
k1 (2 − k1 )
c3 b2
(5.17)
Отметим, что наиболее широкая область допустимых значений параметров
c1 b 1
системы будет достигаться при k1 =
, но при этом должно выполняться
a1 b 2
c1 b1
< 2, что не всегда так. Выясним, как выбирать k1 , если не выпол0<
a1 b 2
c1 b1
няется условие 0 <
< 2. Для этого исследуем левую часть неравенства
a1 b2
(5.17). Имеем
2 0
c1 b 1
c1 b 1 2
c21 b21
c21 b21
2
k1 −
−
2
+
2
2k
k
k
−
2
1
1
a1 b 2
a1 b2 1
a21 b22
a21 b22
= 0,
=
k1 (2 − k1 )
k12 (2 − k1 )2
k1
c1 b1
k1 −
a1 b 2
c1 b 1
c1 b 1
1−
k1 +
= 0.
a1 b2
a1 b 2
Нетрудно проверить, что если
c1 b 1
c1 b 1
c1 b1
< 0 или
> 2, то 0 <
<2и
a1 b 2
a1 b 2
c1 b1 − a1 b2
c1 b 1
доставляет минимум левой части неравенства (5.17). Исходя
c1 b1 − a1 b2
из этого, сформулируем теорему.
Теорема. Система
k1 =
a1 0 0
0 0 b1
ẋ(t) = c1 a2 0 f (x(t)) + 0 0 b2 f (x(t − τ ))
0 0 0
0 c 2 a3
диагонально устойчива тогда и только тогда, когда выполняется одно из условий
1)
0 < c1 b1 < 2a1 b2 ,
31
a32 a23 > c22 b22 ;
2)
0 < c1 b1 < 2c1 b1 − 2a1 b2 ,
a22 a23
c1 b1 (2a1 b2 − c1 b1 )2
< 2 2 − 1.
a21 b22 (c1 b1 − 2a1 b2 )
c2 b2
Теперь рассмотрим систему с матрицами
a1 0 0
0 0 b1
A = 0 a2 0 , B = 0 0 b2 ,
c 1 c 2 a3
0 0 0
где a1 , a2 , a3 < 0.
e + Q,
e где
В данном случае матрицу R можно представить в форме R = R
b1 b2 p1 p2
p21 b21
2p
a
+
1
1
q3
q3
e = b 1 b 2 p1 p2
p22 b22
R
2p
a
+
2 2
q3
q3
p 3 c1
p3 c2
p3 c1
p3 c2
q1 0 0
e=
, Q
0 q 2 0 .
0 0 0
2p3 a3 + q3
Получаем, что в данном случае система (5.1) диагонально устойчива тогда
и только тогда, когда существуют диагональная положительно-определённая
e отрицаматрица P = diag{p1 , p2 , p3 } и число q3 > 0, при которых матрица R
тельно определена.
e АнаРассмотрим выполнение условий критерия Сильвестра для матрицы R.
логично предыдущему случаю можем перейти от условия ∆3 < 0 к условию
∆3
=
p23
2 2
b
b
p
p
a3 q 3
p
b
1
2
1
2
2 + 2 ∆2 − c22 ∆1 + 2c1 c2
− c21 2p2 a2 +
< 0,
p3 p3
q3
q3
a3 q 3
+ . Для того чтобы область допустимых значеp3 p23
ний параметров системы была наиболее широкой, следует выбирать p3 таким
образом, чтобы минимизировать ∆3 /p23 . Откуда следует, что
где от p3 зависит только 2
p3 = −
32
q3
.
a3
Тогда, в данном случае, система (5.1) диагонально устойчива тогда и только
тогда, когда существуют положительные числа p1 , p2 и q3 , при которых отрицательно определена матрица
p21 b21
2p1 a1 + q
3
b 1 b 2 p 1 p2
q3
c1
b1 b2 p1 p2
q3
2p2 a2 +
p22 b22
q3
c2
c1
c2
.
2
a3
−
q3
Рассмотрим условие
p21 b21
∆1 = 2p1 a1 +
< 0.
q3
Неравенство (5.18) выполняется если
p1 = −k1
a1
q3 ,
b21
0 < k1 < 2.
Тогда
∆1 = −k1 (2 − k1 )
a21
q3 .
b21
Теперь рассмотрим неравенство
p22 b22
b21 b22 p21 p22
∆2 = (2p2 a2 +
)∆1 −
> 0.
q3
q32
Подставляя в ∆2 выражения, найденные для ∆1 и p1 , получаем
−2a2 k1 (2 − k1 )
a21 b22
a21
q
−
2k
3
1 2 p2 > 0,
b21
b1
p2 < −(2 − k1 )
p2 = −k2 (2 − k1 )
a2
q3 ,
b22
33
a2
q3 ,
b22
0 < k2 < 1.
(5.18)
Подставляя выражение для p2 в ∆2 , получим
a21 a2 2
∆2 = k1 k2 (2 − k1 ) (2 + 2k2 − k1 k2 − k1 ) 2 q3 .
b1
2
Теперь рассмотрим условие
∆3 =
−k12 (2
2 2 2
a21 c22
a21 a22 a23
2 a1 a2 a3
2 2
q3 + k1 (2 − k1 ) 2 q3 + k1 k2 (2 − k1 )
q3 +
− k1 )
b21 b22
b1
b21 b22
2
a1 a2 c1 c2
a2 c 2
q3 + k1 (2 − k1 ) 22 1 q3 < 0,
b1 b2
b2
что эквивалентно соотношению
+2k1 k2 (2 − k1 )
2 2
a21 c22
a1 a2 c1 c2 a22 c21
2 a1 a2
+ 2k2
+ 2 2 < k1 (2 − k1 ) 1 − k2 2 2 .
b21 a23
b1 b2 a23
b 2 a3
b1 b2
Правая часть этого неравенства положительна, а, значит, наиболее широкая
область допустимых значений параметров будет получена при k1 = 1. Откуда
следует
a2 a2
a1 a2 c1 c2 a21 c22 a22 c21
a21 a22
+
k22 12 22 + 2k2
+
<
.
(5.19)
b1 b2
b1 b2 a23
b21 a23 b22 a23
b21 b22
Для того чтобы область допустимых значений параметров была наиболее широкой следует выбирать k2 таким образом, чтобы минимизировать левую часть
неравенства (5.19). Нетрудно видеть, что минимум левой части неравенства
b1 b2 c1 c2
. Но, кроме того, должно выполняться
(5.19) достигается при k2 = −
a1 a2 a23
b 1 b 2 c1 c2
0 < k2 < 1. Поэтому если −
< 0, то для того чтобы минимизировать
a1 a2 a23
левую часть неравенства (5.19), следует брать k2 достаточно близким к 0. Если
b1 b2 c1 c2
> 1, то k2 следует брать достаточно близким к 1. Таким образом
же −
a1 a2 a23
сформулируем теорему.
Теорема. Система
a1 0 0
0 0 b1
ẋ(t) = 0 a2 0 f (x(t)) + 0 0 b2 f (x(t − τ ))
c 1 c 2 a3
0 0 0
34
диагонально устойчива тогда и только тогда, когда выполняется одно из условий
1)
−b1 b2 c1 c2 < 0,
a21 b22 c21 + a22 b21 c22 < a21 a22 a23 ;
2)
0 < −b1 b2 c1 c2 < a1 a2 a23 ,
a21 b22 c21
+
a22 b21 c22
<
a21 a22 a23
b21 b22 c21 c22
;
+
a23
3)
−b1 b2 c1 c2 > a1 a2 a23 ,
(a1 b2 c1 + a2 b1 c2 )2 < a21 a22 a23 .
35
6. Заключение
В данной работе были получены необходимые условия абсолютной устойчивости и достаточные условия диагональной устойчивости для одного класса
нелинейных систем с переключениями.
Также была рассмотрена нелинейная система с запаздыванием с якобиевой
матрицей коэффициентов при нелинейностях, для которой были установлены достаточные условия диагональной устойчивости. Кроме того, для частного случая такой системы был получен критерий диагональной устойчивости.
Помимо этого, была рассмотрена система, в которой, наряду с постоянным запаздыванием, имеют место допустимые переключения. Для этого случая были
получены достаточные условия диагональной устойчивости.
Были также рассмотрены два специальных случая систем третьего порядка
с запаздываниями и найдены критерии их диагональной устойчивости.
Все полученные условия являются конструктивно проверяемыми и их выполнение гарантирует абсолютную устойчивость рассматриваемым системам.
Часть результатов выпускной квалификационной работы была представлена на международной конференции "Процессы управления и устойчивость.
Control Processes and Stability. CPS’15". По результатам конференции была
опубликована статья [19]. Также часть результатов была представлена на международной конференции "Процессы управления и устойчивость. Control Processes and Stability. CPS’16". Материалы доклада приняты к печати.
36
Список литературы
[1] Александров А. Ю., Платонов А. В. Метод сравнения и устойчивости
движений нелинейных систем. СПб.: Издательство Санкт-Петербургского
униветситета, 2012. 263 с.
[2] Gu K., Kharitonov V. L., Chen J. Stability of Time-Delay Systems. Boston,
MA: Birkhauser, 2003. 343 p.
[3] Свирежев Ю. М., Логофет Д. О.
сообществ. М.: Наука, 1978. 352 c.
Устойчивость
биологических
[4] Александров А. Ю., Александрова Е. Б., Екимов А. В., Смирнов Н. В.
Сборник задач и упражнений по теории устойчивости. СПб.: ООО ”СОЛО”, 2003. 162 с.
[5] Александров А. Ю., Платонов А. В. Об абсолютной устойчивости одного
класса нелинейных систем с переключениями //Автоматика и телемеханика. 2008. № 7. C. 3–18 .
[6] Mason O. Diagonal Riccati stability and positive time-delay systems //
Systems and Control Letters. 2012. Vol. 61, No 1. P. 6–10.
[7] Aleksandrov A. Yu., Mason O. Diagonal Riccati stability and applications //
Linear Algebra and its Applications. 2016. Vol. 492. P. 38–51.
[8] Мишина А. П., Проскуряков И. В. Высшая алгебра. М.: Наука, 1965. 300 c.
[9] Александров А. Ю., Платонов А. В., Старков В. Н., Степенко Н. А. Математическое моделирование и исследование устойчивости биологических
сообществ. СПб.: ООО ”СОЛО”, 2006. 186 c.
[10] Liberzon D., Morse A. S. Basic Problems in Stability and Design of Switched
Systems. // IEEE Control Systems Magazine. 1999. P. 1–19.
[11] Воронов А. А. Основы теории автоматического управления: Автомат. регулирование непрерыв. линейн. систем. М.: Энергия, 1980. 309 с.
37
[12] Кановей Г. В., Логофет Д. О. D-устойчивость матриц 4 × 4 // Журнал
вычислительной математики и математической физики. 1998. том 38. № 9.
C. 1429–1435.
[13] Руш Н., Абетс П., Лалуа М. Прямой метод Ляпунова в теории устойчивости. М.: Мир, 1980. 300 с.
[14] Барбашин Е. А. Функции Ляпунова. М.: Наука, 1970. 240 с.
[15] Kaszkurewicz E., Bhaya A. Matrix Diagonal Stability in Systems and
Computation. Boston, Basel, Berlin: Birkhauser, 1999. 267 p.
[16] Liao X., Yu P. Absolute Stability of Nonlinear Control Systems. New York,
Heidelberg: Springer, 2008. 384 p.
[17] Hale J. K. Theory of Functional Differential Equations. New York, Heidelberg,
Berlin: Springer-Verlag, 1977. 366 p.
[18] Aleksandrov A. Yu., Mason O. Diagonal Lyapunov – Krasovskii functionals for
discrete-time positive systems with delay // Systems and Control Letters.
2014. Vol. 63, No 1. P. 63–67.
[19] Александров А. Ю., Воробьёва А. А. Построение функционалов Ляпунова – Красовского для одного класса нелинейных систем с запаздыванием
// Процессы управления и устойчивость. 2015. T. 2. № 1 С. 17–22.
38
Отзывы:
Авторизуйтесь, чтобы оставить отзыв