Научная статья на тему 'Стабилизация однозвенного манипулятора при неполном измерении состояния: обратная связь по угловой координате звена манипулятора'

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

CC BY
376
48
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
УПРАВЛЕНИЕ / НЕЛИНЕЙНЫЕ ДИНАМИЧЕСКИЕ СИСТЕМЫ / ПРИНЦИП РАЗДЕЛЕНИЯ / ОБРАТНАЯ СВЯЗЬ ПО ВЫХОДУ / НАБЛЮДАТЕЛИ / ОБХОД ИНТЕГРАТОРА

Аннотация научной статьи по математике, автор научной работы — Голубев А. Е.

В настоящей работе представлено решение задачи стабилизации заданного углового положения однозвенного манипулятора при неполном измерении вектора состояния. Рассматривался случай, когда измерениям доступна только угловая координата звена манипулятора. Синтез управления осуществлен при помощи раздельного построения стабилизирующей обратной связи по состоянию и наблюдателя с последующей подстановкой оценки состояния системы наблюдателем в обратную связь, а также с использованием метода обхода интегратора в наблюдателе. По результатам численного моделирования можно сделать вывод о приблизительно одинаковом при рассмотренных начальных данных и использовании одного и того же наблюдателя качестве переходных процессов системы с управлением, найденным при помощи метода обхода интегратора в наблюдателе, и управлением, основанном на принципе разделения и методе линеаризации обратной связью по состоянию. Возможной областью применения полученных в работе результатов является решение задач управления техническими системами в условиях неполноты измеряемой информации о состоянии системы.

i Надоели баннеры? Вы всегда можете отключить рекламу.
iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

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

НАУЧНОЕ ИЗДАНИЕ МГТУ ИМ. Н. Э. БАУМАНА

НАУКА и ОБРАЗОВАНИЕ

Эл № ФС77 - 48211. Государственная регистрация №0421200025. ISSN 1994-0408

электронный научно-технический журнал

Стабилизация одиозвениого манипулятора при неполном

измерении состояния: обратная связь по угловой координате

звена манипулятора

# 11, ноябрь 2012

DOI: 10.7463/1112.0500549

Голубев А. Е.

УДК 519.71

Россия, МГТУ им. Н.Э. Баумана mathmod@bmstu.ru

1. Введение и постановка задачи

Многозвенные манипуляторы представляют собой класс технических систем, активно используемых в промышленности. При решении задач управления такими объектами одной из существенных проблем может являться отсутствие полной информации о состоянии системы. Измеряются, как правило, значения только части переменных, описывающих состояние манипулятора. Причины неполного измерения вектора состояния могут быть различные: высокая стоимость установки датчиков, технологические ограничения и т.п.

При синтезе алгоритмов управления многозвенными манипуляторами важную роль играет решение задач управления для отдельных звеньев.

В настоящей работе рассматривается задача стабилизации заданного углового положения однозвенного манипулятора, уравнения движения которого имеют вид

x 1 = x2,

x2 = M1 sin x1 — k1(x1 — x3),

x3 = x4,

x4 = —b1 x4 + k2(x1 — x3) + u/J, где x1, x2 — угловая координата и угловая скорость звена манипулятора соответственно; x3, x4 — угловая координата и угловая скорость вала двигателя;

u — управляющий момент, создаваемый двигателем; Mgl sin x1 — момент силы тяжести, действующий на звено манипулятора. Константы Mi, b1, k1, k2, J положительны, притом M1 = Mgl/I, k1 = k/I, k2 = k/J, b1 = d/J, где I, J — моменты инерции звена манипулятора и ротора двигателя соответственно, k — жесткость передаточного механизма, d — коэффициент демпфирования, M — масса звена манипулятора.

В качестве стабилизируемого углового положения манипулятора без ограничения общности рассмотрим положение, в котором x1 = 0, x3 = 0. Дляреше-ния задачи управления требуется построить закон управления в виде обратной связи, использующей значения только измеряемого выхода системы, глобально стабилизирующий положение равновесия x = 0, u = 0 системы (1). Здесь x = (x1,x2,x3, Х4)т Е R4 — вектор состояния системы (1), u Е R — управление.

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

