Стабилизация двумерных линейных управляемых систем обратной связью с запаздыванием по Пирагосу

Шумафов Магомет Мишаустович

МАТЕМАТИКА MATHEMATICS

УДК 517.977.1 ББК 22.19 Ш 96

Шумафов М.М.

Доктор физико-матаматических наук, профессор, зав. кафедрой математического анализа и .методики преподавания .математики факультета .математики и компьютерных наук Адыгейского государственного университета, Майкоп, тел. (8772) 59-39-05, e-mail: shumaf@mail.ru

Стабилизация двумерных линейных управляемых систем обратной связью с запаздыванием по Пирагосу

(Рецензирована)

Аннотация. Рассматривается задача стабилизации по выходу двумерных линейных управляемых систем обратной связью с запаздыванием по Пирагосу. Даны необходимые u/imu достаточные условия стабилизируемости двумерных систем. Показаны возможности стабилизации линейных управляемых систем второго порядка обратной связью с запаздыванием по Пирагосу. Полученные результаты сравниваются с необходимыми и достаточными условиями стабилизации обратной связью без запаздывания.

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

Shumafov M.M.

Doctor of Physics and Mathematics, Professor, Head of Department of Mathematical Analysis and Methodology of Teaching Mathematics of Mathematics and Computer Science Faculty, Adyghe State University, Maikop, ph. (8772) 59-39-05, e-mail: shumaf@mail.ru

Stabilization of two-dimensional linear controllable systems by Pyragas' delay feedback

Abstract. In the paper the problem of stabilization of two-dimensional linear time-invariant controllable systems by means of Pyragas ' static time-invariant delayed output feedback is considered. Necessary and/or sufficient conditions of stabilizability of two-dimensional systems are given. The potential of Pyragas ' delayed feedback is shown for stabilization of two-dimensional systems. The obtained results are compared with necessary and sufficient stabilization conditions given by feedback without delay.

Keywords: linear controllable system, delay feedback by Pyragas ', quasipolynom, asymptotic stability, D-decomposition.

1. Введение

Управление динамическими системами является одной из классических задач в прикладных науках, а стабилизация систем - одна из важнейших проблем в теории управления. Различные вопросы, касающиеся проблемы стабилизации динамических систем, интенсивно изучались в последние четыре десятилетия и в настоящее время эти вопросы остаются в центре внимания многочисленных исследователей (см. библиографию в [1], а также обзорные статьи [2, 3]). Повышенный интерес к проблемам стабилизации мотивируется прежде всего запросами практики управления. Наиболее эффективные алгоритмы стабилизации разработаны для линейных стационарных систем. Следует отметить, что в линейной теории управления остается еще много нерешенных трудных задач (см. [3]).

Возникает вопрос: каковы возможности линейной стационарной обратной связи с запаздыванием для стабилизации линейных неустойчивых стационарных систем?

Хорошо известно [4], что для достаточно малых и достаточно больших запаздываний такая стабилизация невозможна. Одной из мотиваций к исследованию задач ста-

билизации путем введения запаздывания в обратную связь явились компьютерные эксперименты К. Пирагоса (K. Pyragas) [5-9] по стабилизации хаоса. В своей пионерской работе [5] Пирагос предложил весьма мощный и эффективный метод стабилизации неустойчивых периодических орбит, встроенных в странные аттракторы динамических систем с хаотическим поведением, в частности, в системах Лоренца и Ресслера. Этот метод состоит во введении в систему обратной связи с запаздыванием специального вида, которая формируется пропорционально разности значений выходного сигнала в текущий момент времени и значения этого сигнала в некоторой предыдущий момент времени. Первоначально примененный для стабилизации неустойчивых периодических орбит этот метод оказался эффективным и для стабилизации неустойчивых состояний равновесия [10-13].

В настоящей статье задача линейной стационарной стабилизации с запаздыванием решается для случая двумерных линейных управляемых систем. Здесь получены необходимые и/или достаточные условия стабилизируемости неустойчивых линейных систем второго порядка с постоянными коэффициентами путем введения в систему обратной связи с запаздыванием по Пирагосу. Эти условия получены на основе D-разбиения пространства параметров рассматриваемой системы. Оказывается, что линейная система в седловом случае не является стабилизируемой. Некоторые частные результаты по стабилизации рассматриваемой системы были анонсированы в [14, 15]. Вопрос о стабилизации рассматриваемой системы с помощью стандартной обратной связи с запаздыванием рассматривался в работе [16].

