СИНТЕЗ МНОГОФУНКЦИОНАЛЬНЫХ УПРАВЛЯЮЩИХ УСТРОЙСТВ
Д.Ю. Муромцев, В.В. Папин
Кафедра «Конструирование радиоэлектронных и микропроцессорных систем», ФГБОУ ВПО «ТГТУ»; crems@crems.jesby.tstu.ru
Ключевые слова и фразы: задачи оптимального управления; задачи синтеза; квадратичный оптимальный регулятор; многофункциональные управляющие устройства; энергосберегающий оптимальный регулятор.
Аннотация: Рассмотрены вопросы синтеза многофункциональных управляющих устройств систем энергосберегающего управления динамическими объектами. Приведен пример сравнения двух оптимальных регуляторов: квадратичного и энергосберегающего.
В общем случае под синтезом понимается процесс создания новой системы управления с оптимизацией ее свойств и показателей. Можно выделить две группы задач синтеза систем энергосберегающего управления (СЭУ). К первой группе относятся задачи синтеза оптимального управления (оу), результатом решения которых являются соотношения для оперативного расчета энергосберегающих управляющих воздействий. Эти соотношения составляют основу алгоритмического обеспечения СЭУ. Ко второй группе относятся задачи синтеза устройств энергосберегающего управления, в том числе выбор стратегии реализации ОУ, структуры СЭУ и т.д.
Алгоритм синтеза оптимального управления, реализуемый управляющим устройством применительно к конкретному объекту, то есть заданной модели динамики, во многом определяется видом минимизируемого функционала и выбранной стратегией.
Определение 1. Управляющее устройство, алгоритмическое обеспечение которого решает задачи синтеза ОУ при одном виде минимизируемого функционала с использованием одной стратегии реализации ОУ, будем называть однофункциональным.
Алгоритмическое обеспечение однофункциональных устройств обычно представляет собой часть результатов полного анализа для конкретной модели задачи оптимального управления (ЗОУ), например, для моделей < М, Э, Пр, О >, <М, Т, Пр, О >, <М, Э, Пз, О > и т.д.
Определение 2. Управляющее устройство, алгоритмическое обеспечение которого позволяет решать задачи при различных функционалах или (и) стратегиях реализации ОУ, будем называть многофункциональным.
Многофункциональные управляющие устройства (МФУУ) являются разновидностью систем с переменной структурой [1]. В МФУУ в зависимости от ситуации может решаться задача синтеза ОУ, например при функционалах затраты
энергии или квадратичном использоваться как позиционная, так и программная стратегии.
Несомненным преимуществом МФУУ является то, что они позволяют объединить достоинства различных однофункциональных управляющих устройств. Рассмотрим это на примере сравнения двух оптимальных регуляторов (ОР): квадратичного (КВОР) и энергосберегающего (ЭОР).
При скалярном управлении в КВОР минимизируется функционал
/кв = xт (*к) N x(tк) + | (xт (t)Gx(t) + qu 2 (/)^/, (1)
*0
здесь х - и-мерный вектор фазовых координат; и - управляющее воздействие; *о, *к - начало и конец временного интервала управления; N, G, q - (и х и)-матрицы и весовой коэффициент, учитывающие важность составляющих квадратичного функционала.
Достоинством ОР с функционалом (1) является быстрое уменьшение ошибки рассогласования.
Если модель динамики объекта управления имеет вид
х = Ах(*) + Ви(*), *е[%*к] , х = г - 2к, (2)
то при функционале I кв обычно и ) ищется в виде
и*({)= К(*)х(*), (3)
где К(*) - (и х1) - матрица, определяемая решением уравнения Риккати и зависящая от А, В, G, q [2].
Реально существующее ограничение на управление
V*е[%*к]: и (* )е[ин,ив], (4)
и другие ограничения при этом не учитываются.
В ЭОР минимизируется функционал
1э =| и 2(* ^ (5)
*о
и управляющие воздействия определяются синтезирующей функцией
и*(*) = S(х(*), *к - *; ^о), (6)
где
Rq =(A, B, uH, uB, x°, хк, t0, tK).
На первый взгляд функционал (5) можно рассматривать как частный случай функционала (1) и поэтому определение синтезирующей функции (6) является простой, хорошо изученной задачей. Однако это не так, ЗОУ с функционалом (5) осложняется тем, что функция (6) может иметь несколько видов. Для определения конкретного вида S (х, т, Rо), зависящего от Rо и обеспечивающего достижение
цели управления х(*к ) = хк, требуется решение ряда задач: определение всех возможных видов функции (6) для исследуемого объекта управления, нахождение границ областей существования видов синтезирующих функций и областей, в которых выполняются ограничения на х и и, разработка алгоритмов поиска вида функции (6) и расчета ее параметров для текущего значения массива R. Неправильное определение вида синтезирующей функции может приводить к срыву процесса регулирования при выполнении задачи с функционалом 1э. Если компоненты
К(І) в КВОР при выбранных G и q зависят только от параметров объекта (матрицы А, В), то вид и параметры синтезирующей функции (6) зависят еще от ин, ив, Х°, Хк, (to - tк).
Достоинством алгоритма (3), по сравнению с (6), является быстрая обработка больших рассогласований (значений х) в задачах стабилизации и слежения [2]. На рис. 1, а приведена качественная картина изменения х^) для п = 1 при использовании КВОР и ЭОР. Вместе с тем использование аппарата синтезирующих переменных [3-5] позволяет обеспечить «приемлемый» профиль траектории х() = (х(ґ)), t є^0, tк ] путем наложения ограничений на некоторые промежуточные значения х(і) (рис. 1, б). Такой ЭОР будем называть модифицированным, он реализует выполнение ограничений к профилю х(-).
Ограничение к профилю траектории х() можно сформулировать, например,
следующим образом: к моменту времени tп = -2^к - tо) рассогласование |х0| (в предположении хк = 0) должно уменьшиться на 2/3, то есть хп = х(/п) = 1 х0.
В общем случае ограничение к профилю траектории х() в задаче с функционалом (5) записывается в виде
х()є X () (7)
или
х(Ц ) = gIАхт, і =1,2, ..., p, (7а)
где X (•) - область допустимых значений х( ); ti, gi - і-й фиксированный момент
времени и коэффициент, задающий х^і), соответственно; Дхт - максимальное рассогласование, обычно |хо - хк|; р - число фиксированных точек требуемой траектории X (• ). В формуле (7 а) предполагается п = 1, но ее не трудно обобщить на случай п > 1.
Наряду с ограничением (7а) в случае необходимости ЭОР должен обеспечивать ограничения:
а) на максимальные значения компонентов вектора х, например, ограничение на скорость (компонента х2) изменения температуры при нагреве тел, то есть
Vt Є[^ ^]: |х2(0| < х2доп, (8)
ускорение движущегося объекта
Vtє[%^]: \х2(0\ < х2доп; (8а)
Рис. 1 Траектории изменения х(0 в задаче стабилизации для КВОР и ЭОР (а) и для модифицированного ЭОР (б)
б) скорость изменения управляющего воздействия
V t е[% tK ]: |U(t)| < г&д0п; (8б)
в) выделенный лимит энергии
|и 2(/) dt < /доп; (8в)
{0
г) значение управления в конечный момент времени
u(tк ) = ик, (8г)
где х2доп, -&2доп, идоп, /доп - допустимые значения Х2, Х2, и и затрат энергии (функционала).
Ограничения на управление и х(/к) в ЗОУ <М,Э, S, О > будем называть основными, а ограничения (7) - (8г) - дополнительными. С учетом дополнительных ограничений размерность массива увеличивается, то есть
Я0 =(А, В, ин, ив, х , х , t0, tк; /п, 8; х2доп, х2доп, идоп, /доп). (9)
В целях повышения надежности функционирования потребуем, чтобы ЭОР
при нарушении обратной связи по х вырабатывал ОУ в соответствии с программной стратегией, а при отсутствии решения ЗОУ находил управление, обеспечивающее наименьшее снижение эффективности системы оптимального управления, например, путем увеличения времени /к, ослабления отдельных дополнительных ограничений на х и и. Придание ЭОР перечисленных функций позволяет отнести его к так называемым интеллектуальным контроллерам [6-8].
В задачах оптимального управления < М,Э, S, О > решаемых МФУУ, можно выделить два вида ограничений по степени их значимости для процесса регулирования.
Определение 3. Ограничения в ЗОУ будем называть сильными (общего характера), если нарушить их физически невозможно или они представляют собой основную цель управления, и частными - если в исключительных случаях для достижения главной цели управления допускается их невыполнение.
Выделение этих двух видов ограничений зависит от особенностей конкретного объекта и условий его реальной работы. Ограничение на величину управляющего воздействия в каждый момент времени (см. (4)) всегда сильное. Сильные ограничения учитываются при построении областей £с и £j. К частным
ограничениям в различных условиях могут быть отнесены следующие: на значения фазовых координат в интервале времени [% /к ], например, характеризующих ускорение или скорость движения объекта (см. (8), (8а)); скорость изменения управления (см. (8б)); лимит используемой энергии (см. (8в); величину временного интервала управления [/о, tк ].
Определение 4. Часть области £с, в которой выполняются частные ограничения конкретной ЗОУ, будем называть допустимой и обозначать £ доп.
Для ряда объектов, например, топливных печей, двигателей внутреннего сгорания и др., накладывается ограничение на скорость изменения управляющего воздействия, то есть
du
dt
А
< идоп = (10)
В качестве примера рассмотрим построение области £18 с £1 для объекта двойного интегрирования. В области £ и*(/)=и (/ )=dо + d1t, следовательно [9]
= KI =
(11)
Из (11) видно, что границы области £18 в пределах L1 определяются отрезками прямых (рис. 2, а):
г0 = г0 2g tK —10
L2 = L1 — ~
3 uB — uH
r0 = r0 , 2g tK —10
L2 = L1 + -----------------------
3 uB — uH
(12)
Для объекта управления с матрицей
0 1
0 a2
уравнения, аналогичные (12), имеют вид:
L2 =
L2 =
Л2ао 1
e 2 — 1
Л2ао 1
e 2 —1
2a2
L1 (1 — «2) — (1 — a2 )e—2a2 )(1 — e2a2)
2a|
(13)
2a0
L1 +~~2 (1 — a2 — (1 + a2 ) e 2«2 X1 — e2«2 ^ a2 = a2 “
tK — “0
2a22
2
Рассмотрим определение границ области £1У, в которой выполняется огра-
1у ,
ничение на ускорение фазовой координаты zl, то есть
<
7 доп 2
(14)
на примере объекта двойного интегрирования. Так как для ОУ первого вида
Z2 = Ьи (/)=Ь^о + dlt), (15)
то границы определяются из условий
|b(d0 + d1t )|<у (16а)
|b(d0 + d1“0^ <Y •
(16б)
Рис. 2. Области L1g (а) и L1y (б)
0
и
В предположении /о = 0 из условия (16а) получаем:
^ = /к — 3 L + /к 2У — Ь(ив + ин) .
1,5/к — 3 1,5/к — 3 Ь(ив — ин)
к к * ч (17)
г"_ /к — 3 т /к 2У — Ь(ив + ин)
L2 =----------^1---------------- -------- ---,
1,5/к — 3 1,5/к — 3 Ь(ив — ин)
а из условия (16б):
L, = 2/к + 3 L + /к 2У — ь(ив + ин) .
1,5/к +3 1,5/к +3 Ь(ив — ин) .1вЛ
(18)
г" = 2/к + 3 т /к 2У — Ь(ив + ин)
L2 =---------------------------- -------- ---.
1,5/к + 3 1,5/к + 3 Ь(ив — ин)
На основе (17), (18) определяются границы области ^, представляющие собой отрезки прямых (рис. 2, б).
В отличие от рассмотренных областей £8, £у аналитическое построение
области £ л с ограниченным лимитом энергии (см. (8в)) достаточно трудоемко. Например, для первого вида функции ОУ в задаче двойного интегратора при /о = о и ив = |ин| граница области £1л имеет нелинейный характер и определяется из соотношения
2 2 4~ /
Ф1(/к)А + ф1,2(/к)L1L2 + Ф2(/к^2 <
(uB uH)
где
Ф1(“к) = i9 — 9“к + 19“к — 6tKФ2(“к) = i9 — 9“к + 14,25“2 — 4,5tK)
А
ф1,2(“к) = (— 18 + 18“к — 33“к2 + 10,5“к )-
Если в ЗОУ присутствуют все три рассмотренных вида ограничений, то допустимая область £доп определяется как пересечение областей £8, £у и £л, то есть
£доп = £8 ^ £у ^ £л.
Утверждение 1. Если Lо е£доп, то можно выбрать такое время квантования, что при отсутствии возмущений в процессе регулирования для всех моментов времени на интервале [/о, /к] ограничения (1о), (14) выполняются, и выделенный лимит энергии не будет перерасходован.
Это следует из того, что при времени квантования 5^ о и*(/) и z*(t)
V / е[% /к ] при программной и позиционной стратегиях совпадают с точностью
до малой величины е. Вместе с тем при Lо е £доп все ограничения для ОУ
с программной стратегией выполняются. Это следует из следующего утверждения.
Утверждение 2. Если Яо е Rj, то есть имеет место программное ОУ и* (•), то задавшись допустимым отклонением 5 и всегда можно найти такие временной шаг квантования 5 /и синтезирующую функцию Sj (х, т; Яо), что при отсутствии возмущений значения управления и(/) V / е[% /к ] при программной (и*пр (/)) и позиционной (и*пр (/)) стратегиях не будут различаться на величину больше 5 и.
Это нетрудно показать, если предположить, что при позиционной стратегии в последовательные моменты времени ^, ^ +51, ^ + 251, ... обеспечивается решение ЗОУ для R при скорректированных значениях х и t. Обозначим скорректированный массив данных ЗОУ в момент ^ + к51 через Rk и через dij ^к) -параметры ОУу-го вида, у = 1, 2,..., рассчитанные по данным Rk.
В качестве основных показателей эффективности работы МФУУ при выборе временного шага квантования 51 можно рассматривать коэффициенты, характеризующие увеличение функционала энергозатрат кф и погрешностей достижения
конечного состояния (фазовых координат zf) кг1, i = 1, п.
Коэффициент kф будем определять по формуле
/„з(50
kф -
(19)
пр
где /пр, Iпз(5t) - значение функционалов при программной и позиционной (с шагом 5/) стратегиях.
Коэффициент кгг соответствует относительной погрешности в конечный момент времени, то есть
kz, - j
.0
z, - z,
-100% или
kz, - ■
z,K - zf(5t)
Uqq %,
(20)
здесь zгк, zi - значения исходных данных ЗОУ; zгmax - максимальное значение Zi на интервале [% ], вторая формула используется в случаях, когда zf = z0 или
разность (г;к - z0) мала.
На рис. 3 приведен качественный характер зависимостей кф(5/) и кг\(5(), на основе которых выбирается шаг квантования 5/.
Таким образом, решение задачи синтеза алгоритмического обеспечения МФУУ связано с выполнением следующих этапов:
1. Построение области допустимых решений £ доп в пространстве синтезирующих переменных для задаваемых входных данных, если £доп = 0, то последовательно ослабляются частные ограничения (см. определение 3).
2. Определение возможных видов синтезирующих функций из области £ доп.
3. Построение области практической устойчивости.
4. Определение временного шага квантования.
5. Имитация работы полученного алгоритма.
При выполнении первого этапа наибольшие трудности могут возникать в связи с частным ограничением (7). В качестве примера рассмотрим построение алгоритма ЭОР, обеспечивающего решение ЗОУ с закрепленным промежуточным значением фазовой координаты (см. рис. 1, б), более подробно.
Определение 5. Алгоритм функционирования ЭОР назовем модифициро-Рис. 3. Качественный характер ванным, если он обеспечивает решение
зависимостей кф (1) и кЛ (2) ЗОУ (2), (5) с дополнительным ограниче-
от шага квантования 5t нием (7а).
тах
z
В задаче закрепления промежуточного значения вектора х = (х^ Х2 )т возможны два случая: 1) в момент времени переключения /п закрепляются все компоненты вектора х; 2) при t = /п закрепляется только первая фазовая координата Х1 (например, температура).
В первом случае алгоритм заключается в следующем:
1) по заданным значениям к/, к^, к2 рассчитываются
= кt(/к -t0), х1 = k1Х1(t0), х2 = к2Х2(/0);
2) определяется и реализуется синтезирующая функция
&(х, X; Rtо = (А,В, «И, «в, х0, хп, /0, /п);
3) в момент времени /п определяется и реализуется синтезирующая функция
&(х, Т; Rtп ), Rtп = (А, В, «н, «в, хп, хк, , /к).
Таким образом, на интервале времени [tо, /к] последовательно используются две синтезирующие функции, то есть
!&(х,т; К. ), t е [%^);
&(х, т; К ) = #
0 [&(х,т; К.п), t е [^,/к].
Основную трудность здесь представляет определение коэффициента к2, который в общем случае зависит от к и к\.
Во втором случае значение х2(/п) определяется из условия минимизации функционала на интервале устранения рассогласования [^, /к ], то есть дополнительно решается задача минимизации суммы функционалов ^ (хп) и 1.п (хп)
по значению хп.
Сопоставление показателей работы КВОР с МФУУ, выполняющим функции ЭОР, показывает, что экономия энергозатрат может достигать 40 % и более.
Заключение. Рассмотренные задачи синтеза оптимального управления и алгоритмического обеспечения СЭУ решаются на основе результатов полного анализа соответствующих ЗОУ. Предложенный математический аппарат позволяет разрабатывать алгоритмы расчета в реальном времени энергосберегающих управляющих воздействий, которые могут реализовываться простыми микропроцессорными устройствами.
Показано, что в процессе эксплуатации объектов могут возникать ситуации, когда требуется изменить вид минимизируемого функционала. Для этих случаев предложен подход построения многофункциональных оптимальных регуляторов.
Список литературы
1. Артемьев, В.В. Теория систем со случайными изменениями структуры /
В.В. Артемьев. - М. : Наука, 1979. - 428 с.
2. Справочник по теории автоматического управления / под ред. А.А. Кра-совского. - М. : Наука, 1987. - 712 с.
3. Муромцев, Ю.Л. Метод синтезирующих переменных при оптимальном управлении линейными объектами / Ю.Л. Муромцев, Л.И. Ляпин, Е.В. Сатина // Изв. вузов. Приборостроение. - 1993. - № 11-12. - С. 19-25.
4. Ляпин, Л.Н. Оптимальный по минимуму затрат энергии регулятор объекта двойного интегрирования / Л.Н. Ляпин, Ю.Л. Муромцев, О.В. Попова // Изв. Акад. наук СССР. Техн. кибернетика. - 1992. - № 3. - С. 39-46.
5. Муромцев, Ю.Л. Теоретические основы энергосберегающего управления : монография / Ю.Л. Муромцев, Д.Ю. Муромцев, В.А. Погонин. - М. : ЯНУС-К, 2010. - 288 с.
6. Муромцев, Ю.Л. Методология «оперативного проектирования» систем оптимального управления промышленными объектами / Ю.Л. Муромцев, И.В. Тюрин // Приборы и системы. Управление, контроль, диагностика. - 2007. - № 7. -
С. 1-3.
7. Зырянов, Ю.Т. Применение энергосберегающих технологий для автоматизированного управления модульными котельными / Ю.Т. Зырянов, А.В. Петров // Энергобезопасность и энергосбережение. - 2011. - № 3 (39). - С. 13-17.
8. Муромцев, Д.Ю. Интеллектуальные управляющие устройства / А.И. Козлов, Д.Ю. Муромцев // Актуальные проблемы информатики и информационных технологий : материалы II Всерос. (VII Тамбов. межвуз.) науч.-практ. конф. (4-5 сент. 2003 г.) / Тамб. гос. ун-т им. Г.Р. Державина. - Тамбов, 2003. - С. 83-84.
9. Ляпин, Л.Н. Анализ и оперативный синтез оптимального управления в задаче двойного интегратора на множестве состояний функционирования / Л.Н. Ляпин, Ю.Л. Муромцев // Изв. Акад. наук СССР. Техн. кибернетика. - 1990. -№ 3. - С. 57-64.
Synthesis of Multifunctional Control Devices D.Yu. Muromtsev, V.V. Papin
Department “Radio-Electronic and Designing of Microprocessor Systems ", TSTU;
crems@crems.jesby. tsty. ru
Key words and phrases: energy saving optimum regulator; multifunctional control devices; problems of optimum control; problems of synthesis; square optimum regulator.
Abstract: The issues of synthesis of multifunctional control devices for energy saving management systems of dynamic objects are considered. The two optimum regulators - square and energy saving ones have been compared.
Synthese der multifunktionalen steuernden Einrichtungen
Zusammenfassung: Es sind die Fragen der Synthese der multifunktionalen steuernden Einrichtungen der Systeme der energiesparenden Steuerung der dynamischen Objekte betrachtet. Es ist das Beispiel des Vergleiches der zwei optimalen Regler angefuhrt: quadratischen und energiesparenden.
Synthese des dispositifs multifonctionnels de commande
Resume: Sont examines les problemes de la synthese des dispositifs multifonctionnels de commande des systemes de la commande conservant de l’energie des objets dynamiques. Est cite l’exemple de la comparaison de deux regulateurs : celui-ci quadratique et celui-la conservant de l’energie.
Авторы: Муромцев Дмитрий Юрьевич - доктор технических наук, профессор, заведующий кафедрой «Конструирование радиоэлектронных и микропроцессорных систем»; Папин Виктор Владимирович - студент, ФГБОУ ВПО «ТГТУ».
Рецензент: Погонин Василий Александрович - доктор технических наук, профессор кафедры «Информационные процессы и управление», ФГБОУ ВПО «ТГТУ».