В настоящей работе рассматривается случай, когда измерениям доступна только угловая координата x1 звена манипулятора, т.е. измеряемый выход системы (1) имеет вид y = x1. Показано, что задача стабилизации заданного углового положения однозвенного манипулятора может быть решена с помощью использования нелинейного принципа разделения [1, 2] и метода обхода интегратора в наблюдателе [3,2].

2. Синтез наблюдателя и обратной связи по состоянию

При построении асимптотического наблюдателя для системы (1) с рассматриваемым выходом y = x1 воспользуемся, например, идеями геометрического метода, изложенного в работах [4, 5].

В переменных х = (х1, Х2, Х3, Х4)т, заданных соотношениями

x1 = Xl,

x2 =Х2 — b1X1, 2 (2) x3 = Хз/к1 — b1X2/k1 + (b1 — k2)x1/k1,

x4 = X4/k1 — b1X3/k1 + (b1 — k2)x2/k1 + (2b1k2 — bf)x1/k1, система (1) с рассматриваемым выходом имеет вид

х1 = Х2 — blXl,

х2 = Хз + M1 sin Х1 — X1(k1 + k2) , Хз = Х4 + b1M1 sin Х1 — b1 k1X1, (3)

Х4 = k2M1 sin Х1 + k1u/J, У = Х1.

Заметим, что отображение x = Ф(х), Ф(0) = 0, определяемое равенствами (2), линейно и является диффеоморфизмом пространств R4 = {х} и R4 = {x}. Асимптотический наблюдатель для системы (3) будем искать в виде

х = АХ + LC (Х — х) + ^(Х1) + Bu, (4)

где х = (хъ х2, хз, х4) , А = (a¿¿), i = 1,4, j = 1,4, — квадратная матрица порядка 4 с элементами aj = 1, если j — i = 1, и aj = 0, если j — i = 1, C = (1,0,0,0), B = (0,0,0,1/J)т, вектор L = (I1,12,1з, ¿4)т коэффициентов усиления задает динамику ошибки e = х — х оценки состояния системы (3),

т