2. Постановка задачи

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

x(t ) = Ax(t)+bu ((), y = cx(t), (1)

где x(t)e R2 - вектор состояния в момент времени t; u(t)e R - управление (вход), y(t)e R - выход; A, b, c - вещественные постоянные матрицы размеров 2 х 2,2х 1 и 1х 2 соответственно.

Системы вида (1) представляют собой важный класс для анализа нелинейных управляемых систем в окрестности состояния равновесия (x = 0, u = 0). К системам вида (1) также сводятся линейные периодические системы, возникающие как системы уравнений в вариациях в окрестности их периодических орбит.

Предположим, что система (1) управляема. Введем в систему (1) обратную связь по Пирагосу [5]

u(t ) = -k (y(t)-y(t -т)), (2)

где k Ф 0 и т > 0 - варьируемые параметры.

Основная задача заключается в том, чтобы найти такие значения параметров k и т, чтобы система (1), замкнутая обратной связью (2), оказалась бы асимптотически устойчивой.

Замкнутая система (1), (2) представляет собой линейное дифференциальное уравнение с запаздывающим аргументом относительно вектор-функции x(t), где k и т -параметры.

3. Формулировка результата

Приведем систему (1) невырожденным линейным преобразованием координат к каноническому виду [17, с. 239]

X1 (t )= X2 (t).

X 2(t) = -«1 Xi(t) - a2 X2(t) - u(t), y(t)= CiX1 (t)+ C2x2 (t),

(3)

где a1, a2, c1, c2 - вещественные параметры. Приведение системы (1) к виду (3) возможно в силу управляемости системы (1).

В работах [14, 15] были даны необходимые и достаточные условия стабилизируе-мости системы (3) в случаях

1) c1 Ф 0, c2 = 0, 2) c1 = 0, c2 Ф 0 и в случае, когда c1 > 0, c2 = 1.

