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

УДК 681.51

DOI: 10.18698/0236-3933-2017-1-67-74

СИНТЕЗ АСТАТИЧЕСКОГО УПРАВЛЕНИЯ ЛИНЕЙНОЙ СИСТЕМОЙ НА ОСНОВЕ ОБОБЩЕННОЙ ФОРМУЛЫ АККЕРМАННА

Н.Е. Зубов1' 2 Nikolay.Zubov@rsce.ru

nezubov@bmstu.ru

Е.А. Микрин1' 2 В.Н. Рябченко1' 2

1 ПАО «Ракетно-космическая корпорация «Энергия» им. С.П. Королёва», Королёв, Московская обл., Российская Федерация

2 МГТУ им. Н.Э. Баумана, Москва, Российская Федерация

Аннотация

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

Ключевые слова

Линейная М1МО-система, формула Аккерманна, астатическое управление, матричные делители нуля

Поступила в редакцию 14.03.2016 © МГТУ им. Н.Э. Баумана, 2017

Исследование выполнено при поддержке гранта Российского научного фонда (проект №14-11-00046)

Введение. Пусть задана линейная многомерная система со многими входами и многими выходами (Multi Input Multi Output System, MIMO-система) [1, 2]

x(t) = Ax(t) + Bu(t) (1)

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

U(t) = Dv(t), (2)

обеспечивающим управление скоростью перемещения исполнительных органов. Здесь x еМ" — вектор состояния системы; u е Мг — вектор входа системы; v е Мг — вектор внешнего управления; М — множество вещественных чисел.

Предполагается, что система (1) с парой матриц (A, B) является полностью управляемой, а матрица D — обратимой.

Требуется определить управляющую функцию

v(t) = -K x(t), (3)

(4)

такую, что замкнутая MIMO-система

x (t) = Ax(t) + Bu(t); U(t) = DK x(t)

имеет заданный (устойчивый) характеристический полином

p(X) = Xn+r + pn+r_1Xn+r_1 +... + p1X + p0, le C,

(5)

где С — поле комплексных чисел.

Обобщение формулы Аккерманна. Для решения поставленной задачи применим известную формулу Аккерманна [3], обобщая ее на случай систем со многими входами и многими выходами по аналогии, как это сделано в работах [4, 5]. Предварительно перепишем систему (1), (2) в следующей блок-матричной форме:

