Например, Бобцов

АДАПТИВНЫЙ РЕГУЛЯТОР СТАБИЛИЗАЦИИ ПРОСТОЙ СТРУКТУРЫ

АДАПТИВНЫЙ РЕГУЛЯТОР СТАБИЛИЗАЦИИ ПРОСТОЙ СТРУКТУРЫ

3 АНАЛИЗ И СИНТЕЗ СЛОЖНЫХ СИСТЕМ

УДК 681.513.6
АДАПТИВНЫЙ РЕГУЛЯТОР СТАБИЛИЗАЦИИ ПРОСТОЙ СТРУКТУРЫ1
В.О. Никифоров, Д.H. Герасимов

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

Введение

В современной теории адаптивного управления предложено решение многих сложных задач: адаптивного управления линейными объектами по выходу с использованием схемы расширенной ошибки [1–3] или алгоритмов адаптации высокого порядка [4–6], адаптивного управления нелинейными объектами с использованием итеративной процедуры синтеза «обратного обхода интеграторов» [3, 7, 8], робастного управления в условиях внешних возмущений [3, 9], адаптивной компенсации заранее неизвестных детерминированных возмущений [6, 10, 11] и т.д.
Одним из основных недостатков большинства предложенных систем адаптивного управления является сложная структура адаптивных регуляторов. Так, динамический порядок адаптивного регулятора
с расширенной ошибкой равен 2n n  m  2 1 (где n – степень знаменателя передаточной функции
объекта управления, а m – степень числителя передаточной функции), а для расчета сигнала управления
требуется выполнение 10n  3 операций умножения и одной операции деления [3]. Очевидно, что реали-
зация на практике регуляторов высокого динамического порядка, содержащих большое число нелинейных элементов (перемножения и деления), весьма затруднительна, несмотря на практические успехи развития микропроцессорных систем управления. В связи с этим одной из актуальных задач современной теории адаптивного управления является синтез достаточно простых и эффективных алгоритмов адаптивного управления.
Для линейного стационарного многомерного объекта с неизвестными параметрами, представленного нижней треугольной канонической формой, предложен регулятор адаптивной стабилизации простой структуры: адаптивный регулятор содержит всего один настраиваемый параметр и один простой алгоритм адаптации интегрального типа. Принцип построения регулятора основан на использовании идеи «большого коэффициента усиления». Представлены также робастные модификации регулятора, обеспечивающие ограниченность траекторий замкнутой системы, действующей в условиях внешних возмущений.

Постановка задачи

Рассматривается линейный объект управления вида x  Ax  enu,

(1)

где x  Rn – вектор состояния, доступный для прямых измерений, ei – n -мерный единичный вектор с единицей на i -м месте, u – сигнал управления, а матрица системы имеет вид

 a11

 

a21

A 

an11

 an1

1
a22 
an 12 an2

0 1
an 13 an3

 0

0

 

 .



1

 

 ann 

(2)

Параметры матрицы ai j считаются неизвестными. Пара A,en  является полностью управляемой,

 а пара e1T , A – полностью наблюдаемой. Требуется синтезировать стабилизирующее управление в фор-

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

limt x(t)  0 .

(3)

1 Работа выполнена при финансовой поддержке ФЦП «Научные и научно-педагогические кадры инновационной России» на 2009–2013 годы (соглашение №14.B37.21.0406)
48 Научно-технический вестник информационных технологий, механики и оптики,
2012, № 5 (81)

В.О. Никифоров, Д.H. Герасимов

Вспомогательный результат

Утверждение 1. Для пары n  n матриц B и G() , где B – произвольная нижняя треугольная

матрица, т.е.

b11

B



b21 

0
b22  

0

0 

  

,

bn1

bn 2



bnn

 

с постоянными коэффициентами bi, j , а G() – сопровождающая матрица характеристического поли-

нома от комплексной переменной p в форме бинома Ньютона

 D( p, )  p   n  pn  Cn1pn1  Cn22 pn2   C1 n1 p  n ,
т.е.

(4)

0



G()



 

0 0

 

n

1 0 0 C1 n1

 0

 

 1

 , 

 Cn1

(5)

существует такое число 0  0 , что для всех   0 матрица S  B  G()

является гурвицевой.

В формулах (4) и (5) Cn j – коэффициент бинома Ньютона степени n , стоящий при j -й степени

переменной  (при этом Cn  1 и C0  1 ). Доказательство. Пусть

R( p)  det  pI  B  pn  rn1 pn1  r1 p  r0

(6)

– характеристический полином матрицы B , а rj – его постоянные коэффициенты. Тогда можно пока-

зать, что характеристический полином N ( p, ) матрицы S можно представить в следующем виде:

N ( p, )  det  pI  S  det  pI  B  G()  pn  Cn1  rn1 pn1 
   Cn22  rn2  n2 (r , 1, C1) pn2  Cn33  rn3  n3 (r , 2 , C2 ) pn3     C1n1  r1  1(r , n2 , Cn2 ) p  n  r0  0 (r , n1, Cn1) ,

