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

DOI: 10.18698/0236-3933-2016-1-51-58

УДК 62-50

ЭКВИВАЛЕНТНОСТЬ УСЛОВИЙ УПРАВЛЯЕМОСТИ ЛИНЕЙНОЙ МНОГОМЕРНОЙ СИСТЕМЫ И РАЗРЕШИМОСТИ ПОЛИНОМИАЛЬНОГО МАТРИЧНОГО УРАВНЕНИЯ СИЛЬВЕСТРА*

Н.Е. Зубов12, Е.А. Микрин12, М.Ш. Мисриханов2, В.Н. Рябченко12

ХМГТУ им. Н.Э. Баумана, Москва, Российская Федерация e-mail: nezubov@bmstu.ru

2Ракетно-космическая корпорация "Энергия" им. С.П. Королева, Королев, Московская обл., Российская Федерация e-mail: Nikolay.Zubov@rsce.ru

Рассмотрены критерии управляемости линейной многомерной системы: ранговый критерий Калмана; модальный (частотный) критерий Попова-Белевича -Хотиса; критерий, когда для управляемости необходимо и достаточно невырожденности ленточной матрицы управляемости. Приведены утверждения, устанавливающие эквивалентность ленточных условий управляемости линейной многомерной системы и условий разрешимости линейного полиномиального матричного уравнения Сильвестра относительно полиномиальной матрицы степени n — 2. Даны дуальные утверждения для критериев наблюдаемости. На практическом примере показано применение разработанного подхода к линейным динамическим системам.

Ключевые слова: критерий управляемости, критерий наблюдаемости, линейные многомерные системы, полиномиальное уравнение Сильвестра.

EQUIVALENCE OF CONTROLLABILITY CONDITIONS OF LINEAR MULTIDIMENSIONAL SYSTEMS AND SOLVABILITY OF THE SILVESTER POLYNOMIAL MATRIX EQUATION

N.E. Zubov12, E.A. Mikrin12, M.Sh. Misrikhanov2, V.N. Ryabchenko12

1Bauman Moscow State Technical University, Moscow, Russian Federation e-mail: nezubov@bmstu.ru

2S.P. Korolev Rocket and Space Corporation Energia, Korolev, Moscow Region, Russian Federation e-mail: Nikolay.Zubov@rsce.ru

In the paper the following controllability criteria of linear multidimensional systems are considered: Kalman rank criterion; Popov-Belevich-Hautus modal (frequency) test, when nondegeneracy of the controllability band matrix is necessary and sufficient for a controllable system. The statements determining the equivalence of band controllability conditions of a linear multidimensional system and the solvability conditions of the Sylvester linear polynomial matrix equation relatively to polynomial matrix ofdegree are presented. The dual statements for the observability criteria are given. The application of the developed method to linear dynamic systems is shown in practice.

Keywords: controllability criterion, observability criterion, linear multidimensional system, Sylvester polynomial equation.

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

Введение. Для анализа управляемости линейной многомерной системы со стационарными параметрами

ах(€) = Ах(*) + ВЦ*), (1)

где х(*) € Мп — вектор состояния; Ц*) € М — скалярное управление; стх(£) = Х(*) — оператор дифференцирования по времени или стх(£) = х(* + 1) — оператор сдвига во времени; М — множество вещественных чисел. Используются два критерия (теста): 1) ранговый критерий Калмана; 2) модальный (частотный) критерий Попова-Белевича - Хотиса.

Согласно критерию Калмана, для управляемости (1) необходимо и достаточно, чтобы матрица

АП-1В

В | АВ |

