Модальная стабилизация некоторого класса нелинейных нестационарных систем Текст научной статьи по специальности «Математика»

УДК 517.938

Вестник СПбГУ. Сер. 1. Т. 2(60). 2015. Вып. 3

МОДАЛЬНАЯ СТАБИЛИЗАЦИЯ НЕКОТОРОГО КЛАССА НЕЛИНЕЙНЫХ НЕСТАЦИОНАРНЫХ СИСТЕМ*

И. Е. Зубер, А. Х. Гелиг

Санкт-Петербургский государственный университет,

Российская Федерация, 199034, Санкт-Петербург, Университетская наб., 7/9 Рассматривается система

— = А(Ь, х)х + Ъ(Ь, х)и, (1)

<и

где А(4, ж) € КпХп, 6(4, ж) = в(4, ж)вь в(4, ж) € К1, в*г = (1, 0,..., 0). Коэффициенты а^(4, ж) матрицы А(4, ж) и в(4, ж) являются ограниченными вещественными непрерывными функциями. Предполагается, что система равномерно управляема:

М |ае^Ап-1(4,ж)6(4,ж),...,А(4,ж)Ь(4,ж),Ь(4,ж)| > 0.

£>¿0,

Стабилизирующее управление определяется формулами

«(4, ж) = «*(4, ж)ж, (2)

«*(*, ж)(А(4, ж) - Л, (4, ж)/)-16(4, ж) = -1, где Л2 (4, ж) —специально выбранный в левой полуплоскости спектр замкнутой системы. Найдены условия, налагаемые на А(4, ж) и в(4, ж), при выполнении которых система (1), (2) глобально экспоненциально устойчива.

В случае, когда а^ (4, ж) и в(4, ж) ограничены лишь при |ж| < Я, построенное управление стабилизирует систему в большом, то есть при |ж(4о)| < Я.

Построено модальное стабилизирующее управление для систем с переменной структурой, у которых

А(4,ж) = Ак(4, ж), в(*,ж)= вк (*,ж) при 4 € (^ ,4к+1),

где 4к+1 > 4к, 4к ^ при к ^ те. Библиогр. 7 назв.

Ключевые слова: нелинейные системы, нестационарные системы, модальная стабилизация, глобальная экспоненциальная устойчивость.

1. Введение. Рассмотрим систему

— = Ах + Ъи, (1)

at

u = s*x, (2)

где A G Rnxn, b G Rnx1, s £ Rnx1, * — знак транспонирования (все величины вещественные). Если A и b являются достаточно гладкими функциями от x, то для определения s(x), при котором система (1), (2) глобально устойчива, можно использовать метод «feedback-linearization» [1]. Однако этот метод применим лишь при условии, что последняя строка матрицы, обратной к матрице управляемости (составленной с помощью скобок Ли), образует потенциальное поле [2]. Относительная степень таких систем n — 1. Стабилизация достаточно гладких нелинейных стационарных систем с относительной степенью, меньшей n — 1, рассматривалась в [3]. В [4] эти результаты были перенесены на достаточно гладкие нестационарные нелинейные системы.

*Работа выполнена при финансовой поддержке РФФИ (проект 14-01-00107) и СПбГУ (тема №6.38.230.2015).

Если в системе (1), (2) А и Ь постоянные и выполнено условие управляемости

...,АЬ,Ь) = о, (3)

то известен метод модальной стабилизации [5], который заключается в определении в из системы алгебраических уравнений

в*с^ = —1, г £ 1, п,

где ¿1 = (А — Л$/)-1 Ь, / — единичная п х п-матрица, Л^ —числа, расположенные в левой полуплоскости и не принадлежащие спектру матрицы А. При этом числа Л^ оказываются спектром матрицы А + Ьв*.

В данной статье метод модальной стабилизации распространен на некоторый класс негладких нестационарных нелинейных систем при условии, что система является равномерно управляемой [6], то есть условие (3) выполнено при всех £ > ¿о,

