УДК 519.7
DOI: 10.17586/0021-3454-2015-58-9-707-712
УПРАВЛЕНИЕ НЕЛИНЕЙНЫМИ ОБЪЕКТАМИ С ЗАПАЗДЫВАНИЕМ НА БАЗЕ МОДИФИЦИРОВАННОГО АЛГОРИТМА БЭКСТЕППИНГА
И. Б. Фуртат1, Е. А. Тупичин2
1Университет ИТМО, 197101, Санкт-Петербург, Россия Институт проблем машиноведения РАН, 199178, Санкт-Петербург, Россия Санкт-Петербургский государственный университет, 199034, Санкт-Петербург, Россия
E-mail: [email protected]
2Университет ИТМО, 197101, Санкт-Петербург, Россия
Предложено решение задачи робастного управления нелинейными объектами с запаздыванием по состоянию. Рассмотрены объекты управления в условиях параметрической неопределенности, действия внешних ограниченных возмущений и измерения только выходной переменной; также изложены условия, при которых нелинейная система с запаздыванием может быть приведена к нормальной форме. Для синтеза алгоритма управления используется модифицированный алгоритм бэкстеппинга (алгоритм обратного обхода интегратора). Процесс синтеза алгоритма управления условно разбивается на n шагов, где n — порядок модели объекта. На каждом шаге синтезируется стабилизирующее управление; на последнем шаге синтезируется основной закон управления, который обеспечивает слежение выходного сигнала объекта управления за гладким эталонным сигналом с требуемой точностью за конечное время. Показано, что для реализации алгоритма достаточно использовать всего один фильтр состояния по сигналу управления и упрощенные законы управления, полученные с использованием реальных дифференцирующих звеньев. Это позволяет существенно упростить расчет и реализацию системы управления. Полученная система управления обладает универсальной структурой по отношению к виду модели объекта управления: линейная, нелинейная, с запаздыванием и без него.
Ключевые слова: робастное управление, метод бэкстеппинга, нелинейная система, запаздывание.
Введение. Рассматривается построение системы управления на базе метода бэкстеппинга. Впервые этот метод был предложен в работе [1] для синтеза адаптивного управления нелинейными объектами по выходу. Использование метода позволяет обеспечить в системе управления параметрическую робастность и возможность учета априорной информации о значениях параметров объекта управления. Последнее свойство наглядно продемонстрировано в работе [2], где представлен эффективный алгоритм бэкстеппинга для параметрически неопределенных объектов с измеряемым скалярным выходом. Другие модификации данного метода рассмотрены, например, в работах [3—5]. Однако предложенные в работах [1—5] методы сложны при аналитическом расчете системы управления и ее технической реализации. Сложность аналитических расчетов заключается в громоздкости вычислений полной производной по времени от стабилизирующих сигналов управления. Трудности, возникающие при технической реализации системы, связаны с большим количеством компонентов и фильтров ее состояния, необходимых для формирования закона управления.
В настоящей статье предложен модифицированный алгоритм бэкстеппинга для роба-стного управления параметрическими неопределенными объектами с запаздыванием. Показано, что в отличие от схем, рассмотренных в работе [6], в предлагаемой системе управления реализуется всего один фильтр размерности, равной порядку модели объекта, а для вычисления производных стабилизирующих сигналов управления используются реальные
дифференцирующие устройства. Данный результат является обобщением алгоритмов, представленных в работах [7—9].
Постановка задачи. Рассмотрим математическую модель объекта управления:
x(t) = f (x(t), x(t-т), У + b(x(t), £)(u(t) + ф(x(t), t)); (1)
y(t) = h (x(t)), x(s) = c(5), ^ е[-т;0],
где x(t) e X e Rn — вектор состояния; u(t) e R — управляющее воздействие; y(t) e 7 e R — регулируемая переменная; t e T e (0, да); fx(t), x(t - т), £), b(x, £) и h(x) — гладкие функции соответствующих размерностей; е S — вектор неизвестных параметров; S — известное ограниченное множество; т > 0 — неизвестное время запаздывания; t) — неизвестная функция; а(5) — кусочно-непрерывная функция на отрезке [-т, 0].
Необходимо синтезировать закон управления, обеспечивающий выполнение целевого условия
\y(t)- Ум(t)| <5 для t > tf, (2)
где 5 > 0, yм(t) — эталонный сигнал, tf > 0 — время переходного процесса. Введем следующие обозначения:
Lfh( x) = f (x, x (t - т), £), Lfh (x) = § gx(t<- ,т) (fh) (x(t -iT), x(t - (i +1) т), £); Lbh(x) = b(x, £), LbLfh(x) = £(Lfh)b(x, £)+ (Lfh)b(x(t-т), £) 5(т);
LbLnf-'h(x) = ^ f . b (x(t-h), £)5(т), Pi = 9 . b (x(t-nr), £).
Здесь Ьу-к( х), Ььк( х) — производные Ли от функции к(х) по направлению векторных полей
_Дх, х(^ - т), £), Ь(х, £); 5(/т) — оператор смещения на величину /т (например, 8(/т)и(7) = и(^ - /т)). Введем вектор г($), определяемый как
-Г
z =
У, У,..., У( n -1)
h(x), Lfh(x),..., Lnflh(x) = 0(x(t), x(t - т),..., x(t - (n-1)т)).
Если выполнены условия
Ьък( х) = ЬЪЬ/И( х) =... = ЬЪЬп/ 2И( х) = 0,
Р(х, х(* - т),..., х(* - (п - 1)т), £) = ЩМх) * 0, то уравнение (1) может быть преобразовано к системе следующего вида:
¿1(*) = 12($);
¿п -1(г) = ¿п (1 X
¿п «) = Ьпгк(х) + ЬьЬпг-1к(х) (и(/) + ф(х(Г), I)).
Предположение 1. Функции /(х(£), х^ - т), £), Ъ(х, £), к(х) — гладкие, и для любых х(^), х(^ - т) и Е, е 2 выполнены следующие условия:
Ьък( х) = ЬЪЬ1/И( х) =... = Ь^}- 2к( х) = 0,
1 дк (х) . . дк (х) . .
где Ь;к(х) = —— /(х, £), Ьък(х) = —— Ъ(х, £) — производная Ли от функции к(х) по на-
дх дх
правлению векторных полей _Дх, £), Ъ(х, £) соответственно.
Производные высших порядков вычисляются по формулам
2 д((А(х)) к д^ед)
Ь2/И(х) = -Ц--(х, & ..., Ьк/И(х) = -11 (X, О.
дх дх
Предположение 2. Функция Ф взаимно-однозначная при х(г - /т) = 0, г = 0, ..., п - 1. Предположение 3. Разностные уравнения ф(х(г), х(г-т),..., х(г-(п-1)т)) = 0,
п-1
^ Ргм (г - /т) = 0 асимптотически устойчивые относительно переменных х(г) и м(г).
г=0
Предположение 4. Функция с(х, х(г-т),..., х(г-(п-1)т), = х) ограниченная
или ограниченная на множестве 2 и липшицева по х(г- /т), г = 0, ..., п - 1.
Предположение 5. Функция ф(х(г), г) ограниченная или ограниченная на множестве Т и липшицева по х(г) е Х.
Предположение 6. Сигналы ум(г), ум (г),..., у^^г) — ограниченные. Аналогичные предположения рассмотрены в работе [10]. Метод решения. Продифференцируем п раз функцию у(г):
рпу(г) = с(х, х(г - т),..., х(г - (п - 1)т), £) + +Р(х, х(г-т),..., х(г-(п-1)т), £) (м(Г) + ф(х, г)), (3)
где р = ё / ёг — оператор дифференцирования.
С учетом выражения (3) запишем уравнение ошибки ^(г) = у(г) - ум(г) в виде
рпв1 (г) = с(х, х(г - т),..., х(г - (п - 1)т), £) + +Р(х, х(г-т),..., х(г-(п-1)т), £)((Г) + ф(х, г))-рпум(Г). (4)
п-1
Введем в рассмотрение оператор 0,п-1(р) = ^ кп-/рг, такой что полином Q(X) = Хп + 0п-1(^)
/=0
гурвицев, где X — комплексная переменная, и перепишем уравнение (4) в виде
Q (р)^(г) = м(г) + у( г), (5)
где
у(г) = с( х, х(г-т),..., х(г-(п-1)т), £) + Р( х, х(г-т),..., х(г-(п-1)т), ^)ф(х, г) +
+(р( х, х(г- т),..., х(г - (п - 1)т), £))))- рпУм (г) - Q"-l(p)e(г). (6)
Рассмотрим фильтр
\>(г) = Лу(г) + (г), (7)
-к1
-к2 1п-1
-кп 0 ... 0
т п-1)х( п-1)
, 1п_1 е ю ' х ' — единичная мат-
где У(г) = [(г), У2(г),..., ^(г)], А =
рица, I = [0, ..., 0, 1]Г.
Перепишем уравнение (5) с учетом выражения (7):
*1(г) = vl(г) + Q-1( р)у(г). (8)
Дифференцируя (8), получаем
¿1(г) = -кЛ(г) + У2(г) + / (г), (9)
где Q(p)/ (г) = ру(г).
Представим процедуру синтеза системы управления следующим алгоритмом.
Шаг 1. Определим у2(г) в виде у2(г) = й^г) и зададим й^г) как
Щ(г) = -а1^-1е1(г) + кЛ(г), (10)
где а1 > 0 и ц > 0 — коэффициенты, выбираемые разработчиком.
Подставив (10) в выражение (9), получим
¿1(г) = -а!|д-1е1(г) + / (г). (11)
Шаг г (2 < г < п - 1). Рассмотрим функцию ошибки ег(г) = уг(г) - й _ 1(г). Взяв производную от ег(г) вдоль траекторий (5), получим
е (г) = -кЛ(г) + ^(г)-йм(г). (12)
Предположим, что функция Vг + 1(г) — сигнал управления в уравнении (12). Пусть V. + 1(г) = й(г), тогда
и (г) = -ае (г) + кл(г) + -1(г), (13)
где аг > 0 — коэффициент, выбираемый разработчиком, й-1(г) — оценка сигнала C/г■_l(г).
Подставив (13) в выражение (12), получим
е (г) = -ае (г) -Л/-1(г), (14)
где Л/-1 (г) = им(г) - й1 -1(г).
Шаг п. Рассмотрим функцию еп(г) = V"(г) - йп- 1(г). Принимая во внимание уравнение (7) и дифференцируя еп(г), получаем
К (г) = -кМг)+и (г)-й^г). (15)
Сформируем закон управления
и (г) = -а^ (г)+к^ ( г) + й^ (г), (16)
где ап > 0 — коэффициент, выбираемый разработчиком, йп-1(г) — оценка функции С/п-1 (г), с учетом которого перепишем выражение (15):
еп ( г) = -апеп ( г) -Лп-1( г), (17)
где Лп-1(г) = й^ г) - йп -1( г).
Для оценки производных й{- (г), / = 2, п, воспользуемся следующими наблюдателями
й/-1(г) = -|-1й/-1(г) + |-1й//-1( г), / = 2^. (18)
Утверждение. Пусть выполнены предположения 1—6. Тогда существуют константы а/ >0, / = 1, п, и ц0 > 0, такие что для ц < ц0 система управления, состоящая из фильтра (7), вспомогательных управляющих воздействий (10), (13), закона управления (16) и наблюдателей (18), обеспечивает выполнение целевого условия (2) для объекта (1).
Доказательство утверждения аналогично приведенному в работах [6—9].
Заключение. Предложен алгоритм слежения выходного сигнала нелинейного объекта управления с запаздыванием за эталонным сигналом в условиях параметрической неопределенности и внешних возмущений. Синтез алгоритма базируется на новой версии метода бэкс-теппинга, предложенного в работе [6]. Полученная система управления содержит всего один фильтр размерности, равной размерности модели объекта, а для реализации производных в законах управления используются реальные дифференцирующие звенья.
Статья подготовлена по результатам работы, выполненной при частичной поддержке Российского фонда фундаментальных исследований (грант 13-08-01014), Министерства образования и науки РФ (проект 14.Z50.31.0031) и Правительства РФ (грант 074-И01); результаты, приведенные в разделе „Метод решения", получены при поддержке Российского научного фонда (проект № 14-29-00142) в ИПМаш РАН.
список литературы
1. Kanellakopoulos I., Kokotovic P. V., Morse A. S. Systematic design of adaptive controllers for feedback linearezable systems // IEEE Trans. Automatic Control. 1991. Vol. 36. P. 1241—1253.
2. Никифоров В. О. Адаптивное и робастное управление с компенсацией возмущений. СПб: Наука, 2003.
3. Khalil H. K. Nonlinear Systems. N. Y.: Prentice Hall, 2002.
4. Zheng Y., Yang Y. Adaptive output feedback control for class of nonlinear systems with unknown virtual control coefficients signs // Adaptive Control and Signal Processing. 2007. Vol. 21, N 1. P. 77—89.
5. Tanner H. G., KyriakopoulosK. J. Backstepping for nonsmooth systems // Automatica. 2003. Vol. 39. P. 1259— 1265.
6. Фуртат И. Б. Модифицированный алгоритм обратного обхода интегратора // Мехатроника, автоматизация, управление. 2009. № 10. С. 2—7.
7. Furtat I. B., Tupichin E. A. Modified simple adaptive-robust backstepping algorithm // Proc. of the 19th Intern. Conf. on Methods and Models in Automation and Robotics (MMAR 2014), Mi^dzyzdroje, Poland. 2014. P. 183—188.
8. Furtat I. B., Tupichin E. A. Modified robust backstepping algorithm for plants with time delay // Proc. of the 6th Intern. Congress on Ultra Modern Telecommunications and Control Systems and Workshops (ICUMT), St. Petersburg, Russia. 2014. P. 541—545.
9. Furtat I. B., Tupichin E. A. Control of nonlinear plant based on modified robust backstepping algorithm // Proc. of IEEE Intern. Conf. on Control Applications (CCA), Antibes, France. 2014. P. 941—946.
10. Цыкунов А. М. Робастное управление с компенсацией возмущений. М.: ФИЗМАТЛИТ, 2012.
Сведения об авторах
Игорь Борисович Фуртат — д-р техн. наук, профессор; Университет ИТМО, кафедра систем
управления и информатики, профессор; ИПМаш РАН, лаборатория управления сложными системами, ведущий научный сотрудник; СПбГУ, кафедра прикладной кибернетики, ведущий научный сотрудник; E-mail: [email protected] Евгений Александрович Тупичин — аспирант; Университет ИТМО, кафедра систем управления и информатики; E-mail: [email protected]
Рекомендована лабораторией Поступила в редакцию
управления сложными системами 22.04.15 г.
ИПМаш РАН
Ссылка для цитирования: Фуртат И. Б., Тупичин Е. А. Управление нелинейными объектами с запаздыванием на базе модифицированного алгоритма бэкстеппинга // Изв. вузов. Приборостроение. 2015. Т. 58, № 9. С. 707—712.
MODIFIED BACKSTEPPING ALGORITHM FOR CONTROL OVER NONLINEAR PLANTS WITH DELAY
I. B. Furtat1, E. A. Tupichin2
1ITMO University, 197101, Saint Petersburg, Russia Institute of Problems of Mechanical Engineering of the Russian Academy of Sciences, 199178, Saint Petersburg, Russia Saint Petersburg State University, 199034, Saint Petersburg, Russia E-mail: [email protected]
2ITMO University, 197101, Saint Petersburg, Russia
A solution to the robust control problem of nonlinear plants with a state delay is proposed. The plants under control are considered to possess unknown parameters; the plants are subject to external disturbances, and only output signal are measurable. The conditions for such a nonlinear system to be reducible to the normal form are formulated. The modified backstepping algorithm is used for the control algorithm synthesis. It is shown that implementation of the algorithm calls for the use of only one filter by the control signal and simplified control laws, and therefore calculations and implementation of the control system may be simplified. The obtained system is shown to possess a universal structure independent of the plant model - linear, nonlinear, with or without delay.
712
H. E. 0ypmam, E. A. Tynmun
Keywords: robust control, backstepping method, nonlinear system.
Data on authors
Igor B. Furtat — Dr. Sci., Professor; ITMO University; Department of Computer Science
and Control Systems; IPME RAS, Laboratory Control of Complex Systems; Saint Petersburg State University, Department of Applied Cybernetics; E-mail: [email protected] Evgeny A. Tupichin — Post-Graduate Student; ITMO University; Department of Computer Science and Control Systems; E-mail: [email protected]
For citation: Furtat I. B., Tupichin E. A. Modified backstepping algorithm for control over nonlinear plants with delay // Izvestiya Vysshikh Uchebnykh Zavedeniy. Priborostroenie. 2015. Vol. 58, N 9. P. 707—712 (in Russian).
DOI: 10.17586/0021-3454-2015-58-9-707-712