АНАЛИЗ И СИНТЕЗ СИСТЕМ УПРАВЛЕНИЯ
т-' I
УДК 62-501.52
ОБ УСЛОВИЯХ ГРУБОСТИ НЕУСТОЙЧИВЫХ НЕАВТВНВМНЫХ ЛИНЕЙНЫХ СИСТЕМ В СМЫСЛЕ СОХРАНЕНИЯ ХАРАКТЕРА УСТОЙЧИВОСТИ
В.П. Жуков
Институт проблем управления им. В.А. Трапезникова РАН, г. Москва
Рассмотрены условия, при которых неустойчивость состояния равновесия неавтономных линейных динамических систем произвольного порядка переходит при определенного класса нелинейных возмущениях их правых частей в неустойчивость состояния равновесия соответствующих нелинейных возмущенных динамических систем (грубость неустойчивых неавтономных линейных систем в смысле сохранения характера устойчивости). Приведены достаточные условия такой грубости относительно некоторых классов нелинейных возмущений.
ВВЕДЕНИЕ. ПОСТАНОВКА ЗАДАЧИ
Будем рассматривать грубость неавтономной линейной динамической системы произвольного порядка п
(1)
/ = і
или в векторном виде
х = А(?)х, х е Rn,
где А(?) — квадратная матрица п х п с непрерывно зависящими от аргумента ? вещественными эле-
ности переменных хр ..., хп, ? вместе со своими первыми частными производными Зт,.(х, ?)/Зх.,
* У
/, j = 1, ..., п (фг(х, ?) е С1, I = 1, ..., п; С1 — наиболее широкий класс непрерывно дифференцируемых по переменным х1, ..., хп функций, включающий в себя как класс аналитических функций, так и класс всех неаналитических функций, непрерывно дифференцируемых по переменным х1, ..., хп не менее одного раза). При указанных условиях на функции фг(х, ?), / = 1, ..., п, для любой точки (х0, ?0) е Gp существует [1] единственное решение системы (2), проходящее через эту точку и максимально продолженное по ? в пределах мно-
ментами <2.(0, /, j = 1, ..., п, которые в общем случае жества Gp. Пусть эти функции фг(х, ?) при любом
отличны от констант. Точка х = 0, очевидно, точка равновесия системы (1).
Наряду с исходной (невозмущенной) системой (1) будем рассматривать возмущенную систему
х = Л(ґ)х + ф(х, ґ), х є Rn, (2)
где возмущение описывается векторной функцией ф(х, ґ), ф(0, ґ) = 0, компоненты ф;.(х, ґ), ф;.(0, ґ) = 0, і = 1, ..., п, которой определены на множестве Gp = G х R+ с Rn + 1 (R+ = [0, то) — положительная
полуось аргумента ґ, G с Rn — область, содержа-
ґ є R удовлетворяют также условиям
дфі (х, о
З*;
= 0, і, 7 = 1,
(3)
х = 0
Условия ф;.(х, ?) е С1, фг(0, ?) = 0, / = 1, ..., п, в совокупности с условиями (3) назовем основными условиями, наложенными на функции фг(х, ?). Основные условия обеспечивают при любом ? е R+ строгую нелинейность каждой функции фг(х, ?), не равной тождественно нулю, по переменным х1, ..., хп,
щая точку х = 0) и непрерывны в Gp по совокуп- т. е. обеспечивают ее нелинейность при отсутствии
п
п
в ней линейной составляющей по этим переменным. Действительно, из условий фг(х, ?) е СХ , / = 1,
..., п, следует, что при любом t е R+ приращение Афг(х, 0 = фг(х, t) — фг(0, t) = фг(х, t) каждой функции ф.(х, t) при переходе из точки (0, 0 е Gp в точку
I р
(х, 0 е Gp представимо в виде [2, 3]
Лфг(х, 0 = фг(х, 0 = І х;дф/
: =і :
+ 0.(х, ґ),
х = 0
ііш °4х 0 -
= 0, і = 1,
п.
При условиях (3) первый член в правой части этого соотношения (главная линейная часть функции ф.(х, 0 = Аф.(х, 0 в точке х = 0) равен нулю, из чего следует, что
ч 0.(х, ґ)
ф.(х, ґ) = 0.(х, ґ), ііт —
= 0, і = 1, ..., п.
Это означает, что функции ф.(х, *) не содержат линейной составляющей и поэтому при любом t е R+ действительно являются строго нелинейными по переменным х1, ..., хп.
Дополнительная характеристика рассматриваемых классов возмущающих функций ф(х, t) будет приведена в § 1 (пп. 1.1 и 1.2, где в виде двух теорем даются соответствующие двум рассматриваемым классам возмущений ф(х, ^ условия грубости систем вида (1)).
Если некоторая невозмущенная динамическая система имеет точку равновесия х = 0 одного из возможных характеров устойчивости (асимптотическая устойчивость, неасимптотическая устойчивость, неустойчивость) и нас интересуют лишь условия, при которых точка равновесия х = 0 возмущенной системы будет иметь такой же характер устойчивости, то возникает понятие грубости невозмущенной системы в смысле сохранения характера устойчивости. При исследовании конкретных динамических систем чаще возникает вопрос именно о сохранении возмущенной системой характера устойчивости, а не о топологической эквивалентности фазовых портретов возмущенной и невозмущенной систем (т. е. о структурной устойчивости невозмущенной системы), требующей существования соответствующего гомеоморфизма, обеспечивающего указанную топологическую эквивалентность.
При исследовании условий грубости невозмущенных систем (1) будем исходить из следующего определения.
Определение. Невозмущенную неавтономную линейную динамическую систему (1) будем называть грубой в смысле сохранения характера устойчивости по отношению к нелинейным возмущениям ф(х, 0 рассматриваемого класса, если характер устойчивости точки равновесия х = 0 возмущенной системы остается таким же, как и у точки равновесия х = 0 невозмущенной системы. ♦
В определении имеется в виду лишь факт сохранения характера устойчивости точки равновесия х = 0 (например, если точка х = 0 системы (1) асимптотически устойчива или неустойчива, то и точка х = 0 системы (2) соответственно асимптотически устойчива или неустойчива). При этом не предполагается существование гомеоморфизма, осуществляющего топологическую эквивалентность фазовых портретов невозмущенной и возмущенной систем, т. е. данное в определении понятие грубости не является определением грубости в смысле структурной устойчивости [4], исследование условий существования которой и их применение связаны с большими трудностями (примером могут служить условия, полученные А.А. Андроновым и Л.С. Понтрягиным для систем второго порядка [5, 6]).
Цель статьи состоит в том, чтобы указать свойства матрицы А(^, обеспечивающие условия грубости (в смысле данного определения) неустойчивых невозмущенных систем (1) по отношению к рассматриваемым классам нелинейных возмущений ф(х, ^.
1. ДОСТАТОЧНЫЕ УСЛОВИЯ ГРУБОСТИ НЕУСТОЙЧИВЫХ ДИНАМИЧЕСКИХ СИСТЕМ (1)
Рассмотрим условия грубости неустойчивых динамических систем (1) по отношению к двум классам возмущающих нелинейных функций ф(х, *).
Пусть оба эти класса (назовем их классами К1 и К2) удовлетворяют описанным во Введении основным условиям (фг(х, 0 е СХ, фг(0, 0 = 0, / = 1, ..., п; функции ф.(х, 0 удовлетворяют соотношению (3)); кроме того, каждый их этих классов удовлетворяет своему дополнительному условию, рассматриваемому далее в одном из пп. 1.1 и 1.2. В этих же подразделах в теоремах 1 и 2 в терминах свойств матрицы А(^ приводятся два достаточных условия грубости неустойчивых систем (1), соответствующие двум рассматриваемым классам возмущений ф(х, ^.
1.1. Сначала укажем дополнительное условие, предъявляемое к нелинейным возмущениям ф(х, 0 класса К1. Оно заключается в том, что для функ-
ции ф(х, t) из класса К1 можно указать такую непрерывную определенно отрицательную функцию ^(х) (^(0) = 0, ^(х) < 0 при х ф 0), что в некоторой
окрестности е точки х = 0 при любом t е R+ выполняется соотношение
d^(x, t) = і дф д x’ t) i w(x).
i = 1
дxl■
(4)
Заметим, что из условия (3) следует div9(x, t)|x=0 =
= 0 при любом t е R+. Это определяет характер функции div9(x, t) в окрестности точки x = 0: при
любом фиксированном значении t е R+ div9(x, t) как непрерывная функция аргумента x проходит в точке x = 0 через нулевое значение. Заметим так-
П
же, что divA(t)x = ^ fl;;(t) = q(t) как след матрицы i = 1
A(t) является непрерывной функцией лишь аргумента t.
Имеет место следующая теорема, дающая достаточные условия грубости неустойчивых неавтономных линейных систем (1) по отношению к нелинейным возмущениям ф(х, t) из класса K1.
Теорема 1. Если матрица A(t), соответствующая системе (1), такова, что divA(t)x = q(t) как непрерывная функция аргумента t е R+ ограничена снизу положительной константой (divA(t)x = q(t) l l C = const > 0), то невозмущенная неавтономная линейная система (1), имея неустойчивую точку равновесия x = 0, является по отношению к возмущениям ф(х, t) из класса K1 грубой в смысле сохранения характера устойчивости. ♦
Доказательство теоремы приведено в Приложении.
Рассмотрим пример применения теоремы 1. Пусть необходимо исследовать грубость в смысле сохранения характера устойчивости точки равновесия x = 0 неавтономной линейной системы
,~h
X1 = —4x1(1 + e-) — Зx2 + 7x3, X2 = (5x1 — 4x2 + x3)(1 + e *), X3 = 8x1 + 7x2 + ^x3(1 + e *).
(5)
Так как для исследуемой системы уравнений
имеем #(^ = divA(t)x = —4(1 + е-) — 4(1 + е-) +
+ 10(1 + е () = 2(1 + е () > 2 > 0, т. е. #(^ > 2 при
любом t е R+, то согласно теореме 1 рассматриваемая система, имея неустойчивую точку равновесия х = 0, является по отношению к возмущениям ф(х, 0 класса К1 грубой в смысле сохранения ха-
рактера устойчивости. Это, очевидно, позволяет утверждать, что точка равновесия х = 0 любой нелинейной системы, полученной возмущением системы (5) нелинейной функцией ф(х, 0, является неустойчивой, если эта функция принадлежит классу К1. В связи с этим рассмотрим нелинейную систему
х 1 = —4х1(1 + е *) — 3х2 + 7х3 — 3х1 х2, х2 = (5х1 — 4х2 + х3)(1 + е *) — 3 х\ х2, (6)
хз = 8х1 + 7х2 + 10х3(1 + е *) — х3 (1 + е %
полученную возмущением системы (5) нелинейной функцией ф(х, 0 с компонентами ф1(х, t) =
= —3х1 х2, ф2(х, t) = —3х1 х2, ф1(х, t) = — х3 (1 + е *). Эта функция удовлетворяет основным условиям. Проверим, принадлежит ли она классу К1. Так как
^ф(х, о = £ =
г = 1 г
= -3х2 - 3х1 - 3х3 (1 + е-'),
то функция ф(х, t) удовлетворяет дополнительному условию: при любом t е R+ divф(x, ^ 1 ^(х) = = —3 х2 — 3 х2 — 6 х3, где ^(х) — определенно отрицательная функция. Поэтому функция ф(х, 0 принадлежит к классу К1. Так как система (5) согласно теореме 1 неустойчива и груба в смысле сохранения характера устойчивости по отношению к возмущениям класса К1, то, очевидно, неустойчивой будет и точка равновесия х = 0 возмущенной системы (6).
1.2. Укажем теперь дополнительное условие, предъявляемое к возмущениям ф(х, 0 класса К2. Оно заключается в том, что для функции ф(х, ^ из класса К2 в некоторой окрестности е точки х = 0
при любом t е R выполняется соотношение div9(x, t) l 0,
(7)
включая тождественное равенство нулю. При этом divф(x, ^|Х = 0 = 0 при любом t е R+, что следует из условия (3).
Следующая теорема дает достаточное условие грубости неустойчивых неавтономных систем (1) по отношению к нелинейным возмущениям ф(х, 0 из класса К2.
Теорема 2. Если матрица А(^, соответствующая системе (1), такова, что divA(t)x = #(^ как непре-
рывная функция аргумента t е R+, не удовлетворяя1 условию теоремы 1 (divA(t)x = q(t) l C = const > 0), удовлетворяет условию
lim fq(T)dT
t — да •*
(8)
то невозмущенная неавтономная линейная система (1), имея неустойчивую точку равновесия x = 0, является по отношению к возмущениям ф(х, t) из класса К2 грубой в смысле сохранения характера устойчивости. ♦
Доказательство приведено в Приложении. Заметим, что при удовлетворении условию теоремы 2 функция q(t) может принимать положительные значения не обязательно во всех точках множества R+, а лишь на некотором его подмножестве M с R+; на множестве R+/M эта функция может принимать либо только нулевые значения, либо отрицательные и нулевые. При этом характер изменения непрерывной функции q(t) при t ^ может быть произвольным: она может как иметь обычный предел lim q(t), так и не иметь его, имея всегда [7]
t —— да
нижний предел
lim q(t) = lim inf q(£)
t —да t — да t <^<да
и верхний предел
lim q(t) = lim sup q(£).
t —да t —да t <^<да
Рассмотрим пример применения теоремы 2. Исследуем грубость в смысле сохранения характера устойчивости точки равновесия x = 0 неавтономной нелинейной системы
x1 =
x + 4x + 5x ,
11 +1 2 3
к~ = (6x, — x, + 8x,) -1- ,
2 v 1 2 3/1 +1
х, = 7x, + 3x + X -------------------.
3 1 2 31 +1
(9)
Для системы (9) имеем #(?) = divA(^)x = —2/(1 + + 0 - 1/(1 + 0 + 7/(1 + 0 = 4/(1 + ?). Так как функция #(?)не удовлетворяет условию теоремы 1
1 В теореме 2 условие q(t) l C = const > 0, из которого, очевидно, соотношение (8) следует, исключено потому, что это условие уже использовано в теореме 1 для получения условия грубости по отношению к возмущениям ф(х, t) класса K1, который шире класса K2, используемого в теореме 2. Сказанное означает, что множества грубых систем (1), соответствующие теоремам 1 и 2, не пересекаются.
(при любом t е R+ q(t) l C = const > 0) и удовлетворяет условию (8), ибо
tt
lim [q(T)dT = lim [t — да ^ t — да ^ 1
оо
+т
■ т = то,
то согласно теореме 2 система (9), имея неустойчивую точку равновесия х = 0, является по отношению к возмущениям ф(х, ?) класса К2 грубой в смысле сохранения характера устойчивости.
ЗАКЛЮЧЕНИЕ
Указаны два непересекающихся класса неустойчивых неавтономных линейных систем (1), которые являются грубыми в смысле сохранения характера устойчивости относительно соответствующих классов К и К2 нелинейных возмущений ф(х, ?).
ПРИЛОЖЕНИЕ
Доказательство теоремы 1. Нужно доказать, что если А(г) — матрица, соответствующая невозмущенной системе (1), а функция q(t) = divA(г)x удовлетворяет условию
q(t) l C = const > О,
(П.1)
то система (1), имея неустойчивую точку равновесия х = 0, является по отношению к возмущениям ф(х, г) класса К1 грубой в смысле сохранения характера устойчивости, т. е. при всех возмущениях ф(х, г) е К1 неустойчивой является также точка равновесия х = 0 возмущенной системы (2).
Условие (П.1) означает [8], что точка равновесия х = 0 системы (1) неустойчива (см. в работе [8] следствие теоремы 2.5 или теорему 2.8). С учетом этого условия, кроме того, имеем для дивергенции от правой части системы (2)
div(A(г)x + ф(х, г)) = divA(г)x + divф(x, г) >
> С + divф(x, г), (П.2)
где ф(х, г) е Согласно дополнительному условию к возмущениям ф(х, г) е К1 в некоторой (соответствующей возмущению ф(х, г)) окрестности е точки х = 0 при любом г е R+ выполняется соотношение (4) divф(x, г) > м?(х), где м?(х) — непрерывная в окрестности е определенно отрицательная функция (м?(х) = 0 при х = 0, м(х) < 0 при х ф 0). Тогда из соотношения (П.2) получаем, что в окрестности е при любом г е R+
div(A(г)x + ф(х, г)) > С + divф(x, г) >
> С + м?(х) = м^(х), (П.3)
где м^(х) = С + м?(х) — непрерывная в окрестности е функция, причем м^(0) = С + м?(0) = С > 0. Значения функции м^(х) могут быть положительны не во всех точках х е е. Но существует такая окрестность е1 с е точки х = 0, что м1(х) > 0 при х е е1, ибо функция м^(х) не-
= оо
о
АНАЛИЗ И СИНТЕЗ СИСТЕМ УПРАВЛЕНИЯ
т-' і
прерывна в окрестности s и м^(0) = C > 0. Соотношение (П.3), будучи справедливым при x е s, очевидно, справедливо и при x е Sj с s и при этом в соотношении (П.3) wj(x) > 0. Таким образом, в окрестности sj точки равновесия x = 0 возмущенной системы (2) дивергенция от правой части этой системы при любом t е R+ не меньше положительной непрерывной функции wj(x). Согласно работе [8] точка равновесия x = 0 этой возмущенной системы неустойчива. Итак, показано, что система (1), имея при условии (П.1) неустойчивую точку равновесия x = 0, является по отношению к возмущениям p(x, t) класса Kj грубой в смысле сохранения характера устойчивости. Теорема 1 доказана.
Доказательство теоремы 2. Необходимо показать, что если функция q(t) = divA(t)x удовлетворяет условию (8) теоремы 2, то невозмущенная неавтономная линейная система (1), имея неустойчивую точку равновесия x = 0, является по отношению к возмущениям p(x, t) класса K грубой в смысле сохранения характера устойчивости: при любом возмущении p(x, t) е K2 неустойчивой является и точка равновесия x = 0 возмущенной системы (2).
Условие (8) означает, что точка равновесия x = 0 невозмущенной системы (1) неустойчива (теорема 2.8 в работе [8]). Покажем, что тогда точка равновесия x = 0 возмущенной системы (2) также будет неустойчивой при любом возмущении p(x, t) е K2. Для этого покажем, что дивергенция от правой части системы (2)
div(A(t)x + p(x, t)) = divA(t)x + divp(x, t) =
= q(t) + divp(x, t)
удовлетворяет следующему соотношению
t
lim f inf (q(s) + divp(x, s))ds = “, (П.4)
t X € B 0
где B — любое компактное множество, которое содержит точку x = 0 и принадлежит при этом окрестности S этой точки (в окрестности s при любом t е R+ выполняется условие (7) divp(x, t) l 0).
Так как при любом t е R+ функция q(t) не зависит от x и divp(x, t)|x = 0 = 0, а в окрестности s divp(x, t) l 0 при
любом t е R+, то тогда
inf (q(t) + divp(x, t)) l inf q(t) = q(t), t е R+
X € B X € B
и, следовательно, учитывая также условие (8) теоремы (2), получаем
tt lim Г inf (q(s) + divp(x, s))ds l lim [q(s)ds = ^.
t ^ X € B t ^ » •*
00
Таким образом, соотношение (П.4) при условии теоремы 2 выполняется, что согласно теореме 2.7 из работы [8] означает неустойчивость точки равновесия x = 0 системы (2). Теорема 2 доказана.
ЛИТЕРАТУРА
1. Понтрягин Л С. Обыкновенные дифференциальные уравнения. — М.: Наука, 1974.
2. Фихтенгольц Г.М. Основы математического анализа. — М.: Наука, 1968. — Т. 1.
3. Шилов Г.Е. Математический анализ (функции нескольких вещественных переменных). — М.: Наука, 1972.
4. Арнольд В.И. Дополнительные главы теории обыкновенных дифференциальных уравнений. — М.: Наука, 1978.
5. Андронов А.А., Витт А.А., Хайкин С.Э. Теория колебаний. — М.: Физматгиз, 1959.
6. Андронов А.А. Собрание сочинений. — М.-Л.: Изд-во АН СССР, 1956.
7. Математическая энциклопедия. — М.: Сов. энциклопедия, 1977. — Т. 1. — С. 670.
8. Жуков В.П. Полевые методы в исследовании нелинейных динамических систем. — М.: Наука, 1992.
® (495) 334-89-61, e-maг7: vpzhukov@ipu.ru
Статья представлена к публикации членом редколлегии
В.Ю. Рутковским. □
УКАЗ
ПРЕЗИДЕНТА РОССИЙСКОЙ ФЕДЕРАЦИИ
О награждении государственными наградами Российской Федерации За заслуги в научной деятельности присвоить почетное звание
«ЗАСЛУЖЕННЫЙ ДЕЯТЕЛЬ НАУКИ РОССИЙСКОЙ ФЕДЕРАЦИИ»
КУЛЬБЕ Владимиру Васильевичу — доктору технических наук, профессору, заведующему лабораторией Института проблем управления им. В.А. Трапезникова Российской академии наук, город Москва.
Президент Российской Федерации
В. Путин
Москва, Кремль 3 апреля 2008 года
Поздравляем члена редколлегии нашего журнала Владимира Васильевича Кульбу с присвоением почетного звания. Желаем ему доброго здоровья и дальнейших творческих успехов.
Редколлегия
Редакция