Арапина-Арапова Елена Сергеевна - ФГБОУ ВПО «Таганрогский государственный педагогический институт им. А.П. Чехова»; e-mail: arapina@mail.ru; 347926, г. Таганрог, ул. Инициативная, 48; тел.: +79185219138; кафедра информатики; к.ф.-м.н.; доцент.
Arapina-Arapova Elena Sergeevna - FGBOU of VPO «Taganrog State Pedagogical Institute of a name A.P. Chekhov»; e-mail: arapina@mail.ru; 48, Initsiativnaya street, Taganrog, 347926, Russia; phone: +79185219138; the department of informatics; cand. of phis.-math. sc.; associate professor.
УДК 681.12
К.В. Бесклубова
ОБЛАСТЬ ПРИТЯЖЕНИЯ ПОЛОЖЕНИЯ РАВНОВЕСИЯ МНОГОМЕРНОГО ОБЪЕКТА С ОГРАНИЧЕННЫМ УПРАВЛЕНИЕМ
Исследуется область притяжения положения равновесия многомерного динамического объекта при ограниченном управлении, формируемом многомерным устройством управления. Построение области производится с применением метода квадратичных функций Ляпунова и инвариантных эллипсоидов. Получены аналитические выражения, определяющие область начальных значений переменных состояния объекта, при которых обеспечивается гарантированная линейность системы управления. Проведено моделирование системы с ограниченным управлением при различных начальных условиях. Результаты моделирования подтверждают справедливость полученных выражений.
Многомерная система управления; ограничение на управление; функция Ляпунова; оценка области устойчивости; переменные состояния.
K.V. Besklubova
STABLE STATE ATTRACTION AREA OF THE MULTIVARIABLE PLANT WITH LIMITED CONTROL
The paper considers the attraction area of stable state of multivariable dynamic plant with limited control formed with the aid of multivariable control unit. The attraction area is formed with application of square Lyapunov's function method and invariant ellipsoids. Analytical expressions, defining area of plant’s initial conditions providing the guaranteed system’s linearity, are received. Simulation of the system with limited control is carried out at the various initial conditions which results prove the specified expressions.
Multivariable control system; limited control; Lyapunov's function; an estimation of linearity area; state variables.
Введение. Управляющие воздействия реальных объектов, как правило, ограничены по величине. В том случае, когда значение управления, формируемого устройством управления, превышает допустимую величину, система может потерять линейность и, как следствие, устойчивость и работоспособность. В связи с этим в данной работе рассматривается задача построения областей притяжения положения равновесия динамической системы с ограничением на управление. Решение задачи получено методом функций Ляпунова с применением инвариантных эллипсоидов.
Постановка задачи. Рассмотрим линейный многомерный объект управления (МОУ), заданный уравнениями
X = Ax + Bu , y = CT x + Du , (1)
где x - n -вектор состояния МОУ; A, B, C, D - матрицы коэффициентов; u, y -q -векторы управлений и выходов МОУ.
Для МОУ (1) методом динамической декомпозиции [5] многомерного объекта синтезировано управление по выходам и воздействиям. При отсутствии задающих воздействий его можно представить в виде функции от переменных состояния замкнутой системы:
f KLш, |K.ю|< и , . —
и; = \ II i = 1, q (2)
[“max,islgn(Ki|KiЮ|> umax,i,
где u , i = 1, q - управление i-того канала замкнутой многомерной системы автоматического управления (МСАУ), Ki - строки матрицы K параметров управления, umax. - допустимые значение управления u . В выражении (2) ю - пзс -вектор переменных состояния замкнутой МСАУ, причем ю = [x z]T , где z e Rr - вектор переменных состояния многомерного устройства управления (МУУ).
Уравнения замкнутой системы с управлением (2) имеют вид
ю = АзСю, y = СЗтСю, (3)
где АЗС , С ЗС - матрицы коэффициентов МСАУ.
Ниже рассматривается приложение метода функций Ляпунова с применением инвариантных эллипсоидов к решению задачи определения диапазонов изменения переменных состояния, в которых сохраняется линейность системы (3).
Решение задачи. Поскольку система (3) является устойчивой, то для неё
можно найти функцию Ляпунова в виде квадратичной формы V(ю) = ютРю.
В этом выражении Р - симметричная, положительно-определенная матрица, которую можно найти путем решения уравнения Ляпунова
АзТС Р + РАзс =- E , (4)
где E - единичная матрица. Производная по времени функции V(ю) = ютРю вдоль траекторий системы определяется выражением V (ю) = -ют Ею . Так как E > 0, то V(ю) является отрицательно-определенной функцией. Поэтому при ограниченной норме ||ю0|| вектора начальных значений ю0 решение ra(t) = ra(t, ю0) систем^1 (3) также является ограниченным при всех t > 0 и стремится к нулю с ростом времени.
В последнее время широкое распространение получила теория инвариантных эллипсоидов [2-4]. Ее методы с применением техники матричных линейных неравенств используются при описании множеств достижимости и притяжения линейных систем, а также для формирования управления в форме линейной статической обратной связи по состоянию [3]. Кроме того, методы теории инвариантных эллипсоидов служат основанием для разработки новых подходов к решению проблемы подавления ограниченных внешних возмущений и, в частности, позволяют проводить синтез оптимальных регуляторов и осуществлять оценку показателей качества замкнутых динамических систем [4].
В данной работе метод инвариантных эллипсоидов используется для оценки допустимого диапазона изменения переменных состоянии системы (3), при котором выполняется неравенство |K. ю| < umaxi.
Приведем определение: эллипсоид
ер(Р) = {юе Rn: ютРю< р2} , Р > 0 (5)
называется инвариантным притягивающим для (3), если для любого ю0 e ер (Р) при всех t > 0 будет выполнено условие
ю(г,ю0) е ер (Р). (6)
Здесь параметр р2 должен быть связан с допустимым значением управления
итах .
Фактически эллипсоид ер представляет собой область в пространстве состояний системы (3), ограниченную поверхностью ют Рю = р2. Если на управление системы не налагается каких-либо условий, т.е. система устойчива при любом векторе начальных значений ю0, все траектории системы с ростом времени стремятся к эллипсоиду ер и попадают в него при некотором г > гр > 0. Однако если при юе 6р(Р) функция V(ю) < р2, то при юг ер(Р) V(ю) = ютРю> р2. В этом и заключается смысл ер (Р).
Согласно лемме 1 из [2], если начальные значения вектора переменных состояния системы (3) удовлетворяют равенству ютРю = р2, то при всех г > 0 выполняется условие
ют (г, ю0)Рю(г, ю0) < р2. (7)
Нужно отметить, что из леммы 1 не следует монотонное затухание ю(г) или
его нормы ||<ю»||. В ней утверждается, по сути, лишь монотонность убывания
функции ют (г)Р ю(г) с ростом г. Причем именно это и требуется для того, чтобы ю(г) е ер (Р), так как ер (Р) определяется условием ютРю<р2. Хотя словами и говорится: если ю0 е ер (Р), то и ю(г, ю0) е ер (Р), на самом деле ограничивается не ю(г, ю0), а функция V (ю) = ют Рю.
При этом, по лемме 2 из [2], обеспечить выполнение неравенства |и| < итх можно при выполнении условий (7) и (8), при У = Кт Р 1
> 0. (8)
u2
л/ max 77
Y У£
Таким образом, задаваясь некоторым значением ишах, по условию (8), можно
ч_/ 2 V-»
наити параметр р , который и определяет размеры эллипсоида єр, являющегося границей области начальных значений х0, при которых обеспечивается неравенство Ц | < м2ахi и линейность системы (3) .
Пример. В рассматриваемой МСАУ (3) q = 3, dim х = 6, dim ю = 30 . Для нее найдена матрица Р из уравнения Ляпунова (4). Соотношения (7), (8) метода инвариантных эллипсоидов позволяют в явном виде построить область притяжения, в
которой выполняется условие |uj < umaxi, umaxi = 10 , i = 1, q . Эта область для объекта (1) с управлением (2) показана на рис. 1. Если начальные значения вектора переменных состояния попадают в область єр, т.е. ю0 є ер (P), то соответствующий график x(t) располагается в этой области, как показано на рис. 1.
Шз0
4
2
®б0
-2
-4
-5 -5
Рис. 1. Проекции траекторий МОУпри Шо є Єр (Р), ш^ РШо < р
-5 -5
Рис. 2. Проекции траекторий МОУ при Шо є Єр (Р), Р1 Хо = р2
б
Рис. 3. Управление МСАУ при Шо є Єр (Р) (а - ш^ РШо < р2, б - л^Рі Ло = р2)
5
5
5
5
10
о
5
5
а
Рис. 4. Проекции траекторий МОУ при Шо £ Єр (Р)
Рис. 5. УправлениеМСАУпри ш0 g ер(P)
На рис. 2 приведены проекции траекторий МОУ для случая, когда начальные значения х0 его переменных состояния расположены на границе области е р ( P).
Как видно, при x^Pjх0 = р2, где р - подматрица P , траектории МОУ также ограничены, однако могут покидать область ер (P). По рис. 3, при ш0 е ер (P) управление не достигает своего максимально допустимого значения итаА!. = 10 и, следовательно, замкнутая система остается линейной.
При дальнейшем увеличении значений х0, когда начальная точка находится за границей области ер (P), т.е. ш0 g ер (P), управляющие воздействия достигают
своего максимально допустимого значения, а система теряет линейность. Это подтверждается изменением переменных состояния объекта, изображенных на рис. 4, и графиками изменения управляющих воздействий на рис. 5.
Заключение. Рассмотренный метод позволяет построить область притяжения положения равновесия многомерного динамического объекта или системы при ограниченном управлении, где гарантируется линейность замкнутой системы.
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1. Гайдук А.Р. Оценки переменных состояния линейных систем. Наука и образование на рубеже тысячелетий: сборник научно-исследовательских работ. Вып. 2. - М.: Учлитвуз, 2011. - С. 23-27.
2. Демидович Б.П. Лекции по математической теории устойчивости. - М.: Наука, 1967.
3. Поляк Б.Т., Щербаков П.С. Множества достижимости и притяжения линейных систем с ограниченным управлением: описание с помощью инвариантных эллипсоидов. Сб. «Стохастическая оптимизация в информатике» / Под ред. О.Н. Граничина. Вып. 4. - СПб.: СПб ГУ, 2008. - С. 3-23.
4. Хлебников М.В. Время установления в линейной динамической системе с ограниченными внешними возмущениями // Автоматика и телемеханика. - 2012. - № 6. - С. 3-17.
5. Гайдук А.Р., Плаксиенко Е.А. Синтез динамических систем по требуемым показателям качества // Мехатроника, автоматизация, управление. - 2008. - № 4. - С. 7-12.
6. Гайдук А.Р. Теория и методы аналитического синтеза систем автоматического управления (Полиномиальный подход). - М.: Физматлит, 2011.
Статью рекомендовал к опубликованию д.т.н., профессор Р. А. Нейдорф.
Бесклубова Ксения Валериевна - Федеральное государственное автономное образовательное учреждение высшего профессионального образования «Южный федеральный университет»; e-mail: kbesklubova@mail.ru; 347904, г. Таганрог, ул. Яблочкина, 8/1, кв. 33; тел.: 88634387349; магистрант.
Besklubova Ksenia Valeryevna - Federal State-Owned Autonomy Educational Establishment of Higher Vocational Education “Southern Federal University”; e-mail: kbesklubova@mail.ru; 8/1, Yablochkina street, app. 33, Taganrog; 347904, Russia; phone: +78634387349; the department of automatic control systems; magister.