где r – множество всех коэффициентов характеристического полинома (6); j – множество всех степе-

ней  j ( j  0, j ); C j – множество всех коэффициентов Cn j ( j  1, n 1); n j (r , j1,C j1) – сумма соот-
ветствующих коэффициентов ri , степеней i и коэффициентов Ci . Так как члены n j содержат настраиваемые коэффициенты  степени j 1 , а стоящие с ними в
скобках при тех же степенях комплексной переменной p коэффициенты гурвицева характеристического

полинома D( p, ) содержат настраиваемые коэффициенты  степени j , то справедливо следующее
равенство: lim N ( p, )  D( p, ) . Последнее означает, что полином N ( p, ) является гурвицевым при бесконечно больших значени-
ях коэффициента  . Однако свойство экспоненциальной устойчивости (задаваемое гурвицевым полиномом D( p, ) ) является грубым, т.е. выполняется при некоторых отклонениях коэффициентов характери-

стического полинома от их номинальных значений. По этой причине полином N ( p, ) будет гурвицевым при конечных значениях достаточно большого коэффициента  . Утверждение доказано.

Структура регулятора

Сформируем адаптивный регулятор в виде

u  n x1  C1 n1x2   Cn22 xn2  Cn1xn ,

(7)

где постоянные коэффициенты C j – коэффициенты характеристического полинома (4) (т.е. коэффициен-

ты бинома Ньютона), а настраиваемый параметр  генерируется алгоритмом адаптации

Научно-технический вестник информационных технологий, механики и оптики, 2012, № 5 (81)

49

АДАПТИВНЫЙ РЕГУЛЯТОР СТАБИЛИЗАЦИИ ПРОСТОЙ СТРУКТУРЫ

  x12 с постоянным положительным коэффициентом усиления  .

(8)

Замечание. Управление (7) может быть записано в «свернутой» форме, минимизирующей при практической реализации число операций взаимного перемножения настраиваемого параметра  :

   u   Cn1xn   Cn2 xn1   Cn3 xn2 x1 .

(9)

Тогда справедливо следующее утверждение. Утверждение 2. Адаптивный регулятор (7), (8) при любых   0 и (0) обеспечивает для объек-

та управления (1) с произвольными начальными условиями x(0) выполнение целевого условия (3).

Доказательство. Нетрудно убедиться, что матрица замкнутой системы, состоящей из объекта (1) и регулятора (7), имеет вид

 a11

 

a21

A



en

 n ,  C1

n1, 

Cn22 ,  Cn1



 



 an 11

1 a22  an 12

0 1
an 13


 

0

0

 

.

1

 

an1  n an2  C1n1 an3  C2n2  ann  Cn1

Иначе говоря, по своей структуре она соответствует матрице S из утверждения 1. В силу полной

 наблюдаемости пары e1T , A измерения переменной x1 , использующегося в алгоритме адаптации, дос-

таточно для установления факта неустойчивости объекта управления. Так как переменная x1 входит в
алгоритм адаптации в четной (второй) степени, то любое ненулевое значение x1 будет приводить к росту
настраиваемого параметра  , который, в соответствии с утверждением 1, будет расти до тех пор, пока характеристический полином замкнутой системы не станет гурвицевым. Объект управления является линейным, поэтому никаких требований на скорость настройки параметров не накладывается, и стабилизация замкнутого объекта управления может быть обеспечена при любых   0 . Утверждение доказано.

Таким образом, предложен достаточно простой регулятор (7) (или (9)) всего с одним настраиваемым параметром  , генерируемым алгоритмом адаптации (8). Другими словами, показано, что для адап-
тивной стабилизации n -мерного динамического объекта, представленного в канонической нижней тре-
угольной форме, можно использовать динамический регулятор первого порядка (в соответствии с порядком алгоритма адаптации (8)), содержащий в своей структуре n операций умножения (для «свернутой»
формы (9)).

Робастная модификация

Рассмотрим объект управления, подверженный влиянию внешнего возмущения:

x  Ax  enu   ,

(10)

где   (t) – заранее неизвестное ограниченное возмущение, а матрица A имеет вид (2). Поставим за-

дачу обеспечения ограниченности всех траекторий объекта (1) и выполнения целевого условия

x(t)   для всех t  T0  0 ,

(11)

где  – произвольное наперед заданное положительное число.

Используя идеи робастного адаптивного управления [3], в случае возмущенного объекта (10) для

достижения целевого условия (11) можно использовать регулятор (7) с модифицированным алгоритмом

адаптации с параметрической обратной связью:

    x12 ,

(12)

где  – произвольное (малое) положительное число (коэффициент параметрической обратной связи).

Можно также показать, что для достижения целевого условия (11) можно использовать статиче-

ский регулятор (7), параметры которого рассчитываются по формуле

  x12 .

Результаты моделирования

