УДК 681.51.013
СИНТЕЗ РЕГУЛЯТОРОВ ПОНИЖЕННОГО ПОРЯДКА ПО ЗАДАННОМУ РАСПОЛОЖЕНИЮ ПОЛЮСОВ ЗАМКНУТОЙ СИСТЕМЫ
О.С. Вадутов
Томский политехнический университет E-mail: vos@ido.tpu.ru
Предлагается алгоритм синтеза регуляторов, обеспечивающих расположение доминирующих полюсов замкнутой системы в заданных точках, а недоминирующих полюсов ~ в заданной области. В основудвухэтапного алгоритма положены метод D-разбие-ния, модифицированный условиями на размещение полюсов системы, и метод поиска наилучшего решения в области параметров, в которой гарантируется желаемое размещение полюсов. Приведен пример.
1. Введение
Актуальность проблемы синтеза динамических регуляторов пониженного порядка объясняется сложностью и практической нецелесообразностью реализации регуляторов полного порядка. Многочисленные примеры показывают, что при помощи регуляторов пониженного порядка можно обеспечить требуемое качество процессов управления [1]. Проблема эта не имеет однозначного решения, и в литературе описаны различные методы синтеза.
В частности, имеется ряд работ, в которых задача синтеза регулятора пониженного порядка сведена к назначению доминирующих полюсов замкнутой системы [2-4]. Число варьируемых параметров динамического регулятора принимается равным числу назначаемых доминирующих полюсов. Значения варьируемых параметров регулятора вычисляются однозначно. Предлагаемые методы опираются на известный факт, что динамические свойства замкнутой системы определяются двумя-тремя полюсами, для которых выполняются условия доминирования. Рекомендации по назначению доминирующих полюсов можно найти, например, в [5]. Возможности этих методов ограничены, поскольку допускается произвольное размещение недоминирующих полюсов системы. Влияние этих полюсов на свойства системы приходится оценивать только на заключительной стадии синтеза регулятора.
В данной статье рассматривается алгоритм синтеза регуляторов пониженного порядка, в основу которого положен метод 1)-разбиения, модифицированный условиями на размещение полюсов системы [6]. Число варьируемых параметров регулятора, в отличие от [2-4], превышает число назначаемых доминирующих полюсов, а на недоминирующие полюсы замкнутой системы наложено ограничение в виде заданной области их расположения. Предлагаемый алгоритм синтеза состоит из двух этапов. На первом этапе при помощи метода 1)-разбиения выделяется та часть пространства параметров регулятора, в которой выполняются условия, наложенные на расположение доминирующих и недоминирующих полюсов системы. Для выбора конкретных значений параметров регулятора в полученной области, осуществляемого на втором этапе, предлагается использовать частотные и инте-
гральные оценки, характеризующие качество процесса управления.
В последнее время метод 1)-разбиения, используемый в данной работе, привлек внимание исследователей. Например, в [1] на его основе был разработан двухэтапный алгоритм синтеза регуляторов. Однако при этом число варьируемых параметров регулятора ограничено двумя, и в плоскости параметров регулятора строится область устойчивости. Оптимальные параметры регулятора определяются по критерию максимальной робастности с помощью численных процедур.
2. Постановка задачи
Рассматривается линейная стационарная система автоматического управления, операторно-структурная схема которой показана на рис. 1.
Л1
m
A(s)
e )—» H(s) u —► B{s)
D(s) A(s)
У
Рис. 1. Операторно-структурная схема системы
Объект управления описывается уравнением A(p)y(t)=B(p)u(t)-C(p)z(t),
где p=d/dt - оператор дифференцирования (формально эквивалентный оператору s преобразования Лапласа); полиномы А(р), В{р) и С{р) являются взаимно простыми, их степени удовлетворяют условиям: degA(p)=n, degB(p)<n, deg С(р)<п.
Уравнение динамического регулятора имеет вид D(p)u(t)=H(p)[g(t)-y(t)l
где àtgH(p)=àtgD(p)=m<n.
Синтез регулятора сводится к выбору его параметров, являющихся коэффициентами полиномов:
D(P) = dmPm + ... + dlP + d0, (1)
H(p)=hmpm + ...+hlP + h0.
(2)
А
А / / /
/ /
/ / / % Яе
/ / У А •х2
а
\
/ / /
Б ;
<
/
•а2
Рис. 2. Варианты размещения полюсов системы
Яе
> Р \ \
| > Яе
Часть параметров регулятора может быть задана заранее или удовлетворять дополнительным условиям. В частности, для статического регулятора целесообразно принять ¿0=1. В случае регулятора с ас-татизмом первого порядка ¿0=0, с11=1 и т. д. С учетом особенностей решаемой задачи могут быть назначены значения и других параметров регулятора. Поэтому число параметров регулятора, подлежащих определению, удовлетворяет условию г<2(т+\).
Система автоматического управления описывается уравнением
[А(р)0(р)+В(р)Н(р)]-у(1) = = В(р)Щр)-8(1)-С(р)0(р)-1(1).
Задачу данной работы сформулируем следующим образом. Необходимо определить значения г варьируемых параметров регулятора, при которых / доминирующих полюсов замкнутой системы принимают предписанные значения Хр у=1,...,/, а остальные п+г—1 недоминирующие полюсы удовлетворяют некоторым условиям в виде неравенств, ограничивающих область Б их размещения на комплексной плоскости. На рис. 2 приведены примеры расположения доминирующих полюсов АьА2Л и возможные границы областей Б размещения недоминирующих полюсов системы.
Поскольку число назначаемых доминирующих полюсов системы меньше числа параметров регулятора, сформулированная задача имеет некоторое множество решений.
3. Вывод основных соотношений
Обозначим произвольным образом коэффициенты полиномов (1) и (2), представляющие собой варьируемые параметры регулятора, через кьк2,...,кг Параметры регулятора входят в характеристическое уравнение системы линейно. Запишем уравнение в следующем виде:
(3)
Варьируемые параметры кьк2,...,кг регулятора разобьем на две группы. К первой группе отнесем параметры, которые назовем свободными. Пусть они образуют вектор к=(кь...,к^, размерность которого равна д=г-1. Во вторую группу включим зависимые варьируемые параметры регулятора, значения которых рассчитываются после выбора сво-
бодных варьируемых параметров из условия, чтобы / полюсов системы приняли предписанные значения. Эти параметры объединим в вектор к=(кч+ь...,к)Т размерностью 1=г-д.
Получим основные соотношения для случая двух свободных варьируемых параметров (д=2). Для этого представим характеристическое уравнение (3) в виде
Xк, ■ (р) + £ VС, (Р) + О0(р) = 0. (4)
1=1 1=3
Подстановкар=кру-1,...,/в (4) дает /уравнений:
1=1 1 = 3
у = 1,...,/.
Эти уравнения связывают варьируемые параме тры кь 1=1,...,г с задаваемыми доминирующими по люсами Хр_/'=1,...,/.
Представим (5) в матричной форме:
(5)
(6)
где
<2п(А) =
■ед) ад)"
<212 (А) =
ад) ад)
ад) - ад) ад) ... ад)
'-аду
и, ад
.-аду
Границу области размещения недоминирующих полюсов системы зададим функцией
Х{]т) = а +5(о) + ](о,8{(о ) = 8 (-а ),со е (-оо,оо).
Для вывода уравнения границы 1)-разбиения на плоскости двух свободных параметров сделаем в (4) подстановкур=а+3(а)+уа) и преобразуем получившееся комплексное уравнение в систему двух вещественных уравнений. В матричной форме она будет иметь вид
021(®)-кс+с>22(®)-кз=:я2(®), (7)
021(а,ю) = (}22(а,ю) =
К2(а,ю)=
где
ЯеСДа,®) КеС2(а,ю) ЬпОДа,®) 1тС2(а,а>) КеС3(а,ю) ••• ЯеСДа,®) 1т03(а,ю) ■■■ 1тСг(а,со) -КеО0(а,ю) -1т С0(а,со)
В результате объединения (6) и (7) получим систему уравнений:
д11(А)кс+д12(А)к3=Н1(А);
021(а,ю) кс+022(а,ю) кз =К2(а,ю). (8)
Систему уравнений (8) можно рассматривать как параметрическое уравнение границы Б-раз-биения плоскости параметров кьк2, образующих вектор кс, когда значение со пробегает границу области Б, при дополнительном условии, что доминирующие полюсы Хр _/'=1,...,/ замкнутой системы принимают предписанные значения. Из первого уравнения системы (8) выразим вектор к, зависимых варьируемых параметров:
кз = <3[2 (А) • (А) - (А) • ; (А) • кс. (9)
Из второго уравнения системы (8) после подстановки (9) получим:
к^а,«))
кс(а,ю) =
к2{а,т)
= [021(а,ю)-д22(а,ю)-д-21(А)-011(А)Г1х х[К2(а,ю)-022(а,ю)-0п1(А)-К1(А)] .
Для каждого конкретного значения на основании (10) рассчитывается числовое значение вектора кс, определяющее одну точку границы 1)-разбиения на плоскости свободных параметров кг и к2 регулятора. Изменяя в (10) со от -со до оо, можно построить всю кривую 1)-разбиения на этой плоскости.
Некоторые значения частоты дают при вычислениях неопределенности. Этим значениям со соответствуют уже не отдельные точки, а так называемые особые прямые. Первая особая прямая соответствует со=0. Для того чтобы получить ее уравнение, в (4) подставимр=а+8(со)+]со\(а=(1=а. Будем иметь
г
к1-С1(а) + к2-С2(а)+^кгС1(а)+С0(а)= 0.
1=3
Запишем это уравнение в векторно-матричной форме:
О21(а)-кс+С>22(а)кз=О0(а), (11)
где
021(а) = [ОДа) С2(а)Ъ 022(а) = [03(а)... Сг(а)].
И, наконец, подставив в (11) выражение (9), получим уравнение первой особой прямой в вектор-но-матричной форме:
ГО21 (а) - <322 (а) • (А) •011(А)]-кс = = С0(а) - ()22(а) ■ (А) • (А).
Вторая особая прямая, соответствующая ю=со, находится путем приравнивания к нулю коэффициента ап при слагаемом со старшей степенью характеристического уравнения (3), если этот коэффициент зависит от варьируемых параметров регулятора.
Штриховка границы 1)-разбиения производится в соответствии со знаком определителя
Д(ю) = С>21 (ю) - С>22 (ю) • <3~2(А) • С>п (А)
системы уравнений (9) на основании известных правил. Штриховка особых прямых, как обычно, выполняется согласно со штриховкой границы 1)-разбиения в окрестности точек сопряжения.
4. Алгоритм синтеза регуляторов
пониженного порядка
На основе полученных соотношений предлагается следующий порядок параметрического синтеза регуляторов пониженного порядка.
Этап 1. Формирование исходных данных для выполнения расчетов по предлагаемому методу синтеза. Исходя из априорных сведений о свойствах объекта управления определяются тип (статический, астатический) регулятора и его порядок. Множество рассчитываемых параметров регулятора, состоящее из коэффициентов передаточной функции, разбивается на свободные и зависимые. В качестве двух свободных параметров, в плоскости которых предполагается строить область 1)-разбиения, выбираются коэффициенты передаточной функции, оказывающие решающее влияние на интересующие свойства системы. Задаются значения Хр/=1,..^/доминирующих полюсов системы. Количество доминирующих полюсов системы должно быть равно числу зависимых параметров регулятора.
Этап 2. Построение области 1>-разбиения в плоскости свободных параметров регулятора. На этом этапе формируются матрицы 0И(А), (^(А), (^(ю), <322(ю), и векторы Ы^А), К2(ю). С помощью формул, полученных в разделе 3, в пространстве свободных параметров строятся граница 1)-разбиения и особые прямые. Для этого предлагается использовать современные специализированные системы программирования, имеющие в своем составе средства решения систем линейных алгебраических уравнений. К таким системам относится, например, МаШСАБ. Далее с помощью правил штриховки выделяется область, в которой обеспечивается заданное расположение доминирующих и других полюсов системы.
Этап 3. Поиск численных значений свободных варьируемых параметров из области 1>-разбиения. На этом этапе решается задача оптимизации по двум свободным параметрам регулятора. В качестве критерия оптимальности могут выбираться различные
показатели качества, характеризующие работу системы в установившихся и переходных режимах. Построенная область 1)-разбиения рассматривается как область допустимых решений задачи. Д ля поиска оптимального решения целесообразно использовать численные методы или методы моделирования с направленным перебором допустимых вариантов.
5. Пример синтеза регулятора
Рассмотрим систему стабилизации, оператор-но-структурная схема которой показана на рис. 3. Объект управления описывается уравнением {Юр +1 )(р +1)2 ХО = 1,2 • и(0 - ОД • гЦ).
В качестве регулятора используется ПИД-регулятор, в котором воздействие по производной формируется с помощью реального дифференцирующего звена непосредственно по выходной координате системы. В отличие от обычно принимаемого предположения, что постоянная времени дифференциатора Тл мала и ею можно пренебречь, будем считать эту постоянную времени четвертым параметром, подлежащим определению в процессе синтеза. Передаточная функция регулятора по управляемой величине у имеет вид
[/(5) Ос +кТ)я2+(к +к Т +к
Ш (сЛ — ^ ' — д п д _ п ид7_
1(5) +5
Найдем значения параметров кю кю кю Тл ПИД-регулятора, которые:
• обеспечивают расположение двух доминирующих полюсов замкнутой системы в точках А[=-0,2+Д 3 и А2=-0,2-ДЗ при условии, что остальные полюсы системы удовлетворяют неравенству
Г-0,2 + 0,5ю при юе (-оо,01,
Яе л\ < ■
[-0,2-0,5ю при юе[0, со);
• максимальное подавление возмущения г при выполнении заданных условий на расположение полюсов замкнутой системы. Качество подавления возмущения будем оценивать с помощью интеграла
со
о
где y¿t) - реакция системы на ступенчатое изменение возмущающего воздействия.
Введем следующие обозначения: кх=ТА, к2=кю к^=кп+кКТю к4=кл+кпТл и запишем характеристическое уравнение системы в виде, удовлетворяющем условиям применения описанного метода:
k1-G1(p) + k2-G2(p) + k,-G,(p)+
+k4-G4(p) + k5-G5(p)+G0(p) = 0,
где
Gl(p)=lOp5 + 21 р4 +12 р3 + р2;
G2(J>)= 1,21;
G,(p)=l,2lp;
G4(p)=l,2lp2;
G0(p)=lOp4 +21 p3 +12 p2 + p.
Параметры регулятора T=k{ и кя=к2 будем считать свободными, акгшк4- зависимыми. На основе полученных выше соотношений и исходных данных для рассматриваемой системы сформированы необходимые векторы и матрицы, проведено D-разбиение по двум свободным параметрам. На рис. 4 построена область в плоскости свободных параметров Тл и кя, соответствующая заданному расположению полюсов замкнутой системы.
В области D-разбиения нанесены линии равного значения критерия подавления возмущающего воздействия, выбранного для определения оптимальных значений параметров регулятора (рис. 4). Минимальное значение критерия достигается при значении Тл=0, то есть для идеализированного ПИД-регулятора. Выберем Гд=0,05 с. При этом максимально возможное значение к =2,2.
После подстановки этих значений в выражение, связывающее свободные и зависимые параметры регулятора, получены следующие значения зависимых параметров: кп=7,68; к=13,43. При найденных значениях параметров регулятора замкнутая система имеет полюсы:
Л =-0,2 + уОД ^=-0,2-70,3;
А, = -0,805 + Д179; Я4 = -0,805 - Д179;
Я5 =-20,09.
Рис. 3. Операторно-структурная схема синтезируемой системы
\
/_________ ----------
-0.2 0 0,2 0,4 0,6 0,8 1,0 Тя
Рис. 4. Область О-разбиения в плоскости свободных параметров регулятора
Заключение
Предложен алгоритм синтеза параметров регулятора пониженного порядка, обеспечивающего
заданное расположение доминирующих и недоминирующих полюсов замкнутой системы. В основу алгоритма положен метод построения границ 1)-разбиения с учетом ограничений на расположение доминирующих полюсов системы в заданных точках комплексной плоскости.
По существу предложенный алгоритм синтеза лежит в русле многокритериальных подходов к конструированию систем автоматического управления, получивших распространение в настоящее время. Он может служить основой для развития человеко-машинных (диалоговых) процедур проектирования систем автоматического управления с использованием, в частности, универсальной среды программирования МаШСАБ, имеющей в своем составе средства решения систем линейных алгебраических уравнений. Приведенный в статье пример подтверждает эффективность предложенного алгоритма.
СПИСОК ЛИТЕРАТУРЫ
1. КиселевО.Н., Поляк Б.Т. Синтез регуляторов низкого порядка по критерию ff° и по критерию максимальной робастности // Автоматика и телемеханика. - 1999. - № 3. - С. 119-130.
2. Скворцов JI.M. Синтез линейных систем методом полиномиальных уравнений // Известия АН СССР. Техническая кибернетика. - 1991. - № 6. - С. 54-59.
3. Скворцов JI.M. Интерполяционный метод решения задачи назначения доминирующих полюсов при синтезе одномерных регуляторов // Известия РАН. Теория и системы управления. -1996. - № 4. - С. 10-13.
4. Lim Choo Min, Flangovan S. Pole Assignment of a Class of System Using Output-Feedback Dynamic Compensators // IEEE Trans, of Circuits and Systems. - 1987. - V. 34. - № 2. - P. 199-201.
5. Райцын Т.М. Синтез систем автоматического управления методом направленных графов. - Л.: Энергия, 1970. - 96 с.
6. Вадутов О.С., Гайворонский С.А. Решение задачи размещения полюсов системы методом Л-разбиения // Известия РАН. Теория и системы управления. - 2004. - № 5. - С. 24-28.
7. Зотов М.Г. Многокритериальное конструирование систем автоматического управления. - М.: БИНОМ. Лаборатория знаний, 2004. - 375 с.
Поступила 28.09.2007г.