CC BY
37
Вычислительные технологии
Том 8, № 3, 2003
СИНТЕЗ И УПРОЩЕНИЕ ТЕХНИЧЕСКОЙ РЕАЛИЗАЦИИ ГИБРИДНОЙ НЕЛИНЕЙНО ПРЕОБРАЗОВАННОЙ СИСТЕМЫ ПРЯМОГО АДАПТИВНОГО УПРАВЛЕНИЯ
Е. Л. ЕРЕМИН, Д. Г. ШЕВКО Амурский государственный университет, Благовещенск, Россия e-mail: [email protected], [email protected]
The method of construction hybrid is considered is nonlinear of a transformed system of direct adaptive management. The possibility of simplification of technical realization of the synthesized system is shown within the framework of computing experiment.
Введение
Существуют различные методы и подходы к решению задач синтеза гибридных систем прямого адаптивного управления (ГСПАУ), прикладное значение которых для управления непрерывными объектами с помощью дискретных управляющих устройств велико в силу широкого применения микропроцессорной техники [1]. Однако проблема синтеза высокоэффективных ГСПАУ по-прежнему остается актуальной. Известно [2], что для непрерывных систем управления, функционирующих в условиях априорной неопределенности, задача построения быстродействующих алгоритмов может быть решена за счет применения метода нелинейного преобразования координат, разработанного Р.У. Брокеттом [3]. Эти результаты можно применить и для ГСПАУ. Действительно, если опираться на результаты синтеза аналоговых алгоритмов, то с помощью метода непрерывных моделей [1] можно построить соответствующие дискретные алгоритмы. Однако недостатком такого подхода к синтезу ГСПАУ является относительно малый шаг дискретизации алгоритмов управления и адаптации.
В свою очередь, при рассмотрении синтеза дискретных адаптивных регуляторов как самостоятельной задачи [4] приходится констатировать тот факт, что непосредственный перенос нелинейных преобразований Р.У. Брокетта на решение задачи синтеза высокоэффективных дискретных алгоритмов для ГСПАУ встречает ряд ограничений. Основные трудности при решении такой задачи возникают уже на начальных этапах синтеза [5]. Поэтому необходимо использовать более гибкие методы, позволяющие преодолеть эти трудности. Такой метод может представлять собой итеративную процедуру синтеза, опирающуюся на использование метода нелинейного преобразования координат [3] и аппарата теории гиперустойчивости и положительности динамических систем [4]. Итеративен Е.Л. Еремин, Д. Г. Шевко, 2003.
ный характер данной процедуры означает, что синтез осуществляется последовательно, в несколько шагов:
— во-первых, представление исходной системы в эквивалентной форме, а именно в виде линейного стационарного блока прямой цепи и нелинейного нестационарного блока обратной связи;
— во-вторых, обеспечение строгой положительности передаточной матрицы линейной части системы;
— в-третьих, нелинейное преобразование над системой;
— в-четвертых, построение дискретных алгоритмов настройки коэффициентов регулятора;
— в-пятых, проверка выполнения целевых условий и свойств адаптивности системы. Однако синтезированные таким образом ГСПАУ обычно сложны с точки зрения технической реализации, поскольку нелинейные преобразования необходимо осуществлять с помощью аналоговых вычислительных устройств (АВУ), а, как известно [1], при построении сложных вычислительных алгоритмов функциональные возможности аналоговой вычислительной техники весьма ограничены. Поэтому требуется обеспечить упрощение структуры такой ГСПАУ, доверив все вычисления управляющему цифровому вычислительному устройству (ЦВУ).
1. Постановка задачи синтеза ГСПАУ
Пусть объект управления описывается уравнением
^ = Ax(t) + bu(t) + f (t), (1)
где x(t) £ Rn — вектор состояния объекта управления (ОУ); u(t) £ R1 — управляющее воздействие; A и b — матрицы заданного размера соответственно состояния и управления; f (t) £ Rn — возмущающее воздействие, удовлетворяющее неравенству
I ||f (t)||2dt< (2)
о
Объект управления (1) функционирует в условиях априорной неопределенности, уровень которой задан в виде
A = A(£), b = b(£), f(t) = f (t), £ £ S, (3)
где £ — набор всех неизвестных параметров; S — известное множество возможных значений £.
Присоединим к объекту управления (1) самонастраивающийся дискретный регулятор, структуру которого зададим выражением
Uk = Xk rk + cT Xk, (4)
где xk £ R1, ck £ Rn — настраиваемые коэффициенты регулятора; rk £ R1 — задающее воздействие; Xk = x(tk), tk = кт — дискретный аналог времени, т = const > 0 — шаг дискретизации, к = 0,1, 2,... — номер шага.
Континуализация управляющего воздействия осуществляется с использованием экст-раполятора нулевого порядка
и(*) = ик при гк < * < гк+1. (5)
Желаемое поведение объекта управления (1) зададим с помощью эталонной модели, описываемой уравнением
dХ(t)
dt
= Am X(t) + Ьм r(t), (6)
где Х(*) £ — вектор состояния модели; Ам и Ьм — постоянные матрица и вектор соответствующих размеров; г(*) = г к при ¿к < * < ¿к+1.
В адаптивной системе управления (1), (4)-(6) требуется определить вид алгоритмов самонастройки коэффициентов Хк и Ск таким образом, чтобы при любых начальных условиях и возмущениях со свойствами (2) было достижимо выполнение следующих целевых условий:
Иш(Х(*) - х(*)) = 0, (7)
t—><x
lim Хк = X* = const, lim ck = c* = const, (8)
k—x k—x
где х* и с* — некоторые скалярная и векторная величины.
2. Синтез нелинейно преобразованной ГСПАУ
Решение поставленной задачи синтеза ГСПАУ осуществим при f (*) = 0, а затем покажем, что синтезированная система адаптивного управления будет работоспособна по отношению к заданным целям и при наличии возмущений со свойствами (2).
2.1. Получение эквивалентного описания системы
Выполняя преобразование результата вычитания уравнения (1) из уравнения (6) с учетом соотношений (4), (5) и обозначения
е(г) = х(*) - х(*), (9)
можно получить следующее эквивалентное математическое описание исследуемой системы:
^ = Ам е(*) + ЬМ*), (10)
dt
V (*) = дТе(*), (11)
V-к = (х* - Хк)гк + (с* - Ск)ТХк, М*) = ¡лк при *к < * < *к+1, (12)
где V(*) £ К1 — обобщенный выход эквивалентной системы; д — постоянный вектор, элементы которого подлежат выбору. Уравнения (10), (11) определяют линейную стационарную часть (ЛСЧ) системы, а уравнения (12) — нелинейную нестационарную часть (ННЧ).
2.2. Определение условий положительности ЛСЧ системы
Проведение синтеза на этой стадии разработки ГСПАУ состоит в решении проблемы положительности относительно линейной части исходной системы управления с эквивалентным математическим описанием вида (10)-(12). Стандартный подход к решению такой задачи состоит в обеспечении свойств вещественности и положительности передаточной функции линейной стационарной части системы:
ТЯП \ Т( п Л дТ(вЕп - Ам) + Ь
W(з)= д - АМ) Ь =-——, (13)
- Ам)
где 5 = ^'и, ] = \/ — 1; и > 0; — единичная матрица порядка п; (зЕга — Ам)+ — матрица, присоединенная к матрице (зЕга — Ам). Известно [2, 6], что для получения W(в) с указанными свойствами необходимо и достаточно вектор д выбрать таким образом, чтобы в условиях априорной неопределенности (3) полином дт(зЕга — Ам)+Ь был гурвицевым степени п — 1 с положительными коэффициентами.
2.3. Этап нелинейных преобразований
Рассмотрим нелинейное преобразование координат исходной системы (10)-(12). Воспользуемся приемом, рассмотренным в работе [7], и осуществим переход от переменных е(£) € Мга к переменным е[2](£) € с помощью выражения
Н (*) = е(^)ет (*), (14)
где Н(¿) — симметричная матрица, элементами которой являются попарные произведения элементов вектора е(£).
Дифференцируя выражение (14) и учитывая уравнение (10), получим
^Н (¿)
<И
АмН(¿) + Н(¿)АМ + ^(¿)(Ьет (¿) + е(£)Ьт).
(15)
Выполняя преобразование выражения (15) по аналогии с [7], уравнение состояния ННЧ нелинейно преобразованной системы запишем как
(16)
а уравнение состояния ЛСЧ нелинейно преобразованной системы запишем следующим образом:
&[2](г)
(Ам )[2]е[2] (*) + Ь[2]^12|(^),
[2] I
(17)
где е[2](£) — вектор размера (^П * 1), = 0.5п(п + 1); ^[2](£) — вектор размера (п * 1); (Ам)[2] и Ь[2] — матрицы размера соответственно (N2*^2) и (N2*п). Для выхода нелинейно преобразованной системы можно записать уравнение
V[2](^) = е(^)ет (*)д = д'[2]е[2](^),
= дт еИ(
(18)
где д|2] — матрица размера (п*0.5п(п +1)), элементы которой удовлетворяют соотношению
д
т д[т2]
т0
т
0д 00
0 0
дт
Таким образом, нелинейно преобразованная система описывается уравнениями (12), (16) - (18).
Теперь введем в рассмотрение функцию z(t), которую определим выражением
z(t) = gT e(t)eT (t)e(t), (20)
тогда с учетом соотношений (16), (18) очевидно следующее равенство:
(^W V(t) = Mt)z (t). (21)
Если аналогично [2] нелинейное преобразование выполнить повторно, т.е. в качестве "новой"исходной системы рассмотреть преобразованную систему (12), (16)-(18), а затем выполнить переход к "новой"преобразованной системе, но уже относительно вектора e[2] (t), то функцию z(t) можно будет представить в виде
z(t) = gTe(t)(eT (t)e(t))2. (22)
По аналогии с [2] при выполнении подобного преобразования неоднократно очевидно, что для исходной системы (10) - (12) функцию z(t) можно будет записать следующим образом:
z(t) = gT e(t)(eT (t)e(t))2q, q = 0,1, 2,... (23)
2.4. Этап синтеза дискретной части ГСПАУ
Согласно [4], для ННЧ исследуемой системы необходимо показать справедливость неравенства
ki
n(0, fci) = — ^ vk > -Y0 = const, vfci > 0, (24)
к=О
где Vk = v (tk).
Аналогично [2] для решения проблемы положительности ННЧ исходной системы (10) -(12) воспользуемся результатами нелинейного преобразования и рассмотрим вместо неравенства (24) неравенство, записанное относительно нелинейно преобразованной системы:
ki
n(0,k1) = — ^^ zk > — y2 = const, Vk1 > 0, (25)
k=О
где Zk = z(tk).
Используя первое из уравнений (12), получим
ki
2
У^ ((Xk — X*)rk + (ck — c*)TXk) Zk > —yO. (26)
k=0
Теперь положим
Xk = Xk-1 + Mzk), (27)
Ck = Ck-1 + 0(zk) (28)
или
k
Xk = ^(zi) + X-1, (29)
1=0
k
Ck = ф(гг) + С—i. (30)
1=0
Тогда получим неравенство
ki / k \ki / k \т
П(0,к0 = ^ Zk(J^ ) + X—i - X* 1 rk + zk(5^ ф(^) + c—i - c* ) xk - -Yo, (31) k=0 \ 1=0 ) k=0 V 1=0 )
которое будет выполняться, если оба члена левой части удовлетворяют неравенству того же типа. Чтобы найти решения для ^ и ф, приводящие к выполнению таких неравенств, применим следующее соотношение [4]:
ki ( k \ ( ki \ 2 ki
]>>k ПС У + C = 0.5 £ yk + C + 0.5^ y2 - 0.5C2 — —0.5C2, C = const. (32) k=0 \1=0 / \k=0 / k=0
Используя (32), получим решения для ^ и ф:
^(zk) = h1zkrk, h1 = const > 0, (33)
ф^) = H2ZkXk, Я2 = diag{h2i}, h2i = const > 0, i = 1,n, (34)
т.е. искомые алгоритмы настройки коэффициентов адаптивного регулятора примут вид
Xk = Xk—1 + h1Zk rk, (35)
Ck = Ck—1 + #2Zk Xk. (36)
2.5. Достижимость поставленных целей управления и адаптации
В силу решения в системе управления (10)-(12) проблем положительности ЛСЧ и ННЧ, причем для любых начальных условий и при наличии априорной неопределенности (3) эту систему, согласно критерию гиперустойчивости [4], следует считать асимптотически гиперустойчивой.
Тогда благодаря выполнению предельного соотношения
lim e(t) = 0 (37)
t^x
цель управления вида (7) также имеет место.
При этом с учетом явного вида алгоритмов самонастройки коэффициентов регулятора Xk и Cfc очевидно, что будут выполнены предельные соотношения
lim xk = const, lim ck = const, (38)
fc^x fc^x
совпадающие с требованиями соответствующих целевых условий (8).
2.6. Работоспособность системы при f (£) = 0
Поскольку эквивалентное математическое описание этой системы включает в себя основной контур
de(t) dt
Am e(t) + b^(t) + f (t), v (t) = gTe(t)
^k = (X* — Xk)rk + (c* — Ck)TXk, Mt) = ^k при tk < t < tk+1
(39)
(40)
и контур адаптации с алгоритмами (35), (36), при затухающих возмущениях (см. свойство (2)) синтезированная система будет сохранять работоспособность, которая, как и ранее, понимается в смысле достижения требований вида (7), (8). Действительно, из выполнения неравенства (2) следует справедливость предельного тождества
lim f (t) = 0
t—
(41)
с учетом которого становятся очевидными как асимптотическая устойчивость эквивалентной системы, так и выполнение целевых условий (7), (8).
Если теперь от математического описания ГСПАУ, представленного в эквивалентном виде, вернуться к исходному, т.е. записать исходные уравнения объекта управления, эталонной модели и адаптивного регулятора, то синтезированная ГСПАУ будет иметь следующий вид:
dt
dx(t) dt
= Ax(t) + bu(t) + f (t), = Am X(t) + Ьм r(t),
e(t) = x(t) — x(t), z(t) = e(t)(eT (t)e(t))2, q = 0,1, 2, Uk = Xk rk + C^Xk, Xk = x(tk), Xk = Xk-1 + hzkrk, h = const > 0, zk = z(tk), Ck = Ck-1 + H2zkXk, Я2 = diag{h2i}, Л2 = const > 0, i = U(t) = Uk, r(t) = rk, tk < t < tk+1.
1, n,
(42)
(43)
(44)
(45)
(46)
(47)
(48)
Таким образом, аналитическая часть проектирования ГСПАУ завершена. Однако ряд числовых значений, в частности значения положительных констант, входящих в алгоритмы контура адаптации, требуют своего задания, что чаще всего осуществляется на этапе имитационного моделирования ГСПАУ вида (42)-(48).
2.7. Пример
Для случая п = 3 рассматривается объект управления (42), где матрица А и вектор Ь имеют известную структуру:
A
0 0
1 0
0 1
Ь
a1 —23 —9
0 0
Ь1
(49)
причем уровень априорной неопределенности задан следующим образом:
0.9 < «1 < 1.1, 0.9 < Ь1 < 1.1,
(50)
а параметры эталонной модели определены в виде
0 1 0
А
м
0 0 1
-15 -23 -9
Ьм
0 0 1
(51)
Требуется в системе (42)-(51) обеспечить выполнение целевых условий вида (7), (8). Для выполнения условий вещественности и положительности линейной части системы (42)-(51), согласно [2, 6], требуется выбрать элементы вектора д таким образом, чтобы полином второго порядка
дт(вЕз — Ам)+Ь = (дзв2 + д2в + д^Ьь дт = ( д2 дз )
(52)
был гурвицевым с положительными коэффициентами.
Результаты имитационного моделирования ГСПАУ (42)-(51), полученные при нулевых возмущениях и начальных условиях и следующих исходных данных:
«1 = 1, Ь1 = 1, д = 0, гк = 0.1 8т(кт ) + 1, т = 8.8,
= 18.47, Н2 = diag{195160; 40000; 8000000},
представлены на рис. 1, причем вектор д был выбран в виде
т
8 6 1 .
(53)
(54)
(55)
Рис. 1. Динамические процессы в системе (42)-(48).
3. Упрощение технической реализации ГСПАУ
Очевидно, что с точки зрения технической реализации синтезированную ГСПАУ можно представить в виде структуры, изображенной на рис. 2, где помимо уже принятых использованы следующие сокращения и условные обозначения: т[к] = ; и[к] = ; х[к] = хк; г[к] = ; АЦП — аналого-цифровой преобразователь; ЦАП — цифроаналоговый преобразователь.
Следует отметить, что ГСПАУ (42) - (48) имеет сложную техническую реализацию. Как уже отмечалось, эта сложность связана с тем, что нелинейные преобразования необходимо осуществлять на аналоговой технике, функциональные возможности которой весьма ограничены [1]. Техническую реализацию ГСПАУ (42)-(48) можно упростить, исключив из контура системы АВУ, доверив при этом все вычисления ЦВУ. Тогда функционирование ЦВУ помимо уравнений (45) - (47) будет описываться дискретными аналогами уравнений (43), (44). Таким образом, очевидно, что ГСПАУ (42)-(48) может иметь
Рис. 2. Структурная схема ГСПАУ (42)-(48).
Рис. 3. Структурная схема ГСПАУ (56)-(62).
Рис. 4. Рассогласования в ГСПАУ (42) - (48) и (56) - (62) по выходу и управлению.
эквивалентную техническую реализацию в виде ГСПАУ, основанной на следующей мате-
матической модели:
dXt) = Ax(t) + bu(t) + f (t), (56) dt
xfc+i = PM xfc + dM rfc, (57)
efc = Xfc - xfc, xfc = x(tfc), zfc = gT efc (e^ efc )2?, q = 0,1, 2,..., (58)
Ufc = Xfc rfc + cT xfc, (59)
Xfc = Xfc-i + hizfcrfc, hi = const > 0, (60)
cfc = cfc-i + H2zfcxfc, Я2 = diag{h2i}, = const > 0, i =1,n, (61)
u(t) = ufc при tfc < t < tfc+i, (62)
где PM = exp(AMt); dM = AMi(PM — )bM. Структура ГСПАУ (56)-(62) представлена на рис. 3.
Проверка эквивалентности ГСПАУ (42) - (48) и ГСПАУ (56) - (62) опирается на результаты вычислительного эксперимента. Вычисления проводились с использованием средства визуального моделирования SIMULINK математического пакета MATLAB. Результаты моделирования, полученные при исходных данных (49), (51), (53)-(55) и нулевых возмущениях и начальных условиях, приведены на рис. 4, где Ai — рассогласование значений элементов вектора состояния объекта управления систем (42)-(48) и (56)-(62), а A2 — рассогласование значений управляющих воздействий.
Заключение
В соответствии с требованиями цели управления, поставленной в условиях априорной неопределенности, решена задача синтеза гибридной системы управления. В рамках вычислительного эксперимента показано, что отличающиеся по структуре ГСПАУ (42)-(48) и ГСПАУ (56) - (62) с точки зрения технической реализации могут считаться эквивалентными. Анализируя полученные результаты имитационного моделирования, можно сделать вывод о целесообразности первоначального синтеза системы управления в виде (42) - (48) с последующим ее упрощением до вида (56)-(62). Иными словами, предложенная новая итеративная процедура синтеза гибридных систем управления может быть дополнена шестым шагом, цель которого — упрощение структуры ГСПАУ путем дискретизации непрерывной части ее контура адаптации.
Список литературы
[1] Цыкунов А. М. Адаптивное управление объектами с последействием. М.: Наука, 1984.
[2] Еремин Е. Л., Цыкунов А. М. Синтез адаптивных систем управления на основе критерия гиперустойчивости. Бишкек: Илим, 1992.
[3] Брокетт Р. У. Алгебры Ли и группы Ли в теории управления // Мат. методы в теории систем / Под ред. Ю. Н. Журавлева. М.: Мир, 1979. С. 174-220.
[4] Landau I.D. Adaptive Control Systems: the Model Reference Approach. N.Y.: Marsel Dekker, 1979.
[5] Баркин А. И., ЗЕЛЕНЦОВСКИй А. Л., Плкшин П. В. Абсолютная устойчивость детерминированных и стохастических систем управления. М.: МАИ, 1992.
[6] ФРАДКОВ А. Л. Адаптивное управление в сложных системах: беспоисковые методы. М.: Наука, 1990.
[7] Баркин А. И. Оценки качества нелинейных систем регулирования. М.: Наука, 1982.
Поступила в редакцию 25 ноября 2002 г.