Рассмотрим неустойчивый объект управления x1  0,1x1  x2 ,
x2  x1  2x2  x3 ,
x3  0, 2x1  0,1x2  x3  u ,

50 Научно-технический вестник информационных технологий, механики и оптики,
2012, № 5 (81)

В.О. Никифоров, Д.H. Герасимов

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

u  3 x1  32 x2  3x3 ,   x12 .

(13)

Результаты моделирования процессов стабилизации при   1 и x(0)  [1, 0, 0]T приведены на

рис. 1 и демонстрируют выполнение целевого условия (3). Пусть теперь на вход объекта в соответствии с формулой (10) вместе с управлением поступает не-

измеряемое возмущение (t)  3sin 2t . Тогда для выполнения целевого условия (11) предлагается регу-

лятор (13) с робастным алгоритмом адаптации (12). Результаты моделирования переходных процессов в

замкнутой системе, подверженной влиянию внешнего возмущения при   1 ,   0,1 и x(0)  [1, 0, 0]T ,

приведены на рис. 2 и демонстрируют небольшую установившуюся ошибку стабилизации, которая мо-

жет быть уменьшена за счет увеличения коэффициента  .

u

Рис. 1. Процессы в адаптивной системе стабилизации выходной переменной объекта (11) с регулятором
(12) и алгоритмом адаптации (13) при   1 и x(0)  [1, 0, 0]T
u

Рис. 2. Процессы в адаптивной системе стабилизации выходной переменной объекта (11) с регулятором
(12) и алгоритмом адаптации (13) при действии возмущения (t)  3sin 2t и   1 и x(0)  [1, 0, 0]T
Заключение
В работе показано, что для адаптивной стабилизации n -мерного динамического объекта, представленного в канонической нижней треугольной форме, можно использовать адаптивный динамический регулятор первого порядка, содержащий в своей структуре n операций умножения. Насколько известно авторам, это наиболее простой (по критерию динамического порядка и числа используемых нелинейных операций) адаптивный регулятор стабилизации. Дальнейшее развитие предложенного подхода должно состоять в его расширении на задачи адаптивного слежения, а также управления динамическими объектами с неизмеряемым вектором состояния.
Литература
1. Narendra K.S., Annaswamy A.M. Stable Adaptive Systems. – Englewood Cliffs. N.J.: Prentice-Hall, 1989. – 495p.
2. Никифоров В.О., Фрадков А.Л. Системы адаптивного управления с расширенной ошибкой // Автоматика и телемеханика. – 1994. – № 9. – С. 3–22.
3. Мирошник И.В., Никифоров В.О., Фрадков А.Л. Нелинейное и адаптивное управление сложными динамическими системами. Серия «Анализ и синтез нелинейных систем» / Под ред. Г.А. Леонова и А.Л. Фрадкова. – СПб: Наука, 2000. – 549 с.
4. Morse A.S., Isidor A., Tarn T.J. (Eds.) Hight-order parameter tuners for adaptive control of nonlinear systems // Systems, Models and Feedback: Theory and Applications. – Birkhauser, 1992. – P. 339–264.

Научно-технический вестник информационных технологий, механики и оптики, 2012, № 5 (81)

51

СТРУКТУРА СОБСТВЕННЫХ ВЕКТОРОВ МАТРИЦ СОСТОЯНИЯ ...

5. Nikiforov V.O. Robust high-order tuner of simplified structure // Automatica. – 1999. – V. 35. – № 8. – P. 1409–1415.
6. Никифоров В.О. Адаптивное и робастное управление с компенсацией возмущений. – СПб: Наука, 2003. – 282 с.
7. Kristic M., Kanellakopoulos I., Kokotovic P.V. Nonlinear and Adaptive control design. – N.Y.: John Wiley & Sons, 1995. – 563 p.
8. Никифоров В.О. Робастное управление линейным объектом по выходу // Автоматика и телемеханика. – 1998. – № 9. – С. 87–99.
9. Ioannou P.A., Kokotovic P.V. Instability analysis and improvement of robustness of adaptive control // Automatica. – 1984. – № 5. – P. 583–594.
10. Никифоров В.О. Адаптивная стабилизация линейного объекта, подверженного внешним детерминированным возмущениям // Известия РАН. Теория и системы управления. – 1997. – № 2. – С.103–106.
11. Никифоров В.О. Адаптивная компенсация внешних детерминированных возмущений // Мехатроника, автоматизация и управление. – 2003. – № 5. – С. 8–12.

Никифоров Владимир Олегович Герасимов Дмитрий Николаевич

– Санкт-Петербургский национальный исследовательский университет информационных технологий, механики и оптики, доктор технических наук, профессор, проректор, nikiforov@mail.ifmo.ru
– Санкт-Петербургский национальный исследовательский университет информационных технологий, механики и оптики, кандидат технических наук, доцент, gerasimovdn@mail.ru

52 Научно-технический вестник информационных технологий, механики и оптики,
2012, № 5 (81)