CC BY
50
yeaeooi i i i a i a o - i i - oaoi e - a n e i a egaai ea
НАУКА и ОБРАЗОВАНИЕ
Эл ^ ФС77 - 30569. Государственная регистрация №0421100025. ISSN 1994-0408
Стабилизация неминимально фазовых аффинных систем с использованием линеаризации по части переменных
77-30569/255087
# 10, октябрь 2011 Ткачев С.Б.
УДК 517.977
МГТУ им. Н.Э. Баумана [email protected]
1. Введение
В теории нелинейных динамических систем с управлением известен ряд важных теоретических результатов, позволяющих для достаточно широкого класса систем решать задачи стабилизации положений равновесия. Основное внимание при этом уделялось аффинным системам — нелинейным системам, линейным по управлению. Среди таких систем были выделены специальные виды, для которых разработаны методы решения задачи стабилизации и с использованием дифференциально-геометрического подхода получены условия эквивалентности аффинных систем и систем специального вида.
Среди систем специального вида выделены канонические и квазиканонические виды стационарных аффинных систем со скалярным и векторным управлением [1-10]. Для аффинных систем, преобразуемых к регулярному каноническому виду [1-6], решение задачи стабилизации положения равновесия известно и один из подходов заключается в преобразовании исходной системы в линейную управляемую систему с помощью нелинейной замены переменных и введения нового управления.
Условия приводимости аффинной системы к каноническому виду выполнены не всегда, и довольно часто система может быть преобразована только к квазиканоническому виду [4] (А.П. Крищенко).
Отметим, что функцию, определяющую преобразование аффинной системы к квазиканоническому виду, часто удобно рассматривать как выход системы. Для стационарных аффинных систем с выходом введены нормальные формы [11-13] (Byrnes C., Isidori A.).
Для решения задачи стабилизации положения равновесия стационарной аффинной системы, преобразованной к квазиканоническому виду или нормальной форме, существенным является наличие свойства минимальной фазово-сти [11,12]. Для минимально фазовых систем известно решение задачи стабилизации положения равновесия статическими и динамическими обратными связями по состоянию [11-17].
В случае, если аффинная система не является минимально фазовой, проблема стабилизации ее положения равновесия оказалась достаточно сложной и подходы к ее решению известны в частных случаях [18-24].
Одним из методов, позволяющих найти стабилизирующую обратную связь для неминимально фазовой системы, является метод виртуальных выходов [25-28]. Применение этого метода связана с анализом нелинейной подсистемы, выделяющейся после преобразования исходной системы к регулярному квазиканоническому виду и линеаризации системы квазиканонического вида обратной связью по части переменных.
Анализ этой подсистемы не всегда приводит к нахождению виртуального выхода с требуемыми свойствами именно из-за существенного влияния нелинейности. Одним из возможных подходов — использование для синтеза стабилизирующего управления линеаризации подсистемы, отвечающих за нулевую динамику.
Исследуем случай, когда соответствующая нелнейная подсистема допускает нетривиальную линеаризацию и опишем метод поиска такого виртуального выхода, которому соответствует нормальная форма с асимптотически устойчивой нулевой динамикой.
В отличие от общего случая поиска виртуального выхода через специально подбираемое виртуальное управление, предлагаемый метод не требует подбора специальных функций, а базируется только на использовании вычислительных процедур.
2. Квазиканонический вид и нормальная форма аффинной системы
Рассмотрим гладкую стационарную аффинную систему со скалярным управлением
х = А(х) + В(х)и, х е Шп, и е К1, А(х) = (а1(х), . . . , ап(х))т, В(х) = (Ъ1(х) аг(х),Ъг(х) е Ож(П), г = 1/п, П С Кп,
('•) = а( ■ )/<И,
где £ — независимое переменное.
Пусть система (1) в некоторой окрестности и0 точки х0 преобразуется к виду
¿1 = ¿2, ... , ¿р-1 = Хр,
¿р = / (г,п) + g(¿,7n)u, (2)
П = я(г, п) + и)(х, п)и,
где г = (г1,..., гр)т, п = (П].,... ,ПтУ, П е Кт, р + т = п, и (г,п) = Ф(х), — соответствующая локальная гладкая невырожденная замена переменных.
Указанный вид называют [4] квазиканоническим видом системы (1). Если коэффициент д(г,п) при управлении в точке (г0,п0) = Ф(х0) отличен от нуля, то квазиканонический вид называют регулярным в этой точке. В этом случае в силу гладкости коэффициент при управлении будет отличен от нуля и в некоторой окрестности указанной точки.
Отметим, что при р = п система (2) содержит управление только в последнем уравнении и такой вид называют каноническим.
Условия, при которых система (1) локально преобразуется к регулярному квазиканоническому виду (2), удобно записать, используя дифференциально-геометрический подход, при котором гладкой аффинной системе (1) на О взаимно-однозначно соответствуют гладкие векторные поля
пп
А = Е а(х) д;-, В = Е ь-(х) (3)
¿=1 ¿=1
Будем использовать обозначение ЬхХ(х) для производной Ли функции \(х) по векторному полю X.
Ъп(х)У,
(1)
Пусть существует достаточно гладкая функция ф(х), удовлетворяющая в некоторой окрестности точки х0 системе нелинейных уравнений в частных производных
ЬеЬАФ(х) = 0, г = 0,р - 2 (4)
для некоторого числа р, где р < п, а соответствующая ей функция 7 (х) = = Ье ЬА_1ф(х) в точке х0 удовлетворяет условию
7 (х0) = Ье ЬА-1ф(х0) = 0. (5)
Приведенные условия (4)-(5) эквивалентны тому, что в окрестности точки х0 производные от функции ф(х) в силу системы (1) до порядка р—1 не содержат управление, а в производную порядка р управление входит с коэффициентом, отличным от нуля в точке х0, и, в силу гладкости, в некоторой окрестности точки х0.
Из условия 7(х0) = 0 вытекает, что В(х0) = 0. Из условия (4) следует, что функции ЬАФ(х), г = 0, р — 2, являются первыми интегралами векторного поля В. Выполнение условия (5) гарантирует локальную функциональную независимость функций г = ЬА—1ф(х), г = 1, р.
К множеству функций г = ЬА—1ф(х), г = 1,р, можно добавить еще т = п — р функций Пк = П^(х) так, чтобы в окрестности точки х0 получить невырожденную замену переменных (г,п) = Ф(х), где г = (г1,..., гр)т, П = (п1,... ,Пт)т, р + т = п. В новых переменных система (1) запишется в виде (2). Поскольку из ЬеЬ^—1ф(х0) = 0 следует, что д(Ф(х°)) = 0, квазиканонический вид будет регулярным в точке Ф(х0).
Заметим, для системы (1) в общем случае могут существовать различные функции ф(х), удовлетворяющие условиям (4), которым могут соответствовать различные значения р от р =1 до р = ртах < п, где ртах есть максимальное значение р, при котором удовлетворяются условия (4), (5). При р =1 имеем Ьеф(х) = 0 и управление входит в первую производную функции ф(х) в силу системы (1).
Вместо системы (4) можно использовать эквавалентную систему уравнений в частных производных первого порядка [1,6]
а^ Вф(х) = 0, г = 0, ...,р — 2,
где а^ В = [А, В] — коммутатор векторных полей А и В,
а^ В = [А, а^11 В], а< В = В.
Квазиканонический вид (2) можно упростить за счет выбора специальной замены переменных. Поскольку В (х0) = 0, в некоторой окрестности точки х0 существуют п — 1 функционально независимых первых интегралов векторного поля В, из которых для построения замены использовались только р — 1 функции = ЬА—1ф(х), г = 1,р — 1. Добавляя к ним функцию ¿р = Ь^"1ф(х) и еще т = п — р первых интегралов Пк = Пк(х) векторного поля В, в окрестности точки х0 можно получить такую невырожденную замену переменных (¿, п) = Ф(х), что в новых переменных система (1) запишется в виде
¿1 = ¿2,..., ¿р—1 = ¿р,
¿р = / (г,П)+ 0(г,пК (6)
П = ^(¿, п).
Такой специальный квазиканонический вид системы (1) более удобен для анализа.
Преобразование системы (1) к некоторому квазиканоническому виду всегда существует [4], но методы стабилизации положения равновесия для таких систем разработаны только в частных случаях.
Приведение аффинной системы к квазиканоническому виду (2) (или (6)) основано на поиске специальной функции ф(х), с использованием которой строится замена по части переменных ¿. Существует подход к преобразованию аффинных систем к специальному виду, используемый в случае, если такая функция уже задана. Эту функцию называют выходом аффинной системы (1) и с ней связывают цель управления — стабилизацию заданного значения выхода или отслеживание заданного изменения выхода как функции времени [13].
В качестве выхода может быть выбрана и произвольная досточно гладкая функция состояния системы — виртуальный выход.
Пусть функция у = ф(х), где ф(х) £ — (виртуальный) выход
аффинной системы (1). Предположим, что существует такое р > 1, что, во-первых, выполнены условия (4) и, во-вторых,
ЬЕ ЬА—1ф(х0) = 0. (7)
Такое число р называют относительной степенью аффинной системы (1) с (виртуальным) выходом у = ф(х) в точке х0 [13].
Если относительная степень р = 1 в точке х0, то Ьвф(х0) = 0. Если р > 1, то первое условие означает, что функция ф(х) в окрестности точки х0 является решением системы уравнений в частных производных (4), причем выполнено условие (7). Приведенные условия есть условия преобразования аффинной системы (1) к регулярному квазаканоническому виду (2) с индексом приводимости р.
Относительная степень в точке х0 может быть неопределена, если не существует такое р < п, что выполнены условия (4) и (7).
Предположим, что для системы (1) задан выход
У = Ф(х), (8)
относительная степень которого в точке х0 равна р. Тогда в силу выполнения условия регулярности (7) система (1) может быть локально преобразована к специальному квазиканоническому виду (6), с которым связана система
¿1 = ¿2, - - - , ¿р-1 = ¿р,
, (9)
¿р = / (^,п) + g(¿,п)u,
П = д^, п), У = ¿1-
Систему (9) называют [11,13] нормальной формой системы (1) с выходом (8).
В случае, если замена переменных (¿, п) = Ф(х), задающая преобразование системы (1) с выходом (8) к нормальной форме (9), определена на всем Кп, и коэффициент при управлении в нормальной форме (9) отличен от нуля при всех (¿, п) е Кп, то такую нормальную форму называют глобально определенной.
В некоторых случаях нахождение таких переменных п, что управление не входит в последние уравнения нормальной формы, является трудоёмким или невозможно в аналитическом виде. Для практических целей часто бывает достаточно записать систему в виде
¿1 = ¿2, . . . , ¿р-1 = ¿р,
(10)
¿р = / (>,п) + g(¿,п)u,
п = д^, п) + п)и, У = ¿1-
В этом случае также говорят, что система записана в нормальной форме [13].
Таким образом, задача преобразования аффинной системы (1) с заданным выходом (8) к нормальной форме (9) является частным случаем задачи о преобразовании аффинной системы к квазиканоническому виду, поскольку в последней задаче требуется найти функцию ф, вводящую новые переменные ¿.
Отметим, что при поиске преобразования к квазиканоническому виду обычно ищут такую функцию ф, которой соответствует максимальное значение р, а при преобразовании к нормальной форме выясняют, определена ли относительная степень системы с заданным выходом в исследуемой точке и какова эта степень.
Любую аффинную систему со скалярным управлением можно рассматривать как систему, записанную в некоторой нормальной форме в окрестности точки х0, где в качестве виртуального выхода выступает такая переменная состояния, что соответствующее этой переменной уравнение системы (производная виртуального выхода в силу системы) содержит управление и коэффициент при управлении отличен от нуля в точке х0. Система с таким виртуальным выходом имеет в точке х0 относительную степень 1.
В дальнейшем будем для системы (1) с некоторым виртуальным выходом использовать терминологию теории нормальной формы.
Поскольку для фиксированной аффинной системы могут рассматриваться несколько виртуальных выходов, удобнее говорить не об относительной степени системы с заданным выходом в некоторой точке, а об относительной степени фиксированного выхода системы.
3. Стабилизация минимально фазовой системы
Не нарушая общности, рассмотрим задачу стабилизации положения равновесия х0 = 0 системы (1). Пусть для нее найден виртуальный выход у = ф(х), ф(0) = 0, относительная степень которого в точке х0 = 0 равна р и р < п. Тогда в окрестности этой точки существует такая гладкая невырожденная замена переменных
г = Ф(х), п = Ф(х), г £ Кр, п £ Кп—р, (11)
где Ф(х) = (ф(х), ЬАф(х),..., ЬА—1ф(х))т, Ф(0) = 0, Ф(0) = 0, что в переменных г, п аффинная система запишется в нормальной форме (9).
При решении задачи стабилизации нулевого положения системы, записанной в нормальной форме, важными являются свойства подсистемы относительно переменных п в (9) при ¿ = 0. Систему уравнений
п = д(0, п)
(12)
называют уравнениями нулевой динамики (нулевой динамикой). Если ее положение равновесия п = 0 асимптотически устойчиво, то аффинную систему (1) с виртуальным выходом У = ф(х) называют минимально фазовой (в точке х = 0) [13].
Введенное для аффинных стационарных систем со скалярным управлением и скалярным выходом понятие нулевой динамики [30] аналогично указанному свойству для линейных систем [33].
Относительная степень выхода [31,32] аффинной системы в случае линейной системы со скалярным управлением и выходом равна разности степеней знаменателя и числителя передаточной функции системы.
Пусть в положении равновесия х0 = 0 система (1) является минимально фазовой. Выбрем для системы (1) управление в виде
р-1
и = ( -ЬАФ(х) - Е СкЬкАф(х)) /ЬвЬА-1ф(х) ,
к=0
(13)
где коэффициенты Ск, к = 0, - - -, р — 1, выбраны так, что корни уравнения Хр + р=0 СкХк = 0 имеют отрицательные действительные части. В переменных ¿, п управление (13) имеет вид
и=
р-1
I^,п) - Е Ск¿к+1) Iд^,п)
Тогда система (9) примет вид
где
А=
0
00
к=0
¿ = A¿, п = п),
0 --0 ---
(14)
(15)
0
0
0
V
С0 -С1 -С2 - - - Ср—2 Ср—1
р-1 у
(16)
1
1
причем матрица А имеет собственные числа со строго отрицательной действительной частью.
Управление (14) стабилизирует систему (9) по части переменных г. При этом стабилизируемость по переменным п зависит от свойств нулевой динамики. При асимптотически устойчивой нулевой динамике положение равновесия (г, п) = (0,0) является асимтотически устойчивым [13].
Действительно, система (9), замкнутая управлением (14), имеет специальный каскадный вид (15). Условия асимптотической устойчивости нулевого положения равновесия таких систем задает следующая теорема.
Теорема 1. [13,29] Система (15) асимптотически устойчива в нулевом положении равновесия (¿,п) = (0,0), если функция д(г,п) непрерывно дифференцируема в окрестности точки (г, п) = (0, 0), д(0,0) = 0, линейная система
г = Аг, (17)
асимптотически устойчива, а система
п = д(0,п), (18)
асимптотически устойчива в точке п = 0.
В случае, если нулевая динамика (18) не является асимптотически устойчивой, проблема стабилизации может быть решена в том случае, если будет найден такой новый (виртуальный) выход системы (1), что соответствующая ему нулевая динамика будет асимптотически устойчива. Тогда для новой нормальной формы, соответствующей этому выходу, можно построить стабилизирующую обратную связь вида (14) и затем переписать ее в исходные переменные.
4. Метод виртуальных выходов
Для аффинной системы (1), считая х0 = 0, фиксируем некоторый выход у = ^(х), ^(0) = 0, при котором относительная степень системы (1) в точке х = 0 равна 1. Запишем систему (1) с этим выходом в соответствующей
нормальной форме
г = / (г,п)+ д (г,пК (19)
п = ^(¿,П ), (20)
У = ¿-
В (19)-(20) ¿ е К1, п = (nl,...,nn-l) е Кп-1, I (0,0) = 0, д(0,0) = 0, д(0,0) = 0, ¿ = Н(х), п = Ф(х), Ф(0) = 0.
Теорема 2. [25] Для того, чтобы аффинная система (1) имела виртуальный выход с относительной степенью р = 1 в точке х = 0 и асимптотически устойчивой нулевой динамикой, необходимо и достаточно, чтобы положение равновесия п = 0 нелинейной системы
п = я (у,п) (21)
с управлением V было стабилизируемо гладкой обратной связью V = V (п). Каждой такой стабилизирующей обратной связи в системе (21) соответствует виртуальный выход ф = ¿ - v(n) = Н(х) - V(Ф(х)) аффинной системы (1) относительной степени р = 1 в точке х = 0 и асимптотически устойчивая нулевая динамика.
Для аффинной системы (1) фиксируем некоторый выход у = Н(х), при котором относительная степень системы (1) в точке х = 0 равна 2. Запишем систему (1) с указанным выходом в соответствующей нормальной форме
¿1 = ¿2, ¿2 = I(¿, п) + д(¿, п)и, (22)
п = я^п)- (23)
В (22)-(23) ¿ = (^,¿2^ е К2, п = (пъ ■ ■ ■ ,пп-22)т е Кп-2, I(0,0) = 0, д(0,0) = 0, я(0, 0) = 0, ¿1 = Н(х), ¿2 = Ькк(х), п = Ъ(х), Ф(0) = 0.
Теорема 3. [25] Пусть система (1) с виртуальным выходом ф, ф = 0,
х=0
имеет в точке х = 0 относительную степень р = 2, а нулевая динамика асимптотически устойчива. Если в переменных ¿,п нормальной формы (22)-(23) ф>г = 0, то существуют функции v1(n), v2(n), Vi(0) = 0, г = 1, 2, стаби-
г=0,п=0
лизирующее положение равновесия п = 0 системы
п = я ^1^2, п) (24)
с управлениями V!, v2, причем
^(п)
¿г
= V2(n). (25)
г]=д(у1(п),у2(п),п)
Теорема 4. [25] Пусть управления v1 = vi(n), v2 = v2(n) стабилизируют положение равновесия п = 0 системы (24) и удовлетворяет условию (25). Если система (22)-(23) с виртуальным выходом 0(z, п) = z1 — v1(n) имеет относительную степень р = 2 в точке (z, п) = 0, то нулевая динамика, соответствующая виртуальному выходу ф, асимптотически устойчива в точке п = 0.
5. Линеаризация нулевой динамики
Успешность применения метода виртуальных выходов к динамической системе с неустойчивой нулевой динамикой зависит от возможности построения виртуальной обратной связи v^), обеспечивающей асимптотическую стабилизацию положения равновесия системы п = Нахождение такого управления часто является достаточно сложной задачей. Например, для системы
г = V + 2 (26)
п = п + sin z + п£ ,
положение равновесия уравнения нулевой динамики п = п неустойчиво, а поиск обратной связи v^), стабилизирующей положение равновесия стационарной системы
п = п + sin v + пv2, (27)
представляет существенную проблему.
Однако для линеаризованного в окрестности нулевого положения равновесия второго уравнения системы (26)
п = п + z
задача построения стабилизирующей обратной связи легко решается. Например, при v(п) = — (c1 + 1)п, c1 > 0, получим т) = —с1п, т. е. указанная обратная связь асимптотически стабилизирует нулевое положение равновесия.
Теорема 5. Если для системы п = положение равновесия п = 0
линеаризованной системы
п = Ап + Bz, (28)
где А = дяЫ
где А = дп
B = dq(z,n) (z,n)=(0,0)' dz
, с виртуальным управлением
(г,п)=(0,0)
г = V стабилизируется линейной обратной связью V = — Кп, то каждой такой
стабилизирующей обратной связи в системе (28) соответствует виртуальный выход у = ¿ + Кп = Н(х) + КФ(х) аффинной системы (1) относительной степени р = 1 в точке х = 0 и асимптотически устойчивая нулевая динамика.
М Пусть обратная связь V = -Кп стабилизирует положение равновесия п = 0 системы (28). Рассмотрим для аффинной системы (1) виртуальный выход, который в переменных нормальной формы (19)-(20) имеет вид ф(¿, п) = ¿ + Кп. В точке (¿, п) = (0, 0) относительная степень этого виртуального выхода равна 1, так как
ф(¿, п) = I(¿, п) + Кд(¿, п) + д(¿, п)и
и в точке (0,0) коэффициент при управлении д(0, 0) = 0. Соотношения
^ = ф^,п) = ¿ + Kn, (29)
п = п,
задают глобальную замену переменных, в которых система (19)-(20) записывается в нормальной форме, соответствующей виртуальному выходу ф
(30)
¿ = / (¿,n) + д (¿,n)u,
п = - Кц,г}),
где _
!(¿^ = ^,п) + КЯ^,п)) \г=г-Кп,п=п , д^^ = д^,п) \г=г-Кп,П=П -Полагая в (30) £ = 0, получаем систему уравнений нулевой динамики
п = Я^Кп,!)- (31)
Рассмотрим линеаризацию уравнений нулевой динамики в окрестности точки п = 0. Получим
л. дя (т,п) п = -
дт
Кп +
дп
• п
т=0,п=0
т=0,п=0
После замены п = п, т = ¿ полученная система совпадает с системой (28), замкнутой стабилизирующей обратной связью ¿ = V = -Кп. Следовательно, положение равновесия п = 0 системы (31) асимптотически устойчиво по первому приближению.
Поскольку замена переменных (29) является линейной, то определена невырожденная замена переменных, задающая преобразование системы (1) к виду (30), который является нормальной формой аффинной системы (1) с выходом ф(г,п) = г + Кп. Этот выход в исходных переменных аффинной системы (1) записывается в виде ф = ^(х) + КФ(х). ►
Одной из особенностей применения метода виртуальных выходов при р = 2 является поиск стабилизирующих обратных связей v1(n) и V2(п), которые связаны соотношением (25). Поиск таких стабилизирующих управлений из-за соотношений (25) является достаточно сложной задачей. Покажем, что для ее решения можно воспользоваться линеаризованными уравнениями.
Теорема 6. Пусть линейное приближение подсистемы п = ^(¿,п) системы (22)-(23) в окрестности точки (г, п) = (0,0) имеет вид
п = Ап + #1*1 + Д^, (32)
где А = Щ
г=0,п=0
в = И
В2 = I1
г=0,п=0
>2 дх2
, #1 = 0, #2 = 0, п £
2, г £ К2, и К = (к1,..., кп—2) £ 2 — такая строка, что система (22)-(23) с виртуальным выходом ф(г, п) = г1 + Кп имеет относительную степень 2 в точке (г, п) = 0, а управления
v1 = —Кп,
—КА + КВ1К (33)
V2 = 1 + КВ2 п, стабилизируют положение равновесия п = 0 системы
пп = Ап + Вт + (34)
Тогда нулевая динамика, соответствующая виртуальному выходу ф, асимптотически устойчива в точке п = 0.
< Покажем, что управление v2 определено корректно и
1 + К#2 = 0. (35)
Поскольку относительная степень виртуального выхода ф(г, п) в точке (г, п) = = 0 равна 2, то значение в этой точке коэффициента
д(г,п) = д(г,п)(1 + К4 (г1, г2, п))
при управлении в
¿2ф(¿, п)
¿г2 (22)-(23)
отлично от нуля. Учитывая, что В2 = д^ ство (35) справедливо. Соотношения
, можно видеть, что неравен-
г=0,п=0
- = Аф(г,Г!) ¿ 2 =
¿ 1 = ф(¿,п) = ¿1 + Кп,
= ¿2 + Кд^, ¿2, п), (36)
(22)-(23)
п = п
задают невырожденную замену переменных в окрестности точки (¿, п) =
д(г1,г2)
1 + Кдг2 (¿l,¿2,n)\
г=0,п=0 1 + КВ2
(0,0), поскольку определитель матрицы Якоби д*1^2] в этой точке равен
и также отличен от нуля.
Пусть обратная для (36) замена имеет вид
¿1 = ¿1 - Кп, ¿2 = ^(¿1 ^Ц), п = п-
В переменных (36) система (22)-(23) с виртуальным выходом ф запишется в нормальной форме
¿1 = ¿2, ¿2 = 7(¿, п) + д(¿, п)и, (37)
п = я(¿1 - Кп, ^2(х 1, £ 2, п), п)- (38)
Полагая — 1 =0, ¿2 = 0 в (38) и сделав замену п = п, получим уравнения соответствующей нулевой динамики
п = Я(-Кп,ы2(0,0,п),п)- (39)
Докажем ее ассимптотическую устойчивость в точке п = 0. Линейное приближение уравнения (39) имеет вид
п =
дд(г, п)
дп
п +
+ дд(г,п) дг2
= Ап — В1Кп + В2
дд(г, п)
дг1
(—Кп) +
г=0,п=0
д^(0,0, п)
г=0,п=0 дп
д^(0,0,0)
п=
(40)
П=0
дп
п.
Поскольку г1 = г2, то г 1 — Кп = 1, г2, п), что при г 1 = 0, г2 = 0 после замены п = п дает соотношение
— Кд(—Кп, эд2(0,0, п), п) = ^2(0,0, п).
(41)
Из (41) получим
д^(0,0, п) дп
_к( дд(г1,г2,п) дд(г1,г2,п) к + \ дп дд(г1, г2, п) ди>2(0,0, п)
дг2
дп
откуда
д^(0,0, п) дп
_К д^ь^п) + к ,п)
дг-1
К
1 + К
дч(г1,г2,п)
дг 2
Х1=~Кп,г2 =W2(0,0,п)
При п = 0 получим
д^2(0, 0,0) —КА + КВ1К
дп
1 + КВ2
Подставляя полученное выражение в (40), получим систему
• , „ ^ ^ —КА + КДК
п = Ап — #1Кп + #2—1 + Кв—п.
которая совпадает с системой (34), замкнутой управлениями
—КА + КВ1К Vl = —Кп, ^ = 1 + КВ п.
(42)
(43)
По условиям теоремы система (43) асимптотически устойчива.
Следовательно, нулевая динамика системы (37)-(38) асимптотически устойчива в точке п = 0 по теореме Ляпунова об асимптотической устойчивости по первому приближению. ►
Отметим, что для нахождения виртуального выхода с устойчивой нулевой динамикой использовалась линеаризация только нулевой динамики исходной системы. Использование полной линеаризации системы (22)-(23) возможно, но в данном случае нецелесообразно, поскольку управление вида (14), построенное для системы (37)-(38) и записанное в переменных ¿, будет стабилизировать нормальную форму (22)-(23 по части переменных ¿ в большей области, чем управление, построенное с использованием полностью линеаризованной системы.
Укажем способ выбора строки К из теоремы 6 в случае п = 4. В этом случае система (32) имеет вид
п = Ап + В^1 + B2V2, (44)
где п ^ А — квадратная матрица второго порядка, а В1 = (Ъ11,Ъ21)т, В2 = (Ъ21М)\ К =(к1,к2).
Характеристический многочлен матрицы
А - В1К + В™
замкнутой системы (44), (33) имеет вид Л2 + &1Л + &0, где
-Ьбз - к2^4 + $0
&0 = & =
кфи + к2Ь22 + 1 '
ЫЬп - $1) + к2(Ь21 - 62) - (й.11 + (122) к1Ъу2 + к2Ь22 + 1
(46)
а
$0 = ¿е! А, 61 = (22Ъ12 - Ъ22а12,
$2 = 111Ъ22 - Ъ121.21, 63 = Ъпа22 - (>^21, $4 = ацЪ21 - Ъ11а2ъ Элементы матрицы К, если они существуют, должны быть выбраны так, чтобы &0 > 0, &1 > 0, что обеспечит гурвицевость матрицы (45). Дополнительное условие 1 + КВ2 = 0 имеет вид к1Ъ21 + к2Ъ22 + 1 = 0. Эти условия на
плоскости Ок1к2 задают область, точку в которой несложно выбрать в случае ее непустоты.
В случае, если значения к1 и к2, удовлетворяющие указанным неравенствам, не существуют, найти виртуальный выход с использованием теоремы 6 не удается.
Рассмотрим задачу стабилизации верхнего неустойчивого положения равновесия перевернутого маятника, установленного на тележке (см. рис. 1). Управляемым параметром является сила ], прикладываемая к тележке. Желательно, чтобы одновременно с маятником стабилизировалось и положение тележки в точке х = 0, х = 0.
Для решения указанной задачи разработаны различные алгоритмы [34]. Ограничимся случаем, когда |$| < |. Пусть тележка имеет массу М, маятник — массу т, а центр масс маятника расположен на расстоянии I от точки подвеса. Для упрощения формул примем момент инерции маятника относительно его центра масс равным 0. В этом случае модель системы имеет вид [34]
6. Стабилизация маятника на тележке
/
о
Рис. 1
тд сое $ Бт $ — т1 Бт $$2 + /
х =
М + т вт2($) (М + т)д Бт $ — т1 сое $ Бт $$2 + / сое $
(47)
(М + т вт2 $)1
Преобразуем систему (47) к нормальной форме относительно виртуального выхода y = Положим zi = z2 = $. Относительно введенных переменных в области < 2 относительная степень выхода равна 2. Найдем первые интегралы векторного поля D при управлении f,
^ 1 д cos $ д
D = —-г^тгг +
(48)
M + m sin2 $dx (M + m sin2 $)/$$'
Решив уравнение Ldф = 0, получим ф\ = $, ф2 = x и ф3 = /$ — xcos $. Примем ni = x, п2 = 1$ — x cos $.
Можно видеть, что замена (z, n) = Ф($, $, x, x)
zi = $, z2 = $, n = x, n2 = 1$ — x cos $
гладкая и невырождена в рассматриваемой области |$| < |, а обратная замена имеет вид
о ó 1z2 — П2
$ = zi, $ = z2, x = ni, x =-.
cos zi
В найденных переменных уравнения для ni, П2 имеют вид
1z2 — П2
ni =-,
cos zi
/ 1z2 — n2 \ •
n2 = g +--z2 sin zi.
cos zi
и дают следующие уравнения нулевой динамики
ni = —n2,
n2 = 0.
Нулевое положение равновесия этой системы не является асимптотически устойчивым.
Будем искать новый виртуальный выход с асимптотически устойчивой нулевой динамикой. Согласно теореме 4, для этого надо найти управления vi (n) и v2(n), стабилизирующие положение равновесия системы
1V2 — n2
ni =-,
cos vi (50)
( , 1v2 — n2 \ .
n2 = g +--z2 sin Vi
V cos vi J
(49)
и связанные соотношением
3vi lV2 - П2 , dvi ( lV2 - П2 \ .
----|- — g +--V2 Sin Vi — V2.
ОЩ cos Vl ОЩ V cos Vl )
Поиск таких управлений представляется достаточно сложной задачей. Для упрощения задачи применим теорему 6.
Линеаризованные в окрестности нулевого положения уравнения (48) после замены ¿1 = v1, ¿2 = v2 имеют вид
т = -V2 +lV2, (51)
П2 — gVi.
Используя соотношения (46), получим
7 —ki д , к2д
do — 1—1—г, di —
kil + V 1 kil + 1'
Из условий d0 > 0, dl > 0, kll + 1 — 0 вытекает, что k2 > 0 и — \ < kl < 0 или
k- < -\ и k2 < 0.
lo > 0, d- > 0, k-l + 1 — 0 вытекает, что k2 > 0 и — j < k-
¿I < — j и k2
Задача нахождения стабилизирующих обратных связей для линеризованных уравнений нулевой динамики решена. Виртуальный выход, обеспечивающий асимптотическую устойчивость нулевой динамики, имеет вид
y — $ + klnl + k2n2 — $ + klx + k2(l$ — x cos $).
Построим теперь стабилизирующее управление для системы (47), для чего найдем нормальную форму относительно найденного виртуального выхода.
Заметим, что для построения стабилизирующего управления достаточно получить нормальную форму в виде заготовки, то есть правую часть уравнений записать в переменных z,rj.
Вычисления удобно проводить с использованием пакета компьютерной алгебры. Ниже приведены с небольшими изменениями результаты, полученные
с использованием пакета Maple:
zi = Z 2,
- fki(ÍZ2 - n2)sin zi k2(lz2 - П2 )z2 sinzA
z2 = -2--1--2- z2 +
\ cos2 z1 cos2 z1 J
+ h(g cos zi + (lz2 - n2)z2)z2 +
+ (i + AL + k2 () sin z,) x (52)
cos zi cos zi
(M + m) sin z1g - mlz| cos z1 sin z1
X
(M + m sin2 z1)l
ki k2z2 sin zi\ ( (lz2 - V2)z2\ .
+--g +--sin zi +
cos zi cos zi cos zi
cos zi + k\l + k2(2lz2 - n2) sin zi — _
+-HTL--u = f + gu
(M + m sin z\)
где zi = y, а
- .j lz2 - П2 . j { . lz2 - П2 \ .
z2 = z2 + ki--+ k2\ g +--zc2 sin zi.
cos zi cos zi
Используя (52), получим стабилизирующее управление в виде
-f - C\z 2 - Coz i
u=
g
Замена переменных (¿, п) = Ф(Щ, Щ, х, Х) позволяет это управление записать в исходных переменных. Поскольку замена задает диффеоморфизм в области Щ < П, то локальная асимптотическая устойчивость сохранится и в исходных переменных.
Численное моделирование выполнялось при М = 2.0, т = 0.1, I = 0.5, Щ(0) = 0, Щ(0) = -0.1, х(0) = 0.0, Х(0) = 0.0, к1 = -1.5, к2 = 0.25, С1 = 0.75, С0 = 0.125.
На рис. 2, рис. 3 приведены графики зависимостей Щ(Щ) и Х(х), полученных интегрированием замкнутой системы при Ь Е [0, Ьк], Ьк = 25 ^
Рис.2
0.5 -0.4 -0.3 -0.2 -0.1 0
X, M
Рис. З
7. Заключение
Поиск виртуального выхода, которому соответствует асимптотически устойчивая нулевая динамика, по-существу является поиском такого многообразия в пространстве переменных состояния системы, при движении по которому траектории динамической системы, замкнутой соответствующим управление, стремятся к положению равновесия. Управление здесь используется только для удержания траектории системы в достаточно малой окрестности указанного многообразия.
Такое управление может быть предпочтительнее управления, полученного с использованием полной линеаризации динамической системы с управлением, поскольку базируется на внутренних свойствах системы.
Работа выполнена при финансовой поддержке гранта РФФИ 10-07-00468 и Программы Президента РФ по государственной поддержке ведущих научных школ (грант НШ-4144.2010.1)
Список литературы
1. Жевнин A.A., Крищенко А.П. Управляемость нелинейных систем и синтез алгоритмов управления // Докл. АН СССР. 1981. Т. 258, № 4. С. 805 - 809.
2. Крищенко А.П. Исследование управляемости и множеств достижимости нелинейных систем управления // Автоматика и телемеханика. 1984. № 6. С. 30-36.
3. Крищенко А.П. Стабилизация программных движений нелинейных систем // Изв. АН СССР. Техническая кибернетика. 1985. № 6. С. 103 -112.
4. Крищенко А.П. Преобразование нелинейных систем и стабилизация программных движений // Труды МВТУ им. Н.Э. Баумана. 1988. № 512. С. 69 -87.
5. Крищенко А.П. Преобразование многомерных аффинных управляемых систем // Управляемые нелинейные системы. 1991. № 2. С. 5-14.
6. Краснощеченко В.И., Крищенко А.П. Нелинейные системы: геометрические методы анализа и синтеза. М.: Изд-во МГТУ им. Н.Э. Баумана, 2005. 520 с.
7. Brockett R.W. Feedback invariants for nonlinear systems // Preprints of VII World Congress IFAC. Oxford: Pergamon Press, 1978. V. 2. P. 1115-1120.
8. Su R. On the linear equivalents of nonlinear systems // Syst.& Cont. Letters. 1982. № 1. P. 48-52.
9. Hunt L.R., Su R., Meyer G. Design for multiinput nonlinear systems // Diff. Geometric Control Theory. Boston: BirkHauser, 1983. P. 24-30.
10. Jacubczyk B., Respondek W. On linearization of control systems // Bul. L'acad Pol. Sciense. 1980. V. 28., № 9-10. P. 517-522.
11. Byrnes C., Isidori A. Local stabilization of minimum phase nonlinear system // Syst. Contr. Lett. 1988. V. 11. P. 9-17.
12. Byrnes C., Isidori A. Asymptotic stabilization of minimum phase nonlinear systems // IEEE Trans. Automat. Contr. 1991. № 10. P. 1122-1137.
13. Isidori A. Nonlinear control systems. London: Springer-Verlag, 1995. 587 p.
14. Byrnes C., Isidori A. New results and examples in nonlinear feedback stabilization // Syst. Contr. Lett. 1989. V. 12. P. 437-442.
15. Isidori A., Moog C., De Luca A. A sufficient condition for full linearization via dynamic feedback // Proc. of the 25th IEEE Conf. Decision Contr. Athlanta, 1986. P. 203-208.
16. Teel A.R. Semi-global stabilization of minimum-phase nonlinear systems in special normal forms // Syst. Cont. Letters. 1992. № 3. P. 187-192.
17. Byrnes C., Isidori A., Willems J. Passivity, feedback equivalence and the global stabilization of minimum phase nonlinear systems // IEEE Trans. Automat. Contr. 1992. V. 37. P. 1004-1017.
18. Allgower F., Doyle III F.J. Approximate input-output linearization of nonlinear systems // Nonlinear model-based process control. Ed. by Barber R., Kravaric C. Dordrecht: Kluwer Academic Publishers, 1998. P. 235-274.
19. Tornambe A. Output feedback stabilization of a class of non-minimum phase nonlinear systems // Syst.& Cont. Letters. 1992. № 3. P. 193-204.
20. Zou Q., Devasia S. Precision preview-based stable-inversion for nonlinear nonminimum-phase systems: The VTOL example // Automatica. 2007. № 1. P. 117-127.
21. Hauser J., Sastry S., Meyer G. Nonlinear control design for slightly nonminimum phase systems: application to V/STOL aircraft // Automatica. 1992. V. 28, №4. P. 665-679.
22. Wright R. A., Kravaris C. Nonminimum-phase compensation for nonlinear processes // AIChE Journal. 1992. V. 38, № 1. P. 26-40.
23. Niemiec P. M., Kravaris C. Nonlinear model-state feedback control for nonminimum-phase processes // Automatica. 2003. V. 39, № 7. P. 1295-1302.
24. Kravaris C., Mousavere D. ISE-optimal nonminimum-phase compensation for nonlinear processes// Journal of Process Control. 2007. V. 17, № 5. P. 453-461.
25. Крищенко А.П., Панфилов Д.Ю., Ткачев С.Б. Построение минимально фазовых аффинных систем // Дифференциальные уравнения. 2002. Т. 38, № 11. С. 1483-1489.
26. Крищенко А.П., Панфилов Д.Ю., Ткачев С.Б. Глобальная стабилизация аффинных систем с помощью виртуальных выходов // Дифференциальные уравнения. 2003. Т. 39, №11. С. 1503-1510.
27. Output maps with associated asymptotically stable zero dynamics / A.P. Krishchenko, D.U. Panfilov, K.E. Starkov, S.B. Tkachev // Nonlinear Control Systems'04: Proc. of VIIFAC Symp. Stuttgart, 2004. V. 1. P. 329-334.
28. Ткачев С.Б. Стабилизация нестационарных аффинных систем методом виртуальных выходов // Дифференциальные уравнения. 2007. Т. 43, № 11. C. 1507-1517.
29. Khalil H.K. Nonlinear systems. 3-d edition. Upper Saddle River: Prentice-Hall, 2002. 750 p.
30. Krener A., Isidori A. Nonlinear zero distributions // Proc. of the 19th IEEE Conf. Decision Contr. Athlanta, 1980. P. 665-668.
31. Byrnes C., Isidori A. A frequency domain philosofy for nonlinear system with application to stabilization and to adaptive control // Proc. of the 23rd IEEE Conf. Decision Contr. Las Vegas, 1985. P. 1031 -1037.
32. Isidori A., Moog C. On the nonlinear equivalent of the notion of transmission zeros // Modeling and Adaptive Control. New York, 1988. P. 203-208.
33. Ким Д.П. Теория автоматического управления. В 2-х томах. Многомерные,
нелинейные и адаптивные системы. М.: ФИЗМАТЛИТ, 2004. Т.2. 464 с.
34. Fantoni I., Lozano R. Non-linear control for underactuated mechanical systems. London: Springer-Verlag, 2002. 295 p.
electronic scientific and t echnical periodical
SCIENCE and EDUCATION
El № FS77 - 30569. №0421100025. ISSN 1994-0408
Stabilization of nonminimum-phase affine systems using a partial linearization
77-30569/255087
# 10, October 2011 S. B. Tkachev
Bauman Moscow State Technical University
For nonlinear dynamical system with scalar control the stabilization problem of an equilibrium point on a basis of the differential-geometric approach is considered. The basic theoretical data on transformation of a dynamical control system in a quasicanonical form and of a dynamic control system with output in a normal form are reviewed. The basic attention is given to the case when zero dynamics of the system isn't asymptotically stable, i.e. the nonlinear system isn't minimumphase. Application for stabilization of such systems of the virtual outputs method not always possible because of difficulties of search of virtual controls with desired properties. It is shown that for construction of virtual outputs with desired properties the linearization of the subsystem defining zero dynamics can be used. For the dynamic system describing the turned pendulum on the cart, the static feedback is designed using of the virtual output method and partial linearization. Designed feedback stabilizes simultaneously the top pendulum equilibrium point and the preset position of the cart.
References
1. Zhevnin A.A., Krishchenko A.P. Upravlyaemost' nelinejnyh sistem i sintez algoritmov upravleniya // Dokl. AN SSSR. 1981. T. 258, N 4. S. 805-809.
2. Krishchenko A.P. Issledovanie upravlyaemosti i mnozhestv dostizhimosti nelinejnyh sistem upravleniya // Avtomatika i telemehanika. 1984. N 6. S. 30-36.
3. Krishchenko A.P. Ctabilizaciya programmnyh dvizhenij nelinejnyh sistem // Izv. AN SSSR. Tehnicheskaya kibernetika. 1985. N 6. S. 103-112.
4. Krishchenko A.P. Preobrazovanie nelinejnyh sistem i stabilizaciya programmnyh dvizhenij // Trudy MVTU im. N.Je. Baumana. 1988. N 512. S. 69-87.
5. Krishchenko A.P. Preobrazovanie mnogomernyh affinnyh upravlyaemyh sistem // Upravlyaemye nelinejnye sistemy. 1991. N 2. C. 5-14.
6. Krasnoshchechenko V.I., Krishchenko A.P. Nelinejnye sistemy: geometrich-eskie metody analiza i sinteza. M.: Izd-vo MGTU im. N.Je. Baumana, 2005. 520 s.
7. Brockett R.W. Feedback invariants for nonlinear systems // Preprints of VII World Congress IF AC. Oxford: Pergamon Press, 1978. V. 2. P. 1115-1120.
8. Su R. On the linear equivalents of nonlinear systems // Syst.& Cont. Letters. 1982. N1. P. 48-52.
9. Hunt L.R., Su R., Meyer G. Design for multiinput nonlinear systems // Diff. Geometric Control Theory. Boston: BirkHauser, 1983. P. 24-30.
10. Jacubczyk B., Respondek W. On linearization of control systems // Bul. L'acad Pol. Sciense. 1980. V. 28., N 9-10. P. 517-522.
11. Byrnes C., Isidori A. Local stabilization of minimum phase nonlinear system // Syst. Contr. Lett. 1988. V. 11. P. 9-17.
12. Byrnes C., Isidori A. Asymptotic stabilization of minimum phase nonlinear systems // IEEE Trans. Automat. Contr. 1991. N 10. P. 1122-1137.
13. Isidori A. Nonlinear control systems. London: Springer-Verlag, 1995. 587 p.
14. Byrnes C., Isidori A. New results and examples in nonlinear feedback stabilization // Syst. Contr. Lett. 1989. V. 12. P. 437-442.
15. Isidori A., Moog C., De Luca A. A sufficient condition for full linearization via dynamic feedback // Proc. of the 25th IEEE Conf. Decision Contr. Athlanta, 1986. P. 203-208.
16. Teel A.R. Semi-global stabilization of minimum-phase nonlinear systems in special normal forms // Syst. Cont. Letters. 1992. N 3. P. 187-192.
17. Byrnes C., Isidori A., Willems J. Passivity, feedback equivalence and the global stabilization of minimum phase nonlinear systems // IEEE Trans. Automat. Contr. 1992. V. 37. P. 1004-1017.
18. Allgower F., Doyle III F.J. Approximate input-output linearization of nonlinear systems // Nonlinear model-based process control. Ed. by Barber R., Kravaric C. Dordrecht: Kluwer Academic Publishers, 1998. P. 235-274.
19. Tornambe A. Output feedback stabilization of a class of non-minimum phase nonlinear systems // Syst.& Cont. Letters. 1992. N 3. P. 193-204.
20. Zou Q., Devasia S. Precision preview-based stable-inversion for nonlinear nonminimum-phase systems: The VTOL example // Automatica. 2007. N 1. P. 117-127.
21. Hauser J., Sastry S., Meyer G. Nonlinear control design for slightly nonminimum phase systems: application to V/STOL aircraft // Automatica. 1992. V. 28, N4. P. 665-679.
22. Wright R. A., Kravaris C. Nonminimum-phase compensation for nonlinear processes // AIChE Journal. 1992. V. 38, N 1. P. 26-40.
23. Niemiec P. M., Kravaris C. Nonlinear model-state feedback control for nonminimum-phase processes // Automatica. 2003. V. 39, N 7. P. 1295-1302.
24. Kravaris C., Mousavere D. ISE-optimal nonminimum-phase compensation for nonlinear processes // Journal of Process Control. 2007. V. 17, N 5. P. 453-461.
25. Krishchenko A.P., Panfilov D.Ju., Tkachev S.B. Postroenie minimal'no fa-zovyh affinnyh sistem. Differencial'nye uravneniya. 2002. T. 38, N 11. S. 1483-1489.
26. Krishchenko A.P., Panfilov D.Ju., Tkachev S.B. Global'naya stabilizaciya affinnyh sistem s pomoshch'ju virtual'nyh vyhodov // Differencial'nye uravneniya. 2003. T. 39, N 11. S. 1503-1510.
27. Output maps with associated asymptotically stable zero dynamics / A.P. Krishchenko, D.U. Panfilov, K.E. Starkov, S.B. Tkachev // Nonlinear Control Systems'04: Proc. of VI IFAC Symp. Stuttgart, 2004. V. 1. P. 329-334.
28. Tkachev S.B. Stabilizaciyanestacionarnyh affinnyh sistem metodom virtual'nyh vyhodov//Differencial'nye uravneniya. 2007. T. 43, N 11. C. 1507-1517.
29. Khalil H.K. Nonlinear systems. 3-d edition. Upper Saddle River: Prentice-Hall, 2002. 750 p.
30. Krener A., Isidori A. Nonlinear zero distributions // Proc. of the 19th IEEE Conf. Decision Contr. Athlanta, 1980. P. 665-668.
31. Byrnes C., Isidori A. A frequency domain philosofy for nonlinear system with application to stabilization and to adaptive control // Proc. of the 23rd IEEE Conf. Decision Contr. Las Vegas, 1985. P. 1031-1037.
32. Isidori A., Moog C. On the nonlinear equivalent of the notion of transmission zeros // Modeling and Adaptive Control. New York, 1988. P. 203-208.
33. Kim D.P. Teoriya avtomaticheskogo upravleniya. V 2-x tomah. Mnogomernye, nelinejnye i adaptivnye sistemy. M.: FIZMATLIT, 2004. T.2. 464 s.
34. Fantoni I., Lozano R. Non-linear control for underactuated mechanical systems. London: Springer-Verlag, 2002. 295 p.