была обратима [1], т.е. гапк ( В | АВ ( В I АВ I • • • I АП-1В ) = 0.

G Rnxn • • I An-1B

(2)

= n, или

В модальном (частотном) критерии Попова-Белевича-Хотиса, в котором для управляемости (1) требуется полнота ранга п матриц (сингулярного пучка матриц) [2] ( А — Лj 1п | В ) € Спх(п+1), где Xj € {Л € С : det (А — Л1п) = 0} — элементы спектра (множества собственных значений матрицы А); С — множество комплексных чисел (комплексная плоскость).

В работе [3] предложен критерий, когда для управляемости (1) необходимо и достаточно невырожденности ленточной матрицы управляемости размером п(п — 1) х п(п — 1):

/в fA 0 0 0 \

Bf BfA 0 0

0 Bf BfA 0

V о 0 0 BfJ

(3)

т.е.

rank

f Bf A 0 0 0 \

BLL BfA 0 0

0 Bf Bf A 0

0 0 0 BfJ

= n(n — 1),

или

det

/BfA 0 0 0 \

Bf в iA 0 0

0 Bf BfA 0

0 0 0 Bi/

= 0.

(4)

Здесь В^ € К(п-1)хп — левый делитель нуля вектора В, удовлетворяющий однородному уравнению [4] В ¿В = 0(п-1)х1, имеющий размеры (п — 1) х п и максимальный ранг п — 1.

Основная цель работы — установление эквивалентности условий управляемости линейной многомерной системы на основе ленточной матрицы (2) и условий разрешимости линейного полиномиального матричного уравнения. В силу дуальности свойств управляемости системы (1) и наблюдаемости системы

ах (г) = Ах (г), у (г) = С х(г), (5)

где у (г) € К — скалярный выход, конечные результаты по аналогии распространяются на задачу анализа наблюдаемости.

Управляемость. Линейным полиномиальным матричным уравнением относительно неизвестной матрицы У (А) называется уравнение [5]

(Гс + АГ 1) У (А) = 0, (6)

где Г г € ктхп, I = 1, 2 — заданные матрицы.

Матрица У (А) представляет собой правое полиномиальное нуль-пространство линейного пучка матриц Гс + АГ1 или, другими словами, правый делитель нуля этой матрицы. Если требуется определить решение матрицы У(А) в виде

У (А) = Ус + АУ1 + ... + А5(у ) У 5(у), (7)

где £(У) — заданное натуральное число (заданная степень), тогда от уравнения (6) переходят к эквивалентной записи [5, 6]

/ F о 0 0 \

F i F0 0

0 F i

F о

0 0 F i /

( Y о \

Yi

\ YS(Y) /

= 0.

Матрицу

S =

/ F о 0 о \

F i F о 0

о F1

F о

о 0 F i /

(8)

называют матрицей Сильвестра.

Справедлива лемма [7].

Лемма. Для разрешимости линейного полиномиального матричного уравнения (6) относительно матрицы Y(Л) с заданной степенью ) необходимо и достаточно неполноты ранга по столбцам матрицы Сильвестра (8).

Следовательно, полнота ранга по столбцам матрицы Сильвестра (8) гарантирует отсутствие решения уравнения (6) относительно матрицы Y(Л) (7) с заданным числом ).

Сравнивая матрицы (3) и (8), нетрудно прийти к выводу, что эти матрицы имеют совпадающую (эквивалентную) структуру. На основании этого факта, а также справедливости приведенной леммы сформулируем теорему.

Теорема 1. Для неуправляемости многомерной линейной системы (1) необходимо и достаточно разрешимости линейного полиномиального матричного уравнения

(В ¿А + ЛВ^) Y(Л) = 0 (9)

относительно матрицы

Y (Л) = Y о + ЛY1 + ... + Лп-2 Yn_2. (10)

Обратная теорема: для управляемости многомерной линейной системы (1) необходимо и достаточно неразрешимости линейного полиномиального матричного уравнения (9) относительно матрицы (10).

< Воспользовавшись сравнением матриц (3) и (8), сопоставим задаче (4) линейное полиномиальное уравнение (9). Ясно, что согласно лемме уравнение (9) будет разрешимо с заданной степенью п — 1, если и только если

' ' ' ' 0

/ВЫ

det

о

V 0

о

о

о

В А

о

о

о

о

= 0,

BtJ

но это соответствует неуправляемости линейной системы (1).

Сделаем замечание: поскольку левый делитель нуля В £ матрицы В определяется с точностью до некоторой обратимой матрицы Т размером (п — 1) х (п — 1) [4], то уравнению (9) может быть поставлено в соответствие множество уравнений вида Т (В£ А + АВ£) У (А) = 0. Это, очевидно, никак не влияет в алгебраическом смысле на условия разрешимости, но может оказать существенное влияние на вычислительные особенности задачи, поскольку в этом случае возникает возможность определенного влияния на обусловленность изучаемых матриц. ►

Наблюдаемость. Одна из модификаций ленточного критерия наблюдаемости многомерной линейной системы (5) имеет вид [3]: для наблюдаемости линейной системы (5) необходимо и достаточно, чтобы

.....0

/ AC

rank

R

0

0

V 0

C

R

AC

R

0

C

R

AC

R

0

0

0

= n(n — 1),

CR

или

det

ACR C R 0 0 \

0 AC R C R 0

0 0 AC R 0

0 0 0 C R )

= 0,