( x(t) Л ( A B

v u(t) J 0rxn 0rxr

( x(t) л "nxm '

u(t)

D

v(t).

(6)

Нетрудно показать, что при полной управляемости пары матриц (А, В) и обратимости матрицы Б пара матриц из (6)

A B

0rxn 0rxr

(0 л

"nxr

D

(7)

также является управляемой.

Для полной управляемости пары (7) наряду с полной управляемостью пары матриц (А, В) необходимо и достаточно, чтобы квадратная матрица Б е Мгхг порядка г имела полный ранг (была обратимой)

rank D = r,

и выполнялось условие [1]

rank(BD | ABD | ••• | An1BD) = n. Действительно, согласно критерию управляемости Калмана, имеем

(8)

(9)

rank

(0 Л

^nxr

( A

B Л

(0 Л

^nxr

A B

Grxn 0rxr

\n+r-1

(0 Л

vnx r

~D

ft

= rank

Gnxr BD An-1BD

D 0rxr 0rxr

(10)

= rank D + rank (BD | ABD | ••• | An-1BD ) = r + n.

Рассмотрим характеристический полином матрицы

A B

0rxn 0rxr

В силу структуры матрицы (11) он равен полиному

f (X) = Xr det (XIn - A ),

(11)

(12)

где

(XIп - А ) = Кп + /п-{Кп-1 +... + /{К + /о Введем разбиение матрицы Б следующего вида:

б = (б|(1 ) = (( аа - аг-1 )емгхг

Б = (((1 (2 - (г-1), ( = (г,

(13)

(14)

и вычислим делитель нуля максимального ранга Б1 [6]. В таком случае его ранг всегда не превосходит 1, т. е. делитель нуля является вектором-строкой. Тогда делитель нуля максимального ранга матрицы

(0 ^

D

v j

(15)

будет равен

In Önxn

Ö1xn D1

(16)

Действительно, вычисляя произведение матриц

Önxn л ( Önxr л ( Önxr

In Önxn

Ö1xn D1

D 0rxr

у v у v j

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

cIn+r

In Önxn

Ö1xn D1

A в

Örxn Örxr

B Л Л ( I

j j

In Ö Л "nxn (<In - A в Л

Ö1xn D1 J Örxn V ccIr J

( <In - A

Ö1x

в Л

(17)

<DD 1

Следовательно, управляющая функция (3), обеспечивающая замкнутой М1МО-системе характеристический полином (5), будет иметь вид

v(t) = -

(к л

k

v j

x(t),

(18)

где

= (©1 |©2)

( <In - A

0

1xn

к = (K1 |K2 ) =

в Л , (19)

= (©1 (coI„ - A)©1 в + <b©2Dx) ;

<D

kT =(kT | kт) = (0 | ••• | 0 | 1)-^-s;

Qc =

((o л

"nx1

d

v v /

A B ^ (o x

D©1 (n - A) ©1B + «©2D 1 d v /

A B

D©1 («In - A) ©1B+ «©2 D1

^+r-1 ( 0 1 ^ "nx1

d

v /

(20)

(21)

A B

D©1 («In - A) ©1B +«©2d 1

" pn+r-1 (

A B

D©1 («In - A) ©1B + «©2 D1

\n+r-1

(22)

-p1

A B

D©1 («In - A) ©1B + «©2 D1

/

\n+r

-p0h

Здесь матрица полного ранга ® = (®1 | ®2) и скаляр ю являются произвольными величинами, кроме того, величина ю не совпадает с каким-либо собственным числом матрицы

(23)

A B

0rxn 0rxr

Формулы (18)-(22) представляют собой аналитический закон астатического управления — обобщенную формулу Аккерманна для астатического управления. Числовой пример. Пусть в уравнениях (1), (2) имеют место матрицы

Г 0 1 0 ^ Г 0 01

A = 0 0 0 , B = 1 0

0 v 0 0 / 0 v 1

, D = 12,

(24)

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

p(X) = Х5 + 5Х4 + 10Х3 + 10Х2 + 5Х +1. (25)

Это соответствует случаю, когда все собственные значения матрицы

(26)

A B

- DK1 - dkT - DK 2 - dk T

(полюсы замкнутой системы) равны -1.

Отметим, что все собственные значения матрицы (23) (полюсы исходной системы) равны 0:

/ (к) = к5. (27)

При этом их алгебраическая кратность равна 6, а геометрическая — 4, т. е. матрица (11) с матрицей А из (24) содержит четыре жордановых клетки [7] (одну размером 2 х 2 и четыре размером 1 х 1).

Введем следующее разбиение матрицы Б (24):

d = (DI d) =

(1

о л 7

(28)

Тогда делитель нуля Б1 можно записать так

Б1=(0^), (29)

при этом матрица (19) будет равна

К = (К |К ) = (ю©11 |ю©12-©11 |ю©13 | ©12 |ю©14-®13 ) , (30)

где ©ц, ю — произвольные элементы с учетом приведенного выше замечания. Для простоты предположим

ю = 1, © = (1 i 1 i 1 i 1), (31)

тогда матрица (30) будет равна

К = ( К1 |К2 ) = (1|0|1|0| -1), (32)

а вектор (20) —

кт = (кт | к2) = (15|1|1б| -25 | 6). (33)

С помощью непосредственных вычислений можно убедиться, что закон управления (18) с матрицами (32), (33)

v(t) = -

(к л

k

v j

x(t) = -

1 0 1 0 -1

15 1 16 -25 6

x(t)

обеспечивает матрице (26) характеристический полином (25).

Отметим, что, если в соответствии с постановкой задачи желательно обеспечить сохранение исходного характеристического полинома системы (27), то в таком случае это будут обеспечивать матрица (32) и

кт =(к[ | к^) = 0,5(-1 | -1 | -1 | 1 | 0),

т. е.

v(t) = -

(к л

k

v j

x(t) = -

1 0 1 0 -1

-1 -1 -1 1 0

x(t).

Осуществить численный синтез астатического закона обратной связи MIMO-системы можно с помощью элементарной модификации функции acker программного продукта Matlab.

Отметим также, что изменение значений в параметрах (31) будет приводить к формированию других астатических законов, которые также обеспечивают замкнутой системе полином (25).

ЛИТЕРАТУРА

1. Kailath T. Linear systems. New Jersey: Prentice Hall, Englewood Cliffs, 1980. 682 p.

2. Zhou K. Essentials of robust control. New Jersey: Prentice Hall, 1998. 425 p.

3. Дорф Р., Бишоп Р. Современные системы управления / пер. с англ. М.: Лаборатория Базовых Знаний, 2004. 832 с.

4. Мисриханов М.Ш., Рябченко В.Н. Модальный синтез управления энергосистемой на основе обобщенной формулы Аккерманна // Автоматика 2008: Доклады XV международной конференции по автоматическому управлению. Т. 1. Одесса: ОНМА, 2008. С. 362-365.

5. Синтез стабилизирующего управления космическим аппаратом на основе обобщенной формулы Аккерманна / Е.А. Воробьева, Н.Е. Зубов, Е.А. Микрин, М.Ш. Мисриханов, В.Н. Рябченко, С.Н. Тимаков // Известия РАН. Теория и системы управления. 2011. № 1. С. 116-126.

6. Терминальное релейно-импульсное управление линейными стационарными динамическими системами / Н.Е. Зубов, Е.А. Микрин, М.Ш. Мисриханов, А.С. Олейник, В.Н. Рябченко // Известия РАН. Теория и системы управления. 2014. № 3. С. 134-149.

7. Гантмахер Ф.Р. Теория матриц. М.: Наука, 1967. 576 с.

Зубов Николай Евгеньевич — д-р техн. наук, заместитель руководителя по науке научно-технического центра ПАО «РКК «Энергия» им. С.П. Королёва» (Российская Федерация, 141070, Королёв, Московская обл., ул. Ленина, д. 4-а), профессор кафедры «Системы автоматического управления» МГТУ им. Н.Э. Баумана (Российская Федерация, 105005, Москва, 2-я Бауманская ул., д. 5).

Микрин Евгений Анатольевич — д-р техн. наук, академик РАН, генеральный конструктор ПАО «РКК «Энергия» им. С.П. Королёва» (Российская Федерация, 141070, Королёв, Московская обл., ул. Ленина, д. 4-а), заведующий кафедрой «Системы автоматического управления» МГТУ им. Н.Э. Баумана (Российская Федерация, 105005, Москва, 2-я Бауманская ул., д. 5).

Рябченко Владимир Николаевич — д-р техн. наук, ведущий научный сотрудник научно-технического центра ПАО «РКК «Энергия» им. С.П. Королёва» (Российская Федерация, 141070, Королёв, Московская обл., ул. Ленина, д. 4-а), профессор кафедры «Системы автоматического управления» МГТУ им. Н.Э. Баумана (Российская Федерация, 105005, Москва, 2-я Бауманская ул., д. 5).

Просьба ссылаться на эту статью следующим образом:

Зубов Н.Е., Микрин Е.А., Рябченко В.Н. Синтез астатического управления линейной системой на основе обобщенной формулы Аккерманна // Вестник МГТУ им. Н.Э. Баумана. Сер. Приборостроение. 2017. № 1. C. 67-74. DOI: 10.18698/0236-3933-2017-1-67-74

SYNTHESIS OF ASTATIC LINEAR SYSTEMS CONTROL BASED ON GENERALIZED ACKERMANN'S FORMULA

N.E. Zubov1, 2 Nikolay.Zubov@rsce.ru

nezubov@bmstu.ru

E.A. Mikrin1, 2 V.N. Ryabchenko1, 2

1 S.P. Korolev Rocket and Space Corporation Energia, Korolev, Moscow Region, Russian Federation

2 Bauman Moscow State Technical University, Moscow, Russian Federation Abstract Keywords

The main purpose of the article was to examine the syn- Linear MIMO-system, Ackermann's thesis of astatic law of linear multiple-input and output formula, astatic control, matrix zero system control. The synthesis algorithm hinges upon the divisors generalized formula of Ackermann's controller for astatic control. Generalization is achieved by using an input matrix special partition and the use of matrix zero divisors technique

REFERENCES

[1] Kailath T. Linear systems. New Jersey, Prentice Hall, Englewood Cliffs, 1980. 682 p.

[2] Zhou K. Essentials of robust control. New Jersey, Prentice Hall, 1998. 425 p.

[3] Dorf R.C., Bishop R.H. Modern control systems. Addison-Wesley (Russ. ed.: Sovremen-nye sistemy upravleniya. Moscow, Laboratoriya Bazovykh Znaniy Publ., 2004. 832 p.).

[4] Misrikhanov M.Sh., Ryabchenko V.N. Modular synthesis of power system control using generalized Ackermann's formula. Avtomatika 2008: doklady XV mezhdunarodnoy konfer-entsii po avtomaticheskomu upravleniyu. T. 1 [Automatics 2008: Proc. XV int. conf. on automatic control. Vol. 1]. Odessa, ONMA Publ., 2008, pp. 362-365 (in Russ.).

[5] Vorob'yeva E.A., Zubov N.E., Mikrin E.A., Misrikhanov M.Sh., Ryabchenko V.N., Tima-kov S.N. Synthesis of stabilizing spacecraft control based on generalized Ackermann's formula. Journal of Computer and Systems Sciences International, 2011, vol. 50, no. 1, pp. 93-103. Available at: http://link.springer.com/article/ 10.1134/S1064230711060189

DOI: 10.1134/S1064230711010199

[6] Zubov N.E., Mikrin E.A., Misrikhanov M.Sh., Oleynik A.S., Ryabchenko V.N. Terminal bang-bang impulsive control of linear time invariant dynamic systems. Journal of Computer and Systems Sciences International, 2014, vol. 53, no. 3, pp. 430-444.

Available at: http://link.springer.com/article/ 10.1134%2FS1064230714030174 DOI: 10.1134/S1064230714030174

[7] Gantmakher F.R. Teoriya matrits [Matrix theory]. Moscow, Nauka Publ., 1967. 576 p.

Zubov N.E. — Dr. Sci. (Eng.), Deputy Director for Science, Research and Development Centre, S.P. Korolev Rocket and Space Corporation Energia (Lenina ul. 4-a, Korolev, Moscow Region, 141070 Russian Federation), Professor of Automatic Control System Department, Bauman Moscow State Technical University (2-ya Baumanskaya ul. 5, Moscow, 105005 Russian Federation).

Mikrin E.A. — Dr. Sci. (Eng.), Academician of the Russian Academy of Sciences, General Designer of S.P. Korolev Rocket and Space Corporation Energia (Lenina ul. 4-a, Korolev, Moscow Region, 141070 Russian Federation), Head of Automatic Control System Department, Bauman Moscow State Technical University (2-ya Baumanskaya ul. 5, Moscow, 105005 Russian Federation).

Ryabchenko V.N. — Dr. Sci. (Eng.), Leading Researcher, Research and Development Centre, S.P. Korolev Rocket and Space Corporation Energia (Lenina ul. 4-a, Korolev, Moscow Region, 141070 Russian Federation), Professor of Automatic Control System Department, Bauman Moscow State Technical University (2-ya Baumanskaya ul. 5, Moscow, 105005 Russian Federation).

Please cite this article in English as:

Zubov N.E., Mikrin E.A., Ryabchenko V.N. Synthesis of Astatic Linear Systems Control Based on Generalized Ackermann's Formula. Vestn. Mosk. Gos. Tekh. Univ. im. N.E. Baumana, Priborostr. [Herald of the Bauman Moscow State Tech. Univ., Instrum. Eng.], 2017, no. 1, pp. 67-74. DOI: 10.18698/0236-3933-2017-1-67-74

В Издательстве МГТУ им. Н.Э. Баумана вышло в свет учебное пособие под редакцией А.И. Николаева

«Радиолокационные системы»

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

По вопросам приобретения обращайтесь:

105005, Москва, 2-я Бауманская ул., д. 5, стр. 1

+7 (499) 263-60-45

press@bmstu.ru

www.baumanpress.ru