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

УДК 007 ББК 22.213

АЛГЕБРАИЧЕСКОЕ ДОСТАТОЧНОЕ УСЛОВИЕ УСТОЙЧИВОСТИ ПЕРИОДИЧЕСКИ НЕСТАЦИОНАРНЫХ СИСТЕМ УПРАВЛЕНИЯ

Чечурин Л. С.1

(Государственное образовательное учреждение высшего профессионального образования «Санкт-Петербургский государственный политехнический университет»)

Предлагаются достаточные алгебраические условия устойчивости систем управления c периодически нестационарным коэффициентом усиления обратной связи. Результат получен на основе теоремы Харитонова для интервальных полиномов и достаточного критерия Бонджиорно для нестационарных систем. Эффективность метода иллюстрируется примером.

Ключевые слова: линейные периодически нестационарные системы, устойчивость, теорема Харитонова

1. Введение

Один из важных разделов теории робастных систем управления относится к анализу линейных стационарных систем, параметры которых неизвестны, но принадлежат определенным интервалам. Если описание такой системы получено в форме характеристического полинома D(p) с интервальными коэффициентами, то вопрос о робастности характеристического полинома имеет решение в форме необходимого и достаточного критерия устойчивости - теоремы Харитонова. Робастность

I Леонид Сергеевич Чечурин, кандидат технических наук, доцент (cepreu4@gmail.com).

24

проверяется гурвицевостью четырех определенным образом составленных полиномов [3]. На комплексной плоскости при частоте о* каждому такому полиному соответствует угловая точка прямоугольника семейства возмущенных годографов D(ja)*). Далее критерий Михайлова гарантирует, что устойчивость всех четырех полиномов будет означать и устойчивость всего семейства. Поскольку устойчивость можно проверить алгебраическими критериями, мы получаем алгебраический критерий робастной устойчивости интервальных стационарных систем. Тот факт, что из четырех годографов ближайшим к началу координат будет только один, позволяет после дополнительных построений и вовсе свести проверку робастности к анализу одного годографа (критерий Цыпкина-Поляка). Отметим также обобщение анализа В. Л. Харитонова на случай полиномов с комплексными интервальными коэффициентами (см., например, [6]). Несмотря на достигнутые значительные результаты, направление исследования интервальных систем продолжает оставаться в центре внимания, поскольку содержит ряд не имеющих простого решения задач. Современный и детальный обзор состояния дел в этой области можно найти в [2].

С инженерной точки зрения кажется полезным иметь возможность оценить робастность в случае, когда у интервальной системы снимается предположение о стационарности. Ограничимся случаем, когда коэффициент усиления обратной связи может изменяться внутри указанного интервала по некоторому неизвестному периодическому закону (см. рис. 1, слева). Системы рассматриваемого класса достаточно широко распространены, кроме того, к такому описанию сводятся и некоторые случаи описания систем в форме дифференциальных уравнений с периодическими нестационарными коэффициентами. Обычно это означает переход к качественно другому классу систем - периодически нестационарным, где анализ устойчивости проводится специфическими и достаточно сложными инструментами [1].

В работе показывается, как на основании достаточного критерия устойчивости Бонджиорно можно построить простой

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

Опубликованная в [5] идея существенно доработана, снабжена комментариями, примером расчета границы устойчивости и сравнением результата с полученным из известного критерия и из численного эксперимента.

2. Построение критерия

Рассмотрим систему управления, описываемую уравнением в операторной форме:

