CC BY
21
электронное научно-техническое издание
НАУКА и ОБРАЗОВАНИЕ
_Эл № ФС77 - 30569. Государственная регистрация ^0421100025. ISSN 1994-0408_
Синтез управления плоской моделью корабля при неполном измерении вектора состояния
77-30569/233615
# 10, октябрь 2011 А. Е. Голубев
УДК 519.71
МГТУ им. Н.Э. Баумана [email protected]
1. Введение
При моделировании технических систем одной из существенных проблем является отсутствие полной информации о системе. Измеряются, как правило, значения только части параметров и переменных, описывающих состояние системы. Причины неполноты измеряемой информации могут быть различные: высокая стоимость установки датчиков, технологические ограничения и т.п.
Одним из способов восстановления полной информации о состоянии на основе имеющихся измерений является построение наблюдателя. Наблюдатель представляет собой специальную динамическую систему, состояние которой с течением времени достаточно быстро, например асимптотически, приближается к состоянию исследуемой системы. Состояние наблюдателя в произвольный момент времени можно рассматривать как оценку состояния системы в этот момент времени.
Оценка состояния системы наблюдателем может быть использована при синтезе стабилизирующих законов управления. Одним из подходов к решению задач стабилизации при неполном измерении состояния является принцип разделения. Основная идея принципа разделения заключается в раздельном построении стабилизирующей обратной связи по состоянию и наблюдателя с последующей подстановкой в обратную связь состояния наблюдателя вместо состояния системы.
Принцип разделения позволяет находить стабилизирующие управления для линейных динамических систем [1, 2]. Однако, для решения задач стабилизации нелинейных динамических систем в общем случае принцип разделения неприменим [3, 4]. Причина этого - возможный неограниченный рост за конечное время решений системы с управлением, полученным заменой в ста-билизирущей обратной связи по состоянию состояния системы на его оценку наблюдателем [3, 4]. Известны только некоторые классы нелинейных систем, для которых принцип разделения справедлив и позволяет находить решения задач стабилизации в условиях доступности измерениям только части вектора состояния [4].
В настоящей работе показано, что задача стабилизации заданного положения корабля на плоскости может быть решена с помощью раздельного построения стабилизирующей обратной связи по состоянию и наблюдателя с последующей подстановкой оценки состояния системы наблюдателем в обратную связь. Установлено, что для рассматриваемой нелинейной модели плоского движения корабля выполнены достаточные условия справедливости принципа разделения.
Рассмотрим уравнения движения корабля на плоскости, имеющие следующий вид [5, 6]:
где п = (х,у,ф)т Е К3, V = (и, v,т)т Е К3 образуют вектор состояния системы; т = (т1,т2,т3)т Е К3 - управление; А = —М—1В, В = М—1 - постоянные матрицы размера 3 х 3, матрица J(п) имеет вид
Здесь х, у - координаты центра масс корабля в фиксированной декартовой системе координат, связанной с Землей; ф - угол рысканья; и и V - соответственно продольная и поперечная составляющая скорости корабля в системе
2. Модель движения и постановка задачи
п = J (п^,
V = Av + Вт,
(1)
координат, жестко связанной с корпусом; г - угловая скорость вращения корабля относительно вертикальной оси; т - вектор управляющих сил и момента, создаваемых гребными винтами в кормовой части корабля и вспомогательными движителями; матрицы М = Мт > 0 и Б > 0 соответственно моментов инерции и коэффициентов демпфирования имеют вид
Шц 0 0 \ / ¿ц 0 0 \
М = 0 Ш22 Ш23 , Б = 0 ¿22 ¿23
V 0 Ш32 Ш33 / 0 ¿32 ¿33 /
Отметим, что математическая модель (1) описывает динамику корабля при малых скоростях, например, при швартовке.
Решим задачу приведения корабля из произвольного положения на плоскости в заданное. Будем предполагать, что измерениям доступны только координаты х, у центра масс корабля и угол рысканья ф, то есть измеряемый выход системы (1) имеет вид у = п. В качестве стабилизируемого положения корабля рассмотрим положение, при котором п = 0. Для решения задачи управления требуется построить закон управления в виде обратной связи, использующей значения только измеряемого выхода системы, глобально асимптотически стабилизирующий положение равновесия п = 0, V = 0, т = 0 системы (1).
3. Синтез обратной связи по состоянию и наблюдателя
Согласно результатам, полученным в работе [5], асимптотический наблюдатель, работающий глобально, для системы (1) с рассматриваемым выходом у = п имеет вид
п = 3 (п)£ + ¿1(п - п^ (2)
V = А> + ¿2 (п - п) + Вт,
где ¿1 = -Р-1Я1, ¿2 = -Р-1 ,т(п)Р1. Здесь Р = рт > 0 и Ql = > 0 -произвольные симметрические положительно определенные матрицы, Р2 = = Р2т > 0 и Q2 = Q2 > 0 - симметрические положительно определенные матрицы, удовлетворяющие уравнению Ляпунова
1 т
-(АтР2 + Р2А) = ^2. (3)
Далее будем предполагать, что матрица А = —М—1В имеет собственные числа только с отрицательными действительными частями. В этом случае при любой фиксированной матрице = ^2 > 0 решение Р2 = Р2 > 0 уравнения Ляпунова (3) существует и единственно.
В случае, когда данное предположение не выполнено, асимптотический наблюдатель, работающий глобально, для системы (1) можно найти, например, с помощью метода, рассмотренного в работе [7].
Отметим, что построение наблюдателя (2) основано на асимптотической стабилизации при помощи статической обратной связи по выходу системы, описывающей динамику ошибки оценки состояния системы (1) и имеющей вид
п= J (п)77 + ¿177,
7 = А7 + Ь2 п,
где п = п — п, 7 = О — V. Согласно работе [5] при указанном выборе матриц Ь1 и Ь2 положение равновесия п = 0, 7 = 0 системы (4) экспоненциально устойчиво в целом.
Найдем закон управления в виде обратной связи по состоянию, глобально асимптотически стабилизирующий положение равновесия п = 0, V = 0, т = 0 системы (1). Используем метод нелинейной стабилизации, рассмотренный в работе [8]. В переменных
п = п, ^ = J(пК (5)
где £ = (£ь £2, £3)т, при управлении
т = 7(п, V) = B—1 Jт(п)(—J(п)£(т^ — J(п)Av — Кщ — ) (6) система (1) без выхода примет вид
Xх = £ъ
у = ^
Ф = ^ (7)
£1 = —к}1х — х^ь
£2 = — К22У — к22£2 ,
£3 = — к13ф —
где к^ > 0, % = 1, 3, 3 = 1, 2, - произвольные положительные постоянные.
В выражении (6) матрицы К1, К2 и Б(т) имеют вид К1 = diag(к11, к22, к^3) К2 = ^^, к22, к33) и
Б (т) =
0 —т 0 т00 000
соответственно.
Положение равновесия п = 0, £ = 0 системы (7) экспоненциально устойчиво в целом. Так как соотношения (5) задают отображение, являющееся диффеоморфизмом пространств К6 = {п, V} и К6 = {п,£}, и справедливо равенство ||(пт, vт)т|| = ||(пт, £т)т||, то положение равновесия п = 0, V = 0 системы (1) с управлением (6) также экспоненциально устойчиво в целом. Тогда согласно работе [9] существует функция Ляпунова V (п, V) такая, что при всех (пт, vт)т Е К6 выполняются неравенства
С11(пт, Vт)т||2 < V(п, V) < С2||(пт, Vт)т||2,
^(п, V)
д(пт, vт)т
< С3|(пт, vт)т|,
дV(п, V)
^(п^) < —с4||(пт, V1')1
(8)
(9)
д(пт, vт)т
где ^(п, V) = ((J(п^)т, (Av + В7(п, V))т)т; с1, с2, с3, с4 - некоторые положи тельные константы. Здесь и далее в работе || • || - евклидова норма.
4. Применение принципа разделения
Рассмотрим нелинейную динамическую систему, имеющую общий вид
X = / (х, и),
у = Мх)
(10)
где х Е - вектор состояния системы, и Е Кт - вход (управление), у Е -измеряемый выход системы, ](•, •) и ) - непрерывно дифференцируемые функции своих аргументов, f (0,0) = 0, ^(0) = 0.
2
Обозначим через хы(£, х0) решение х = х(£) системы (10) с произвольной непрерывной и ограниченной на интервале [0, функцией и = и(£) на входе, удовлетворяющее начальному условию х(0) = х0. Здесь хи(0, х0) = х0
Предположение 1. 1) Существует непрерывно дифференцируемая обратная связь и = к(х), к(0) = 0, такая, что система (10) при управлении и = к(х) асимптотически устойчива в целом в точке х = 0; 2) V(х) > 0 -функция Ляпунова для замкнутой обратной связью и = к(х) системы (10), удовлетворяющая при всех х £ неравенствам
^(ЦхЦ) < V(х) < а2(||х||), (11)
^^/(х, к(х)) <-«3(|х|), (12)
где г = 1, 2, - некоторые функции класса а а3( ) - функция класса К.
Предположение 2. 1) Система (10) допускает построение наблюдателя
х = #(х,^(х),и), #(0, 0,0) = 0, (13)
уравнение ошибки е = х - х оценки которым состояния системы имеет вид
е = р(х„(^,х0) + е, й(х„(£,х0)),и(£)) - /(х„(£,х0),и(£)) =
= ^(хи(*,х0),е,и(*)), ^(хи(^,х0), 0,и(£)) = 0 У > 0; (14)
2) положение равновесия е = 0 системы (14) с произвольной непрерывной и ограниченной на интервале [0, функцией и(£) при произвольном решении хы(£, х0) системы (10) с и = и(£) асимптотически устойчиво в целом.
На основе результатов, приведенных в работе [4], сформулируем следующую теорему.
Теорема 1. Пусть для системы (10) выполнены предположения 1, 2 и при всех х £ е £ справедливы неравенства
(х - е, к(х)) - I(х, к(х))) < ^1(|х|) + ^2(|е|), (15) ^^^(х - е, е, к(х)) < ф1(||х||) + ф2(||е||), (16)
где ^ (•) и ф (•), % = 1, 2, - некоторые функции класса . Тогда при управлении и = к(х) система
X = f (х, и), , ч
х V (17)
хх = #(х, ^(х), и),
составленная из уравнений системы (10) и уравнений наблюдателя (13), асимптотически устойчива в целом в точке х = 0, х = 0, если ¿(й) = а3(й) — — — ф1(й), й Е [0, является функцией класса
Отметим, что для системы (1) с выходом у = п выполняются предположения 1 и 2. Покажем, что для системы (1) и наблюдателя (2) справедливы
. т л т \ т ^ / т т \ т
неравенства (15) и (16). Действительно, при всех (п , /> ) Е К6 и (п , 7 ) Е К6 имеем
дv(п,>) / J(п — п)^ — 7) — J(п)v д (пт, />т)тД —А 7 <
< С3|(пт,/>т)т| ^(п — п)(/^7) — ^ | <
< С311 (пт, >т)т 11| | J(п — п)(> — 7) — J(>) > 11 + С311 (пт, >т)т 11| | А711 < < С3||(пт,/>т)т|и(п — п) — J(п)||||>|| + С3||(пт,>т)т|и(п — 7)||||7||+
+ С3|(пт,>т)т||А||7||< 2^3с3|(пт,>т)т|2 + ^3с3|(пт, ¿УИИ^, 7т)т|| +
+ С3 ||АПН(пт,>т)тНП(пт,7т)тн =
= 2-У3с3||(пт,/>т)т||2 + (^3С3 + С3|А|)|(пт,>т)т||(пт, 7
Воспользовавшись неравенством Юнга
1 1 р
ху < -хр + £р-1 ур-1, Ух > 0, Уу > 0, Ур > 1, Уе > 0, (18)
£
получим
дv(п,/) / J(п — п)(/> — 7) — J(п)>
д(пт, >>т)М —А7
<
< 2^3с3|(пт, >>т)т|2 + 1(^3С3 + С3|А|)|(пт, >т)т|2+
£
+ £(У3С3 + С3|А|)|(пт, 7т)т||2 = й(||(п)т,/)т)т||) + 02(|(пт,7т)т|)
Здесь 01 (в) = (2^3сз + -(-У3сз + сз||А||))52, 02(5) = е(^3сз + сз||А||)52, 5 е
£
е [0, £ > 0 - произвольная положительная константа.
Далее при всех (>, )т е К6 и (пт, от)т е К6 справедлива оценка
дУ(п,>) / J(п — п)о + п \ д(>т, + Ь2п У"
^ II /Л ^т,, / J(п — п)г> — Р1—1^1п А ^
< Сз|(п,г/)| ( АО — Р2—^(п> — п^п] < < сз||(пт^т)т|||/(п—п)г>—-Р-^пН + СзНОЖЛНИР — -¿"^й —«Р^Н <
<У3сз|(пт ,>т)т ||||(пт йт)т || + cзНP-1ННQlНН(чт,i>т)т ||||(пт ,от)т ||+ + ^Зсз | -2—1 |||^1ПП(пт У )т ||||(пт йУ || +
+ сз!И||||(пт ,>т)т ||||(пт ,от)т || = (^3сз + cзНPl-1ННQl Н+
+ ^Зсз НР—1ННР1Н + сзНА|)Н(>т ,>т)т НН(пт ^
Воспользовавшись неравенством Юнга (18), получим
дУ(п,>) / J(п — п)о + п \ д(пт, Vх)т ^ + Ь2п )<
< ^(^Зсз + сзНР1—1НН^1Н + -^Зсз^ННАН + сзНАН)Н(пт,>т)т||2+
£
+ £(л/3сз + сз|Р—11 ^л/3сзУРГ1УУР1У + сз||А||)||(пт,от)т||2 =
>т)т ||)+ ф2(|(пт ,От)т ||).
Здесь ф1(в) = -(^Зсз + сзНРТ1НН^1Н + ^3сзНР2—1ННР1Н + сз|А|)в2,&(«) =
£
= £^\/3сз + СзНР—1 ННН + ^/Зсз^ННАН + сз||А||)в2, 5 е [0, + ю), £ > 0 произвольная положительная константа. При с4 > 2л/3сз и
> (2У3сз + сзНР—ЧМ + УЗсз^ЩДН +2сзНАН)
£ >-т=-
С4 — 2у3з
при всех в > 0 справедливо неравенство
а(в) = аз(в) — 01(5) — ф1(в) = С4в2 — (2^3сз + -(^3сз + сз|А|))в2—
£
-(^3сз + сзНР^НШ + ^3сзНР2—'ННРхН + сз|А|)52 > 0,
£
то есть функция а( ) является функцией класса Кж. Следовательно, при с4 > 2л/3сз согласно теореме 1 система, составленная из уравнений системы (1) и уравнений наблюдателя (2) с управлением т = &(?), 7), асимптотически устойчива в целом в точке п = 0, V = 0, П = 0, 7 = 0.
Рис. 1. Переходные процессы системы (сплошная линия) и наблюдателя (пунктир)
Результаты численного моделирования системы (1) и наблюдателя (2) при управлении т = &(?), 7) представлены на рис. 1 при следующих значениях параметров и начальных данных рассматриваемой системы и наблюдателя: шц = 1.1274, Ш22 = 1.8902, т2з = тз2 = -0.0744, тзз = 0.1278, = 0.0358, ¿22 = 0.1183, ¿2з = -0.0124, ¿з2 = -0.0041, ¿зз = 0.0308 (собственные числа матрицы А = -Мимеют вид А1 = -0.2428, Л2 = -0.0627, Аз = = -0.0318), Р1 = ^(5, 5, 5), = ^(5, 5, 5), = ^(1,1,1), К = = ^(0.0225, 0.0225,0.0225), К = diag(0.3,0.3,0.3),
(х(0), у(0), ^(0), и(0), ^(0), г(0))т = (3, 3,3.14,0.05, 0.05,0.01)т, (£(0),у(0),^(0),м(0),£(0),г(0))т = (3.5,3.5,3,0.2,0.2,0.1)т.
Отметим, что рассматриваемые значения параметров т^ и , г = 1,3, ] = 1, 3, соответствуют с учетом масштабирования параметрам математической модели реального судна [5, 6].
5. Заключение
В настоящей работе представлено решение задачи стабилизации заданного положения корабля на плоскости на основе принципа разделения. Установлено, что для рассматриваемой нелинейной модели плоского движения корабля выполнены достаточные условия справедливости принципа разделения. Синтез управления осуществлен при помощи раздельного построения стабилизирующей обратной связи по состоянию и наблюдателя с последующей подстановкой оценки состояния системы наблюдателем в обратную связь. Полученные теоретические результаты подтверждены результатами численного моделирования.
Работа выполнена при поддержке гранта РФФИ №11-01-00733 и Программы Президента РФ по государственной поддержке ведущих научных школ (грант №НШ-4144.2010.1).
Список литературы
1. Уонэм М. Линейные многомерные системы управления: Геометрический подход. М.: Наука, 1980. 376 с.
2. Машиностроение. Энциклопедия / Ред. совет: К. В. Фролов (пред.) и др. Автоматическое управление. Теория. Т. I-4 / Е. А. Федосов, А. А. Красовский, Е. П. Попов и др. Под общ. ред. Е. А. Федосова. М.: Машиностроение, 2000. 688 с.
3. Freeman R. Global internal stabilizability does not imply global external stabilizability for small sensor disturbances // IEEE Trans. on Autom. Control. 1995. V. 40, № 12. P. 2119-2122.
4. Голубев А.Е., Крищенко А.П., Ткачев С.Б. Стабилизация нелинейных динамических систем с использованием оценки состояния системы асимптотическим наблюдателем (обзор) // Автоматика и телемеханика. 2005. № 7. С. 3-42.
5. Fossen T.I., Grovlen A. Nonlinear output feedback control of dynamically positioned ships using vectorial observer backstepping // IEEE Trans. on Contr. Syst. Technol. 1998. V. 6, № 1. P. 121 -128.
6. Fossen T.I. Guidance and Control of Ocean Vehicles. New York: Wiley, 1994. 300 p.
7. Robertsson A., Johansson R. Comments on "Nonlinear output feedback control of dynamically positioned ships using vectorial observer backstepping" // IEEE Trans. on Contr. Syst. Technol. 1998. V. 6, № 3. P. 439-441.
8. Крищенко А.П. Стабилизация программных движений нелинейных систем // Известия АН СССР. Техническая кибернетика. 1985. № 6. С. 103 -112.
9. Красовский Н.Н. Некоторые задачи теории устойчивости движения. М.: Физматгиз, 1959. 212 с.
electronic scientific and t echnical periodical
SCIENCE and EDUCATION
El № FS77 - 30569. №0421100025. ISSN 1994-0408
Output Feedback Control of Ships Planar Motion 77-30569/233615
# 10, October 2011 A.E. Golubev
Bauman Moscow State Technical University
This note deals with a ship positioning control problem. The dynamics considered are nonlinear ship dynamics on the plane. Only the center of mass coordinates and yaw angle are supposed to be physically measured.
The main result states that separation principle holds for the nonlinear model in question. Specifically, the stabilizing control law is synthesized in two steps. First, the state feedback control law is found using the feedback linearization technique. Then, the global asymptotical observer is designed and the observer's state is put into the state feedback control instead of the system's state.
The resulting system in closed-loop with the observer state feedback is shown to be globally asymptotically stable. The proved theoretical results are also illustrated through simulation.
References
1. Uonjem M. Linejnye mnogomernye sistemy upravlenija: Geometricheskij podhod. M.: Nauka, 1980. 376 s.
2. Mashinostroenie. Jenciklopedija / Red. sovet: K. V. Frolov (pred.) i dr. Avtomaticheskoe upravlenie. Teorija. T. I-4 /E. A. Fedosov, A. A. Krasovskij, E. P. Popov i dr. Pod obw. red. E. A. Fedosova. M.: Mashinostroenie, 2000. 688 s.
3. Freeman R. Global internal stabilizability does not imply global external stabilizability for small sensor disturbances // IEEE Trans. on Autom. Control. 1995. V. 40, № 12. P. 2119-2122.
4. Golubev A.E., Krishchenko A.P., Tkachev S.B. Stabilization of nonlinear dynamic systems using the system state estimates made by the asymptotic observer // Automation and Remote Control. 2005. V. 66, № 7. P. 1021-1058. Translated from Avtomatika i Telemekhanika. 2005. № 7. P. 3-42.
5. Fossen T.I., Grovlen A. Nonlinear output feedback control of dynamically positioned ships using vectorial observer backstepping // IEEE Trans. on Contr. Syst. Technol. 1998. V. 6, № 1. P. 121 -128.
6. Fossen T.I. Guidance and Control of Ocean Vehicles. New York: Wiley, 1994. 300 p.
7. Robertsson A., Johansson R. Comments on "Nonlinear output feedback control of dynamically positioned ships using vectorial observer backstepping" // IEEE Trans. on Contr. Syst. Technol. 1998. V. 6, № 3. P. 439-441.
8. Krishchenko A.P. Stabilizacija programmnyh dvizhenij nelinejnyh sistem // Izvestija AN SSSR. Tehnicheskaja kibernetika. 1985. № 6. S. 103-112.
9. Krasovskij N.N. Nekotorye zadachi teorii ustojchivosti dvizhenija. M.: Fizmatgiz, 1959.212 s.