= (—blХl,Ml sinх1 — х1 (k1 + k2),b1M1 sinх1 — blklХl,k2Ml sinх^ .

Для дальнейших построений можно использовать линейную технику. Система, описывающая динамику ошибки e = х — х оценки наблюдателем (4) состояния системы (3) при одинаковом управлении в системах (3) и (4), имеет вид

e = (А + LC )e, (5)

где вектор L коэффициентов усиления наблюдателя выбирается таким образом, что матрица А + LC имеет собственные числа только с отрицательными действительными частями. Указанный выбор матрицы L возможен, так как пара

(А, С) наблюдаема [6]. Следовательно, положение равновесия е = 0 системы (5) экспоненциально устойчиво в целом и ошибка оценки состояния не зависит от управления и экспоненциально стремится к нулю.

Отметим, что функцией Ляпунова для системы (5) является положительно

решение которого существует и единственно в силу указанного выше выбора спектра матрицы А + ЬС. Здесь Ц = Цт > 0 — произвольная симметрическая положительно определенная матрица. Для производной по времени функции W(е) в силу системы (5) при всех е Е справедлива следующая оценка:

где Ат1п(Ц) — наименьшее собственное значение матрицы Ц, || • || — евклидова норма в К4.

Далее найдем закон управления в виде непрерывно дифференцируемой обратной связи по состоянию, глобально экспоненциально стабилизирующий положение равновесия х = 0, и = 0 системы (3).

Положение равновесия х = 0, и = 0 аффинной системы (3) (без выхода) во всем пространстве состояний можно стабилизировать методом нелинейной стабилизации, предложенным в работе [7], поскольку эта система во всем пространстве состояний эквивалентна регулярной системе канонического вида, тоже определенной во всем ее пространстве состояний. Преобразование аффинной системы (3) к каноническому виду определяется функцией ф(х) = Хь которую можно найти, следуя работе [8]. Дифференцируя эту функцию в силу аффинной системы (3), находим новые переменные для записи системы канонического вида

определенная квадратичная форма W(е) = етРе, где матрица Р = Рт > 0 удовлетворяет уравнению Ляпунова

(А + ЬС)тР + Р(А + ЬС) = -Ц, Ц = Цт > 0,

(6)

Ж(е) = -етЦе < -Ат1П(Ц)|е

2 ^

£1 = Хъ 6 = £1 = Х2 - ЬхХъ

£з = £2 = Х3 + М1 втХ1 - Х1(к1 + к2) - Ь1Х2 + Ь2lХl,

£4 = £3 = Х4 + М1Х2 СОв Х1 - Ь1М1Х1 сое Х1 - (к1 + к2)Х2 +

(7)

+Ь1к2Х1 + Ь1(к1 + к2)Х1 - Ь1Хз + Ь1Х2 - ь3Х1.

В переменных ¿ система (3) без выхода имеет канонический вид

á = б,

¿2 = 6, (8) 6 = ¿4,

¿4 = f(¿) + kiu/J,

где ¿ = (¿1, ¿2, ¿з, ¿4)т — вектор состояния системы, f(¿) = k2M1 sin¿1 + +Mi cos ¿1 (¿3 + bi¿2) - Mi¿22 sin¿1 - (ki + k2)¿3 - biki¿2 - bi¿4.

Соотношение ¿ = Ф-1(х), Ф-1 (0) = 0, имеющее вид (7), действительно представляет собой замену переменных, так как оно разрешимо относительно X, х = Ф(£). Отметим, что отображение Ф является диффеоморфизмом пространств R4 = {¿} и R4 = {х}, притом функции ¿ = Ф-1(х) и х = Ф(£) таковы, что при всех ¿ G Rn и х G Rn справедливы неравенства

ll¿II = 1|Ф-1(х)11 < 7Ф-1 Ы, 11х11 = ||Ф(е)ll < 7*ll¿ll, (9)

где 7ф, 7^-1 — некоторые положительные константы. Поэтому задача глобальной экспоненциальной стабилизации положения равновесия х = 0, и = 0 системы (3) без выхода эквивалентна аналогичной задаче для положения равновесия ¿ = 0, и = 0 системы (8).

Непрерывно дифференцируемая обратная связь по состоянию, глобально экспоненциально стабилизирующая положение равновесия ¿ = 0, и = 0 системы (8), имеет вид

з

U = k(¿) = J (-/(¿) - ¿ Ki¿i+i) , k(0) = 0. (10)

Замкнутая этим управлением система (8) запишется следующим образом:

¿1 = ¿2,

¿2 = ¿3, (11) ¿3 = ¿4,

¿4 = -K0¿1 - Ki¿2 - K2¿3 - K3¿4.

Здесь постоянные к > 0, г = 0,3, выбираются таким образом, что матрица

А4 = (а4-), г = 1,4, ] = 1,4, с элементами а4- = 1, если ] - г = 1, а4- = -если г = 4, и а4 = 0, если ] - г = 1 и г = 4, имеет собственные числа только с отрицательными действительными частями. Тогда из экспоненциаль-

ной устойчивости в целом положения равновесия £ = 0 системы (8), замкнутой управлением и = ), имеющим вид (10), в силу неравенств (9) следует экспоненциальная устойчивость в целом положения равновесия х = 0 системы (3), замкнутой управлением и = к(Ф-1(х)).

3. Применение нелинейного принципа разделения

Рассмотрим динамическую систему

х = Ах + ?Кх1) + ви,

п ( )

У = <^хъ

где х = (х1,..., Хп)т е А е Кпхп, С е К1хп — постоянные матрицы, пара (А, С) наблюдаема; отображение ^ : К ^ непрерывно-дифференцируемо,