(1) Мр)х(0 + М(р)[а(0х(0] = 0,

где Мр), М(р) - полиномы оператора дифференцирования р, а(0 - периодически нестационарный параметр, т.е.

М(р) = Ьрп + Ъ„_1РП-1 + ... + Ьр + Ьо,

(2) М(р)=еп_1рп 1 + СП_2РП 2 + ... + С1Р + Со, а(0 = а0 + а~(0,

где а0 - постоянная составляющая (среднее значение) параметра, а~(0 - ограниченная центрированная периодическая составляющая, т.е. | а~(0 | < а.

Достаточные условия устойчивости системы (1) представленной в структурной форме с обратной связью на рис. 1 слева, где Щр) = М(р)/Мр), дает частотный критерий Бонджиорно [7], имеющий ясный геометрический смысл. Для устойчивости (1) достаточно того, чтобы при устойчивости стационарной системы (устойчивости Щ(/'о)/(1 + а0Щ(/о)) годограф обратной частотной амплитудно-фазовой характеристики Щ1(/ю) не попадал внутрь окружности радиуса а с центром в точке (-а0, 0) (см. рис. 1, справа). Аналитическое условие этого, очевидно, имеет вид

(3) | а0 + Щ'Цо) | = | а0 + [Що)/М^о)] | > а.

Рис. 1. Структурная схема нестационарной системы управления (слева) и иллюстрация критерия Бонджиорно (справа)

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

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

(4) Щ ~1(/'ю) = -(а0 + ае^) = -(а0 + аео8 р +jasinф).

Это равносильно тому, что с учетом обозначений (2) годограф D(/'m) = N(/0) + М(/0)[а0 + асо$р+/а8тр] при некотором р пересекает начало координат. В таком случае по крайней мере один из четырех годографов с интервальными комплексными коэффициентами

D(ja') = N(/0) + М(/о)[а0 ± а ± /а] не охватывает начало координат (обходит начало координат слева). Для данного случая мы можем связать этот факт с потерей устойчивости одного из полиномов семейства

(5) D(p) = Фпрп + <3п-1рп- + . + йхр + ф = 0, коэффициенты которого имеют вид ф = Ь + (а0 + асоър + +/аътрс, и их действительные Яе Ф и мнимые 1т ф части ограничены интервалами [а, Д] и [1, т] соответственно, причем

а, = Ь, + (а0 + а)сг,

Дг = Ьг + (а0 - а)с,,

1 = а, т = -а.

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

Утверждение. Для устойчивости колебаний периодически нестационарной системы (1) достаточно робастной параметрической (интервальной) устойчивости по р системы, описываемой характеристическим уравнением (5), устойчивость которой проверяется гурвицевостью восьми полиномов

А+(р) = (а + А) + (а + /т)р + (Д + /'т)р2 + (Д + А)р3 +

+ (а + /14)р4 + (а5 + /т)р5 + .,

(6) А+(р) = (а + /т) + (Д + /т)р + (Д + /1)р2 + (аз + /1з)р3 +

+ («4 + />4)р4 + (Д + /т)р5 + .,

А+(р) = (Д + Л)) + (а1 + /1)р + (а2 + /т2)р2 + (Дз + /тз)р3 +

+ (Д4 + /14)р4 + (а5 + /1^)р5 + .,

А+р = (Д + /т) + (Д1 + /1)р + а + Л)р2 + (аз + /тз)р3 +

+ (Д4 + />4)р4 + (Д5 + /'^)р5 + . ,

А_(р) = (а> + А) + (Д1 + Л)р + (Д2 + /т)р2 + («з + /тз)р3 +

+ (а4 + /14)р4 + (Д5 + /1^)р5 + .,

А”(р) = (а> + /т) + (а1 + /1)р + (Д + /1)р2 + (Дз + /тз)р3 +

+ (а + /т)р4 + (аз + А)р5 + .,

А“(р) = (Д0 + Л)) + (Д + /А)р + (а2 + /тт)р2 + (аз + /1з)р3 +

+ (Д4 + /14)р4 + (Д5 + /тЗ)р5 + .,

А(р) = (Д0 + /т) + (а + /т)р + (а2 + Л)р2 + (Дз + /1з)р3 +

+ (Д4 + /'Цл)р4 + (аз + /тз)р5 + . .

3. Обсуждение

Геометрический смысл сформулированного утверждения лучше всего пояснить на частном случае, когда М(р) = 1. Тогда можно прямо указать возмущения стационарного характеристи-

ческого полинома, вызванные снятием условия стационарности одного из его коэффициентов. Пусть у полинома В(р) свободный член является интервальным стационарным

(7) В(р, а)=Фпрп + ... + Ф1р + а = 0,

а < а < а.

Устойчивость интервального полинома определятся по теореме Харитонова, однако из-за того, что неопределенный параметр только один, возможно непосредственно построить семейство годографов Михайлова. Для частоты о сечение семейства годографов есть отрезок между Б(/о, а) и Б(/о, а), параллельный действительной оси. Эта ситуация приведена на рис. 2, где а0 - центр интервала [а, а].

Рис. 2. Построение семейства возмущенных годографов Михайлова для одного интервального коэффициента

Предположим теперь, что а = а(0 = а0 + а~(0. Отнесем а0 к стационарной части полинома, обозначим его 0(р,а0) и применим критерий Бонджиорно в такой формулировке: для устойчивости (7) достаточно, чтобы на произвольной частоте о окружность с центром 0(/о, а0) и радиусом а не охватывала начало координат (рис. 2). Если вместо окружности использовать описанный около нее квадрат с параллельными осям сторонами, то условия критерия усиливаются, однако мы получаем возможность оценить робастность по теореме Харитонова. Нетрудно видеть, что в данном случае (один интервальный комплексный

коэффициент), условия устойчивости сводятся к гурвицевости не восьми, а лишь четырех полиномов А_4(р) (6). Рисунок явным образом иллюстрирует принципиальный факт: если интер-вальность коэффициента характеристического полинома для систем такого вида вызывает возмущения годографа вдоль одной лишь оси, то допущение к этому его нестационарности -уже вдоль двух. Иными словами, вопрос достаточной устойчивости системы, описываемой линейным дифференциальным уравнением с переменным коэффициентом, может быть сведен к вопросу устойчивости стационарного характеристического полинома с (неопределенным по фазе) комплексным коэффициентом.

В случае, когда М(р) является полиномом, критерий остается работоспособным, однако переход к анализу устойчивости с характеристическим полиномом О(р) уже не будет так нагляден, т.к. переменный коэффициент а(0 и его производные войдут в него сложным образом.

Приведенные математические манипуляции дают возможность лучше понять эффект от снятия предположения о стационарности интервального коэффициента усиления обратной связи с точки зрения теории управления. Если линейная динамическая система может быть представлена в форме системы с обратной связью, условие ее устойчивости может быть установлено критерием Найквиста, проверяющим наличие в системе баланса фаз и амплитуд (случай пересечения точки (-1, 0) на комплексной плоскости годографа Найквиста). Интервальность коэффициента усиления обратной связи приведет возмущению годографа разомкнутой системы, проявляющемуся в неравномерном переносе вдоль действительной оси. Нестационарность же этого коэффициента вызовет возмущение годографа уже в обоих направлениях, т.е. вносит в разомкнутую систему дополнительный сдвиг по фазе.

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

системы часто предсказуемы. Во-первых, это резонансные частоты В(/т, а0), точнее, резонансная частота, соответствующая преобладающему пику амплитудно-частотной характеристики \№(1а>)1(\ + а0Ж(/ш)\, или частота, определяющая \Щ^(/ю)І(1 + а0Ж(/ о))\\¥. Во-вторых, если частота изменения параметра а(ґ) задана и равна О, то это частоты первого ОІ2, второго О, третьего 3ОІ2 и т. д. параметрических резонансов системы. В таком случае имеет смысл привести формулировку упрощенного «точечного» критерия устойчивости, сводящегося к проверке гурвицевости четырех полиномов с действительными коэффициентами.

Пусть левая часть (3) достигает минимума на частоте ор. Что означает, что на этой частоте годограф а0 + ^г(/О) ближе всего к началу координат, или, что то же самое, на годографе W(/ю)І(1 + а0Щ/'ю)) различим резонанс. Таким образом, для выполнения (3) достаточно его справедливости для частоты сор.

Рассмотрим колебания на интересующей нас частоте ор . Используем, следуя [4], замену р = /ор и / = рІор , тогда равенство (4) переписывается в виде

W _1(р) + [а0 + асо8р + (арІюр)8Іи р] = 0, откуда с учетом обозначений (2) следует характеристическое уравнение

(8) О(р) = Ы(р) + М(р)[а0 + асо$р+ (арІюр)8Іпр] =

=dn рп + dn_\pn 1 + ... + d\p + do = 0.

Коэффициенты характеристического полинома, в отличие от (5), действительны и определяются для і = 1, ., п

(9) di = Ьі + (а0 + асо$,р)еі + (Сі-1аІюр)8Іпр =

2 2 2 1І2 = Ьі + асСі + а(Сі + с- Ор ) 8Іп(р + аг^с/с^ Ор)

и d0 = Ь0 + (а0 + асо%р)с0. Таким образом, коэффициенты di ограничены и имеют максимальные di и минимальные ^ значения, причем для і = 1, ., п

(10) di = Ьі + а0Сі + а(с2 + Сi-12Юр-2)1І2,

& = Ьі + а0Сі - а(Сі2 + С_12Юр-2)1/2,

а для і = 0

d0 = Ь0 + а0с0 - ас0, й0 = Ь0 + а0с0 + ас0.

Фазу р можно считать неопределенностью системы. Таким образом, для устойчивости колебаний периодически нестационарной системы (1) достаточно робастной параметрической (интервальной) устойчивости по р системы, описываемой характеристическим уравнением (8), устойчивость которой проверяется гурвицевостью четырех полиномов

М(р) = ^0 + ^1р + $2р + d3p3 + ^4р4 + й?5р5 + . ,

(11) 02(р) = ^0 + ^р + $2р + ^3р3 + ^4р4 + ^5р5 + . ,

А(р) = ^0 + ^р + $_2р + d3p3 + d4p4 + ^5р5 + ... ,

А(р) = d0 + d1p + d2p2 + d3p3 + d4p4 + d5p5 + ... .

4. Пример

Проведем анализ устойчивости динамической системы, описанной характеристическим полиномом Б(р) = Ь2р2 + Ь1р + Ь, где Ь принадлежит интервалу [Ь, Ь]. Пусть Ь0 - центр этого интервала, а Л - его полудлина, т.е. Ь е [Ь0 - Л, Ь0 + Л].

Задача эквивалентна анализу устойчивости колебательного звена с обратной связью, имеющей интервальный коэффициент усиления рис. 3, где а е [-Л, Л].

Рис. 3. Структурная схема интервальной системы с обратной связью

Если исследуемая система стационарна, то интервальность коэффициента породит необходимость рассмотреть семейство характеристических полиномов или семейство годографов

Михайлова. В данном случае они строятся сдвигом невозмущенного годографа А(/ю, Ь0) вдоль действительной оси на значения от -Л до Л. Определим устойчивость системы. Формальное применение для этого теоремы Харитонова в данном случае потребует проверки гурвицевости лишь двух простых полиномов

Вз(р) = Ь0 + Л + Ьр + Ь2р2.

Теперь снимем условие стационарности коэффициента Ь, оставив прежними границы его изменения. В соответствии с критерием Бонджиорно в форме (3), для устойчивости достаточно выполнения на любой частоте условия

I А(/ю, Ь0) I = I Ь + Ь/а + Ь2(/'ю)21 > Л.

Минимум левой части достигается на (резонансной) частоте

С другой стороны, разрешимость (соответствующего случаю равенства) уравнения границы устойчивости, очевидно, могла бы быть истолкована как наличие такого р, при котором выполняется

А(/ю, р) = Ь0 + Ь/а + Ь2(/ю)2 + Лсо$р+/Л$тр= 0.

Иными словами, речь идет об устойчивости полинома с комплексным коэффициентом А(р, р) = Ьр2 + Ьр + Ь,

Ь = Ь0 + Лс08р + /Л8тр

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

А(р) = Ьо - А + Ьр + Ьр2

и

его значение для Л:

I А(/юр, Ь0) I дает достаточное граничное значение

зоваться теоремой Харитонова для интервальных полиномов с комплексными коэффициентами. Соответственно, для устойчивости достаточно гурвицевости четырех полиномов А(р) = Ь - Л -/Л + Ь1р + Ьр2,

А2(р) = Ь0 - Л + /Л + Ь1р + Ь^р2,

А3(р) = Ь + Л - /Л + Ь1р + Ьр2,

А4(р) = Ь + Л + /Л + Ьр + Ьр2.

Для рассматриваемого случая граница устойчивости по Л определится лишь одним полиномом А(р), из которого можно получить напрямую алгебраическое условие устойчивости

(13) -Ь1 ± Яе (Ь12 - 4Ь0Ь2 + 4ЛЬ2 - /4Ь2Л)1/2 < 0.

Наконец, получим «точечную» оценку для частоты юр. При этом годограф возмущенного полинома имеет вид (8):

А/ю, р) = Ь0 + Лоо8р + (Ь1 + Л 8тр/юр)/ю + Ь2(/ю)2 = 0. Последнее означает, что на некоторой частоте а возмущения годографа А(/а, р) пересекают начало координат. Для того, чтобы этого не произошло, в соответствии с критерием Михайлова, достаточно, чтобы при любом р был устойчив полином А(р, р) = d2p2 + dp + dо,

где

d2 = Ь2,

d1 = Ь1 + Л 8тр®р, d0 = Ь0 + Л 008р.

Коэффициенты d1 и do стационарны, принадлежат интервалам d0 е [Ь0 - Л, Ь0 + Л], d1 е [ Ь1 - Л/юр, Ь1 + Л/юр] и связаны через р. Для применения теоремы Харитонова будем считать их независимыми и проверим (более) достаточное условие устойчивости нестационарной системы, сформулированное, однако, в алгебраической форме: требуется гурвицевость четырех полиномов

ААр) = Ь - Л + (Ь1 - Л/ар)р + Ьр2,

А(р) = Ь - Л + (Ь1 + Л/Юр)р + Ь^2,

А3(р) = Ь0 + Л + (Ь1 + Л/ар)р + Ь^2,

А4(р) = Ь0 + Л + (Ь - Л/ар)р + Ь^2.

В данном случае допустимая достаточная граница для А находится из этого условия (также по одному полиному) немедленно:

(14) Л < Ьюр = 2Ь_^4Ь2Ь - Ь2 .

2Ь2

Сравнение в явном виде с (12) показывает степень консервативности новой оценки.

Проведем примерный расчет, задавшись следующими значениями параметров: Ь2 = 1, Ь1 = 0,1, Ь е [1 - А, 1 + А]. При стационарном Ь, очевидно, полином устойчив для любого Л < 1. При Ь = Ь(0 имеем

М - 0,12

юр = » 0,99875 (рад/е),

и достаточную оценку по (12):

А < I А(/ар, Ь0) I = 2Ь-^4Ь2Ь0 - 4Ь2 » 0,09991.

Достаточная оценка по (13):

А < 0,0944 Достаточная оценка по (14):

А < -Ь-^4Ь2Ь0 - Ь12 » 0,09987.

2Ь2

Численный эксперимент, проведенный в пакете МаАаЬ, показывает, что устойчивость данной нестационарной системы теряется в форме возбуждения первого параметрического резонанса лишь при А > 0, 16 для наихудшего в смысле устойчивости закона изменения а(0 - меандра с частотой 2юр. Разъяснение этого эффекта и методику уточнения границы устойчивости на основе метода гармонической стационаризации дано в [5].

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

а + ЛрЯ < 0,

где а - действительная часть ближайшего к мнимой оси корня А(/а, Ь0); р - модуль наибольшего корня А(ю, Ь0); Я - сумма модулей вычетов Ж(р) = 1/А(р, Ь0). В нашем примере а=-0,05, р= 1, Я » 1,0012, соответственно, приближенно получаем, что для устойчивости системы достаточно выполнения Л < 0,049,

что в практически в два раза осторожнее предлагаемых в статье оценок.

5. Выводы

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

Автор благодарит А. А. Трембу за ряд замечаний и советов.

Литература

1. Нестационарные системы автоматического управления: анализ, синтез и оптимизация / Под ред. К. А. Пупкова,

Н. Д. Егупова. - М.: Изд-во МГТУ им. Н. Э. Баумана, 2007.

- 631 с.

2. ПОЛЯК Б. Т., ЩЕРБАКОВ П. С. Робастная устойчивость и управление. М.: Наука, 2002. - 303 с.

3. ХАРИТОНОВ В. Л. Асимптотическая устойчивость семейства систем линейных дифференциальных уравнений // Дифференциальные уравнения. - 2002. - Т. 1, вып. 11. - С. 2086-2088.

4. ХЛЫПАЛО Е. И. Расчет и проектирование нелинейных корректирующих устройств в автоматических системах. Л.: Энергоиздат, 1982. - 272 с.

5. ЧЕЧУРИН Л. С. Устойчивость нестационарных систем управления с периодически изменяющимися коэффициентами // Изв. вузов. Приборостроение. - 2007. - Т. 50, №. 10.

- С. 23-25.

6. BHATTACHARYYA S. P., CHAPELLAT H., KEEL L. H.

Robust Control: the Parametric Approach. - Upper Saddle River: Prentice Hall PTR, 1995. - 625 p.

7. BONGIORNO, J. J., JR. An extension of the Nyquist-Barkhausen stability criterion to linear lumped parameter systems with time-varying elements, IEEE Trans. - 1963. - AC-8, P. 166-172

ALGEBRAIC SUFFICIENT STABILITY CONDITION FOR PERIODICAL TIME-VARIANT CONTROL SYSTEMS

Leonid Chechurin, St.Petersburg State Polytechnical University, St.Petersburg, Cand.Sc., Dept. Head (cepreu4@gmail.com).

Abstract: Algebraic sufficient stability condition for periodical time-variant control systems is given. The result is based on Kharitonov Theorem for interval systems and Bongiorno sufficient condition for time-variant system stability. The efficiency of the method is illustrated by an example.

Keywords: linear periodical time-variant systems, stability, Kharitonov theorem.

Статья представлена к публикации членом редакционной коллегии П. С. Щербаковым