х е К".

2. Модальная стабилизация. Рассмотрим систему

— = Аи, х)х + Ъи, х)и, (4)

от

и = в* (£, х)х, (5)

где t > ¿о, х е К", Ь(£, х) = х)в1, е* = (1,0,..., 0). Предполагается, что элементы ац (¿,х) матрицы А(£, х) и скалярная функция в(£,х) непрерывны и удовлетворяют следующим условиям:

0 < в1 < в(£,х) < в < ^ при £ > ¿0, х е К", (6)

эир х)\ = < оо приг(Е1,п, (7)

х) при г £ 2,п, j £ 1,п, j ^ г,

_ (8)

—х) + ж) при = г £ 2, п,

где £ > 0,

inf а(£,х) = ао > пк, (9)

эир |ау(£,х)| < при г € 2, п, ] £ 1, п. (10)

Здесь и далее к — положительные константы. Рассмотрим определитель управляемости

Д„(г, х) = ае1(А"-1(£, х)Ь(£, х),..., А(£, х)Ь(£, х), Ь(£, ж)).

Очевидно представление

Д„(£,х) = в"(*,х)Д(*,х), (11)

где Д(£, х) = det(A"-1(t, х)е1,..., А(£, х)е1, е1). В силу свойства (8) Д(£, х) зависит от £.

Справедлива формула

Д(£,х)= £0,5"("-1)До (¿,х), (12)

где До (£, х) от е не зависит. Чтобы убедиться в этом, надо разложить Д(£, х) по элементам последнего столбца. В полученном определителе (п — 1)-го порядка последний столбец имеет вид е(«21(^, х),..., ап1 (£, х))*. Элементы к-го столбца, отсчитываемого справа, являются полиномами относительно е степени к, причем они обращаются в ноль при е = 0. Умножив последний столбец на а(£, х) — ац(4,х) и прибавив его к предпоследнему столбцу, мы ликвидируем в предпоследнем столбце члены с е в первой степени. Продолжая этот процесс, мы ликвидируем в к-м столбце члены с множителями е® (г < к) и приходим к соотношению (12). Из формул (11), (12) и свойства (6) вытекает оценка

|ДИ(М)| > вГе°'5"(п-1)|До(М)|.

Таким образом, условие равномерной управляемости сводится к требованию

Ы |Д°(¿,х)| = > 0, (13)

и для Д„(г, х) получается оценка

|ДИ(М)| > КЗе°'5п(п-1), где КЗ = вП¿0. (14)

Зададим спектр замкнутой системы в виде

Xj(t,x) = — х) — ] € 1,п, (15)

где ф^- —различные положительные числа, и определим х) в формуле (5) из алгебраических уравнений

в*(г,х)<1з(г,х) = -1, (16)

где ^(£, х) = (А(£, х) — А^ (£, х)/)-16(£, х). Выберем теперь числа ф^ в формуле (15) таким образом, чтобы числа (15) были отделимы от спектра матрицы А(£, х). Обозначим через х) собственные числа матрицы А(£, х). Согласно кругам Гершгорина [7] и свойству (10) при г > 1 справедлива оценка

|м®(^,х) — (—«(¿,х) + еай(г,х))| < ^ (¿,х)|е < К2е.

гфз^. 1,п

Отсюда следуют неравенства

п

|м®(^,х) — (—«(¿,х))| < ^¿(¿,х) — (—«(¿, х)) — еай(г, х)| + е|ай(г,х)| < ^ (¿, х)|е.

5=1

Поэтому в силу (10) имеет место оценка

х) — (—«(¿,х))| < ек2.

Пусть фг > £К2 (г €= 1,"-)- Тогда все числа \(1,,х) отделимы от х), ..., х). Согласно кругам Гершгорина

(£, х) = ац(£, х) + V(£, х),

где

"

И*,х)| <Х>и (¿,х)|.

Ц=2