^7(0) = о.

Асимптотическим наблюдателем для системы (12) является система

х = Ах + Ь с(х - х) + ^(х1) + в и, (13)

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

Теорема 1 ([1]). Пусть: 1) вектор-функция ^(хО глобально липшицева; 2) существует непрерывно дифференцируемая обратная связь и = к1(х), к1(0) = 0, по состоянию, глобально экспоненциально стабилизирующая положение равновесия х = 0, и = 0 системы (12). Тогда при управлении и = к1 (х) система

х = Ах + ^(х1) + в u,

(14)

х = ах + Ь с(х - х) + ^7(х1) + в и,

составленная из уравнений системы (12) и уравнений наблюдателя (13), экспоненциально устойчива в целом в точке х = 0, х = 0.

Заметим, что для системы (3) выполнены условия теоремы 1. Следовательно, благодаря линейности замены переменных х = Ф(х), заданной соотношениями (2), система, составленная из уравнений исходной системы (1) и уравнений наблюдателя (4), записанного в переменных X = Ф(х), при упра-

влении и = к(Ф 1 (Ф 1 (X)) экспоненциально устойчива в целом в точке х = 0, X = 0.

Результаты численного моделирования системы (1) и наблюдателя (4), записанного в переменных X = Ф(Х), при управлении и = к(Ф-1(Ф-1(х)) представлены на рис. 1 при следующих значениях параметров и начальных данных рассматриваемой системы и наблюдателя: М = 0,21 кг, I = 0,0093 кг • м2, J = 0,0037 кг • м2, к = 0,18 —, ё = 0,046 ^^, / = 0,15 м, д = 10 м,

рад с с2

ко = 561,1, К1 = 461,9, К2 = 142,4, кз = 19,5, /1 = -17,5, /2 = -114,7, /з = -333,6, /4 = -363,4, (х1(0), Х2(0), хз(0), х4(0)) = (3,14, 0, 3,14, 0), (Х1(0), Х2(0), хз(0), Х4(0)) = (2, 1,2, 3, 1,4).

012345 015345

с

Рис. 1. Переходные процессы системы (сплошная линия) и наблюдателя (пунктир)

при управлении и = к(Ф-1(Ф-1 (X))

4. Метод обхода интегратора в наблюдателе

Рассмотрим систему

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

(15)

х1 = х2 — e2 + ^Х^

х2 = хз + te + ^2 (Х1),

хз = х + to + ^х^

ki

Х4 = ¿4e1 + ^4(Х1) + JJu,

e = (А + LC )e, У = Хl,

где х = (хъ /^2, х^з, е К4 — вектор состояния наблюдателя (4); х

= (х1, х2, х3, х4)т е К4 — вектор состояния системы (3); е = (е1, е2, е3, е4)т = (х1- Xl, Ь - X2, хз- хз, х4- х4)т;

^(Х1) =

( — Ь1Х1 Х

Ml sin Х1 — Х1 (kl + k2)

V

blMl sin Х1 — blklХl

h¿Ml sin Х1

Заметим, что динамическая система (15) в силу линейного соотношения

/

e = х — х эквивалентна системе

х = АХ + ^(Х1) + Bu, Х = АХ + LC (Х — х) + ^(Х1) + Bu,

(16)

состоящей из уравнений системы (3) и наблюдателя (4).

На основе метода обхода интегратора в наблюдателе [2, 3] можно предложить следующий алгоритм построения управления в виде динамической обратной связи и = 7(х, у) по выходу, глобально асимптотически стабилизируещего положение равновесия х = 0, е = 0, и = 0 системы (15). Отметим, что в силу соотношения х = х — е найденная обратная связь будет также и решением задачи стабилизации положения равновесия х = 0, х = 0, и = 0 системы (16).

Шаг 1. Рассмотрим функцию

^(¿1, е) = 1 ¿2 + ^(е),

где = Х1, ^(е) — функция Ляпунова для системы (5). Для удобства используем далее также обозначение ¿2 = Х2 - а1(Х1), где а1( ) — некоторая гладкая функция. Для производной по времени функции V в силу системы (15) справедлива следующая оценка:

V = ад + Ж(е) < ^1(Х2 + ^Ы - е2) - А||е||2 <

< ;1(;2 + а1(Х1) + ^(хО - е2) - Ае2.

Выбрав функцию

а1(Х1) = -(^1 + ё1);1 - ^(хО = -^¿1 - + 61Х1,

где с1 > 0, ё1 > 0 — произвольные положительные константы, получим

V < ^(-ад + ;2 - ¿1^1 - е2) - Ае2 =

О /О О \ о

= -с^ + ^1^2 - (¿1^1 + ;1е2 + Ае2) = -с^ + - 51,

где

12 1 1 * = + 2ЖеУ + (А - Ш4 - 0 при А>

2д/ё1 / V 4^7 2" 4ёГ

Шаг 2. Рассмотрим функцию

V(¿1, ¿2, е) = ^(¿1, е) + 2¿22 + W(е) > 0.

Для удобства используем далее обозначение ;з = Хз - а2(Хъ Х2, хО, где а2( ) — некоторая гладкая функция своих аргументов. Для производной по времени функции V в силу системы (15) справедлива следующая оценка:

1>2 = V + ¿2 ¿2 + <

< -ад2 + ¿1 ¿2 - + ад + < -ад2 - 51 +

+ ¿1 + & + 12(Х1 - Х1) + ^2(Х1) - дХ1 (Х2 + ^1(Х1) - е2^ - Ае2 =

= -ад2 - + ¿2^¿1 + ¿з + о^ъ^хО + /2(:х1 - Х1) + ^2(Х1) -

да1 (х2 + ^1(Х1) - е2^ - Ае2.

дХ1

a2(:xl,xi) = c2z2- zi - ¿2(xi - xi) - ^2(x0 +

(dai \2

+ ^(» + WxO)- d<ä^) Z2 =

= -C2Z2 - zi - ¿2(xi - Xi) - Mi sin Xi + Xi(ki + k2) +

dai (dai\2

+ (X2 - blXl) - 4d^)Z2-

где c2 > 0, d2 > 0 — произвольные положительные константы, получим

da А 2 2 dai Z2 dXi = -CiZ2 - C2Z| + Z2Z3 - Si - S2,

V < -Ciz2 - C2z| + Z2Z3 - (d^ z| - ^ Z262 + Ae2) - Si =

где

да 1 2 1 1

52 = а§22 — 2^/32е2) + (А — 44)е2 -0 при Л> 4^.

Шаг 3. Рассмотрим функцию

^з(21, 22, 2з, е) = ^2(21, 22, е) + 12з2 + W(е).

Для удобства используем далее обозначение 24 = х4 — аз(хъ х2, хз, хО, где аз( ) — некоторая гладкая функция своих аргументов. Для производной по времени функции Уз в силу системы (15) справедлива следующая оценка:

= У2 + 2з2з + W <

< —С122 — С222 + 222з — 51 — 52 + 2з2з + W < < —С122 — С22| — 51 — 52 + 2з ( 22 + + ¿з^ — х0 + ^з(х0 —

/ - / / ч ч da2 ^ da2 ^ \ > 2

■(Х2 + ^1(Х1) - e2) - dXT-Xi - дхх-^J - Ae2 =

dxi dxi

= -ciz2 - c2z2 - Si - S2 + zm z2 + z4 + о^хъЬ^хО + ¿3(:xi - xi) +

, , , da2, Л , Л Л da2 ^ da2 ^ \ Л 2 + ^3(Xl) - ^+ WxO - e2) - dXX2Xi - ^J - Ae2.

а3(Хъ Х2, ^ Xi) = -C3z3 - z2 - ¿3(X1 - Xi) - ^a(xi) +

доад , , , д«2 х , д«2 х , (д«2\2

+ дх;+ WxO) + д-хi + дх"х- <ЦдХг) Z3 =

= -C3Z3 - Z2 - ^(xi - Xi) - biMi sin Xi + bikiXi +

, даад , , , да2х д«2 х , (д«2\2 + ад(Х2 - biXi) + адXi + адX2 - d3 ад z3,

5xi dxi дХ2 v<9xiу

где c3 > 0, d3 > 0 — произвольные положительные константы, получим

• 2 2 2 ( (да2 \2 2 да2 Л 2\

V3 < -ciZi - C2Z2 - C3Z3 + Z3Z4 - I д^) Z3 - e2 + Ae^ -

- Si - S2 = -CiZ2 - C2z| - C3Z3 + Z3Z4 - Si - S2 - S3,

где

(д 1 \2 / 1 \ 1 ^ дХ2з - еу ПА - ^е2 - 0 при А> ^ •

Шаг 4. Рассмотрим функцию

V:(¿1, ¿2, ¿з, ¿4, е) = Уз(^1, ¿2, ¿з, е) + 2¿42 + W(е) > 0.

Для производной по времени функции V в силу системы (15) справедлива оценка

зз

^ = + ад + < - CiZ? * + ад + ад + <

¿=1 ¿=1

з з / к

< - ^Cгz1 * + ¿з +МхЬ- Х1) + ^4(Х1) + 4и -

¿=1 ¿=1 ^

з о \ 4

(Х2 + ^i(Xi) - e2) - ¿ X i) - A ¿ e2.

Тогда при выборе закона управления

u = а4(х, Xi) = J I C4Z4 - z3 - ¿4(;^i - Xi) - ^4(Xi) +

+ дх;(X2 + 'Ata)) + L дх3Xi- <4241 =

= J ^-C4Z4 - Z3 - ^(Xl - Xi) - k2Mi Sin Xi +

, A 5a3 х да3 х Э«3 х , (да3)2 \

+дх-ь -bi ы+д- * i+дх* 2+дьх 3 - Ч дХ3)24 J, (17)

где c4 > 0, d4 > 0 — произвольные положительные константы, для производной по времени функции V4 в силу замкнутой системы (15) справедлива оценка

V < - E сггг2 - ¿ Si - Ae2 - Ae2 - Ae2

i=1 i=i

CiZ / Si Ae Aeq Ae4

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

= -£^ - + e2)2 - E(^¡-Tz - e2

4 1

i=1

При

- E (a - ¿>2 - Ael - Ae2 - Ae4. (18)

A > maxj-1-, -1-, -1-, -M (19)

14di 4d2 4d3 4d4 J

4 4dr' 4d2' 4d3' 4d4

справедливо неравенство

44

~2 A ^e2

V4 Ciz2 - A ^>2, (20)

i=i i=i

4 ( 1

где A = mm{a, E (a - 4d^)} >0

¿=1

Отметим, что выполнения условия (19) можно добиться, например, зафик-

сировав коэффициенты 3 > 0, % = 1,4, и подобрав матрицу Ц = Ц > 0 в уравнении Ляпунова (6), удовлетворяющую данному условию.

2

Соотношения

¿1 = Уъ

¿2 = Х2 - а1(У1^

= Хз - а2(у1,у2,Хl),

¿4 = Х4 - аз(у1,у2

представляют собой гладкую замену переменных, определенную глобально. В

переменных ¿¿, г = 1,4 и е система (15) без выхода, замкнутая управлением (17), примет вид

¿1 = -ад + ¿2 - ^¿1 - е2,

7 (да1 \2 да1

¿2 = £^¿2 - ¿1 + ¿з - «2 -ТГ" ¿2 + -ТГ" е2,

Ч^/ дУ1

/да2\ 2 д«2

¿3 = -cз¿з - ¿2 + ¿4 - «з ¿3 + ТГ"е2, (21)

Ч^У дХ1 'даз\2 даз

ТТ" ¿4 +

^У^ дХ1

е = (А + ЬС )е,

/даз\2 даз

¿4 = ^¿4 - ¿з - «4 -— ¿4 + -— е2,

Чду^ ду1

где ¿ = (¿1, ¿2, ¿з, ¿4)т е Ж4.

Положение равновесия ¿ = 0, е = 0 системы (21) экспоненциально устойчиво в целом. Тогда, так как а^(0) = 0, г = 1,3, положение равновесия у = 0, е = 0 системы (15), замкнутой управлением (17) асимптотически устойчиво в целом. В силу соотношения у = у - е положение равновесия у = 0, у = 0 системы (16) при управлении (17) также асимптотически устойчиво в целом.

Далее, так как отображение Ф, заданное соотношениями (2), является диффеоморфизмом пространств Ж4 = {у} и Ж4 = {у}, из асимптотической устойчивости в целом положения равновесия у = 0, у = 0 системы (16), составленной из уравнений системы (3) и уравнений наблюдателя (4) при управлении и = а4(у, у0, следует асимптотическая устойчивость в целом в точке х = 0, X = 0 системы, состоящей из уравнений исходной системы (1) с рассматриваемым выходом и уравнений наблюдателя (4), записанного в переменных X = Ф(у), при управлении и = а4(Ф-1 (х),х1).

Результаты численного моделирования системы (1) и наблюдателя (4), записанного в переменных £ = Ф(Х), при управлении и = а4(Ф"1(ж),ж1) представлены на рис. 2 при следующих значениях параметров и начальных данных

рассматриваемой системы и наблюдателя: М = 0,21 кг, I = 0,0093 кг • м

^, ё = 0,046 кЛ "

рад с

" = 0,0037 кг • м2, к = 0,18 —, ё = 0,046 / = 0,15 м, д = 10 м, С1 = 2,

С2 = 2, сз = 2, С4 = 2, ¿1 = 0,1, ¿2 = 0,1, ёз = 10~з, ёА = 10~з, /1 = -17,5, /2 = -114,7, /з = -333,6, /4 = -363,4, (Х1(0), Х2(0), Хз(0), Х4(0)) = = (3,14, 0, 3,14, 0), (Х1(0), Х2(0), Хз(0), Х4(0)) = (2, 1,2, 3, 1,4).

2

Рис. 2. Переходные процессы системы (сплошная линия) и наблюдателя (пунктир) при управлении и = а4(Ф-1(£),£1)

5. Заключение

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

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

Работа выполнена при поддержке гранта РФФИ №11-01-00733, №12-07-329 и Программы Президента РФ по государственной поддержке ведущих научных школ (грант №НШ-3659.2012.1).

Список литературы

1. Голубев А.Е., Крищенко А.П., Ткачев С.Б. Принцип разделения для аффинных систем // Дифференциальные уравнения. 2001. T. 37, № 11. C. 1468-1475.

2. Голубев А.Е., Крищенко А.П., Ткачев С.Б. Стабилизация нелинейных динамических систем с использованием оценки состояния системы асимптотическим наблюдателем (обзор) // Автоматика и телемеханика. 2005. № 7. С.3-42.

3. Krstic M., Kanellakopoulos I., Kokotovic P.V. Nonlinear and adaptive control design. New York: John Wiley and Sons, 1995. 563 p.

4. Krener A.J., Respondek W. Nonlinear observers with linearizable error dynamics // SIAM J. Control and Optimization. 1985. V. 23, no 2. P. 197-216.

5. Крищенко А.П., Ткачев С.Б. Нелинейные k(x)-двойственные системы и синтез наблюдателей // Дифференциальные уравнения. 1999. T. 35, № 5. C. 648-663.

6. Уонэм М. Линейные многомерные системы управления: Геометрический подход. М.: Наука, 1980. 376 с.

7. Крищенко А.П. Стабилизация программных движений нелинейных систем // Известия АН СССР. Техническая кибернетика. 1985. № 6. С. 103-112.

8. Жевнин А.А., Крищенко А.П. Управляемость нелинейных систем и синтез алгоритмов управления // Доклады АН СССР. 1981. Т. 258, № 4. С. 805-809.

SCIENTIFIC PERIODICAL OF THE BAUMAN MSTU

SCIENCE and EDUCATION

EL № FS77 - 48211. №0421200025. ISSN 1994-0408

electronic scientific and technical journal

Single-Link Manipulator Output Feedback Control: Manipulator Link Angular Coordinate Feedback # 11, November 2012 DOI: 10.7463/1112.0500549 Golubev A. E.

Russia, Bauman Moscow State Technical University

mathmod@bmstu.ru

In this paper the author considers solving the problem of stabilization of a set angular position of a single-link manipulator when the measurement can only be done by the angular coordinate of the motor shaft. It was sown that synthesis of stabilizing control laws, as well as in the case when measurements can only be done by the angular coordinate of the link manipulator, can be carried out by the principle of separation and bypass integrator in the observer. According to the results of numerical simulation one can draw a conclusion about approximately the same (with the considered initial data and use of the same observer) quality of transient processes of the system with the control found by the bypass integrator in the observer and the control based on the principle of separation and linearization technique using feedback according to the state. The possibility of applying the bypass method to the problem of stabilization can solve this problem also in the case of system disturbances and uncertainties. Possible range of application of the results obtained in the work is solving solve problems of control of technical systems with incomplete information about the state of the measured system.

References

1. Golubev A.E., Krishchenko A.P., Tkachev S.B. Printsip razdeleniia dlia affinnykh sistem [Principle of distribution for affine systems]. Differentsial'nye uravneniia, 2001, vol. 37, no. 11, pp. 1468-1475.

2. Golubev A.E., Krishchenko A.P., Tkachev S.B. Stabilizatsiia nelineinykh di-namicheskikh sistem s ispol'zovaniem otsenki sostoianiia sistemy asimptotich-eskim nabliudatelem (obzor) [Stabilization of nonlinear dynamic systems with the use of the assessment of the status of the system asymptotic observer (review)]. Avtomatika i telemekhanika, 2005, no. 7, pp. 3-42.

3. Krstic M., Kanellakopoulos I., Kokotovic P.V. Nonlinear and adaptive control design. New York, John Wiley and Sons, 1995. 563 p.

4. Krener A.J., Respondek W. Nonlinear observers with linearizable error dynamics. SIAMJ. Control and Optimization, 1985, vol. 23, no 2, pp. 197-216.

5. Krishchenko A.P., Tkachev S.B. Nelineinye k(x)-dvoistvennye sistemy i sintez nabliudatelei [Nonlinear k(x) dual systems and synthesis of observers]. Differ-entsial'nye uravneniia, 1999, vol. 35, no. 5, pp. 648-663.

6. Wonham W.M. Linear Multivariable Control: A Geometric Approach. New York, Springer, 1979. (Russ. ed.: Uonem U.M. Lineinye mnogomernye sistemy up-ravleniia: Geometricheskiipodkhod. Moscow, Nauka, 1980. 376 p.).

7. Krishchenko A.P. Stabilizatsiia programmnykh dvizhenii nelineinykh sistem [Stabilization of programmed motions of non-linear systems]. Izvestiia AN SSSR. Tekhnicheskaia kibernetika [Proceedings of Academy of Sciences of the USSR. Technical Cybernetics], 1985, no. 6, pp. 103-112.

8. Zhevnin A.A., Krishchenko A.P. Upravliaemost' nelineinykh sistem i sintez algo-ritmov upravleniia [Controllability of nonlinear systems and synthesis of control algorithms]. Doklady AN SSSR [Reports of Academy of Sciences of the USSR], 1981, vol. 258, no. 4, pp. 805-809.

i Надоели баннеры? Вы всегда можете отключить рекламу.