Параметрический синтез линейного регулятора интервальной системы с гарантированными корневыми показателями качества Текст научной статьи по специальности «Математика»

УДК 681.5

ПАРАМЕТРИЧЕСКИЙ СИНТЕЗ ЛИНЕЙНОГО РЕГУЛЯТОРА ИНТЕРВАЛЬНОЙ СИСТЕМЫ С ГАРАНТИРОВАННЫМИ КОРНЕВЫМИ ПОКАЗАТЕЛЯМИ КАЧЕСТВА

М.С. Суходоев, С.А. Гайворонский, C.B. Замятин

Томский политехнический университет E-mail: smike@aics.ru

Рассматривается система автоматического управления, содержащая пропорционально-интегральный регулятор и объект управления, который имеет интервально-заданные параметры. Используя робастное расширение метода корневого годографа, разработана методика синтеза параметров пропорционально-интегрального регулятора, гарантирующего минимально допустимую степень устойчивости и максимально допустимую колебательность системы. В основу методики положен вершинный анализ корневых показателей качества с применением уравнений Теодорчика-Эванса. Приводится числовой пример.

Введение

В реальных системах автоматического управления возможны случаи, когда некоторые их параметры не известны точно, либо меняются в процессе эксплуатации системы по заранее неизвестным законам, причем их значения не могут быть доступны измерению. Если известны диапазоны возможных значений постоянных параметров или пределы изменения нестабильных параметров, то говорят о параметрической интервальной неопределенности. Системы, имеющие интервально-нео-пределенные параметры, получили название интервальных систем.

При проектировании интервальной системы основная задача состоит в обеспечении желаемого качества её функционирования при любых возможных значениях интервально-неопределенных параметров. Введем в рассмотрение корневые показатели качества системы: степень устойчивости а и колебательность ср. Очевидно, что при нестабильности параметров системы эти показатели качества могут изменяться. Поэтому представляет интерес задача обеспечения в интервальной системе гарантированной минимальной степени устойчивости и максимальной колебательности.

Для задания желаемого качества системы, соответствующего этим корневым показателям, может быть использована ломаная линия АВСБ (рис. 1), задающая границу области локализации корней Г.

В

/

/

Re

Im —►

1. Постановка задачи

Рассмотрим систему автоматического управления, рис. 2.

2) * Линейный регулятор Объект управления

г

Рис. 2. Схема системы автоматического управления

Пусть линейный регулятор имеет передаточную функцию пропорционально-интегрального (ПИ) регулятора:

1¥Р(з) = Кп'° + К\ (1)

где Кп, Ки - настройки регулятора, а передаточная функция объекта управления имеет вид:

.Ж*)

Woy(s) =

B(s)

(2)

где A(s) = ]Г а. • s\ B(s) = X V sJ> bj ~ bj ~ bj-

1=0 j=0

Тогда интервальный характеристический полином системы может быть записан в виде:

+ = (3)

Область возможных значений интервально-нео-пределенных параметров системы (многогранник Р„ являющийся прямоугольным гиперпараллелепипедом) отображается на комплексную плоскость корней в виде областей локализации комплексно-сопряженных корней и отрезков вещественной оси, где локализуются вещественные корни (рис. 3).

Рис. 1. Область желаемого расположения корней

А

^^^ /

--

/ / ^^jlí Re

1 W

< а >

г

У —

V

Рис. 3. Отображение параметрического многогранника Р3

Необходимо определить настройки ПИ-регуля-тора, обеспечивающие расположение корневых областей интервального характеристического полинома в заданном секторе Г при любых значениях интервальных параметров объекта управления.

2. Задание минимальной степени устойчивости интервальной системы

Приведем интервальный полином (3) к виду

О, (4)

1=0

где коэффициенты с,- могут являться функциями интервальных параметров объекта управления и настроек регулятора. Интервальные коэффициенты с,, образуют параметрический многогранник, вершины которого определяются их крайними значениями.

