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

СЛЕДЯЩАЯ СИСТЕМА КОМБИНИРОВАННОГО УПРАВЛЕНИЯ

В.О. Никифоров, Г.В. Лукьянова

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

УДК 681.5.01
СЛЕДЯЩАЯ СИСТЕМА КОМБИНИРОВАННОГО УПРАВЛЕНИЯ
В.О. Никифоров, Г.В. Лукьянова

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

Введение

В подавляющем большинстве практических случаев построение высокоточных систем автоматического управления не может быть осуществлено без учета влияния внешних воздействий – сигналов задания и (или) внешних возмущений. Один из подходов к синтезу систем управления, работающих при внешних воздействиях, состоит в использовании принципа внутренней модели [1–6]. В соответствии с данным принципом, внешнее воздействие рассматривается в качестве выхода автономной динамической системы (командного генератора или генератора возмущения), а для асимптотического воспроизведения (или полной компенсации) внешнего воздействия модель генератора должна быть соответствующим образом воспроизведена в структуре замкнутой системы управления.
Принцип внутренней модели является хорошо разработанным для задач асимптотического слежения и предполагает построение систем комбинированного управления, содержащих в своей структуре обратные связи по ошибке слежения и прямые связи по сигналу задания [4–8]. Однако большинство представленных в литературе решений ориентировано на использование методов пространства состояний, что, в частности, требует измерения координат векторов состояния объекта управления и командного генератора.
В работе предложен метод синтеза, позволяющий избежать необходимости измерения координат векторов состояния и предполагающий использование специальных наблюдателей старших производных сигнала задания.

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

Рассматривается линейный объект управления

y  k ( p) u , ( p)

(1)

где y  y(t) – регулируемая переменная; u  u(t) – сигнал управления; p  d / dt – оператор дифферен-

цирования; ( p)  pn  n1pn1  n2 pn2  0 и ( p)  pm  m1 pm1  m2 pm2  0 – нормированные полиномы с известными постоянными коэффициентами i и  j ; k – «высокочастот-

ный» коэффициент усиления. Дополнительно полагается, что полином ( p) является гурвицевым и

nm.

Будем считать, что сигнал задания g  g(t) представим в виде выхода линейной модели (команд-

ного генератора) конечной размерности:

x  x ,

(2)

g  cT x .

(3)

В модели (2)–(3) x – q -мерный вектор состояния командного генератора с начальным значением

x(0) ;  – известная q  q матрица постоянных коэффициентов; c – вектор постоянных коэффициентов

соответствующей размерности. Без потери общности будем считать, что пара (cT , ) является полно-

стью наблюдаемой. Будем также полагать, что модель (2)–(3) является виртуальной, вектор состояния командного генератора x прямым измерениям недоступен, а измеряется зашумленный сигнал задания вида

  g ,

(4)

где   (t) – шум измерений (заранее не известная функция времени).

Рассматриваемая задача состоит в обеспечении условия limt (g(t)  y(t))  0 при (t)  0 ,

(5)

Научно-технический вестник Санкт-Петербургского государственного университета информационных технологий, механики и оптики, 2011, № 6 (76)

39

СЛЕДЯЩАЯ СИСТЕМА КОМБИНИРОВАННОГО УПРАВЛЕНИЯ

или в обеспечении малой ошибки слежения e  g  y в случае малой нерегулярной составляющей  .
Замечание 1. Здесь малость сигнала (т.е. малость функции времени) понимается в смысле какойлибо нормы амплитуды данного сигнала.

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

Для синтеза наблюдателя производных сигнала задания используем следующий результат. Лемма 1 [9, 10]. Сигнал задания g может быть представлен в виде выхода модели
  G  l g ,

(6)

g  T  ,

(7)

где G – произвольная q  q гурвицева матрица постоянных коэффициентов, образующая с q 1 векто-

ром постоянных коэффициентов l полностью управляемую пару; T  cT P1 – 1 q вектор постоян-

ных коэффициентов; вектор состояния  связан с вектором состояния x командного генератора (3)–

(4) соотношением подобия

  Px ,

(8)

а q  q матрица P является решением матричного уравнения

P  GP  lcT .

(9)

Тогда легко показать, что i -ая производная сигнала g может быть представлена в виде [10]

g (i)  T (G  lT )i  ,

i  1, 2, 3,

Замечание 2. В силу гурвицевости матрицы G начальные условия модели (6) могут быть выбраны

произвольными. Однако если x(0) и (0) не согласованы (например, если (0)  0 ), то выражения (7) и

(8) будут содержать дополнительное экспоненциально затухающее слагаемое, которым в дальнейшем рассмотрении будем пренебрегать.
При построении замкнутой системы автоматического управления будем использовать следующий результат.
Лемма 2 [11–13]. Пусть заданы произвольные нормированные полиномы

A( p)  pM  aM 1pM 1  aM 2 pM 2  a0 ,

