CC BY
78
А
нализ и синтез систем управления
УДК 517.977
Н -УПРАВЛЕНИЕ ЭНЕРГЕТИЧЕСКОЙ СИСТЕМОЙ
¥
В АВАРИЙНОМ РЕЖИМЕ.
Ч. 1. Теоретические иснпны синтеза робастных Н¥-регуляторон
А.П. Курдюков, В.Н. Тимин
На основе метода размещения полюсов замкнутой системы в заданной области на комплексной плоскости и Нш-метода формирования контура с помощью аппарата линейных матричных неравенств решена задача подавления внешних возмущений в энергетической системе.
Ключевые слова: управление энергетической системой, Ню-теория управления, линейные матричные неравенства, подавление возмущений, Ню-метод формирования контура, расположение полюсов замкнутой системы в заданной области на комплексной плоскости.
ВВЕДЕНИЕ
Современная энергетическая система (ЭС) представляет собой сложную многоуровневую и многосвязную динамическую систему. Задачи управления и защиты в нештатных ситуациях решаются с помощью элементов силовой электроники и микропроцессорной техники. Создание управляемых гибких линий электропередачи переменного тока (Flexible AC Transmission Systems—FACTS) [1—4] на основе элементов силовой электроники позволило существенно улучшить управление потоками мощности в электрических сетях в установившихся и переходных режимах работы. Устройства FACTS позволяют увеличить пропускную способность линии электропередачи и управлять параметрами, определяющими функционирование ЭС, включая реактивное сопротивление, ток, напряжение, углы фазовых сдвигов в узлах сети, затухание колебаний на различных частотах и др. В процессе функционирования ЭС при коммутациях в сети и изменении нагрузки возникают электромеханические колебания, компенсация которых обычно осуществляется автоматическими регуляторами возбуждения и системными стабилизаторами. Применение устройств FACTS совместно с системными стабилизаторами позволяет дополни-
тельными управляющими воздействиями улучшить качество функционирования ЭС. Синтез законов управления в ЭС с целью повышения робастной устойчивости и качества переходных процессов, демпфирования колебаний перетоков мощности и улучшения других характеристик, представляет собой достаточно сложную и актуальную задачу. В целях ее решения в данной работе исследуется возможность применения Нт-ро-бастных регуляторов для управления ЭС в аварийном режиме.
В первой части статьи кратко изложены теоретические основы синтеза Н^-робастных регуляторов с помощью метода размещения полюсов замкнутой системы в заданной области на комплексной плоскости и Н^-метода формирования контура. Во второй части будет дана постановка и предложен метод решения задачи синтеза Нт-регулятора для подавления возмущений, возникающих в энергетической системе.
1. СОВРЕМЕННЫЕ КОРНЕВОЙ И ЧАСТОТНЫЙ МЕТОДЫ СИНТЕЗА СИСТЕМ
С инженерной точки зрения под синтезом системы управления можно понимать построение системы управления, обеспечивающей выполнение
всех предъявляемых к ней технических требований. Проектируемая система управления должна обеспечивать приемлемое качество переходных процессов, необходимую точность и грубость (робастность) к неопределенности в модели объекта.
Требования к проектируемой системе могут формулироваться как во временной, так и в частотной областях [5]. Качество переходных процессов в смысле запаса устойчивости (удаленности от границы устойчивости) и быстродействия (быстроты реагирования на управляющие и возмущающие воздействия) может характеризоваться во времени видом переходного процесса на типовое воздействие, расположением полюсов и нулей передаточной функции замкнутой системы и др.
Во временной области требования к запасу устойчивости по виду кривой переходного процесса при единичном воздействии могут задаваться в виде, например, максимального значения регулируемой величины, значения перерегулирования, допустимого числа колебаний и др. Быстродействие системы может определяться по длительности переходного процесса.
Расположение и вид корней характеристического полинома (полюсов передаточной функции замкнутой системы) также определяют характер переходного процесса, так как они однозначно определяют решение однородного дифференциального уравнения, описывающего модель объекта. Поэтому требования к запасу устойчивости и быстродействию, не рассматривая кривые переходного процесса, можно сформулировать, накладывая ограничения на корни характеристического полинома.
Следовательно, задавая ограничения на расположение и вид корней характеристического полинома, мы задаем ограничения во временной области на вид переходного процесса и, очевидно, на максимальные значения регулируемой величины, значения перерегулирования и времени переходного процесса и др. Поэтому корневой метод размещения плюсов передаточной функции замкнутой системы в заданной области на комплексной плоскости достаточно удобен для задания требований к качеству переходных процессов.
Таким образом, требования к расположению полюсов передаточной функции замкнутой системы в некоторой заданной области позволяет задавать ограничения на:
• быстродействие, которое может быть оценено по степени устойчивости (абсолютной величине действительной части ближайшего к мнимой оси корня);
• запас устойчивости, который может быть оценен по колебательности (максимальному отношению мнимой части корня к его действитель-
ной части) или параметру затухания (коэффициенту демпфирования);
• максимальную угловую частоту собственных колебаний системы (максимальную абсолютную величину мнимой части корня).
В теории управления широко используется аппарат линейных матричных неравенств (ЛМН) [6—10]. Он хорошо подходит для синтеза регулятора в обратной связи, когда различные требования к замкнутой системе можно задавать в виде линейных матричных неравенств.
Концепция расположения полюсов в выпуклой области, которая может быть описана линейными матричными неравенствами, предложена в работе [6]. Класс областей, описываемый ЛМН (см. п. 2.2), представляет собой достаточно широкое множество областей для размещения полюсов, включая полуплоскости, круги, секторы, вертикальные и горизонтальные полосы, а также любые пересечения этих областей. В соответствии с данной концепцией процедура синтеза регулятора сводится к минимизации какого-либо критерия оптимальности на данной ЛМН-области. Так, в работах [6, 7] задача синтеза Н^-оптимального регулятора с расположением полюсов передаточной функции замкнутой системы в ограниченной выпуклой области на комплексной плоскости сводится к решению ЛМН.
Вопросы чувствительности к возмущениям и шумам, а также точностные характеристики при таком подходе считаются менее сложными и не рассматриваются, однако они требуют дополнительного решения.
Альтернативный путь дает синтез системы управления в частотной области. Он основан на частотных свойствах объекта, частотных критериях устойчивости и качества переходных процессов, частотных методах анализа и синтеза систем [5].
Одним из широко распространенных современных методов синтеза многомерных систем в частотной области является Н^-метод формирования контура, предложенный в работе [17]. Под методом формирования контура замкнутой системы в частотной области понимают метод синтеза, при котором разработчик непосредственно определяет форму амплитудно-частотной характеристики разомкнутой или замкнутой систем. В общем случае передаточная функция многомерной системы матричная. Элементы матрицы представляют собой передаточные функции от каждого входа к каждому выходу системы. Началом построения более строгой классической теории управления в частотной области, распространенной на многомерные системы, послужила статья [11], посвященная роли сингулярных чисел в анализе и синтезе многомерных систем управления. Зависимость сингулярных чисел матричной передаточной функции
от частоты представляет собой обобщение амплитудной частотной характеристики. В системах с одним входом и одним выходом эти понятия тождественны.
Сингулярные числа комплекснозначной матрицы Л определяются как квадратные корни из собственных чисел матрицы Л*Л, где Л* — эрмитово сопряженная матрица. Они характеризуют по соответствующему выходу многомерной системы коэффициент усиления на заданной частоте [11]. Для квадратной матрицы передаточных функций число сингулярных чисел равно числу выходов системы. На заданной частоте минимальный (максимальный) коэффициент усиления многомерной системы характеризуется соответственно минимальным (максимальным) сингулярным числом.
Для многомерных систем требования к форме амплитудно-частотных характеристик замкнутой системы формулируются в терминах ограничений типа неравенств на сингулярные числа многомерных передаточных функций замкнутой [11—15] или разомкнутой [16, 17] систем.
Требования к замкнутой системе могут быть сформулированы в частотной области в виде требований к максимальным и минимальным значениям сингулярных чисел многомерной передаточной функции разомкнутой системы [16]. Принципиальная идея Н^-метода формирования контура в данном случае в многомерных системах заключается в том, что максимальные сингулярные числа матричных передаточных функций замкнутой системы, характеризующие коэффициент усиления по соответствующему выходу, могут быть непосредственно аппроксимированы в соответствующих частотных диапазонах максимальными или минимальными сингулярными числами соответствующих матричных передаточных функций разомкнутой системы.
Таким образам, основные требования состоят в:
• высоком коэффициенте усиления (ограничение снизу на минимальные значения сингулярных чисел (обозначение ст(*)) на низких частотах, позволяющем уменьшить ошибку слежения и влияние возмущений на входе и выходе объекта;
• малом коэффициенте усиления (ограничение сверху на максимальные значения сингулярных
чисел (обозначение а (•)) на высоких частотах для парирования шума и повышения робастности (грубости) к неопределенности в модели объекта;
• плавном переходе в области частоты среза контура.
В области частоты среза коэффициент усиления контура не должен убывать быстрее, чем —20 дБ на декаду, чтобы были желаемые робастная устой-
чивость, запасы по амплитуде и фазе, перерегулирование и параметр затухания [18, 19].
Рассмотрим замкнутую систему управления с объектом управления G и регулятором К, изображенную на рис. 1, где г — командный сигнал, е — ошибка слежения, и — управление, у — выход объекта, d — возмущение, п — шум измерения. Передаточная функция разомкнутой системы Ь = GK. Передаточная функция замкнутой системы с отрицательной обратной связью от командного сигнала г к ошибке слежения е называется функцией
чувствительности ^ = (I + GK) 1, а передаточная функция от командного сигнала г к выходу у объекта называется функцией дополнительной чувствительности Т = GK(I + GK) 1, где I — единичная матрица соответствующей размерности. В принятых обозначениях на рис. 2 приведена графическая иллюстрация основных требований к формированию контура разомкнутой системы. Под желаемой формой контура понимается желаемое поведение графиков сингулярных чисел от частоты для многомерной передаточной функции разомкнутой системы. На рисунке по оси абсцисс отложена угловая частота ю в логарифмическом масштабе, а по оси ординат А — значения сингулярных чисел передаточных функций в децибелах.
Желаемая форма контура передаточной функции разомкнутой системы определяется тем, что на графики минимального а( GK) и максимального а (GK) сингулярных чисел передаточной функции разомкнутой системы накладываются частотные ограничения (рис. 2). В области низких частот ю < ю1 график минимального сингулярного числа расположен выше заданного проектировщиком ограничения Ах(ю). В области высоких ю > юи график максимального сингулярного числа расположен ниже заданного проектировщиком ограничения Ац(ю). Ограничения в области средних частот не налагаются, в чем состоит недостаток метода.
При выполнении частотных ограничений на сингулярные числа передаточной функции разомкнутой системы гарантируется, что графики сингулярных чисел передаточной функции замкнутой системы также удовлетворяют некоторым ограничениям. В области низких частот ю < ю1
Рис. 1. Замкнутая система управления
Рис. 2. Требования к сингулярным числам разомкнутой системы
график обратного значения максимального сингулярного числа 1/ст ) функции чувствительности
^ = (I — ОК )-1 расположен выше графика минимального сингулярного числа передаточной функции разомкнутой системы ст(ОК). В области высоких частот ю > юи график максимального сингулярного числа ст (Т) функции дополнительной
чувствительности Т = ОК(1 — ОК) 1 расположен ниже графика максимального сингулярного числа
ст (ОК) передаточной функции разомкнутой системы. Более подробно указанные свойства передаточных функций рассматриваются в п. 2.6.
В работе [17] предложен оригинальный Нт-ме-тод формирования контура для многомерных систем (см. п. 2.6). В ней впервые выдвинута идея применения нормализованной взаимно простой факторизации [20] модели объекта для обеспечения робастной устойчивости на основе решения задачи Н^-оптимизации (см. п. 2.4). Центральный элемент в предложенном методе заключается в решении задачи робастной стабилизации [21].
Задача робастной стабилизации состоит в нахождении максимальной неопределенности (робастного запаса устойчивости), при которой единственный регулятор К обеспечит устойчивость замкнутой системы для семейства возмущенных объектов с ограниченными по Н^-норме возмущениями (см. п. 2.5).
В методе Н^-формирования контура задача робастной стабилизации решается не для исходного, а для нового обобщенного объекта, образованного с помощью формирующих функций (префильтра и постфильтра) — более подробно см. п. 2.6. Требуемые характеристики качества и робастности замкнутой системы обеспечиваются благодаря выбору соответствующих формирующих функций, а
решение задачи робастной стабилизации гарантирует устойчивость замкнутой системы. Отметим, что в данном методе не требуется информации о фазовой характеристике номинального объекта О. Кроме того, значение запаса робастной устойчивости является единственным индикатором, характеризующим степень достижения задаваемых требований к замкнутой системе.
Каждый из кратко описанных методов (корневой и частотный) синтеза регулятора обладает своими достоинствами и недостатками. Корневой метод по расположению плюсов передаточной функции замкнутой системы в заданной области на комплексной плоскости очень удобен для задания требований к качеству переходных процессов (к быстродействию и запасу устойчивости). Кроме того, он позволяет в явном виде задавать требования к замкнутой системе в виде ограничения на параметр затухания (коэффициент демпфирования). Однако вопросы чувствительности к возмущениям и шумам, а также точностные характеристики требуют дополнительного решения.
Н^-метод формирования контура позволяет в частотной области задавать требования к быстродействию, точности слежения, уровню подавления шумов и возмущений. Требование к качеству переходных процессов в виде ограничения на колебательность (параметр затухания) в явном виде не задается. В частотной области апериодичность переходного процесса определяется наклоном в —20 дБ на декаду графиков сингулярных чисел (в логарифмическом масштабе) в области частоты среза. В Н^-методе формирования контура удовлетворение ограничения на колебательность (близость к наклону в —20 дБ на декаду) может быть обеспечено только с помощью формирующих функций (см. п. 2.6), выбор которых представляет собой достаточно сложную задачу (особенно в классе формирующих функций низкого порядка).
Комбинированный подход, объединяющий Н^-метод формирования контура [17] и расположение полюсов замкнутой системы в заданной ЛМН-области [6], предложен в работе [22] для решения энергетической задачи. В ней с помощью аппарата ЛМН решается классическая задача робастной стабилизации для объекта, представленного в виде нормализованной левой взаимно простой факторизации.
2. ИСПОЛЬЗУЕМЫЕ В РАБОТЕ МАТЕМАТИЧЕСКИЕ ОПРЕДЕЛЕНИЯ И МЕТОДЫ
2.1. Стандартный объект
Произвольная многомерная замкнутая система с объектом управления ) и регулятором К(я) в положительной обратной связи (см. рис. 1) может
быть представлена в виде структурной схемы, изображенной на рис. 3. Объект управления (7(5) называется номинальным, а объект Р(5) — стандартным объектом.
Стандартный объект имеет два входа, а именно: обобщенный вектор всех входов м>, включающий в себя командный сигнал, внешние возмущения и шумы, и вектор управления и.
Выходами стандартного объекта является обобщенный вектор управляемых выходов г, обычно включающий в себя сигнал ошибки, и вектор измеряемых выходов у.
Стандартный объект задается уравнениями вида
Z( s)
_Y( s )
Pll( s ) Pi2(s)
P21( s 1 P22( s )
W( s) U( s)_
(2.1)
где 2(5), 7(5), Щя) и и(5) — преобразования Лапласа сигналов соответственно г, у, w и и, s — переменная оператора Лапласа.
Стандартный объект Р(5) обязательно включает в себя номинальный объект G(s) и всевозможные связи, а также, при необходимости, произвольные передаточные функции (см., например, работу [12]). Произвольные передаточные функции (например, префильтр и постфильтр) используются для формирования требований к замкнутой системе в частотной области. Обычно блок Р22(я) является номинальным объектом 6(5).
Замкнутая система управления определяется уравнениями вида
Д5) = РП(5)Ж(5) + Ри№(5),
7(5) = Р^)^) + Р22(*)и(5),
и(5) =К(я)7».
Передаточная функция замкнутой системы Тжг(Р, К) от обобщенного вектора входов w к обобщенному вектору контролируемых выходов г определяется выражением
ТШ(Р, К) = Рп + Рпщ + Р22к)-1Р21. (2.2)
Стандартный объект управления Р(5) описывается в пространстве состояний системой линейных обыкновенных дифференциальных уравнений
х = Ах + Biww + ^ z = C1x + D11w + D12u, y = C2x + D21w + D22u,
(2.3)
где х — вектор внутреннего состояния системы.
Связь между матричными передаточными функциями Р,-(5) стандартного объекта и представле-
]
нием в пространстве состояний стандартного объекта определяется уравнениями вида
Ру = С(5/ - Л)~1В] + Бу, /, ] = 1,2.
2.2. ЛМН-область
Описание выпуклой области на комплексной плоскости линейными матричными неравенствами приведено в работе [6].
Рассмотрим область Б левой комплексной полуплоскости С , которую можно определить следующим образом:
Б = {г е С- : Ь + гМ + 1ыт}, т
где Ь = Ь и М — действительные матрицы, г —
комплексная переменная, а г — комплексно сопряженная. Область Б, описываемая линейными матричными неравенствами, называется ЛМН-об-ластью.
_ т
Матричнозначная функция /д(г) = Ь + гМ + гМ называется характеристической функцией множества Б [6]. Значения матриц Ь и М определяют конкретный вид выпуклой области на комплексной плоскости.
В теории управления в аспекте расположения полюсов передаточной функции замкнутой системы интересны следующие виды элементарных ЛМН-областей.
• Б1(г) — полуплоскость Re(г) < — п:
/d1 (z) = z + z + 2n <
0.
Матрицы, задающие вид области, имеют вид Ь = 2п, М = 1. Характеристическая функция множества /д (г) определяет все точки г на комплексной плоскости, расположенные левее прямой Re(г)< —п.
• Б2(г) — сектор, образованный двумя лучами, выходящими из начала координат и симметрично расположенными относительно вещественной оси под некоторым углом ф:
/d2 (z) =
sinф(z + z) cosф(г- z) < о.
- cos ф( z - z) sin ф( z + z)_ Матрицы, задающие вид области, имеют вид:
L = 0, M =
sin ф cos ф -cos ф sin ф
Характеристическая функция множества fD (г)
определяет все точки г множества на комплексной плоскости, расположенные между двумя лучами.
• Б3(г) — круг с центром в точке (—q, 0) и радиусом г
fD3 (г) =
< 0.
-г q + г
_q + г -Г _
Матрицы, задающие вид области, имеют вид.
ь = -г q , м = 0 1
_q -г 0 1
Характеристическая функция множества /0 (г)
определяет все точки г множества на комплексной плоскости, расположенные внутри окружности с центром в точке (—q, 0) и радиусом г.
Основные свойства ЛМН-областей [6, 7].
• пересечение ЛМН-областей также может быть описано линейными матричными неравенствами, т. е. является ЛМН-областью;
• любая выпуклая область на комплексной плоскости, симметричная относительно действительной оси, может быть аппроксимирована ЛМН-областью с любой заданной точностью. Важной с точки зрения разработчика систем
управления представляет изображенная на рис. 4 ЛМН-область, образованная пересечением элементарных областей п п Б3, с характе-
ристической функцией вида
fDl 0 0
0 /Ъ2 0
0 0 fD„
(2.4)
Размещение полюсов замкнутой системы в указанной области гарантирует:
— затухание колебаний переходных процессов не менее п;
— колебательность не более ц = tgф или параметр затухания (коэффициент демпфирования) не менее £ = cosф;
— максимальную угловую частоту собственных незатухающих колебаний не более ю0 = а^ф.
Напомним, что задача размещения полюсов в заданной области может рассматриваться в совокупности с другими целями управления, такими как минимизация нормы передаточной функции замкнутой системы [6, 7].
2.3. Н.-норма передаточной функции
Для устойчивой (не имеющей полюсов при Re(s) 1 0) и правильной (порядок полинома числителя не более порядка полинома знаменателя) передаточной функции Ду) Д^-норма определяется как [23]
11Д»|и = sup{ а (Д»), ю е [0, да)},
где
ст (ю) = стп
» = тах{ 12 (П»Д»)} —
Рис. 4. Заданная ЛМН-область
максимальное сингулярное значение матрицы Д у'ю), X. — собственные числа произведения матриц, «звездочка» означает комплексное сопряжение.
Для скалярной функции Д^-норма равна максимальному значению амплитудной частотной характеристики. Класс функций с конечной Д^-нормой состоит из устойчивых физически реализуемых передаточных функций. Использование Д^-нор-мы в качестве критерия оптимальности при синтезе многомерных систем основано на том факте, что она может служить мерой усиления системы. Иначе, Д^-норма передаточной функции есть энергия выхода системы при подаче на вход сигнала с единичной энергией. Если выходом является ошибка, а входом — возмущение, то минимизируя Д^-норму передаточной функции, мы минимизируем энергию ошибки для наихудшего случая входного возмущения.
2.4. Задача ^-оптимизации
Задача нахождения Д^-оптимального регулятора формулируется следующим образом [23]:
для объекта управления G найти регулятор К такой, чтобы:
— замкнутая система была внутренне устойчива (все состояния системы асимптотически стремятся к нулевому из произвольного начального состояния при нулевых внешних входах ^);
— для Д^-нормы передаточной функции замкнутой системы (2.2) выполнялось равенство
||Р (Р *)||„ = уп
(2.5)
где уш1п — минимально возможное значение нормы для заданного объекта Р.
В работе [23] был предложен метод решения задачи Н^-оптимизации, так называемый «2-Рикка-ти подход». Его основная идея состоит в том, что задача нахождения регулятора, минимизирующего Н^-норму передаточной функции замкнутой системы от внешнего входа w к контролируемому выходу г, сводится к задаче поиска регулятора, который обеспечивает заданное значение у ее Н^-нормы, т. е. решается следующая субоптималь-ная задача:
найти регулятор К такой, чтобы выполнялось неравенство
1Тг(Р *)Ц
<
Y,
(2.6)
где у > уш1п и К — допустимый регулятор.
Эта задача решается в пространстве состояний. Вычисление Н^-субоптимального регулятора основывается на решении двух обобщенных алгебраических уравнений Риккати. Для решения задачи минимизации (2.5) применяется итерационная процедура по у. На каждом шаге решается субоп-тимальная задача (2.6). На первом шаге у задается заведомо больше у^п. Затем значение у уменьшается, и субоптимальная задача решается до тех пор, пока существуют неотрицательно определенные решения алгебраических уравнений Рикка-ти—Лурье. Полученное в результате итерационной процедуры минимальное значение у, «близкое» к уш1п с заданной степенью точности, а также решения обобщенных алгебраических уравнений Рик-кати—Лурье однозначно определяют робастный Н^-субоптимальный регулятор.
В работе [6] представлен алгоритм построения Н^-субоптимальных регуляторов, альтернативный «2-Риккати подходу». Обобщенные алгебраические уравнения Риккати—Лурье по управлению и фильтрации [23] заменяются неравенствами Риккати— Лурье, и множество решений данных неравенств используется для параметризации всех Н^-субоптимальных регуляторов, включая регуляторы пониженного порядка. Синтез регулятора производится при минимизации значения у на ос-
нове множества пар симметричных матриц Я и ^, удовлетворяющих системе матричных неравенств
AR + RAJ RC
С1R В,
1
—/1
т
11
В
Б
11
Б
< 0, (2.7)
г _ т Ат5 + 5А в 1 ст -
N 0 0
0 I вТ? БТ1 0 I
_ С1 Біі 1
< 0, (2.8)
(2.9)
где 0 и I соответственно нулевая и единичная матрицы соответствующей размерности, а матрицы и определяют базис нуль-пространств (ядро) т т
соответственно (В2 Б12) и (С2 В21).
Матрицы А, В1, В2, С1, С2, Б11, Б12, Б21 и Б22 определяются (см. п. 2.1) стандартным объектом (2.3).
2.5. Задача робастной стабилизации
Более общим описанием неопределенности, которое позволяет как нулям, так и полюсам передаточной функции находиться в правой полуплоскости, является взаимно простое описание неопределенности [20]. Нормализованная левая взаимно простая факторизация [21] передаточной функции 6(5) модели объекта означает представление в виде 6(5) = М(5) 1Л(^), если матричные передаточные функции (факторы) М(5) и N(5) принадлежат множеству устойчивых действительных рациональных матричных передаточных функций и М(я)М *(5) + Л(5)Л *(5) = I для всех 5, находящихся на мнимой оси.
Задача робастного проектирования состоит в том, чтобы одним регулятором стабилизировать, т. е. обеспечить устойчивость замкнутой системы, не только номинальный объект 6(5), но и все семейство возмущенных объектов с ограниченными по Н^-норме возмущениями.
Множество возмущенных объектов с ограниченным по Н^-норме возмущениями АМ(я) и АЛ(я) задается следующим образом:
ПЕ = {(М + ДМ)
-1
(Л + ДЛ) : ||[ДМ, ДЛ]|1 < в}.
Таким образом, задача робастной стабилизации состоит в том, чтобы найти регулятор К в обратной связи, который стабилизирует для заданного значения неопределенности е все возмущенные объекты из множества возмущенных объектов О .
Следовательно, задача состоит в нахождении максимальной неопределенности, определяемой значением ешзх, при которой регулятор К обеспечит устойчивость замкнутой системы для всех возмущенных объектов. В рамках ^-оптимизации данная задача формулируется следующим образом: Найти стабилизирующий регулятор К такой, что
(I - вК) 1М 1 К( I - вК)-1 M 1
-1
(2.10)
где К выбирается из множества регуляторов, стабилизирующих объект управления <7.
В работе [21] показано, что для случая нормализованной левой взаимно простой факторизации модели объекта выполняется равенство Н^-норм двух передаточных функций, а именно:
(I - вК)1М 1 К( I - вК)-1М 1
(I - вК) 1 (I - вК) 1 в
К( I - вК) 1 К( I - вК)-1в
(2.11)
Из равенства норм вытекает принципиально важный вывод: Н^-оптимальный регулятор, полученный при решении задачи робастной стабилизации (2.10), также является Н^-оптимальным при минимизации Н^-нормы передаточной функции замкнутой системы вида
ТДв, К) =
(I - вК) 1 (I - вК) 1 в
К( I - вК) 1 К( I - вК)-1в
. (2.12)
Данное свойство используется в Н^-методе формирования контура, предложенном в работе [17].
2.6. Н¥-метод формирования контура
Метод формирования контура позволяет достичь компромисса между качеством переходных процессов и робастной устойчивостью замкнутой системы.
Кратко поясним его основные идеи. Рассмотрим многомерную замкнутую систему управления с объектом в и регулятором К = Кт, изображенную на рис. 5. Пусть передаточные функции ^ ^ и
W2 = !2, где ^ и !2 — единичные матрицы соответствующих размерностей. Тогда на замкнутую систему действуют два внешних входа: йх на входе объекта в и й2 на его выходе. Матричная передаточная функция замкнутой системы Т от вектора
Рис. 5. Замкнутая система управления с формирующими функциями
внешних входов w = [й1, йу к вектору выходов г = [г1, г2]т имеет вид (2.12).
Рассмотрим матричную передаточную функцию Т . Малые значения максимальных сингу-
wZ *
лярных чисел передаточных функций (I — БК) 1 и (I — БК )16 обеспечивают требуемое качество функционирования замкнутой системы (подавление возмущений, уменьшение ошибки слежения). Для обеспечения робастности системы необходимо обеспечивать малые значения максимальных сингулярных чисел передаточных функций К(1 — СК) 1
и 6К(1 — (Ж)16 (см. работу [17]). Как уже отмечалось, принципиальный момент Н^-метода формирования контура заключается в том, что максимальные сингулярные числа матричных передаточных функций замкнутой системы могут быть аппроксимированы в соответствующих частотных диапазонах максимальными или минимальными сингулярными числами соответствующих матричных передаточных функций разомкнутой системы и регулятора (см. рис. 2).
На низких частотах ю < юх, при которых ст( 6К) . 1 имеем следующие аппроксимации
4-1
а ((I - вК) 1) * 1/а(вК),
ст ((I - 6К) хв) * 1/ст(К). (2.13)
Для случая высоких частот ю > юи, при которых ст (6К) П 1, имеем аппроксимации вида
а (Щ - вК) 1) * а (вК)
а (Щ -
вК) 1в)
а (к).
(2.14)
Обеспечение требований к качеству переходных процессов (а(вК) и а(К) большие) на низких частотах вступает в конфликт с требованием робастной устойчивости (а (К) и а (вК) малые) на высоких частотах. Компромисс достигается благодаря тому, что требования к качеству переходных процессов (а( вК) и а(К) большие) существенны на низких частотах, а требования робастной устой-
чивости (ст (К) и ст (СХ) малые) важны на высоких частотах.
Далее рассмотрим замкнутую систему (см. рис. 5) при произвольных передаточных функциях ^ и №2. Аппроксимации максимальных сингулярных чисел передаточных функций (2.13) и (2.14) будут соответственно иметь вид
ч-1ч
а ((I - W2GW1KJ-1) * l/<5(W2GWlK!т), а ((I - W2GWlK„)-1W2GWl) * 1/а(Х),
а (Кю(і - ж2єжхКт)~х) * а ^Є^Х),
а (Х(1 - Ж2ЄЖХК^)~ХЖ2ЄЖХ) * а (К) (2.15)
Таким образом, требования к максимальным сингулярным числам замкнутой системы могут быть сформулированы в терминах сингулярных чисел разомкнутой системы.
В соответствии с этим методом разработчик формирует новый обобщенный объект с префиль-тром W1(s) и постфильтром W2(s), как показано на рис. 5. Это позволяет задать желаемое поведение сингулярных чисел разомкнутой системы в соответствии с основными требованиями формирования разомкнутой системы (высокий коэффициент усиления на низких частотах, малый коэффициент усиления на высоких частотах и плавный переход в области частоты среза контура). Требования к сингулярным числам разомкнутой системы задаются в виде минимальных и максимальных сингулярных чисел обобщенного объекта на всем диапазоне частот. Таким образом, графики сингулярных чисел обобщенного объекта задают желаемую форму разомкнутой системы.
Заданная желаемая форма контура достигается, если Но-норма матрицы передаточных функций замкнутой системы от внешних входов йх и ё2 к выходам z1 и і2 минимизируется по всем стабилизирующим регуляторам Ко так, чтобы получить минимальное значение у:
У
Т Л'
^2
(I- GSK„ )-1 (I- GsК» )-1 Gs
Кт(I- Є5Ка})-1 Ксю(I- GsКт)-Ч
і ■
'орі’
где Є = Ж2ЄЖ1 — обобщенный («взвешенный») объект, Ко — регулятор. Величина, обратная у, называемая запасом робастной устойчивости є, принимает значения между 0 и 1. Чем больше є, тем больше запас робастной устойчивости. В работе
[17] показано, что величина ешзх также является индикатором совместимости требований к замкнутой системе.
Если значение ешзх мало, то это значит, что требование к качеству переходных процессов и требование робастной устойчивости не могут выполняться одновременно. В этом случае разработчик должен выбрать новые префильтр №^) и постфильтр №2^) для задания новой формы обобщенного объекта.
Процедура формирования контура состоит из следующих шагов.
Шаг 1. Формирование контура. Выбираются префильтр №^) и (или) постфильтр №2^) номинального объекта С для формирования требований к сингулярным значениям желаемой разомкнутой системы. Новый обобщенный («формирующий») объект представляется в виде О = Ж2С^Ж1. Предполагается, что W1 и W1 такие, что не содержит скрытых неустойчивых мод.
Шаг 2. Оценка совместимости требований к формированию контура. Решается задача робастной стабилизации для объекта О и вычисляется максимальный запас робастной устойчивости ешзх. Для этого минимизируется Н^-норма передаточной функции замкнутой системы 7 по всем стабилизирующим регуляторам Кт для получения оптимального значения уор1 в виде
Уорі єтах
(I- ад-1 (I- )-Ч
кт(I- о,Кт)-1 Кт(I- )-1 Є
Далее:
— если є тах П 1 (обычно єтах < 0,25), то возвращаемся к шагу 1 и выбираем новые W1(s) и Ж2(я);
— иначе выбираем є < єтах; при выбранном є синтезируем регулятор Ко, который робастно стабилизирует объект Є,, с областью робастной устойчивости є, применяя подход, описанный в работе [21].
Шаг 3. Синтез искомого регулятора. Окончательный регулятор К представляется комбинацией регулятора Ко с префильтром Ж1(я) и постфильтром Ж2(я) в виде К = Ж1Ко Ж2.
Теоретический базис Но-метода формирования контура состоит в том, что регулятор Ко не модифицирует желаемую форму контура существенным образом на низких и высоких частотах, если достигнутое значение єтах слитком мало [24].
30
о
Решающую роль в синтезе играет процедура нахождения формирующих функций (префильтра №^) и постфильтра №2^)). Обобщенный объект 6,, = Ж2^Ж1 задает желаемое поведение сингулярных чисел (требуемое разработчиком) разомкнутой системы во всем диапазоне частот. Значение уор1 характеризует меру близости (точность приближения) между сингулярным числом желаемой разомкнутой системы 6 = Ж26Ж1 (в виде обобщенного объекта) и полученной в результате синтеза Кт регулятора передаточной функции разомкнутой системы Ь = Ж26Ж1Кю Соотношение между ними имеет вид
°(ВД = а^С^К.) = а^б^) ± Уор( (дБ).
Таким образом, на графиках соответствующие сингулярные значения синтезированной разомкнутой системы находятся в некоторой трубке диаметром 2уор1 относительно соответствующих сингулярных значений желаемого обобщенного объекта.
ЗАКЛЮЧЕНИЕ
В первой части статьи приведены основные теоретические сведения по синтезу робастного Н^-регулятора с применением Н^-метода формирования контура и метода расположения полюсов замкнутой системы в заданной выпуклой области на комплексной плоскости. Оба метода основаны на аппарате линейных матричных неравенств. Во второй части будет приведена постановка и предложен алгоритм решения задачи синтеза регулятора для подавления возмущений в энергетической системе в аварийном режиме.
ЛИТЕРАТУРА
1. Кочкин В.И., Нечаев О.П. Применение статических компенсаторов реактивной мощности в электрических сетях энергосистем и предприятий. — М.: Изд-во НЦ ЭНАС, 2000.
2. FACTS: Modelling and Simulation in Power Networks / Tacha E., et al. — N.-Y.: John Wiley and Sons, 2004.
3. Flexible AC Transmission Systems (FACTS) / Ed. Yong Hua Song and Allan T. Johns. — London, IEE, 1999.
4. Hingorani N.G., Gyugyi L. Understanding FACTS: Concepts and Technology of Flexible AC Transmission Systems. — N.-Y.: IEEE Press, 2000.
5. Бесекерский В.А., Попов Е.П. Теория систем автоматического регулирования. — М.: Наука, 1972.
6. Chilali M., and Gahinet P. H Design with Pole Placement Constraints: An LMI Approach // IEEE Trans. AC. — 1996. — Vol. 41, N 3. — P. 358—367.
7. Chilali M., Gahinet P., Apkarian P. Robust Pole Placement in LMI Regions // IEEE Trans. AC. — 1999. — Vol. 44, N 12. — P. 2257—2270.
8. Linear Matrix Inequalities in Systems and Control Theory / Boyd, S., et al. — Philadelphia: SIAM books, 1994.
9. Баландин Д.В., Коган М.М. Синтез законов управления на основе линейных матричных неравенств. — М.: Физмат-лит, 2007.
10. Scherer C., Gahinet P., Chilali M. Multiobjective Output-Feed-back Control via LMI Optimization // IEEE Trans. AC. — 1997. — Vol. 42, N 7. — P. 896—911.
11. Zames G. Feedback and optimal sensitivity: model reference transformations, multiplicative seminorms, and approximate inverses // IEEE Trans. AC. — 1981. — Vol. AC-26, N 2. — P. 301—320.
12. Francis B.A., A Course in H Control Theory, Lecture Notes in Control and Information Sciences. — N.-Y.: Springer-Verlag. 1987.
13. Essays on Control: Perspectives in the Theory and its Applications. Progress in System and Control Theory / Ed. H.L. Tren-telman, J.C. Willems. — Birkhauser, 1993. — Vol. 14.
14. Postletweit L, Tsai M., Gu D. Weighting function selecting in H design // 11th IFAC World congress, Tallinn, Estonia, August 13—17. — 1990. — Vol. 5. — P. 104—109.
15. Gjerrit Meinsma. H control with Unstable and Nonproper Weights // Proc/ of 3rd European Control Conference, Roma, Italy, September 1995.
16. Doyle J.C, Stein G., Multivariable Feedback Design: Concepts for a Classical / Modern Synthesis // IEEE Trans. AC. — 1981. — Vol. AC-26, N 1. — P. 4—16
17. McFarlane and Glover K. A Loop Shaping Design Procedure Using H Synthesis // IEEE Trans. AC. — 37(6). — P. 759—769.
18. Papageorgiou G. Robust Control System Design: H Loop Shaping and Aerospace Applications // PhD thesis, University of Cambridge, Department of Engineering, July 1998.
19. Skogestad S. and Postlethwaite I. Multivariable Feedback Control. — London: John Wiley and Sons, 2001.
20. Vidyasagar M. Control System Synthesis: A Coprime Factorization Approach. — USA, Cambridge: MIT Press, 1985.
21. McFarlane D.C. and Glover K. Robust Controller Design Using Normalized Coprime Factorization Description. — N.-Y.: Springer-Verlag, 1990.
22. Bikash Pal, Balarko Chaudhuri, Robust Control in Power Systems. — N.-Y.: Springer Science + Business Media, Inc., 2005. — P. 140—150.
23. State-space solutions to standard H2 and H control problems / J.C. Doyle et al. // IEEE Trans. AC. — Vol.°°34, N 8. — 1989. — P. 831—847.
24. Zhou K., Doyle J.C, Glover K. Robust and Optimal Control. — New Jersey: Prentice-Hall, 1996.
Статья представлена к публикации членом редколлегии
С.Д. Земляковым.
Курдюков Александр Петрович — д-р техн. наук, зав. лабораторией, S (495) 334-92-61, e-mail: [email protected],
Тимин Виктор Николаевич — канд. техн. наук, ст. науч. сотрудник, S (495) 334-92-41, e-mail: [email protected],
Институт проблем управления им. В.А. Трапезникова, г. Москва.