В настоящей статье рассматривается случай произвольных c1 и c2 таких, что c1 Ф 0, c2 Ф1. В этом случае здесь нами получено достаточное условие стабилизируе-мости системы (3) обратной связью (2). Это условие расширяет область стабилизации {a2 > cxj, полученной в [14, 15].

Введем в рассмотрение числа

sin а

. ( sin( m0 = min cos(H--

ce[0,2п]^ ( J

m

■ min -

ee[0,2п] (

t

m„

-1,0419, т ~ -0,2172). Имеет место следующая

Теорема. Пусть в системе (3) с1 Ф 0, с2 Ф 0. Тогда система (3) стабилизируема обратной связью (2), если выполнено хотя бы одно из условий

а) схс2 > 0, аг > 0 ; Ь) с1с2 < 0, аг > 0, а2 > с1/с2 ;

с1 п с) с1с2 < 0, a1 > — a2 +

Г Л 2

4

V C2 J

a2

- те < a2 <

2с2 тф2

"2 '

/

0,

2с1

тс

2J

т$- 2сх/ с2

Замечание. В обоих случаях схс2 > 0 и схс2 < 0 условие а1 > 0 является необходимым для стабилизируемости системы (3) обратной связью по Пирагосу. Таким образом, в случае схс2 > 0 условие а) теоремы является необходимым и достаточным.

4. Доказательство теоремы

Замкнутая система (3), (2) имеет вид

x(t )= (A - kbc )x(t) + kbcx (t -t),

где

Г 0 1 > 0

A = , b = -1

V-a1 -a2

C = (c1 C2 )•

(4)

(5)

Характеристическое уравнение системы (4), (5) представляет собой квазиполином ^(; а1, а2; с1, с2; к, т) = г2 +(а2 -кс2) + (а1 - ксг)+ке~т (с2г + с1). (6)

Как хорошо известно [4, с. 118], для асимптотической устойчивости линейной системы (4), (5) необходимо и достаточно, чтобы все корни характеристического полинома (6) лежали в левой полуплоскости Иег > 0 плоскости комплексного переменного г е С.

Далее доказательство теоремы аналогично доказательству теоремы 1 из [16] или теоремы 3 из [15]. Оно основано на применении метода О-разбиения пространства

2

с

1

R2 = (а15 а2} коэффициентов квазиполинома (6) с последующим использованием теоремы Руше о нулях аналитической функции [18, с. 160].

Определим границу Гкх D-разбиения пространства коэффициентов R 2а квазиполинома F из (6). Она представляет собой образ мнимой оси z = iy плоскости комплексного переменного z в R 2а при отображении z ^ а = (а1, а2), определяемом уравнением F(z;а;c)= 0 (c = (c1,c2) - фиксировано). Имеем F(iy;а;c)= 0, или после разделения вещественной и мнимой частей

[а1 = y2 - kc2y sin (ту)- Ác^cos (y)+kc1, \а2 = kc1sin(ry )/ y - kc2cos(ry)+kc2,

при y Ф 0 и а1 = 0 при y = 0.

Таким образом, граница Гк Т состоит из особой прямой K = {а1 = 0} и кривой Lk т, определяемой параметрическими уравнениями (7): Гкт= K u LkT. В силу четности функций а1 (y) и а2 (y) достаточно рассмотреть (7) для y > 0.

Кривая Lk т, выходя (при y = 0) из точки (а1 = 0, а2 = kTc1), при достаточно малых Т колебательным образом поступательно навивается на ось абсцисс а1, а при больших значениях Т она движется возвратно-поступательно, т.е. возникают точки самопересечения и, следовательно, кривая Lk Т содержит петли. Поэтому при малых т кривая (7)

устроена проще. С помощью теоремы Руше будет показано, что выше кривой (7) располагается область устойчивости S . Затем подберем значения k, т такие, чтобы точка (а1 > 0, а2) оказалась в области S. Проще всего это сделать для малых т . Ясно, что условием отсутствия петель является неравенство а'(y )> 0.

Обозначим через DkT область в плоскости R 2 параметров а1 и а2, точкам которых соответствуют квазиполиномы вида (6) с одинаковым числом p корней в правой полуплоскости Rez > 0 плоскости комплексного переменного z .

Как и в доказательстве теоремы 3 из [15], применяя теорему Руше к функциям

f (z;а1,а2) = z2 + (а2 -kc2)z + (а1 -kc1)

и

(p(z) = k (c2 z + c1 )exp(- Tz),

устанавливаем, что единственной областью DlT (p = 0) будет область, лежащая выше кривой LkT в правой полуплоскости а1 > 0. Возможны четыре случая:

A) c1 > 0, c2 > 0; B) c1 < 0, c2 < 0; C) c1 < 0, c2 > 0; D) c1 > 0, c2 < 0. Случай A: c1 > 0, c2 > 0.

Пусть (а1, а2) - произвольная точка с а1 > 0. Будем считать k < 0. Кривая LkT начинается (при y = 0) в точке (0, kTc1) на отрицательной полуоси а2. Возьмем отрицательное число в < а2. Выберем k и Т такими, чтобы kTc1 = в. При таком выборе параметров k и Т производная

а'(y) = (o¡Tcl )c1 + (- в)2 (cos а + (1 -Tc^lc2 )sin (J¡а)], ( = Ty), (8)

заведомо положительна на интервале (0, п/ 2т) , а производная

a2(y) = (вт/j2hcosa+(c2a2¡щ -1)sin j] (9)

отрицательна на интервале (0, 70 (т)/т) при достаточно малых т, где 70 (т) - наименьший положительный корень уравнения

ctgj = ( - c2 j2 )тс1<. (10)

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

(Здесь 70(т)е (я/2, я); 70(т) —я при т —0.)

Заметим, что при т — 0

da2 /da1 (0; т) = вФ^ - Tc1 )3(2c1 - 2в + pTcx) — -0 и da2¡da1 < 0 на том участке кривой LkT, который соответствует промежутку 0 < y < я/ 2т.

Возьмем прямую l, проходящую через точку (0, в) оси a2 и точку (a1 (я/т), a2 (я/т)) на кривой Lk т. Так как

a1 (я/т) = (я2 + 2рт)1 Т, a2 (я/т) = 2в Щ,

то при достаточно малом т угловой коэффициент прямой l будет отрицательным и, следовательно, точка (a1, a2) с a1 > 0 и a2 > в окажется выше прямой l, а значит, и выше кривой LkT, то есть (a1, a2 )е DkT при достаточно малом т и соответствующем k = P/tc1. Таким образом, стабилизация системы (6) обратной связью (2) возможна при любом a1 > 0 и любом a2 .