Поэтому справедлива оценка

"

М£,х)| |а1Ц(£,х)1, ц=1

откуда в силу (7) вытекает неравенство

|м(г,х)| < пк1.

Из этого неравенства и свойства (9) следует, что все х) отделимы и от х).

Поскольку х) являются собственными числами матрицы Р(£, х) = А(£, х) + Ь(£, х)в*(£, х), можно рассмотреть спектральное разложение этой матрицы

Р (¿,х) = ^ ^(¿,х)М^,х), (17)

¿=1

где

^¿(¿,х) = ¿¿(¿,х)^(^х), (18)

¿¿(¿, х) и х) —собственные векторы матриц Р(£, х) и Р*(£, х) соответственно, обладающие свойством

{1 при г = 7,

(19)

0 при г = 7.

Построим матрицы С(£, х) = (¿1 (£, х),..., ¿"(¿, х)) и С(£,х) = (01 (£, х), ..., х)). Из свойства (19) вытекает равенство х)С(¿,х) = /. Поэтому

С(£,х) = (С* (¿,х))-1. (120)

Для замкнутой системы

^=Р(г,х)х (20)

рассмотрим функцию Ляпунова V(х) = х*/х. Производная ¿У/А, взятая в силу системы (20), с учетом формулы (17) имеет вид

- = ж I) I, ж ж,

А V . у .

где ^(¿,х) = Л¿(t, x)(M¿(t, х) + М*(£,х)). В силу (15) справедливо представление

""

£(г,х) = — а(£,х) ^(^¿(¿,х) + М*(£,х)) + £ ^ х) + М*(£,х)).

¿=1 ¿=1

Имеет место равенство

" ¿=1

в чем легко убедиться, умножив его слева на С(£, х) и воспользовавшись свойством (19). Поэтому в силу условия (9) справедлива оценка

D(t,x) < -2а0/ + eñ(t,x),

(21)

где

|R(t,x)| < 2V>o¿ |Mi(t,x)|,

(22)

фо = тах фг, | • |—евклидова норма матрицы. Нашей дальнейшей целью является

гЕ 1,п

оценка х). Из (18) вытекает неравенство

|Мг(£, х)| < ^(¿,х)| • х)|. (23)

Очевидно, что ¿¿(¿, х) имеет вид

^¿(¿,х)

1

di(t, x) =

Ai(í,x)

(24)

где Д¿(t, х) = det(A(t, х) — Л¿(t,ж)/), х) —миноры первой строки матрицы

(А(£, х) — Л¿/). Очевидна оценка

|V>fc¿(í,x)| < х4е™ (A; G 1,п).

(25)

Здесь и в дальнейшем к — абсолютные константы, не зависящие от £. Из (24) вытекает оценка

|di(t,x)| < Ввиду (24) матрица C*(t, x) имеет вид

Кб£

1-1

|Ai(í,x)| '

(26)

C *(t,x)

/ ipii(t,x) Ai (4,ж)

1p2l(t,x)

A2{t,x)

lpnl{t,x)

\ а

j>ln{t,x) \ Ai (i,®)

1p2n(t,x)

A 2(M)

Ipnnjt, x) An(t,x) /

Пусть g*(t, x) = (7¿1 (t,x),... ,7¿„(t, x)). Очевидно, что 7¿i(í, x),..., 7¿n(t,x) являются алгебраическими дополнениями г-й строки матрицы C* (t,x), деленными на det C(t, x). Поэтому

7¿fc (t,x)

wifc(í,x)

где в силу (26) справедливо неравенство

|w¿fc(í, x)| < Кб£

(n-1)2

Отсюда вытекает оценка

к7е(п-1)2

Для det С(£,х) имеет место следующая элегантная формула:

п

det С (*,х) = ПД-1(*,х) П (А® — Ад- )Д„(*,х). (28)

®=1 1<®<5<п

п

