УДК 681.3
ЧИСЛЕННО-АНАЛИТИЧЕСКАЯ ПРОЦЕДУРА ОПТИМАЛЬНОГО СИНТЕЗА ГИБРИДНЫХ СИСТЕМ*
Г.Н. Рогачев, В.А. Егоров
Самарский государственный технический университет 443100, г. Самара, ул. Молодогвардейская, 244 E-mail: grogachev@mail.ru
Статья посвящена способу решения задачи оптимального синтеза гибридных систем управления, содержащих непрерывный объект управления и дискретный по времени регулятор. Рассмотрен численно-аналитический алгоритм оптимального синтеза гибридных систем, включающий численную процедуру поиска оптимальных параметров регулятора и моментов его срабатывания и базирующуюся на достаточных условиях оптимальности аналитическую процедуру поиска оптимального позиционного закона управления с полной обратной связью. Приводится пример возможного использования численно-аналитической процедуры для решения задачи оптимального синтеза гибридной системы.
Ключевые слова: оптимальный синтез, гибридная система, закон управления с полной обратной связью.
В современных условиях, когда большинство систем управления техническими средствами строится на базе компьютерной техники, приоритетным направлением теории управления является исследование гибридных систем. В таких системах непрерывно изменяющиеся компоненты отражают физические законы, технологические или технические принципы, а дискретно меняющиеся моделируют работу устройств управления. Важнейшим этапом проектирования таких гибридных систем является разработка цифровых законов управления непрерывными объектами. Известно несколько подходов к решению этой задачи [1, 2]:
- переоборудование регулятора - синтез непрерывного регулятора и последующая его замена (аппроксимация) дискретной моделью;
- дискретизация объекта - построение дискретной модели непрерывного объекта и последующий синтез дискретного регулятора методами теории дискретных систем;
- прямой синтез цифрового регулятора для непрерывного объекта без каких-либо упрощений и аппроксимаций.
Первые два подхода являются приближенными и фактически означают замену одной задачи другой с целью применить известные результаты теории непрерывных или дискретных систем. В первом случае игнорируется наличие цифровой части (импульсного элемента, дискретного регулятора и экстраполятора). Однако дискретизация полученного аналогового регулятора иногда не позволяет добиться желаемого эффекта. При использовании второго подхода не учитывается поведение системы в промежутках между моментами квантования, что может дать принципи-
*Работа выполнена при финансовой поддержке грантов РФФИ №09-08-00297-а, №10-08-00754-а; ФЦП «Научные и научно-педагогические кадры инновационной России на 20092013 годы» (проекты НК 66П/11, 2010-1.3.1-230-009/8); АВЦП «Развитие научного потенциала высшей школы» (проект №2.1.2/4236).
Геннадий Николаевич Рогачев - к.т.н., доцент.
Владимир Анатольевич Егоров - аспирант.
ально неверные результаты, например вследствие скрытых колебаний. Поскольку приближенные методы проектирования могут приводить к неработоспособным решениям, на современном этапе в теории систем управления основное внимание уделяется прямым методам анализа и синтеза [3, 4]. В качестве одной из основных при разработке прямых методов синтеза оптимальных цифровых регуляторов можно выделить проблему создания методов, применимых к широкому классу задач и обладающих вычислительной надежностью.
Для решения этой задачи в настоящей работе предлагается подход, основанный на совместном использовании методов численной оптимизации для поиска оптимальных параметров регулятора и моментов его срабатывания и базирующейся на достаточных условиях оптимальности аналитической процедуры поиска оптимального закона управления с полной обратной связью. Это позволяет рассмотреть ряд новых задач, а уже известные решать в рамках единого подхода.
Рассмотрим численно-аналитическую процедуру синтеза, оптимального по квадратичному критерию качества. Примем в качестве математической модели гибридной системы совокупность линейных дифференциальных и разностных уравнений. Динамическая часть гибридной системы, задающая движение объекта управления, описывается дифференциальными уравнениями, а логическая (автоматная) часть, моделирующая работу устройства управления, - рекуррентными включениями или уравнениями. Такая модель применима для описания широкого класса систем автоматического управления техническими комплексами и технологическими процессами. Изменение состояний дискретной части может происходить в заранее заданные (тактовые) или в произвольные моменты времени. Более того, выбор тактовых моментов может считаться ресурсом управления и подлежать оптимизации так же, как параметры регулятора и оптимальный закон управления.
Пусть поведение системы управления описывается совокупностью линейных дифференциальных и разностных уравнений:
x(t) = A(t)x(t) + B(t)u(t) + F(t)y(tk )Л ^ t < tk+i>
y(tk+i) = Ck+iy(tk) + Dk+iv(tk+i) + Gk+ix{tk+i), k = 0Д,...,N-1, (1)
где x, у - векторы состояния непрерывной и дискретной частей системы соответственно, x е Rn, y е Rm ; u, v - векторы управления непрерывной и дискретной частями
соответственно; u eU œ Rq, v eV œ Rs, U и V - заданные множества допустимых значений управления; t - время, t eT = [t0, t^ ), T - промежуток времени функционирования системы, на котором выделены моменты tk,k = 0,i,...N—i, разбивающие множество Т на непересекающиеся подинтервалы Tk = [tk,tk+i),k = 0,i,...N—i;
f (t, x, y, u): T x Rn x Rm xU ^ Rn, x(/n ) = x(/n — 0) . Входящие в (1) матрицы A(t), B(t), F(t) имеют соответственно размер (n x n), (n x q), (n x m) , их элементы непрерывны; Ck+i,Dk+i, Gk+i — матрицы размера (m x q), (m x s), (m x n) соответственно;
на управление ограничений не наложено, u е U = Rq, v eV = Rs. Начальное состояние системы (1) задано:
x(t0) = x0 е Rn, y(t0) = y0 е Rm. (2)
Конечное состояние (x(tN), y(tN)) системы произвольно. Предполагается, что при управлении используется информация о текущем времени и о векторе состояния
(х,у). Множество допустимых управлений с полной обратной связью ип образуют
функции и(/, х, у): Т х Яп х Я1” ^ и и v(t, х, у):(^,..., х Яп х Ят ^ V, которые на траекториях системы (1) для различных начальных состояний (2) порождают допустимые процессы ё = (х *(•), у *(•), и *(•), V *(•)) е О(/0, х0, Уо). На множестве О(/0, х, Уо ) задан функционал качества управления:
I =
N—1 ^к+1 р 1 *
- хТ (/ )5 (/) х(/) + - ит ( / )<2^ )и(/) + хТ (/ )Ф(/) у^к)
N —1 к+1
2« I
1 г
& + 2 У (/к)Як+1у(/к) +
/к
1 Т/ \ г> / \ Т/ / \ 1 Т/
+ - v (/к+1)рк+^(/к+1) + х (/к+№+1У(/к) + - х (/к+1)Лк+1х(/к+1)} + (3)
1 т 1 Т Т
+ ~ х (/К )Л№+1х(/№ ) + ^ У (/К )Я№+1у(/№ ) + х (/К У^К+1у(/К X
где 5(/), Як+1, Лк+1 - неотрицательно определенные симметрические матрицы размера (п х п), (т х т) , (п х п) соответственно; Q(t), р+1 - положительно определенные матрицы размера (д х д), (5 х 5) соответственно; Ф(/), ^+1 - матрицы размера (п х т) . Определению подлежат матрицы С к+1, Ок+1, +1, моменты /к, к = 0,1,...N—1
разбиения промежутка времени функционирования системы Т и оптимальное управление (и *(/, х, у), V *(/, х, у)) еип с полной обратной связью, минимизирующие функционал (3).
Поставленная задача естественным образом разбивается на две подзадачи. Для определения параметров регулятора и режима его работы (матриц С к+1, Ок+1, Ок+1 и моментов ^, к = 0,1,...N—1) целесообразно использовать один из численных методов глобальной оптимизации (генетическое программирование, метод отжига). Эта процедура является внешней по отношению к процессу вычисления предельно возможного значения критерия качества оптимального управления с полной обратной связью, который должен выполняться на каждом шаге внешней процедуры при фиксированных значениях параметров регулятора. Внутренняя аналитическая процедура, подробно описанная в [5], базируется на достаточных условиях оптимальности позиционного закона управления с полной обратной связью и состоит в решении на каждом из временных подинтервалов задачи Коши для уравнения с частными производными первого порядка относительно функции Беллмана. Так, в линейном случае функция Беллмана ф(/, х, у) ищется в виде
ф/,х,у)=\хК) ^ + щ) хТу+2 ум(/) уТ,
где К(/), Ь(/),М(/) - неизвестные матрицы размера (п х п), (п х т), (т х т) соответственно. Для нахождения матриц К(0, Ь(0, М(1) используются уравнения
К + АТК + КА + KBQ ~1БТК — 5 = 0,
Ь + АТЬ + КБ + KBQ ~1БТЬ — Ф = 0, < / < /*+и (4)
М + БТЬ + ЬТ Б + Ь BQ ~1БТ Ь = 0,
условия перехода к следующему полуинтервалу (условия скачка)
34
К (к+1 — 0) = К + ЬО + ОТЬТ + ОТ МО — Л + [Ь + ОТМ ]О[Р — ОМО]~1О [ЬТ + МО],
Ь(/к+1 — 0) = ЬС + ОТ МС — ¥ + [Ь + ОТМ ]О[Р — ОМБ]~1 ОТМС, (5)
М (гк+1 — 0) = СТМС — Я + СТМОР — ОТМО]—1 ОТМС и конечные условия
К(tN ) = —ЛN+1, Ь(/¥ ) = —^¥+1, М(tN ) = —^N+1. (6)
В правых частях выражений (5) для вычисления К (/к+1 — 0), Ь(/к+1 — 0), М (/к+1 — 0) матрицы К (/), Ь(/), М (/) берутся в момент / = /¿+1, а матрицы Р, Я, ¥, Л, С, О, О имеют индексы к+1.
Оптимальное управление на каждом интервале вычисляется так:
и *(/, х, у(/к)) = Q“1(t)BT(/)[К(/)х + Ь(/)у(/к)], /к < / < /к+1,
V * (^1, х(/к+1), у(гк)) = [Р+1 — О+М С/*+1) Ок+1]—1ОТ+1 х {[ ЬТ (/а+1) + М (/к+1)Ок+1]х(/к+1) +
+М(/к+1)Ск+1у(/к)), к = 0,1,...,N — 1. (7)
Алгоритм вычисления оптимального управления состоит из следующей последовательности шагов:
1. Из конечного условия (6) определяются матрицы К(/к+1), Ь(/к+1), М(/к+1) при к = N-1.
2. По формуле (5) вычисляются матрицы К(1к+1 — 0), Ь(1к+1 — 0), М(1к+1 — 0) .
3. Решается система уравнений (4) справа налево от момента / = — 0 до / = ^
с конечными условиями, найденными в п. 2. В результате находятся матрицы К(0, Щ, М(г).
4. Выполняются пп. 2, 3 для к = N — 2, N — 3,..., 1,0. В результате находятся матрицы К(0, Ь(0, М(1) на всех полуинтервалах Тк.
5. Подставляя найденные матрицы К((),Ь((),М(0 в (7), получим искомое оптимальное управление с обратной связью.
Минимальное значение функционала (3) находится по формуле:
тп 1(ё) = —„(0, х0, у)) = —1 хТ К (^К — хТ Ь(/0)у0 —1 у^М^^. (8)
ёеО(/0,х0,у0) 2 2
Выделим основные достоинства численно-аналитической процедуры:
— используемые численные методы оптимизации позволяют выявлять глобально-оптимальные настройки и режим работы регулятора;
— критерий качества (3) весьма общего вида позволяет решать самые разнообразные задачи;
— с помощью такой процедуры определяется точное решение задачи оптимального управления с полной обратной связью;
— значение критерия качества может быть рассчитано по (8) без определения оптимального управления, что позволяет снизить объем вычислений на промежуточных этапах поиска оптимальных настроек и режимов работы регулятора.
Помимо этого, такой подход позволяет поставить и решить ряд новых задач, а уже известные рассматривать с единых вычислительных позиций. Например, отличительной чертой цифрового регулятора является то, что при подаче на вход единичного дискретного импульса или ступенчатого сигнала переходный процесс полностью заканчивается за конечное число тактов, т.е. за конечное время. Это явление не имеет аналога для непрерывных систем. Такое управление называют апериодическим (в англоязычной литературе - deadbeat control) [1]. Рассмотрим в качестве примера задачу по переводу объекта управления в виде двойного интегратора из произвольного начального состояния в нулевое конечное за минимальное число шагов срабатывания цифрового регулятора. Подобная проблема исследуется в [4], там она решена как проблема размещения полюсов передаточной функции, определен цифровой регулятор второго порядка с дискретной передаточной функцией
2.5 - 1.5z
-1
— . Включаясь в работу, он переводит объект управления в нулевое
Ж =
1762 1 + 0.75 г_
конечное состояние за три такта. Сформулируем эту же задачу как задачу оптимального синтеза гибридной системы управления. В этом случае
"0 1" "0" О О О О
(( N 0 0 (( N 1 0 1 (( N 1 0 II t) © II t) to і 0 0 1
"0 0" "0 0" "1 0"
Rk -^ = , к = 1,2,., N +1 , Лк - ,к -1,2лЛк -
0 0 0 0 0 1
, t єТ = [t0, tN ),
к - 3,4,..., N +1.
Из-за требования положительной определенности Q(t) и р вынужденно примем Q(t) = 0.01,ї еТ = [ї0,) и р = 0.01,к = 1,2,...,N+1. Считаем, что структура регулятора остается неизменной в процессе управления, т.е. матрицы Ск+1, Ок+1, Ок+1 постоянны, Ск+1 = С,Ок+1 = О,Ок+1 = О,к = 0,1,...,N—1. Матрицы С,О,О подлежат определению. Если будет выбран аритмический режим работы регулятора, то дополнительно будет необходимо определить длительности подинтервалов или, что эквивалентно, моменты їк, к = 0,1,...N—1 разбиения промежутка времени Т функционирования системы. С целью предотвращения скрытых колебаний в системе можно
Г1 0“
изменить параметры критерия качества (3), выбрав £ (ї) =
0 1
t ЄТ - [t0, tN ) .
Таким образом, задача синтеза апериодического регулятора может быть поставлена и решена как задача оптимального синтеза гибридной системы управления.
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1. Острем К., ВиттенмаркБ. Системы управления с ЭВМ. - М.: Мир, 1987. - 480 с.
2. Джури Э.И. Импульсные системы автоматического регулирования. - М.: Физматгиз, 1963. - 456 с.
3. Розенвассер Е.Н. Линейная теория цифрового управления в непрерывном времени. - М.: Наука, 1994. - 455 с.
4. Поляков К.Ю. Основы теории цифровых систем управления / СПбГМТУ. - СПб., 2006. - 161 с.
5. ПантелеевА.В. Теория управления в примерах и задачах. - М.: Высшая школа, 2003. - 583 с.
Статья поступила в редакцию 2 сентября 2010 г.
UDC 681.3
NUMERICAL ANALYTIC PROCEDURE OF THE OPTIMAL SYNTHESIS OF HYBRID SYSTEMS
G.N. Rogachev, V.A. Egorov
Samara State Technical University
244, Molodogvardeyskaya st., Samara, 443100
E-mail: grogachev@mail.ru
The article is devoted to the method of solving the problem of optimal synthesis of hybrid control systems that contain continuous control object and discrete-time controller. A numerical analytic algorithm for optimal synthesis of hybrid systems, which includes a numerical procedure for searching the optimal parameters of the controller and the moments of its release and based on sufficient optimality conditions analytic procedure for searching the optimal position control law with the complete feedback. Examples of use of numerical analytic procedure for the solving some problems of optimum synthesis of hybrid systems are given.
Keywords: Optimal Synthesis, Hybrid System, Control Law with Full Feedback.
Gennadiy N. Rogachev - Candidate of Technical Sciences, Associate professor. Vladimir A. Egorov - Postgraduate student.