где матрица называется ленточной матрицей наблюдаемости и имеет размер п(п — 1) х п(п — 1); С^ — правый делитель нуля матрицы С, удовлетворяющий уравнению СС^ = 01х(п-1) и имеющий размер п х х (п — 1) и максимальный ранг п — 1 [4].

В силу дуальности задач управляемости и наблюдаемости [1, 2] приведем без доказательства еще одну теорему.

Теорема 2. Для ненаблюдаемости многомерной линейной системы (5) необходимо и достаточно разрешимости линейного полиномиального матричного уравнения

X (А) (АС£ + АС^) =0 (11)

(12)

относительно матрицы

X (А) = X о + АХ 1 + ... + АП-2Х п-

п-2.

Обратная теорема: для наблюдаемости многомерной линейной системы необходимо и достаточно неразрешимости линейного полиномиального матричного уравнения (11) относительно матрицы (12).

0

0

0 1 1

0 0 0

1 0 0

Пример. Пусть задана линейная система (1) с матрицами

А = ( 0 | 0 | 0 | , В = ( "Г

Вычислим левый делитель нуля матрицы В из (13): В£ = Теперь составим уравнение (9), в результате получим

(В ¿А + АВ ¿) У (А) =

0 0

1 0 0

10 0 0

1 1

0 0 Л

Л 0 0

+

где матрица У (А) имеет степень п — 2 = 1, т.е.

*01 \ / ^11 \ / *01 + А^11 + А

Y (Л) = 0, (14)

Y(Л) = I У

02

Y1

12

= I Y02 + ЛУ1

12

(15)

^03 / \ ^13 / \ ^03 + А^13

Уравнение (15) неразрешимо, поскольку ленточная матрица

/

B fA 0

Bf BfA

0 Bi

1 0 0 0

0 1 1

0 0 1 1 0 0

1 0 0 0 1 1

0 0 0 1

1 0 0

/

является невырожденной, что, в свою очередь, соответствует условию полной управляемости линейной системы с матрицами (13).

Отметим, что уравнение (14) оказывается разрешимым, но уже относительно матрицы У(А) степени п — 1 = 2:

_0_

Y(Л) = |

+ Л

т

о

0

Л

+ Л2

-1 I = | 1- Л2

0

что в таком случае является недостаточным.

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

ЛИТЕРАТУРА

1. Калман Р., Фалб П., Арбиб М. Очерки по математической теории систем. М.: Мир, 1971.

2. Zhou K. Essentials of robust control. Upper Saddle River, New Jersey: Prentice Hall, 1998.

3. Мисриханов М.Ш., Рябченко В.Н.Ленточная формула решения задачи Крылова //АиТ. № 12. 2011. С. 53-69.

4. МисрихановМ.Ш., Рябченко В.Н.Алгебраические и матричные методы в теории линейных MIMO-систем // Вестник ИГЭУ 2005. Вып. 5. С. 196-240.

5. Antsaklis P.J., Gao Z. Polynomial and rational matrix interpolation: theory and control applications // Int. J. Contr. 1993. Vol. 58 (2). P. 349-404.

6. Kailath T. Linear Systems. New Jersey: Prentice Hall. Englewood Cliffs. 1980.

7. Henrion D., Sebek MEfficient algorithms for polynomial matrix equation triangularization // Proc. of the IEEE Mediterranean Conference on Contr. and Automat, 1998. Alghero, Sardinia, Italy.

REFERENCES

[1] Kalman R.E., Falb P.L., Arbib M.A. Topics in mathematical system theory. McGraw-Hill, 1969.

[2] Zhou K. Essentials of robust control. Upper Saddle River. N.J., Prentice Hall, 1998.

[3] Misrikhanov M.Sh., Ryabchenko V.N. The band formula for A.N. Krylov's problem. Automation and Remote Control, 2007, vol. 68, no. 12, pp. 2142-2157.

[4] Misrikhanov M.Sh., Ryabchenko V.N. Algebraic and matrix methods in the theory of linear MIMO-systems. VestnikIGEU, 2005, vol. 5, pp. 196-240 (in Russ.).

[5] Antsaklis P.J., Gao Z. Polynomial and rational matrix interpolation: theory and control applications. Int. J. Contr, 1993, vol. 58 (2), pp. 349-404.

[6] Kailath T. Linear Systems. N.J., Prentice Hall. Englewood Cliffs, 1980.

[7] Henrion D., Sebek M. Efficient algorithms for polynomial matrix equation triangularization. Proc. of the IEEE Mediterranean Conference on Contr. and Automat. Alghero, Sardinia, Italy, 1998.

Статья поступила в редакцию 22.12.2014

Зубов Николай Евгеньевич — д-р техн. наук, заместитель руководителя по науке научно-технического центра РКК "Энергия" им. С.П. Королeва, профессор кафедры "Системы автоматического управления" МГТУ им. Н.Э. Баумана. МГТУ им. Н.Э. Баумана, Российская Федерация, 105005, Москва, 2-я Бауманская ул., д. 5.

РКК "Энергия" им. С.П. Королeва, Российская Федерация, 141070, Московская обл., Королeв, ул. Ленина, д. 4-а.

Zubov N.E. — Dr. Sci. (Eng.), Deputy Director for Science, Research and Development Center, S.P. Korolev Rocket and Space Corporation Energia; Professor of Automatic Control Systems department, Bauman Moscow State Technical University. Bauman Moscow State Technical University, 2-ya Baumanskaya ul. 5, Moscow, 105005 Russian Federation.

S.P. Korolev Rocket and Space Corporation Energia, ul. Lenina 4-a, Korolev, Moscow Region, 141070 Russian Federation.

Микрин Евгений Анатольевич — д-р техн. наук, первый заместитель генерального конструктора РКК "Энергия" им. С.П. Королeва, заведующий кафедрой "Системы автоматического управления" МГТУ им. Н.Э. Баумана.

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

РКК "Энергия" им. С.П. Королeва, Российская Федерация, 141070, Московская обл., Королeв, ул. Ленина, д. 4-а.

Mikrin E.A. — Dr. Sci. (Eng.), First Deputy of Chief Designer, Research and Development Center, S.P. Korolev Rocket and Space Corporation Energia, Head of Automatic Control Systems department, Bauman Moscow State Technical University. Bauman Moscow State Technical University, 2-ya Baumanskaya ul. 5, Moscow, 105005 Russian Federation.

S.P. Korolev Rocket and Space Corporation Energia, ul. Lenina 4-a, Korolev, Moscow Region, 141070 Russian Federation.

Мисриханов Мисрихан Шапиевич — д-р техн. наук, ведущий научный сотрудник научно-технического центра РКК "Энергия" им. С.П. Королева. РКК "Энергия" им. С.П. Королeва, Российская Федерация, 141070, Московская обл., Королeв, ул. Ленина, д. 4-а.

Misrikhanov M.Sh. — Dr. Sci. (Eng.), Senior Staff Scientist, Research and Development Center, S.P. Korolev Rocket and Space Corporation Energia.

S.P. Korolev Rocket and Space Corporation Energia, ul. Lenina 4-a, Korolev, Moscow Region, 141070 Russian Federation.

Рябченко Владимир Николаевич — д-р техн. наук, ведущий научный сотрудник научно-технического центра РКК "Энергия" им. С.П. Королeва, профессор кафедры "Системы автоматического управления" МГТУ им. Н.Э. Баумана. МГТУ им. Н.Э. Баумана, Российская Федерация, 105005, Москва, 2-я Бауманская ул., д. 5.

РКК "Энергия" им. С.П. Королeва, Российская Федерация, 141070, Московская обл., Королeв, ул. Ленина, д. 4-а.

Ryabchenko V.N. — Dr. Sci. (Eng.), Leading Researcher, Research and Development Center, S.P. Korolev Rocket and Space Corporation Energia, Professor of Automatic Control Systems department, Bauman Moscow State Technical University. Bauman Moscow State Technical University, 2-ya Baumanskaya ul. 5, Moscow, 105005 Russian Federation.

S.P. Korolev Rocket and Space Corporation Energia, ul. Lenina 4-a, Korolev, Moscow Region, 141070 Russian Federation.

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

Зубов Н.Е., Микрин Е.А., Мисриханов М.Ш., Рябченко В.Н. Эквивалентность условий управляемости линейной многомерной системы и разрешимости полиномиального матричного уравнения Cильвестра // Вестник МГТУ им. Н.Э. Баумана. Сер. Приборостроение. 2016. № 1. C. 51-58. DOI: 10.18698/0236-3933-2016-1-51-58

Please cite this article in English as:

Zubov N.E., Mikrin E.A., Misrikhanov M.Sh., Ryabchenko V.N. Equivalence of controllability conditions of linear multidimensional systems and solvability of the silvester polynomial matrix equation. Vestn. Mosk. Gos. Tekh. Univ. im. N.E. Baumana, Priborostr. [Herald of the Bauman Moscow State Tech. Univ., Instrum. Eng.], 2016, no. 1, pp. 51-58. DOI: 10.18698/0236-3933-2016-1-51-58