B( p)  pN  bN 1pN 1  bN 2 pN 2 b0 , где M  N . Тогда существуют единственные полиномы

R( p)  rN 1pN 1  rN 2 pN 2  r0 ,

D( p)  pM N  dM N 1pM N 2  d0 , удовлетворяющие уравнению

A( p)  B( p)D( p)  R( p) .

Основной результат

Сформируем управление в виде

u



1 k

d( p) D( p)(

p)



g



R( p) d( p)

y 

,

(10)

где d( p) – произвольный нормированный гурвицев полином степени n 1 , нормированный полином

D( p) степени n  m 1 и полином R( p) степени n 1 являются решениями уравнения

d( p)( p)  ( p)D( p)  R( p)

(11)

с произвольным гурвицевым полиномом ( p)  pnm  nm1pnm1  nm2 pnm2  0 . При этом сигнал g  g(t) формируется по правилу

g  gnm  nm1gnm1  nm2gnm2  0g ,

(12)

где для формирования оценок gi i -ых производных сигнала задания g используется специальный наблюдатель, состоящий из вспомогательного фильтра

  G   l 

(13)

и блока расчета производных

gi  T (G  lT )i  , i  1, 2, n  m .

(14)

40 Научно-технический вестник Санкт-Петербургского государственного университета
информационных технологий, механики и оптики, 2011, № 6 (76)

В.О. Никифоров, Г.В. Лукьянова

Размерность, свойства и правила выбора (расчета) матрицы G и векторов l и  определены в лемме 1.
Утверждение. Управление, формируемое в соответствии с выражениями (10), (12)–(14), обеспечивает для объекта управления (1) выполнение целевого условия (5), а также малую ошибку слежения e  g  y в случае малости нерегулярной составляющей  .
Доказательство. С учетом подстановок A( p)  d( p)( p) , B( p)  ( p) , N  n и M  2n 1 из
леммы 2 следует существование полиномов D( p) и R( p) , удовлетворяющих уравнению (11). Подстав-
ляя далее управление (10) в (1), после элементарных преобразований (при этом происходит сокращение полиномов ( p) и d( p) ) получаем уравнение замкнутой системы:

y 1 g. ( p)

(15)

При (t)  0 из леммы 1 следует, что уравнения (13) и (15) определяют старшие производные сиг-

нала задания и, таким образом (см. (12)), g  ( p)g . Подставляя последнее равенство в (14), получаем

y  g , т.е. выполнение целевого условия (5). На самом деле имеет место не абсолютное равенство, а

асимптотическое стремление y к g , так как мы пренебрегли экспоненциально затухающим слагаемым,