Случай B: c1 < 0, c2 < 0.

Пусть (a1, a2) - произвольная точка с a1 > 0. Будем считать k > 0. Выбираем значения в, k и т, как и выше для случая A. Тогда совершенно аналогично, как и в случае A, получаем, что при достаточно малых т точка (a1, a2 )e D^, т.е. стабилизация возможна при любом a1 > 0 и любом a2 .

Случай C: c1 < 0, c2 > 0.

Пусть (a1, a2) - произвольная точка на плоскости R га с a1 > 0 и a2 > c1¡c2. Пусть k > 0. Возьмем число в, удовлетворяющим условиям: c1¡ c2 <в< a2, если a2 < 0, и в < 0 , если a2 > 0. Выберем k и т такими, чтобы krc1 = в. Тогда из (8) и (9) следует, что при достаточно малых т соответствующая функция a1 (y) будет возрастающей на всей полуоси [0; (в силу неравенства в > c1¡c2), а функция a2(y) - возрастающей на промежутке [0;<0(т)/т) и убывающей на интервале (<0(т)/т;2я/я), где <г0(т) -наименьший положительный корень уравнения (10). (Здесь /<<0(т) < 3я/ 2, 7°(т) —я при т — 0.) Поскольку функции (jcosj-sinj))<7 и sinj/j ограничены на полупрямой [0, (их пределы при 7 — +0 существуют), то производная

da2 _Т¡тccl^c2)• (jcosj-sinj)/j3 + sin jj] << = т) da1 2c1/c2в - [cos 7 + (1 -тcl|c2 )sin j j]