Для ее доказательства надо матрицу С(£, х) умножить слева на ^ (А(£, х)—А® (£, х)1) и

®=1

воспользоваться перестановочностью матриц (А(£, х) — А® (£, х)1) и (А(£, х) — А^ (£, х)1). Поскольку А® — А^ = е(ф^ — ф®), из (28) и (14) вытекает оценка

п

^С(¿,х)| > к8Ц ^х) — еп(п-1).

®=1

Отсюда в силу (27) получаем соотношение

|Л(*,х)| < кд|Д®(£, х)|е(1-п), из которого ввиду (26), (23) следует оценка

|М®(^х)| < К1°. Поэтому согласно (22) имеет место свойство

|Д(г,х)| < кц

и, следовательно, спектральная норма

||Д(*,х)|| = ^/^¡Щад

имеет такую же оценку:

||Д(*,х)|| < кц.

Отсюда ввиду симметричности матрицы Д(£, х) вытекает соотношение

|Ашах|Д(^, х)}| < Кц . Поэтому согласно неравенству Рэлея

Д(£, х) < кцI и, следовательно, в силу (21) имеет место неравенство

^(¿,х) < —(2а° — ек11)1. Таким образом, система (20) глобально экспоненциально устойчива при условии

£ < (29)

К11

Сформулируем полученный результат.

Теорема 1. Пусть выполнены условия (6)-(10), (13), (29) и вектор s(t,x) определяется из уравнений (15), (16). Тогда система (4), (5) глобально экспоненциально устойчива.

3. Стабилизация в большом. Требование равномерной ограниченности во всем пространстве коэффициентов системы (4) исключает из рассмотрения системы, у которых коэффициенты являются полиномами относительно составляющих вектора x. Однако предложенный метод модальной стабилизации применим и к системам с неограниченно растущими коэффициентами, если начальные условия принадлежат некоторому шару M = {x : |x|2 < Д2}. Действительно, если в условии теоремы 1 sup и inf брать при t > to, |x| < Д, то все решения с начальными условиями, принадлежащими шару M, будут при t — стремиться к устойчивому по Ляпунову состоянию равновесия x = 0.

4. Системы с переменной структурой. Предположим, что в системе (4) при t G (tfc, tfc+i) A(t, x) = Afc(t, x), e(t, x) = вк(t, x) и tfc —у то при k — то.

Теорема 2. Если Ak (t,x) и вк (t, x) удовлетворяют условиям теоремы 1 при всех k, то система с переменной структурой будет глобально экспоненциально устойчива, если при tk < t < tk+1 управление определяется формулой

u == Sk(t, x)x, в которой sk(t, x) находится из системы

s*k{t,x){Ak{t,x) - Xkj(t,x))-1ßk(t,x)e1 = -1 (j G T~rä),

где Дк (t,x) = —ak(t,x) — еф.

Доказательство очевидно, поскольку функция Ляпунова и оценка ее производной не зависят от k. Замечание об устойчивости в большом справедливо и для системы с переменной структурой.

5. Заключение. Рассмотрен некоторый класс негладких нестационарных нелинейных равномерно управляемых систем, к которым применим метод модальной стабилизации, заключающийся в специальном выборе спектра замкнутой системы и определении коэффициентов обратной связи в виде решения линейной системы алгебраических уравнений. При этом замкнутая система становится глобально экспоненциально устойчивой. Показано, что этот метод применим для стабилизации некоторого класса нелинейных нестационарных систем с переменной структурой, а также для стабилизации в большом, когда коэффициенты ограничены лишь в заданном шаре пространства состояний и начальные условия берутся из этого же шара.

Литература

1. Isidory A. Nonlinear Control Systems. II. London: Springer, 1999. P. 520.

2. Zak S. H. Systems and Control. Oxford: Oxford Univ. Press, 2002. P. 490.

3. Крищенко А. П., Панфилов Д. Ю., Ткачев С. Б. Глобальная стабилизация афинных систем с помощью виртуальных выходов // Дифференциальные уравнения. 2003. Т. 39, №11. С. 1503—1510.