вызванным несогласованностью начальных условий командного генератора и наблюдателя производных (см. замечание (2)). Так как для формирования оценок производных сигнала задания используется асимптотически устойчивый динамический фильтр (13) и операции умножения на постоянные коэффициенты (14), то можно показать, что присутствие в сигнале измерений  малой (в смысле амплитудной

нормы) шумовой составляющей (t) приведет к малым ошибкам слежения (в отличие от попыток ис-

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

управления происходит сокращение двух полиномов – ( p) и d( p) . По этой причине они должны быть

гурвицевыми. Полином d( p) назначается разработчиком и может быть интерпретирован как характери-

стический полином наблюдателя вектора состояния объекта управления. Полином ( p) является числи-

телем передаточной функции объекта управления, и его гурвицевость означает требование минимальной

фазовости объекта управления. Как видно из выражения (15), полином ( p) определяет динамику за-

мкнутой системы и поэтому может быть интерпретирован как желаемый характеристический полином замкнутой системы. Известно, что минимальный динамический порядок замкнутой системы может быть обеспечен равным относительной степени исходного объекта управления, т.е. равным n  m .
Предложенный алгоритм управления (10) является комбинированным, так как содержит обратную связь по регулируемой переменной y , а также прямую связь g по сигналу задания, для формирования

которой используется наблюдатель производных сигнала задания (13), (14).

Пример

Рассмотрим неустойчивый объект управления

y



3

p3



p2 2p2 

p

u 1

.

Пусть цель управления состоит в асимптотическом слежении за сигналом задания

g(t)  5sin2t  0,41, который можно представить в качестве выхода командного генератора вида

 x1 x2

  



0  4

1 x1 

0

x2

 

,

g  x1 , x1(0)  2 , x2 (0)  9,17 .

Выберем ( p)  p2  5 p  6 , d( p)  p2  3p  2 . Тогда, решая уравнение (11), находим

D( p)  p  6 и R( p)  12 p2  33p  6 . Пусть

G



0  6

1  5 ,

l



0 6

.

Решая уравнение (9), находим

P



0,115 1,153

 0,288

0,115

 

и T  0,334 0,834 .

Тогда наблюдатель первой и второй производной сигнала задания будет описываться выражениями

Научно-технический вестник Санкт-Петербургского государственного университета информационных технологий, механики и оптики, 2011, № 6 (76)

41

СЛЕДЯЩАЯ СИСТЕМА КОМБИНИРОВАННОГО УПРАВЛЕНИЯ

1  2 ,

 2  61  52  6 , g1  3,331  0,332 ,

g2  1,331  3,332 ,

где   g   является зашумленным сигналом задания, а  – белый шум определенной интенсивности.

Выберем 1(0)  0 , 2(0)  0 . Окончательно искомое управление принимает вид

u



1 3

p2 p2

 3p  2  8 p 12

 

g

 12 p2  33 p  6 p2 3p  2

y

 

,

где сигнал g формируется по правилу

g  g2  5g1  6 . 

Рис. 1. Зашумленный сигнал задания, доступный измерениям

Рис. 2. Оценка первой производной ( g – первая производная сигнала задания, g1 – оценка первой производной сигнала задания g )
Рис. 3. Оценка второй производной ( g – вторая производная сигнала задания, g2 – оценка второй производной сигнала задания g )
Результаты моделирования переходных процессов в замкнутой системе приведены на рис. 1–4. На рис. 1 представлен зашумленный сигнал задания   g   . На рис. 2, 3 представлены первая и вторая производные сигнала задания и их оценки. Как видно из рис. 2, 3, несмотря на существенную зашумлен-
42 Научно-технический вестник Санкт-Петербургского государственного университета
информационных технологий, механики и оптики, 2011, № 6 (76)

В.О. Никифоров, Г.В. Лукьянова
ность сигнала задания, представленный наблюдатель обеспечивает точное восстановление его производных. На рис. 4 приведена ошибка слежения, демонстрирующая выполнение цели управления, т.е. высокую точность слежения за командным сигналом g(t) .

Рис. 4. Ошибка слежения
Заключение
Представленная процедура синтеза позволяет построить для объекта «вход–выход» комбинированную систему слежения за командным генератором сигнала задания. Для построения управления не требуется измерения координат вектора состояния объекта или командного генератора. Предложенное решение основано на использовании специального наблюдателя производных сигнала задания. Как показано в работе и проиллюстрировано результатами моделирования, даже существенная зашумленность сигнала задания не приводит к потере работоспособности предложенной системы.

Литература

1. Johnson C.D. Accommodation of external disturbances in linear regulator and servomechanism problems // IEEE Trans. on Automatic Control. – 1971. – V. 16. – № 6. – Р. 635–644.
2. Francis B.A. and W.M. Wonham The internal model principle for linear multivariable regulators // Appl. Mathematics and Optimization. – 1975. – № 2. – Р. 170–194.
3. Davison E.J. The robust control of a servomechanism problem for linear time-invariant multivariable systems // IEEE Trans. on Automatic Control. – 1976. – V. 21. – № 1. – Р. 25–34.
4. Уонем М. Линейные многомерные системы управления: Геометрический подход. – М.: Наука, 1980. – 376 с.
5. Синтез дискретных регуляторов при помощи ЭВМ / В.В. Григорьев, В.Н. Дроздов, В.В. Лаврентьев, А.В. Ушаков. – Л.: Машиностроение, 1983. – 245 с.
6. Дроздов В.Н., Мирошник И.В., Скорубский И.В. Системы автоматического управления с микроЭВМ. – Л.: Машиностроение, 1989. – 284 с.
7. Акунов Т.А., Алишеров С., Оморов Р.О., Ушаков А.В. Матричные уравнения в задачах управления и наблюдения непрерывных объектов. – Бишкек: Илим, 1991. – 61 с.
8. Мирошник И.В., Никифоров В.О. Синтез линейных систем автоматического управления. – СПб: СПбГУ ИТМО, 2000. – 54 с.
9. Nikiforov V.O. Adaptive non-linear tracking with complete compensation of unknown disturbances // European Journal of Control. – 1998. – V. 4. – № 2. – Р. 132–139.
10. Никифоров В.О. Адаптивное и робастное управление с компенсацией возмущений. – СПб: Наука, 2003. – 282 с.
11. Острем К., Виттенмарк Б. Системы управления с ЭВМ. – М.: Мир, 1987. – 480 с. 12. Лукьянова Г.В., Никифоров В.О. Алгоритм компенсации внешних детерминированных возмущений:
операторный метод синтеза // Научно-технический вестник СПбГУ ИТМО. – 2003. – № 10. – С. 5–10. 13. Бобцов А.А., Лукьянова Г.В., Никифоров В.О. Алгоритм компенсации внешнего гармонического воз-
мущения неизвестной частоты для систем активной виброзащиты // Научно-технический вестник СПбГУ ИТМО. – 2007. – № 11. – С. 39–43.

Никифоров Владимир Олегович Лукьянова Галина Владимировна

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

Научно-технический вестник Санкт-Петербургского государственного университета информационных технологий, механики и оптики, 2011, № 6 (76)

43