Из теории корневого годографа известно, что если корень лежит на вещественной оси, то угол выхода корневого годографа из этого корня составляет 0° или 180°. Очевидно, что для того, чтобы вещественный корень интервального полинома двигался влево при изменении любых интервальных коэффициентов, необходимо, чтобы углы выхода по всем этим коэффициентам составляли 180°. Набор интервальных коэффициентов, который обеспечивает данное требование, определяется на основе следующего утверждения.

Утверждение

Если интервальные коэффициенты полинома заданы чередующимися пределами с^с^с^..., начиная с максимального для с0, то данный набор коэффициентов определяет вещественный корень, углы выхода корневых годографов из которого по всем интервальным коэффициентам составляют 180°.

Доказательство

Из основных фазовых соотношений теории корневого годографа известно, что угол выхода из корня при увеличении с,- можно найти по формуле:

®f= 180°-¿0^+200, (5)

р=1

где &р и 0О - углы между вещественной осью и векторами, направленными из корня соответственно к р-ому полюсу и к i нулям функции

Пусть полином (4) имеет правый вещественный корень S*. Тогда для любой пары комплексно-со-пряженных корней и для любого вещественного корня, лежащих левее S*, на основе свойств корневого годографа можно заключить, что 0Л+0Д=36О°, 0Я=О°, 0о=18О0 (рис. 4).

Следовательно, угол выхода из корня S* при увеличении с,- можно найти по формуле: 0;®=18О°+г©о, а при уменьшении, 0;?=/0о.

Таким образом, если интервальный полином будет иметь чередующиеся пределы коэффициентов с^с?!то углы выхода из образа соответ-

ствующей этому набору вершины параметрического многогранника Рл будут равны 180°.

полинома

На основании данного утверждения можно сделать вывод, что минимальная степень устойчивости интервальной системы может быть задана вертикальной прямой, проходящей через вещественный корень, соответствующий чередующимся пределам коэффициентов интервального полинома.

Пусть минимальная степень устойчивости задается корнем s*=a*. Введем информацию об этом корне в полином (4), задавая в нем чередующиеся значения Ь( и значение s *=a *. Выразив далее ^через Кп, получим выражение:

Ки=/(Кп,а*,Ъ]), (6)

где bj - граничные значения коэффициентов полинома B(s), обеспечивающие чередование пределов с;.

Для дальнейшего синтеза ПИ-регулятора получим интервальный характеристический полином с одним варьируемым коэффициентом

(Ки (Kn,a*,b])+Kn-s)-A(s) + s- B(s) = 0.

3. Локализация корней интервального

характеристического полинома в заданной области

на основе уравнения Теодорчика-Эванса

Известно [1], что в случае интервальной неопределенности характеристического полинома минимальная степень устойчивости и максимальная колебательность определяются образами некоторых вершин многогранника интервальных коэффициентов, то для решения поставленной выше задачи необходимо определить эти вершины. Для этого предлагается воспользоваться методикой, разработанной в [2].

После нахождения набора проверочных вершин необходимо на основе уравнений корневого годографа Теодорчика-Эванса [3] определить для каждой из них значения настраиваемого параметра регулятора Кп, при котором происходит пересечение корневого годографа с заданной на рис. 1 границей области локализации корней.

Представим полученный полином (9) в виде: D<j>) = 0(p)+Kn4>(p)=0.

Введем следующие обозначения:

Е(а,ю) = Яе(Ф(р)), Р(а,т)=ЯеС¥(р)),

Согласно [3], для определения искомых значений Кп, необходимо подставить уравнение границы области Г в (10), (11) и решить уравнение

^(ю)Р(ю)-£(ю)Д(ю) = 0.

(7)

Полученные действительные корни ю,- следует подставить в выражение:

Е(р)Р(р)+Р(со)К(со)

кп=-

(8)

Проделав данную процедуру для каждой проверочной вершины, можно получить соответствующие интервалы значений Кп, при которых корни интервального полинома лежат в заданной области Г. Для нахождения интервала значений Кп, удовлетворяющего всем вершинам, следует найти пересечение всех интервалов.

Для окончательного определения настроек ПИ-регулятора необходимо выбрать значение Кп из полученного интервала и подставить это значение в выражение (6).

