Абсолютная устойчивость одной системы прямого управления
Автор: Поливенко Виталий Кузьмич
Журнал: Математическая физика и компьютерное моделирование @mpcm-jvolsu
Рубрика: Математика
Статья в выпуске: 10, 2006 года.
Бесплатный доступ
В работе решается задача абсолютной устойчивости системы регулирования на плоскости с прямым управлением. Найдены области изменения значений параметров управлений и свойства двух исполнительных органов (нелинейных функций), при которых система абсолютно устойчива.
Короткий адрес: https://sciup.org/14968587
IDR: 14968587
Текст научной статьи Абсолютная устойчивость одной системы прямого управления
В работе решается задача абсолютной устойчивости системы регулирования на плоскости с прямым управлением. Найдены области изменения значений параметров управлений и свойства двух исполнительных органов (нелинейных функций), при которых система абсолютно устойчива.
Рассмотрим математическую модель объекта управления в виде системы
Xk = fM\ А(0) = 0, £ = 1,2, (А)
где {Ti(t), x2(t^ — вектор состояния объекта, (иДя?!,^), ui^i,^)} — вектор управления.
Будем предполагать:
I) движения Xk(t, t0, х°\ к = 1,2,х° = {ж®,^} Е R2 = Oxi.x2 наблюдаются на промежутках /(t0) = [i0, 4-оо)\тД при управлении вида:
ui = axi + bx2, . - l л
(В)
△ = аа — ос = 0; и2 = сху + ах2;
-
II) в В? функции fk непрерывны и таковы, что система (А):
-
1. Не имеет линейной части, то есть матрица А = О — нулевая.
-
2. Имеет единственное решение задачи Коши х^ = Xk^to^Yk = 1,2, каковы бы ни были допустимые управления Ик.
В теории автоматического управления (регулирования) (ТАУ) (А), (В) — называют системой прямого управления (регулирования) с двумя исполнительными органами /i(«i) и /2(^2)-
(с) В.К. Поливенко, 2006
Для системы (А) с (В) решаем задачу анализа ТАУ: найти область S изменения значений параметров a,b,c,d в управлениях Uk, для которых система возмущения (А) асимптотически устойчива при любых начальных возмущениях х° и любом выборе Д, удовлетворяющих обобщенным условиям Рауса — Гурвица, то есть система (А) абсолютно устойчива.
Если множество S не пусто, то управления (В) называют законом обратной связи или законом регулирования. Его реализация осуществляется с помощью совокупности измерительных приборов, усилителей, преобразователей и исполнительных устройств, называемых регулятором.
В силу предположения II. 1 мы не можем использовать методику исследования абсолютной устойчивости работы [1, гл. III, IV]. Поэтому определение области S осуществляем на основании прямого метода Ляпунова и некоторых качественных приемов исследования динамических систем на плоскости, которые приведены в работе [2] и других.
В системе (А) перейдем от координат (xi,^) к координатам («1,^2), преобразование (В), осуществляющее этот переход, согласно Д ^0 является невырожденным. Новая система будет иметь вид
< Й1 = af^u^ + bf2(u2), ^2 = с/1(гц) + d/2(u2).
Система (1) является системой непрямого регулирования.
Пусть всюду в дальнейшем функции hk^k) = М^1 при Uk ^ о, k = 1, 2, удовлетворяют обобщенным условиям Рауса — Гурвица системы (1):
' ah-ituy) + dh2(u2) < 0,
(ad — be) x /ii(ui)/i2(u2) = Д x /ii(ui)/i2(u2) > .0
Будем считать, что △ = ad — be > 0, тогда из (3) следует, что и пусть ради определенности:
Mui) > о, /г2(г12) > 0,
Ui + 0, и2 ^ 0
Jк^к^^к > 0, к — 1, 2.
Тогда вследствие (2) и (4) случай а > 0, с? > 0 невозможен.
Поэтому возможны следующие случаи:
-
1) ad < 0, при этом из Д > 0 следует, что Ьс < 0;
-
2) а < 0 и d < О, в этом случае произведение Ьс знаконеопределено.
В случае 1), не умаляя общности, будем считать:
а > О, d < О
(если а < 0, d > 0, то, переименовав переменные Ui на и2 и и2 на Ui, придем к тому же случаю.)
В рассматриваемом случае 1) из (2) следует, что функция h^uiYui 7^ 0 ограничена сверху, а функция h,2(u2\u2 ^ 0 — снизу.
Положим:
ci = sup Мщ^; с*к = inf hk(ukY к = 1,2, Ui/0
тогда вследствие (2)-(4) при Uiit2 ^ 0 имеют место неравенства:
оКх^ + dc2 < О, ^ < /12(^2); (51)
hx (ui) < Ci, асх + dh2M < 0, (52)
причем в каждом из условий (51) и (5г) хотя бы одно из неравенств является строгим, где ci > 0, с^ > 0.
Таким образом, в рассматриваемом случае 1):
а > 0, d < 0, Ci > О, С2 > О
Введем в рассмотрение условия:
° UW1) - A(Mi)l + dc2(ux - и") > О их< и", /2(^2) — С2М2 — /2(^2) — С2«2 U1 < м2"
Теорема 1. Пусть b < 0 и имеет место случай!): (5*). Если выполнены условия (2), (4) и (6), (7), то для абсолютной устойчивости системы (А) необходимо и достаточно, чтобы выполнялись условия:
U1
lim
U1—* + 00
У fx^dux - (aMui) + dc*^ux ,0
—^— lim \acx + dh2('U2^U2 = +00, (8) Z\C2 U2—*+OO
lim
Uj—>—00
Ui
У /i(ui)cZui + (a/zi(ui) + dc2)ux
-O
lim \acx4-dh2(u2^ti2 = +00. (9) AC^ U2->-OO
Доказательство. Достаточность. Покажем, что неограниченные траектории системы (1), если они существуют, при возрастании t образуют раскручивающиеся против часовой стрелки спирали.
Рассмотрим функцию
U1
Д Г 1
V(^i;m2) = ^ №x)dux + -(dux - bu2f.
О
В силу (51) либо a/ii(ui) + dc2 < 0,ui ^ 0; либо /12(^2) > ^,«2 + 0. Поэтому производную функции y(ui,u2) в силу системы (1) можно представить в виде
V(ui,u2) = 4 (hx(ux^ahx(ux)-V dc*^ + bhx^tix)Vh2^^^^ - c^]uiM2) -C2
Вследствие A > 0, (5i), c2 > 0 и b < 0 (из be < 0 —> c > 0)
V(ui,m2)<0, wiu2 > 0, то есть в первой и третьей четвертях плоскости 0ttxtt2, причем V(0, ц2) = 0 (на оси 0и2).
Если
±00
j fi^dtl^ = 4-00, (11)
О то функция y(u!,u2) является бесконечно большой и ее линии уровня V(ui,w2) = С в первой и третьей четвертях пересекают оси координат фазовой плоскости Ои^г.
Поэтому неограниченная положительная полутраектория f+(p,t),p Е I четверти должна покинуть эту четверть. Из
^2|u2=o ch^uy > 0,щ /* 0, zzi|U1=0 6/z2zz2 < (Ь/2 > О и рассмотрения поля направлений системы (1) ясно, что это может произойти только через положительную полуось ординат 0и2.
Пусть условие (11) нарушается и пусть ради определенности сходится интеграл
+оо j f^u^dui. о
Тогда вследствие (8) f+(p, t) не уходит в бесконечность в I четверти. Покажем это.
Заметим, что в рассматриваемом случае линии уровня V(u1,u2) = С функции (10) при достаточно больших С состоят из двух ветвей, расположенных по обе стороны прямой du^ — Ьк2 = 0 и лежащих внутри полосы вида
\du! - М < С.(12)
Очевидно, полутраектория f+(p,t) может уходить в бесконечность только внутри этой полосы.
В силу (8) либо lim [a/ii(ui) + dc^tii = —сю(13)
U1—»+оо либо lim [d/z2(u2) + aci]u2 = —oo.(14)
U2—*+OO
Если выполняется (13), то вдоль f+(p,t) = {ui(<),u2(i)},t > 0 при достаточно больших Их iii = [a/zi(ui) + dc^i + ^[^2(^2) — c^u2 — c^du^ — bu2) <,
< [a/ii(ui) + dc^tti + c^chi — 6u2| < [a/ii(ui) + dc^ui + C < 0, то есть ux(t) убывает вдоль f+(p, t).
Кроме того,
Й2|и2=о = ск^и^гц > 0 с > О, Ai(mi) > 0, мт ^ О, то есть f+(p,t) пересекает ось Omi в первой четверти снизу вверх. Пусть теперь (13) не выполняется. Это возможно только при асг 4- dc^ = 0.
Поэтому при больших «2 > 0 в полосе (12) в силу b < 0, d < 0 и (5г) имеем
Mi = (ahi + dc^tii + -\dh2 — bc^W — ^(dui — 6m2) < d
b
< (ahi + dc^tii + ^[dd2 + ac^u2 + c^dui — bu2\ < "[d/i2 + ac]tt2 + c^C < 0.
d
Из этих рассуждений следует, что траектория f+(p,t) не может уходить в бесконечность в I четверти вдоль прямой du,i—bu2 = 0. Следовательно, траектория f+(p,t) должна покинуть I четверть только через полуось 0м2,м2 > 0.
Во II четверти M2(t) вдоль траектории f+(p,t) убывает, поскольку в силу (4) с > 0 и d < 0,
1 du^ 2~dT
= chi(ui)tiiti2 + dh2(ti2)ti2 < 0
Mi < 0.
При дальнейшем возрастании t f+(p,t) пересекает отрицательную полуось 0мг. Покажем это.
С этой целью введем в рассмотрение функцию
Vi(mi,M2) = |М1| - ~|м2| = MiSgHMi + 3M2SgnM2 М1М1 + 0. d
В силу системы (1)
■ (be \
Vi(ui,tt2) = hiUi asgnMi + —sgnM2 + b/z2M2(sgnMi + sgnM2).
\ “/
При M1M2 < 0 имеем sgnM2 = —sgnMi и
Vi = —/iiMisgnMi < 0 при mi ^ 0.(15)
Поскольку Vi(mi,m2) бесконечно большая функция, то ее линии уровня Vi(mi,m2) = С пересекают ось Омь Поэтому во II четверти вследствие (15) f+(p, t) пересекает ось 0mi,Mi < 0.
В нижней полуплоскости рассуждения аналогичны. Отсюда видно, что уходящая в бесконечность (по предположению) положительная полутраектория f+(p, t) должна образовывать раскручивающуюся против часовой стрелки спираль.
Далее, дословно повторяя соответствующие рассуждения из доказательства теоремы 3.2 гл. I работы [2], можно показать, что вследствие (6) и (7) ни одна траектория системы (1) при возрастании времени t не может раскручиваться против часовой стрелки. Полученное противоречие показывает, что в условиях теоремы все положительные полутраектории f+(p,t) ограничены.
Затем, ссылаясь на лемму 3.1 из [2], завершаем доказательство достаточности условий (8) и (9).
Необходимость этих условий доказывается от противного, как в теореме 3.2. [2]. Теорема доказана.
Замечание. Если aci + dc2< 0, то условия (8) и (9) выполняются автоматически.
Это неравенство следует из неравенства (2) в том случае, когда точные границы Ci и с2 функций /ii(ui) и /12(^2) достигаются соответственно при каких-то u'f 7^ 0 и и2 7^ 0, то есть
Су = SUp /11(^1) = /11(^1), ui/o
= inf Tti(u2) = /t2(u2). U2#0
Для случая 2): а < 0, d < 0 имеет место следующая теорема.
Теорема 2. Пусть выполняются условия (2), (4) и имеет место случай 2). Тогда для абсолютной устойчивости системы (А) достаточно выполнения неравенства inf h(u2) = с2 > 0, be < —ad
112 7^0 ~
Доказательство. Заметим прежде всего, что в случае 2) при выполнении (2) и (3) имеем Д = ad — be > 0. Введем в рассмотрение определенно положительную функцию
Ui U2
V2(U1,U2) = |С1 У fl^dU! + |6| j /2(u2)du2.
о о
В силу системы (1)
Ы*С, и^ = a\c\f2(u^ + [Ь|с| + |6|c]/i(ui)/2(tt2) + d|6|/a (^2) < о, (17)
согласно критерию Сильвестра:
А = а|с| < 0,
В = Ь|с| + |6|с,
К = d\b\ < 0
и
D = В2 - 4АК = 2|Ьс|[|Ьс| + Ьс - 2nd] =
—4ЬсД < 0
4bcad < 0
Ьс > О, Ьс<0 '
то есть V2 = W(/i,f2) — определенно отрицательная квадратичная форма относи тельно Jk, k = 1,2.
Поэтому, если функция V2 является бесконечно большой, то вследствие I теоремы Барбашина — Красовского [1] нулевое решение и^= «2 = 0 системы (1) устойчиво в целом, а система (А) абсолютно устойчива.
Пусть РДм^иг) не является бесконечно большой. Вследствие (17) для любой ТОЧКИ р фаЗОВОЙ ПЛОСКОСТИ 0«1«2
У2и+(р,ед < VM
Выберем С > 0 так, чтобы было ^(р) < С. Тогда согласно (17) f+(p, t) не покидает области V^-u^u^ < С и при достаточно больших |«2| содержится внутри некоторой полосы |«i| < ki, а при достаточно больших |«i| внутри некоторой полосы «2| < к2.
Но при |«i| < ку, вследствие (16), d < 0 и достаточно больших |м2| имеет место неравенство т-7- < с/1(«1)«2 + dc2u2 < 0.
-
2 at
Поэтому вдоль f+(p,t) |«2(<)1 ~ ограничен.
Вдоль положительной полу траектории ограничен также и |«i(t)|. Покажем это.
В зависимости от знака параметра b дальнейшее доказательство подразделяется на два подслучая.
Подслучай А: b < 0. При \и2\ < к2 и достаточно больших |«i| имеем
1 du^ а
«1«2 > 0,
—— = апуш + оп^иуи^ < 0
2 dt ’ '
причем знак «=» только при «1 = 0, а функция 14 («1, «г) имеет полную производную (15) в силу системы (1)
Vi(ui,«2),) < 0
Mj«2 < 0.
Вследствие (15), (17) и (18) линии уровня Vi(«i,«2) = С, V^Ux,^ = С и прямые х = С пересекаются траекториями системы (1) в сторону начала координат.
Поэтому в этом подслучае |«i(i)| вдоль f+(p,t) также ограничен.
Таким образом, ограниченность f+(p,t) доказана. Далее, рассуждая как в разделе 1.2 главы 1 работы [2], нетрудно показать, что любая полутраектория f+(p,t) —> (0,0) и состояние равновесия «х = и2 = 0 системы (1) устойчиво по Ляпунову.
Подслучай Б: b > О. Рассмотрим функцию Vi(«i,«2), которая в этом подслучае имеет вид
Vl(«i,«2) = «lSgn«i - -M2SgnM2, d а ее полная производная в силу системы (1) и того, что b > О, Д > 0, (4) и (16)
у _ 1 7^iWiSgn«i < 0 «1«2 > 0.
) ^^^«xSgnUi — 2b/z2«2Sgn«2 < 0 «1«2 < 0.
Поскольку в рассматриваемом подслучае функция Vi(ui,U2) является бесконечно большой, то есть 14 (гл,г<2) —* +оо при Itti^l ~* +оо. а ее производная 14 < О при всех Ui и и2, то по I теореме Барбашина — Красовского нулевое решение системы (1) устойчиво в целом.
Поэтому для этого подслучая в теореме введено дополнительное условие (16): ad + bc > 0. Теорема доказана.
Заключение. Итак, нами показано, что в случае 1):
-
△ = ad — be > 0, ad < О
при выполнении условий (2), (3), (4), (6) и (7) для абсолютной устойчивости системы (А), (В) необходимо и достаточно выполнения условий (8) и (9);
в случае 2): а < 0, d < О, Д > 0 при выполнении обобщенных условий Рауса — Гурвица (2), (4) и inf A2(rz2) = С2 > 0, ad + be > О
U2^O . '
достаточно для абсолютной устойчивости системы (А), (В).
Список литературы Абсолютная устойчивость одной системы прямого управления
- Барбашин Е.А. Функции Ляпунова. М.: Наука, 1970. 240 с.
- Поливенко В.К. Глобальное поведение решений систем дифференциальных уравнений на плоскости. Волгоград, 1997. 236 с.