4. Ткачев С. Б. Стабилизация нестационарных афинных систем с помощью виртуальных выходов // Дифференциальные уравнения. 2007. Т. 43, №11. С. 1507—1517.

5. Гелиг А. Х., Леонов Г. А., Якубович В. А. Устойчивость нелинейных систем с неединственным состоянием равновесия. М.: Наука, 1978. 400 с.

6. Зубер И. Е., Якубович В. А., Гелиг А. Х. Стабилизация некоторого класса равномерно управляемых систем // Доклады Академии Наук. 2012. Т. 446, №6. С. 1—3.

7. Го,нт.мо,хер Ф. Р. Теория матриц. М.: Наука, 1967.

Статья поступила в редакцию 26 марта 2015 г.

Сведения об авторах

Зубер Ирина Ефремовна —доктор технических наук, ведущий научный сотрудник; zuber-yanikum@mail.ru

Гелиг Аркадий Хаимович — доктор физико-математических наук, профессор; agelig@yandex.ru

MODAL STABILIZATION OF A CERTAIN CLASS OF NONLINEAR NONSTATIONARY SYSTEMS

Irina E. Zuber, Arkadiy Kh. Gelig

St.Petersburg State University, Universitetskaya nab., 7/9, St.Petersburg, 199034, Russian Federation; zuber-yanikum@mail.ru, agelig@yandex.ru The system

— = A(t, x)x + b(t, x)u (1)

dt

is considered. Here A(t, x) € RnXn, b(t, x) = e(t, x)eb e(t,x) € R1, e\ = (1, 0,..., 0). Matrix coefficients aij (t, x) for A(t, x) and e(t, x) are real bounded continuous functions. It is supposed that system (1) is uniformly controllable:

inf |det(An-1(t,x)b(t,x),...,A(t,x)b(t,x),b(t,x)| > 0.

t>tо, x£Rn

Stabilizing control is defined by formula

«(t,x) = s*(t, x)x, (2)

s*(t, x)(A(t, x) - Xj(t, x)J)-1b(t, x) = -1,

where Xj (t, x) are special choised functions of spectrum for matrix of closed-loop system. Such conditions are formulated for A(t, x) and e(t, x) so that fulfilment these conditions provides global exponential stability for closed-loop system (1), (2).

In the case then aij(t, x) and e(t, x) are bounded only for |x| < R the constructed modal control stabilizes the closed-loop system in the large, i.e. |x(to)| < R.

The modal stabilizing control is constructed for systems with varying structure where

A(t,x) = Ak(t,x), e(t,x) = вк(t,x) by t € (tfc,tfc+i).

Here tk+i > tk, tk ^ for к ^ <x. Refs 7.

Keywords: nonlinear system, nonstationary system, modal stabilization, globally exponentially stability.

References

1. Isidory A., "Nonlinear Control Systems" II (Springer, London, 1999).

2. Zak S.H., "Systems and Control" (Oxford Univ. Press, Oxford, 2002).

3. Krishchenko A. P., Panfilov D. V., Tkachev S. B., "Global stabilization of affine systems via virtual outputs", Differential Equations 39(11), 1585-1592 (2003).

4. Tkachev S. B., "Stabilization of Nonstationary Affine Systems by the Virtual Output Method", Differential Equations 43(11), 1546-1557 (2007).

5. Yakubovich V. A., Leonov G.A., Gelig A. Kh., "Stability of Stationary Sets in Control Systems with Discontinuous Nonlinearities" (World Scientific, Singapore, 2004).

6. Zuber I.E., Yakubovich V. A., Gelig A.Kh., "Stabilization of a Certain Class of Uniformly Contrallable Systems", Doclady Mathematics 86(2), 1-3 (2012).

7. Gantmacher F.R., The Theory of Matrices (Chelsea, 1959).