по абсолютной величине будет сколь угодно малой для всех je [0; + при достаточно малых т . Поэтому при некотором достаточно малом т > 0 кривая Lk т окажется ниже точки (a1, a2), и, следовательно, (a1, a2 )e D\T при достаточно малом т и соответствующем k = в|ccí . Значит, стабилизация системы (6) возможна при a1 > 0 и a2 > c1¡c2.

Пусть теперь k < 0. Тогда в> 0, где в = Ыс1. Из условия а'(у )> 0 отсутствия петель у кривой Lk т имеем из (8):

со87+(1 -т^/с2)п7/7> 2с,/сгр. (11)

Неравенство (11) будет выполнено, если

m - m1тc1/ с2 > 2с1/ с2в. (12)

Будем брать в из интервала (0, 2с 1/mc2). Из неравенства (12) получаем промежуток изменения т> 0, для каждого значения которого соответствующая функция а1( у) - возрастающая на всей полуоси [0, . Этот промежуток зависит от в:

0<т<т (в), т (в) = (mc2в-2с)/^св, в е (0,2с,/mc2). Далее, производная а'2(у) из (9) отрицательна на промежутке [0, 70(т)), положительна на (<т°(т)/т,2п/т] и а'(<т°(т)/т) = 0, где <7°(т) - наименьший положительный корень уравнения (10) (7° (т)^п при т^ 0). Поэтому точка Q(a1 (70 /т), а2 (70 /т)), 7° = 7° (т), на кривой Lk т будет самой нижней на том участке этой кривой, который

соответствует значениям у е (0,2п/т). При изменении т в промежутке (0, т*(в)] точка Q будет описывать нижнюю границу области, где возможна стабилизация. Чтобы получить явные алгебраические уравнения границы Г подобласти О в области стабилизации возьмем «соседнюю» с точкой Q точку P на кривой Lk т, соответствующую значению у = л/т, и найдем ее координаты. Имеем

\а1 = +п2/т2,

I — 2kc2,

(13)

где k = в/щ , те (0, т (в)]. Из (13), учитывая, что

— < k < Г (в), г (в)= mlв

mß — 2 c1/ c2

получаем уравнение границы Гв

c1 п

a1 — —1L a2 +--

1 c2 2 4

2 ( Л2 2 о /32

a2 ___ ___ 2mic2ß

c1

V c2 У

— ^ < a2 <

в2 2 mв-2сх!сг

где 0<в< 2с1/mc2. Соответствующая подобласть Ов задается неравенствами

c, п

a, > —1L a2 +--

1 c2 2 4

2 ( Л2 2 О п2

a2 ___ ^ ___ 2mic2ß

c1 V c2 У

— те < a2 <

ß2 mß — 2 cj c

2

Объединение всех областей О в (0 <в< 2 с,/ mc2) будет областью, где возможна стабилизация системы (3).

Случай Б: с1 > 0, с2 < 0.

Пусть (а1, а2) - произвольная точка с а1 > 0 и а2 > с1/с2. Пусть k < 0. Выбрав k и т так, чтобы kтc1 = в, где число в удовлетворяет неравенствам с1/с2 <в<а2, если а2 < 0 и неравенству в < 0, если а2 > 0, как и выше в случае С, получим, что стабилизация системы (3) возможна при указанных выше значениях параметров а1 и а2.

2

Пусть теперь k > 0. Тогда совершенно аналогично, как и в предыдущем случае С, для k < 0 придем к заключению, что область, определяемая неравенствами из условия с) теоремы, является подобластью О области стабилизации Теорема доказана полностью.

5. Заключение

В статье получены необходимые и/или достаточные условия стабилизируемости двумерных линейных стационарных систем обратной связью с запаздыванием по Пира-госу. Выявлены возможности стабилизации двумерных систем обратной связью по Пи-рагосу. В случае с1с2 > 0 область стабилизации по Пирагосу уже по сравнению с областью стабилизации (К1 ), даваемой обратной связью без запаздывания и() = £у(). В

случае с1с2 < 0 условие с) теоремы расширяет область стабилизации, даваемая теоремой 3 из [15].

Примечания:

1. Леонов ГА., Шумафов М.М. Методы стабилизации линейных управляемых систем. СПб.: Изд-во СПбГУ, 2005. 420 с.

2. Static output feedback: a survey / V.L. Syrmos, C.T. Abdallah, P. Dorato, K. Grigoriadis // Automática. 1997. Vol. 33, No. 2. P. 125-137.

3. Поляк Б.Т., Щербаков П.С. Трудные задачи линейной теории управления. Некоторые под// .

2005. № 5. С. 7-46.

4. . ., . . -рию дифференциальных уравнений с отклоняющимся аргументом. М.: Наука, 1971. 296 с.

5. Pyragas K. Continuous control of chaos by self-controlling feedback // Phys. Lett. A. 1992. Vol. 170. P. 421-428.

6. Pyragas K. Control of chaos via extended delay feedback // Phys. Lett. A. 1995. Vol. 206. P. 323330.

7. Nomajunas A., Pyragas K., Tamasevicius A. Stabilization of an unstable steady state in a Mackey-Class system // Phys. Lett. A. 1995. Vol. 204. P. 255-262.

8. Pyragas V., Pyragas K. Delayed feedback control of the Lorenz system: An analytical treatment at a subcritical Hopf bifurcation // Phys. Lett. E. 2006. Vol. 73. 036215. P. 1-10.

9. Delayed feedback control of periodic orbits without torsion in nonautonomous systems: Theory and experiment / A. Tamasevicius, G. Mycolaitis, V. Pyragas, K. Pyragas // Phys. Rev. E. 2007. Vol. 76. 026203. P. 1-6.

10. Hövel P., Schöll E. Control of unstable steady states by time-delayed feedback methods // Phys. Rev. E. 2006. Vol. 72. 046203.

11. Control of unstable steady states by long delay feedback / S. Yanchuk, M. Wolfrum, P. Hövel, E. Schöll // Phys. Rev. E. 2006. Vol. 74. 026201.

12. Dahms T., Hövel P., Schöll E. Stabilization of fixed points by extended time-delayed feedback control // Phys. Rev. E. 2007. Vol. 76. 056213.

13. Huijberts H., Michiels W., Hijmeijer H. Stabi-

References:

1. Leonov G.A., Shumafov M.M. Methods of stabilization of the linear operated systems. SPb.: SpbSU Publishing House, 2005. 420 pp.

2. Static output feedback: a survey / V.L. Syrmos, C.T. Abdallah, P. Dorato, K. Grigoriadis // Automatica. 1997. Vol. 33, No. 2. P. 125-137.

3. Polyak B.T., Shcherbakov P.S. Difficult problems of the linear theory of control. Some approaches to solution // Automatic equipment and telemechanics. 2005. No. 5. P. 7-46.

4. Elsgolts L.E., Norkin S.B. Introduction to the theory of differential equations with deviating argument. M.: Nauka, 1971. 296 pp.