УДК 621.8
ПРОДУКЦИОННЫЙ МЕТОД СИНТЕЗА РЕГУЛЯТОРОВ СИСТЕМ УПРАВЛЕНИЯ НЕПРЕРЫВНО-ДИСКРЕТНЫМИ ОБЪЕКТАМИ
Г.Н. Рогачев
Самарский государственный технический университет 443100, г. Самара, ул. Молодогвардейская, 244
E-mail: grogachev@mail.ru
Решена задача синтеза совокупности правил работы регулятора системы управления непрерывно-дискретным объектом. Для этой цели использовалась численноаналитическая процедура, основанная на совместном применении методов конечномерной оптимизации для определения условной части каждого из правил как логического выражения темпоральной логики и базирующаяся на достаточных условиях оптимальности аналитической процедуры поиска на каждом из временных интервалов закона управления с полной обратной связью.
Ключевые слова: гибридная система, непрерывно-дискретная система управления, совокупность правил работы регулятора.
Одним из важных направлений современного развития теории управления является исследование гибридных непрерывно-дискретных систем управления (ГНДСУ). В таких системах компоненты с непрерывными сигналами отражают физические законы, технологические или технические принципы, которым подчинено функционирование объектов управления, а дискретные элементы моделируют работу цифровых устройств управления. Различные классы гибридных систем изучались в работах Емельянова С.В., Уткина В.И. (системы с переменной структурой); Васильева С.Н., Федосова Е.А. (логико-динамические системы); Куржанского А.Б., Цыпки-на Я.З., Li Z., Silva G.N., Soh Y., Vinter R.V., Wen C. (импульсные системы); Аграче-ва А.А., Савкина А.В., Antsaklis P.J., Brockett R.W., Evans R.J., Hedlund S., Liberzon D., Rantzer A., Rischel H. (переключательные системы); Hespanha J.P., Johansson K.H., Lygeros J., Sastry S.S., Tomlin C.J. (автономные и неавтономные гибридные системы).
Задачи управления применительно к гибридным системам оказываются качественно более сложными по сравнению с аналогичными задачами управления непрерывными или дискретными системами ввиду целого ряда принципиальных особенностей этих систем, связанных с наличием в них существенным образом взаимодействующих непрерывных и дискретных подсистем, сочетанием непрерывной динамики управляемых процессов с дискретными процедурами получения информации о параметрах объекта, вычисления по этим данным управляющих воздействий и передачи их на исполнительные устройства. Ситуация усугубляется необходимостью учета эффекта квантования сигналов в гибридных системах управления как по времени, так и по уровню, влияния различных возмущающих факторов на качество их работы. Так, в реальных системах управления передача управляющего воздействия может происходить с меняющимися по случайному закону временными задержками и потерей части информации. На входные и выходные сигналы накладывается шу-
Геннадий Николаевич Рогачев (к.т.н., доц.), доцент кафедры «Автоматика и управление в технических системах».
22
мовая составляющая. Эти эффекты особенно присущи распределенным сетевым системам управления.
К гибридным относится широкий круг управляемых объектов, охватывающий как традиционные, так и новейшие технологии в самых различных областях техники. Практическая реализация таких технологий с требуемыми качественными показателями невозможна без построения соответствующих алгоритмов автоматического управления, что и явилось главным стимулом к созданию теории и методов управления гибридными системами. Программа фундаментальных научных исследований Российской академии наук относит к числу первоочередных исследований в области теории систем вопрос создания общей теории управления сложными техническими и другими динамическими системами, в том числе единой теории управления, вычислений и сетевых связей. Научный семинар ИПУ РАН «К единой теории управления, вычислений и связи» (3 апреля 2008 г.) констатировал, что движение к единой теории управления, вычислений и связи (Control + Communication + Computation = C3) - важная тенденция, обусловленная запросами практики и осознанная мировым научным сообществом. Группа ведущих ученых в коммюнике [1] «Future Directions in Control in an Information-Rich World» («Грядущие изменения в управлении в информационном обществе») определяет стратегию развития теории и практики управления техническими системами. Среди основных задач - создание технических устройств, самостоятельно генерирующих программное обеспечение своих систем управления и при необходимости модернизирующих его, автоматический синтез алгоритмов управления с одновременными их верификацией и валидацией.
Указанные проблемы могут быть решены путем разработки новых эффективных и реализуемых на практике методов решения задач управления как рассматриваемых в качестве задачи организации на цифровом регуляторе необходимых вычислительных процедур с учетом особенностей, привносимых цифровым характером устройств управления и сетевым способом обмена информацией между ними. В том числе оказываются необходимыми новые способы синтеза оптимальных по принятым критериям эффективности алгоритмов и программ работы цифровых устройств управления для непрерывно-дискретных систем, учета различных возмущающих факторов и ресурсных ограничений, характерных для цифровых и сетевых систем управления. Такие методы позволят повысить эффективность процесса разработки автоматических регуляторов за счет исключения промежуточных этапов и гарантировать техническую реализуемость получаемого в виде псевдокода закона управления.
В работе [2] рассмотрена численно-аналитическая процедура оптимального синтеза процедурной модели регулятора в виде совокупности правил его работы, основанная на совместном использовании методов конечномерной оптимизации для определения количества правил и поиска условной части каждого из них как логического выражения темпоральной логики и базирующейся на достаточных условиях оптимальности аналитической процедуры поиска на каждом из временных интервалов закона управления с полной обратной связью (части «действие» каждого из правил).
В настоящей работе эта процедура используется для решения ряда новых задач управления совокупностью линейных дифференциальных и разностных уравнений вида
x(t) = A(t)x(t) + B(t)u(t) + F(t)y(tk X tk < t < tk+1, y(tk+1) = Ck+1 y(tk) + Dk+1v(tk+1) + Gk+1 x(tk+l), k = 0,1,...,N -1,
где х, у - векторы состояния непрерывной и дискретнои частей системы соответственно, х е Я",у е Ят ;
и, V - векторы управления непрерывной и дискретной частями соответственно, и еи с Яя, V е¥ с Я;
и и V - заданные множества допустимых значений управления; t - время, t е Т = [t0,^);
Т - промежуток времени функционирования системы, на котором выделены моменты tk,к = 0,1,...N -1, разбивающие множество Т на непересекающиеся подин-тервалы
Тк = Ш, tк+1), к = 0,1,...N -1.
Так, рассматривалась задача оптимального управления ГНДСУ вида
) = У^к) + u(t), t е Тк, х(0) =1 У(/к+1) = х(*к+1) + ^к+1) , Я0) = 1, к = 0,1
х е Я;и е Я,у е Я; V е Я^ е Т = [0,2);
1 ^ ^ (*гк+1 2 2 1 2 2 1 критерий качества I = — > I и (г)Л + V (ґк+1) + — [л (2) + у (2)] ^ шт
2 к=0 - к -I 2
Известно [3] решение этой задачи при наличии одной степени свободы регулятора, когда условием г є Т = [0,2), Т0 = [0,1), Т = [1,2), N = 2, г0 = 0, ^ = 1, Ґ2 = 2 заданы количество правил системы продукций и часть «если» каждого правила. Это - оптимальное управление (и *(г, л, у), V *(г, л, у)) с полной обратной связью. Система продукций регулятора в этом случае (вариант № 1) имеет вид
, гт 3л - 3 уг + 3у 0
«если г є Т0, то и =----------------, V = -0.75л »;
0 3г - 5
. т 3л - 3уг + 6у
«если г є л , то и =----------------, V = -0.5л ».
1 3г - 8
Возможность использовать вторую степень свободы регулятора позволяет изменять часть «если» правила его работы. Этим в рассмотрение включаются варианты неравномерного во времени срабатывания регулятора (аритмия). Найдено решение задачи при наличии двух степеней свободы регулятора. Установлено, что такая логика его работы приводит к улучшению качества управления в широком диапазоне начальных условий. Количество правил задано и по-прежнему равно двум, часть «если» каждого правила определяет неравномерность срабатывания регулятора, далее часть «то» каждого правила подлежит определению посредством описанной ранее процедуры.
Вариант № 2: г є Т = [0,2),Т0 = [0,1.5), Т = [1.5,2),N = 2,г0 = 0,г1 = 1.5,г2 = 2 . Система продукций регулятора имеет вид
гт 0.85лг - 1.8л - 0.85лг2 + 3.1>>г - 2.7 у
«если г є Т0, то и = 260000-------------------------------------7-^-ч---, V = -0.529л »,
0 (560г - 1200)(400г - 840)
. т 3л - 3уг + 6у
«если г є Т , то и =----------------, V = -0.5л ».
1 3г - 8
Вариант № 3: г є Т = [0,2),Т0 = [0,0.5),Т1 = [0.5,2),N = 2,г0 = 0,г1 = 0.5,г2 = 2 . Система продукций регулятора имеет вид
«если і є Т0, то и = 84000
- 0.38хі + 0.46 х - 0.38хі + 0.65 уі - 0.23 у
(220? - 270)(і40і -170)
, V = -0.85х »,
«если і є Т1, то и =
3х - 3 уі + 6 у 3І-8
V = -0.5х ».
Результаты анализа качества работы регуляторов с полученными системами продукций представлены на рис. 1. На этих графиках более светлые участки отображают область изменения начальных условий, при которых указанный вариант алгоритма превосходит другие. Как видно из приведенных результатов, наименьшую область предпочтительности имеет вариант № 1, характеризующийся равномерным режимом срабатываний регулятора, наибольшую - вариант № 3.
-10 в;
Вариант № 1 Вариант № 2 Вариант № 3
Рис. 1. Область предпочтительности вариантов № 1 , № 2, № 3
Далее рассматривалась задача оптимального синтеза цифрового программного регулятора линейным непрерывным объектом второго порядка. Требуется перевести ГНДСУ вида
X = - х1 + х2;
х2 = - х2 + у ); у(ік+1)=+0,
где х1 є Я, х2 є Я, и є Я, у є Я; V є Я, і є Т = [0,3), Т = [і -1, і), і = 1, 2, ..., 6,
Х1 (0) = х10, х2 (0) = х20, Я°) = 0 ,
в малую окрестность начала координат, при этом эвклидова норма вектора управляющих воздействий должна быть минимальной,
5 г -I 1 5
1 2 ,
■ тіл .
1 = 2 ^[х12(ік+1) + х2 (ік+1) ]+ 2 2 700 v2 (ік+1) 2 к=2 2 к=0100
Сравнивались два варианта решения: с одной степенью свободы регулятора при равномерном шаге и с двумя степями свободы регулятора при аритмическом регуляторе с неравномерным и определяемым в процессе поиска шаге его срабатывания. Оптимальное значение функционала равно
ОПТ1
и їґ
(0-
= (1.45х20 + 7.07х10х20 + 8.91х20 )/10000
= (0.772х20 + 3.44х10х20 + 4.11х20 1/100 00 соответственно
,)/
1ОПТ2 - \у-' '^Л10 "г-’-‘т‘тл10л20 "Г‘т-Ііл20/І Анализ полученных решений показал, что квадратичная фор-
2
■ 7.07х10х20 + 8.91х20 )/10000 - 1.878(0.772х20 + 3.44х10х20 + 4.11x2
)/10000
ма ^1.45х10 + 7.и/хюх20 + 0.71x2^/ 10000 — 1.0/ оуи. / /2Л]о + 3.44x10x20 + 4.11x20) является положительно полуопределенной, т. е. качество управления по второму сценарию, даже умноженное на коэффициент 1.070, не хуже качества управления по
первому сценарию при том же самом времени достижения малой окрестности начала координат (3 с). Так, при нулевом значении х20 = 0 качество управления по второму сценарию лучше в 1.47/.772 = 1.904 раза.
Кроме того, был выявлен эффект уменьшения необходимого количества правил работы регулятора при его аритмической работе. При наличии двух степеней свободы регулятора поставленная задача фактически решается за два такта срабатывания регулятора вместо трех (рис. 2, 3).
і
Рис. 2. Поведение ГНДСУ с ритмическим регулятором
І
Рис. 3. Поведение ГНДСУ с аритмическим регулятором
Таким образом, рассмотренная численно-аналитическая процедура оптимального синтеза процедурной модели регулятора позволяет решать задачи синтеза закона управления линейными объектами, оптимального по квадратичному критерию качества. Однако попытка использовать его для решения более широкого круга задач сопряжена с рядом проблем, подчас трудноразрешимых. Так, затруднен синтез системы продукций регулятора, управляющего нелинейным объектом. Причина - невозможность установления в этом случае аналитического выражения для функции Беллмана. Определенные трудности представляет синтез системы продукций регулятора при наличии ограничений на фазовые переменные и управляющие воздействия. Частично снять эти проблемы можно при использовании альтернативных вариантов, например эволюционных вычислений как стохастического метода глобальной оптимизации [4].
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1. Murray R.M. [et al.] Future Directions in Control in an Information-Rich World // IEEE Control Systems Magazine. - 2003. - Vol. 23, No. 2. - PP. 20-33.
2. Рогачев Г.Н., Егоров В.А. Численно-аналитическая процедура оптимального синтеза гибридных систем // Вестник Самарского государственного технического университета. Серия: Технические науки, 2010. - № 7. - С. 32-38.
3. Пантелеев А.В., Бортаковский А. С. Теория управления в примерах и задачах. - М.: Высшая школа, 2003. - 583 с.
4. Рогачев Г.Н. Использование генетического алгоритма с отсечением по времени в задаче синтеза программного регулятора для машины Дубинса // Вестник Самарского государственного технического университета. Серия: Технические науки, 2011. - № 3. - С. 27-33.
Статья поступила в редакцию 20 декабря 2012 г.
PRODUCTION SYNTHESIS METHOD OF CONTROLLERS FOR CONTINUOUS-DISCRETE OBJECTS
G.N. Rogachev
Samara State Technical University
244, Molodogvardeyskaya st., Samara, 443100
The problem of controller design for continuous-discrete object is solved. A set of rules of the regulator is synthesized. Numerical-analytical procedure is used for this purpose. It is based on the joint application offinite-dimensional optimization methods for the determination of the conditional part of each of the rules as a logical expression of temporal logic and on sufficient optimality conditions for the analytical procedure of searching for each of time intervals of the control law with full feedback.
Keywords: hybrid system, continuous-discrete control system, a set of the regulator’s rules.
Gennady N. Rogachev (Ph.D. (Techn.)), Associate Professor.