4. Методика синтеза ПИ-регулятора

На основании проведенных исследований разработана методика параметрического синтеза ПИ-регулятора, включающая:

1. Задание требуемых показателей качества (минимально допустимой степени устойчивости и максимально допустимой колебательности интервальной системы).

2. Задание в характеристическом полиноме 5*=а * и пределов интервальных параметров объекта

управления, соответствующих вершине с координатами сё С[

3. Получение зависимости (6) и приведение характеристического полинома с двумя варьируемыми параметрами ПИ-регулятора к полиному с одним варьируемым параметром.

4. Нахождение возможных вершин многогранника Р„ отображающихся на границу области Г.

5. Определение для каждой найденной граничной вершины на основе выражений (7), (8) интервала значений варьируемого параметра Кп, при которых ветви корневого годографа находятся в области Г

6. Определение пересечения найденных интервалов Кп, при котором корни полинома располагаются в заданной области Г для всех граничных вершин.

7. Выбор значения Кп из области пересечения интервалов и определение значения Ки на основе выражения (6).

5. Пример синтеза

Пусть задана передаточная функция разомкнутой системы с единичной обратной связью:

К-{Т-5 +1)

кп

где Кр, Тр - настраиваемые параметры регулятора, Кц - постоянный коэффициент передачи объекта управления, 02, аь щ - интервальные параметры объекта управления, (я0е[О,О7;О,О8]; [0,3;0,4]; а0е[2;3]).

Необходимо определить настройки ПИ-регуля-тора, гарантирующие а*= 1 и <р=10°.

Приведем характеристическое уравнение данной системы к виду:

с3-53+с2-52 + с1-5+с0= 0, (9)

Область желаемой локализации корней Г

Ьп

Ие

-1,5

Рис. 5. Области локализации полюсов системы с синтезированным ПИ-регулятором

где с}=а2, сг=аъ с^щ+К^-Тр, с0=

Подставив в (9) соответствующие вершине адс^з... граничные значения коэффициентов аь получим зависимость:

т а2^ъ+а1-А'2+а0-А' + К0-Кр

На основе (10), характеристическое уравнение системы приведем к виду:

к0кр

а2а *3 +арс *2 *+КдК

К0Кра*

*+К0Кр= 0.

После преобразований, получим: а2 -53 + ах -52 +а0 - 5- (а2-а *2+а1- а * 2+а „)• 5 +

'К» "^-1 = 0.

(И)

Составим на основе (11) уравнения Теодорчи-ка-Эванса вида (7), (8). В соответствии с методикой [2] для шести вершин многогранника интервальных коэффициентов найден общий интервал Кг, при котором корни полинома располагаются внутри области Г: ^е[2,4733;5,0160].

Выберем значение Кг из найденного интервала, ^,=5,0160, и на основании выражения (15) получим значение второго искомого параметра Гр=0,4458.

Области локализации корней интервального характеристического полинома с найденными настройками ПИ-регулятора представлены на рис. 5.

Из рис. 5 видно, что корневой годограф интервального характеристического полинома располагается в области Г, ограниченной минимальной степенью устойчивости и максимальной колебательностью, следовательно, синтезированные параметры линейного регулятора гарантируют заданные корневые показатели качества системы.

СПИСОК ЛИТЕРАТУРЫ

1. Гусев Ю.М., Ефанов В.Н., Крымский В.Г., Рутковский В.Ю. Анализ и синтез линейных интервальных динамических систем (состояние проблемы). Анализ с использованием интервальных характеристических полиномов // Известия АН СССР. Техническая кибернетика. - 1991. - № 1. - С. 3-23.

2. Гайворонский С.А. Вершинный анализ корневых показателей качества системы с интервальными параметрами // Известия

Томского политехнического университета. - 2006. - Т. 309. -№ 7. - С. 6-9.

3. Римский Г.В. Корневой метод решения задач устойчивости интервальных систем // Весщ АН Беларуси. Сер. физ.-техн. на-вук. - 1994. - № 4. - С. 53-61.

Поступила 28.09.2007г.