Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

МОТС ИТМО

.pdf
Скачиваний:
77
Добавлен:
14.04.2015
Размер:
11.04 Mб
Скачать

Утверждение 12.1 (У12.1). Введем в рассмотрение обобщенное уравнение Риккати

AT P PA PBVBT P Q ,

(12.1)

где V V T , V 0.

 

 

Тогда уравнение Ляпунова

 

 

F T P PF Q

(12.2)

с матрицами F A BK и

F T AT K T BT , где K R 1BT P ,

совпадает с матричным уравнением Риккати (12.1) при V 2R 1 .

Доказательство. Для доказательства У12.1 достаточно в (12.1)

положить V 2R 1 , в (12.2)

подставить F A BR 1BT P

и

F T AT PBR 1BT , и убедиться в совпадении левых и правых частей

(12.1) и (12.2).

 

Утверждение 12.2 (У12.2). Если уравнение Риккати используется для синтеза закона управления, доставляющего системе качественную экспоненциальную устойчивость, так, что оно записывается в виде

AT P PA PBVBT P P ,

 

 

(12.3)

то оно может быть сведено к уравнению Ляпунова вида

 

( A

1

1I )T P 1 P 1 ( AT

1

1I ) BVBT .

 

(12.4)

 

 

 

2

 

2

 

 

 

Доказательство. Умножим уравнение Риккати (12.3) слева и

справа на матрицу P 1 , тогда получим

 

 

P 1 AT AP 1 BVBT P 1

 

 

В полученном

матричном

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

P 1 все члены ее содержащие разместим в левой части уравнения, а не

содержащие ее – в правой, тогда получим (12.4).

 

Утверждение 12.3 (У9.3).

Матрицы P и

Q : P PT ,

P 0;

Q QT , Q 0 вида

 

 

 

 

 

P (M 1 )T M 1 , Q (M 1)T (Γ ΓT )M 1 ,

 

(12.5)

где матрицы M

и Γ удовлетворяют матричному уравнению

Сильвестра

 

 

 

 

 

M Γ AM BH

 

 

 

 

(12.6)

удовлетворяют

матричному уравнению

Ляпунова

(12.2).

 

 

 

 

 

Доказательство. Пусть матрица K обратной связи найдена в

результате решения задачи модального управления так, что K HM 1 , где M – решение уравнения Сильвестра (12.6).

Тогда подстановка H KM в (12.6) приводит к уравнению подобия

M Γ FM ,

(12.7)

которое позволяет записать

 

219

F M M 1 , FT (M 1)T ΓT M T .

(12.8)

Подстановка в (12.2) матриц P и Q вида (12.5), а также F и F T

вида (12.8) дает

 

FT P PF (M 1 )T ΓT M T(M 1 )T M 1 (M 1 )T M 1M M 1

(12.9)

,

(M 1)T (Γ ΓT )M 1 Q

 

что доказывает справедливость утверждения 12.3.

Примечание 12.1 (П12.1). Необходимо заметить, что матрица

Γ ΓT , в которой матрица Γ имеют спектр собственных значений

i : Re λi 0;i

 

с отрицательными вещественными частями, не в

1,n

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

Эта матрица является гарантировано положительно определенной,

если матрица

Γ

является диагональной

Γ Λ или

блочно –

диагональной

Γ

~

 

 

 

Λ .Действительно в этом случае имеет место

представление Γ ΓT diag 2Re i ;i

 

.

 

1,n

 

Утверждение 12.4 (У12.4). Матрицы P и Q вида

 

P (M 1 )T M 1 , Q (M 1 )T ( T ) H T BR 1BT H M 1, (12.10)

где матрицы , M , H и R удовлетворяют модифицированному

УС

 

 

 

 

 

M AM BR 1BT H , HM T I ,

 

(12.11)

вкотором ( A, B) – управляема, ( , R 1BT H ) – наблюдаема, а также

A Ø, удовлетворяют обобщенному уравнению Риккати

(12.1) при V R 1 . □

Доказательство. Для доказательства У12.4 сконструируем на основании (12.11) представления матриц A и AT

A M M 1 BR 1BT (M 1)T M 1,

 

(12.12)

AT (M 1)T ΓT M T (M 1)T M 1BR 1BT .

 

 

Тогда подстановка в (12.1) V R 1 , а также представлений

A и

AT (12.12) дает

 

 

 

 

AT P PA PBR 1BT P (M 1)T ΓT M T (M 1)T M 1BR 1BT (M 1)T M 1

(M 1)T M 1 M M 1 BR 1BT (M 1)T M 1

 

 

(M 1)T M 1BR 1BT (M 1)T M 1 (M 1)T (Γ ΓT ) H T BP 1BH M 1 Q,

что доказывает справедливость сформулированного утверждения. ■

Примечание 12.2

(П12.2).

Управление

u K x ,

где

K R 1BT HM 1 является

модальным

управлением,

оптимальным в

смысле квадратичного функционала качества (15.17). Это управление доставляет матрице F A BK состояния системы структуру мод,

220

носителем которой является матрица Γ , и минимизирует функционал (15.17) с матрицами P и Q вида (12.10) на траекториях системы.

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

уравнения вида

 

MS RM U

(12.13)

или

 

MS RM VW ,

(12.14)

в которых соблюдаются соотношения размерностей матричных

компонентов

 

dimM dimS dimR dimU dim VW n n .

(12.15)

12.2 Методы прямого решения линейных матричных уравнений

Под прямым решением уравнений (12.13) и (12.14) понимается

решение этих матричных уравнений относительно

матрицы M

соответственно в формах

 

 

M M S, R,U ,

M M S,W , R,V .

(12.16)

Первым аналитическим способом решения матричного уравнения вида (12.13) является способ, основанный на сведении матричного уравнения к векторно–матричному с использованием

кронекеровской суммы матриц R и S T , в результате чего становится

справедливой запись

 

 

 

 

 

 

 

 

 

 

~ ~

~

 

 

 

 

 

 

 

 

(12.17)

 

 

TM U ,

 

 

 

 

 

 

 

 

 

 

~

 

T

R I I

S

T

,

 

 

(12.18)

 

 

где T R S

 

 

 

 

 

 

 

 

, - символы кронекеровских суммы и произведения матриц,

n

2

 

 

 

 

~

~

 

 

 

 

 

 

 

 

мерные векторы M и U имеют структуру

 

 

 

 

 

~

T

 

T

T T

~

 

T

T

T T

(12.19)

 

 

M M1

, M 2

,...,M n

, U

U1

,U2 ,...,Un ,

 

 

 

 

 

 

M и U, i

 

.

 

 

Mi ,Ui

i

столбец соответственно

матриц

1, n

Решение уравнения (12.17) находится стандартным способом в форме

~

~ 1 ~

 

 

 

(12.20)

M

T U .

~

 

 

Из

столбца

формируется

искомая

матрица

M

M M1 M 2 M n .

 

 

 

 

Примечание 12.3 (П12.3). Ключевым моментом реализуемости

предлагаемого способа

решения является обратимость

~

матрицы T

221

размерности n2 n2 , которая, являясь кронекеровской суммой матриц R и ST , имеет своими собственными значениями попарные суммы

собственных

 

значений

матриц

R и S T .Таким образом обратимость

матрицы

~

будет

обеспечена, если алгебраические спектры

T

 

собственных значений матриц S

и R не будут пересекаться, то есть

будет выполняться соотношение S R Ø.

Следует также заметить, что предложенный способ не гарантирует

невырожденности матрицы M , так как она «собирается» из элементов

~

вектора M .

В связи с отмеченным выше прежде, чем продолжить рассмотрение способов прямого решения матричных уравнений вида (12.13) и (12.14) сформулируем условие существования единственного невырожденного решения этих уравнений, опираясь на представление (12.14) с декомпозированной правой частью, ограничившись случаем матриц S и R простой структуры.

Утверждение 12.5 (У9.5). Линейное матричное уравнение вида (12.14) имеет единственное невырожденное решение, если матричные компоненты этого уравнения удовлетворяют условиям:

непересекаемость алгебраических спектров собственных

значений матриц S

и R , то есть

выполнение соотношения

S R Ø;

 

 

 

управляемость пары матриц R ,V ;

 

 

наблюдаемость пары матриц S ,W .

 

 

Доказательство.

Необходимость

выполнения

условия

непересекаемости алгебраических спектров собственных значений матриц S и R доказана при рассмотрении первого способа решения матричных уравнений вида (12.14).Тем не менее, осуществим еще раз доказательство необходимости выполнения этого условия в контексте подготовки новых способов решения матричного уравнения (12.14).

Для

этих

целей

предположим, что матрица

S

задана в

 

 

S S diag si ; i

 

так,

 

 

диагональной

форме

1,n

что

уравнение

(12.14)

принимает

вид M S RM VW .Запишем

последнее

матричное уравнение в виде системы векторно – матричных уравнений вида

 

 

 

 

 

 

 

 

M Si RMi

VWi ; i 1,n.

 

 

 

(12.21)

Столбец Si

имеет

вид Si 0Ti 1 | si | 0Tn i T , подстановка

которого в

 

векторно–матричное уравнение

(12.21) дает

 

 

 

 

 

 

 

( si R)Mi VWi , i 1, n ,

 

 

 

 

откуда для матрицы M получим представление

 

M row M

 

 

 

R 1VW ; i

 

.

(12.22)

i

si

1,n

 

 

 

 

i

 

222

Решение матричного уравнения (12.14) в форме (12.22) возможно только при обратимости матрицы ( si R) , что возможно только при

выполнении условия S R Ø.

Покажем теперь справедливость требования управляемости

пары матриц

R ,V . Для этих

целей

воспользуемся

разложением

Д.Фаддеева матрицы (

si

R) 1 ,

в соответствии с которым получаем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

развернутое по степеням матрицы R представление

 

 

 

 

 

Ι R 1

1

D

(

 

)Ι D

(

 

)R D

 

(

 

)R2 D (

 

)Rn 1 ,(12.23)

si

 

si

si

 

si

si

 

 

D( si )

n 1

 

 

 

n 2

 

n 3

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где

D(

si

) det(

si

Ι R) n

a n 1

a

 

n 2 a

 

 

si

a

n

;

 

 

 

 

 

 

 

 

 

 

 

 

 

si

1 si

2

 

si

n 1

 

 

 

 

 

 

 

 

 

 

n 1

 

 

 

 

n 2

an 2 si an 1;

 

 

 

 

 

 

Dn 1 ( si ) si

 

a1 si

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n 2

 

 

 

 

n 3

an 3 si an 2 ;

 

 

 

 

 

 

 

Dn 2 ( si ) si

 

a1 si

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

D1 ( si ) si a1;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

D (

si

) 1.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Введем обозначения

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

D

j

(

si

)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(

 

)

 

; j 0, n 1.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

D

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

j

si

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

D( si )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

В

 

 

результате

 

 

 

для

(

si

R) 1

получаем

 

 

 

si Ι R 1

| R | R2 | | Rn 1 col

 

j ( si ); j

 

.

 

 

 

 

 

 

n 1,0

 

 

 

 

 

D

 

 

Подстановка (12.26) в (12.22) дает для матрицы M

M V | RV | | Rn 1V row col Dj ( si )Wi ; j n 1,0 ; i 1, n ..

(12.24)

(12.25)

запись

(12.26)

(12.27)

Из (12.27) следует, что матрица M оказывается невырожденной,

то есть M 1 , если

rang

V | RV | | Rn 1V n ,

то

есть

если пара

матриц R , V управляема.

 

 

 

 

 

 

 

 

 

 

Для целей дальнейших исследований транспонируем матричное

уравнение (12.14), тогда получим

 

 

 

M T RT ST M T W TV T .

 

 

 

 

 

 

 

 

(12.28)

Теперь сделаем предположение, что матрица R задана в

 

 

R R diag Ri ; i

 

 

 

 

 

диагональной

форме

1,n

так,

что

уравнение

(12.28) принимает вид M T

R

ST M T W TV T .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Запишем последнее матричное уравнение в виде системы

векторно – матричных уравнений вида

 

 

 

M T

 

S T M T W TV T ;i

 

.

 

 

(12.29)

Ri

1, n

 

 

 

i

 

 

 

i

 

 

 

223

 

 

 

Столбец

 

Ri

 

 

 

 

имеет

 

 

вид

 

 

 

 

 

 

Ri

0Ti 1 |

Ri | 0Tn i T , подстановка

которого

 

 

в

 

 

векторно

 

 

 

 

 

матричное

 

уравнение

 

 

 

(12.29) дает

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

M T получим

(

Ri

ST )M

i

W TV T , i

1, n

,откуда

 

 

 

 

 

 

 

для

 

 

 

матрицы

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

M T row M T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i 1, n .

 

 

 

 

 

 

 

 

 

 

 

(12.30)

 

 

 

Ri

ST

W TV T ;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Решение матричного уравнения (12.14) в форме (12.30) возможно

только при обратимости матрицы (

Ri

S T ) , что возможно только при

выполнении условия S R

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ø.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Покажем теперь справедливость требования наблюдаемости

пары матриц

 

S ,W . Для этих целей

 

 

воспользуемся

разложением

Д.Фаддеева матрицы (

Ri

ST ) 1 ,

в соответствии с которым получаем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

развернутое по степеням матрицы S T

 

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ι S T 1

 

1

 

 

D

 

 

(

 

 

)Ι D

 

 

 

 

(

 

 

 

)S T D

 

 

 

 

(

 

)S T 2

D (

 

)S T n 1 ,

 

 

 

Ri

 

 

 

 

 

 

n 1

Ri

n 2

Ri

n 3

Ri

Ri

 

 

 

 

 

 

 

 

 

 

 

 

D( Ri )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

D(

 

 

) det(

 

 

Ι S T ) n

 

a n 1

a

 

n 2

 

a

n 1

 

 

 

 

a

 

;

 

 

 

 

 

 

 

 

Ri

 

 

 

 

 

 

Ri

 

 

 

 

 

 

 

 

 

 

 

Ri

 

 

1

 

Ri

 

 

 

 

 

 

 

2

Ri

 

 

 

 

 

 

 

 

 

Ri

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n 1

 

 

 

 

 

n 2

an 2 Ri an 1;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Dn 1 ( Ri ) Ri

 

a1 Ri

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n 2

 

 

 

 

 

 

n 3

an 3 Ri an 2 ;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Dn 2 ( Ri ) Ri

 

 

a1 Ri

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.(12.31)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

D1 ( Ri ) Ri a1;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

D (

Ri

) 1.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Введем обозначения

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

D

j

(

Ri

)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(

 

)

 

 

 

 

 

 

 

0, n 1.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

D

 

 

 

 

 

 

 

 

 

 

 

;

 

 

 

j

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(12.32)

 

 

 

 

 

 

 

 

 

 

 

 

 

j

Ri

D( Ri )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

В

 

 

 

 

результате

 

 

 

 

для

 

 

 

 

 

(

Ri

ST ) 1

 

 

получаем

 

 

запись

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ri

Ι ST

 

 

 

| ST

 

| ST 2 | | ST n 1 col D

j

(

Ri

);

 

 

j n 1,0

 

 

 

 

(12.33)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Подстановка (12.33) в (12.30) дает для матрицы M T

 

 

 

 

 

 

 

 

 

M T W T | STW T | | ST n 1 W T row col

 

 

 

 

 

 

 

 

)V T ; j

 

; i

 

..

 

 

 

 

 

 

 

(

 

 

 

 

n 1,0

1,n

 

 

(12.34)

D

j

Ri

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Из (12.34) следует, что матрица M оказывается невырожденной,

то есть

 

 

M 1 , если

rang

W T | STW T | | ST n 1 W T n ,

то есть если

пара матриц S , W наблюдаема.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Примечание 12.4 (П12.4).В процессе доказательства утверждения 12.5 сформированы два аналитических способа решения матричного

224

уравнения (12.14), представленных соотношениями (12.22) и (12.30), которые соответственно названы вторым и третьим.

Четвертый (поэлементный) способ решения уравнений (12.13), (12.14).

При этом способе решения уравнения (12.13) строятся n2 скалярных соотношений вида

M i S

j

Ri M

j

U

ij

,

 

 

 

(12.35)

 

 

 

 

 

 

 

 

 

где M i , Ri i

строки матрицы M и R, S

j

,M

j

j -е столбцы

 

 

 

 

 

 

 

 

 

 

матриц S и M , Uij (i, j) -ый элемент матрицы U . В результате приведения подобных членов получается система содержащая n2

условий для n2

неизвестных

M ij

матрицы M, которая

решается

стандартными методами линейной алгебры.

 

Пятый способ (с использованием итеративной процедуры)

Рассмотрим

матрицу

Z (t) ,

определяемую

решением

дифференциального матричного уравнения

(12.36)

Z (t) RZ ZS U ,

 

Z (0) .

 

 

 

 

 

 

 

 

 

Прямой подстановкой в выражение (12.36) нетрудно убедиться,

 

 

 

 

 

t

 

что решение уравнения имеет вид Z (t) eRt Z (0)eSt eR UeS d .

 

 

 

 

 

0

 

Если lim eRt Z (0)eSt

0, то очевидно, что

 

t

 

 

 

 

 

 

 

 

 

 

 

 

 

M lim Z (t) eR UeS d .

 

(12.37)

t

0

 

 

 

 

 

 

 

 

 

 

 

Построим суммарное приближение интегрального выражения

(12.37) в виде

 

 

 

 

 

 

 

 

 

 

 

 

 

eR UeS d lim h

eRhkUeShk .

 

(12.38)

0

h 0

k 0

 

 

 

 

 

 

 

Если для матричных экспонент использовать аппроксимацию Е.

Девисона:

 

 

 

 

 

 

eRh ГR [12I 6hR h2R2 ] 1[12I 6hR h2R2 ]

 

eSh ГS [12I 6hS h2S 2 ] 1[12I 6hS h2S 2 ],

 

то получим

 

 

 

 

 

 

eRhk ГK , eShk Гk .

 

 

 

R

 

S

 

 

 

Для суммарного представления интеграла можно записать

 

 

 

 

 

 

 

eR UeS d lim hГkRUГkS .

 

 

 

0

h 0 k 0

 

 

 

 

Введем в рассмотрение частичную сумму

225

M ( ) hГkRUГkS h[U ГRUГS ГRUГS ]

k 0

На базе частичной суммы нетрудно построить рекуррентную процедуру

M ( 1) M ( ) Г 1(hU 1,

M (0) hU

(12.39)

 

R

S

 

 

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

перехода

 

 

 

 

M lim

M ( )

 

 

 

 

 

 

 

 

Алгоритм 12.1 (A12.1)

решения уравнения (12.13) пятым способом принимает следующий вид:

1.Выбор шага дискретизации h из условия

h

 

, max max( i , j ),

 

 

 

 

| max |

 

 

i, j

 

 

 

 

 

где

i – собственные значения матрицы R;

j – собственные

значения матрицы S; – константа;

2.Вычисление матриц ГR , ГS ;

3.Вычисление матрицы M (0) hU ;

4.Для шагов рекуррентной процедуры 1,2, вычисление

матрицы

M( 1) M ( ) ГRR (hU S S ;

5.Останов рекуррентной процедуры на основании удовлетворения

неравенству

( )M ( 1) M ( )

где рекомендуется выбирать в пределах от 10-7 до 10-5; в качестве правила останова можно использовать условие

m

tr ( ) ii ( ) .

i 1

Рассмотренные способы решения матричных уравнений (12.13), (12.14) не исчерпывают весь банк способов, с некоторыми из них можно ознакомиться в приводимой литературе.

Следует также добавить, что в случае, если правые части матричных уравнений (12.13) и (12.14) оказываются нулевыми, то они превращаются в однородные матричные уравнения тем самым становятся матричными условиями подобия матриц S и R , поэтому

алгебраическим условием нетривиального решения этих уравнений в оговоренном случае является совпадение алгебраических спектров собственных значений матриц S и R .

Также необходимо отметить проблему вычислительной устойчивости рассмотренных матричных уравнений, которую можно проконтролировать с помощью C МУ – числа обусловленности

226

«матричного уравнения» вида (12.13), (12.14), вычисляемого в силу соотношения

C МУ max

S

 

R

 

min

 

S

 

 

 

1

(12.40)

j

 

j

R

,

i, j

i

 

 

 

i

 

 

 

 

 

 

 

 

 

i, j

 

 

 

 

 

 

 

где - сингулярное число матрицы .

 

 

 

 

 

В заключение следует сказать, что в настоящее время самым

эффективным способом

решения

матричных

уравнений

(12.13)

(Ляпунова) и (12.14) (Сильвестра) является способ, использующий оператор «lyap» в программной оболочке MatLab, который для уравнения (12.13) имеет запись M lyap S, R,U .

12.3 Инверсное решение линейных матричных уравнений

В основу инверсного способа решения матричных уравнений вида (12.13),(12.14) положено то соображение, что матрица М есть

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

4.

Под инверсным решением уравнений (12.13) и (12.14) понимается

решение этих матричных уравнений относительно матрицы U

в форме

U U S, R, M

(12.41)

применительно к уравнению (12.13) и в формах

 

W W S, R,V , M ,

(12.42)

V V S, R,W , M

(12.43)

применительно к матричному уравнению (12.14). Причем инверсное решение матричного уравнения (12.14) в форме (12.42) используется в основном в задачах управления, а в форме (12.43) – в

задачах наблюдения.

Инверсный способ решения матричных уравнений вида (12.13),(12.14) не всегда реализуем. Сформулируем условия реализуемости инверсного способа решения матричных уравнений вида (12.13),(12.14).

Применительно к матричному уравнению (12.13), которое в задачах исследования устойчивости и динамических свойств систем при стохастических экзогенных воздействиях стационарных в широком смысле принимает вид матричного уравнения Ляпунова, эти условия принимают вид.

Условие 12.1 (УС12.1). Условие разрешимости матричного уравнения (12.13), используемого в задачах исследования устойчивости непрерывных систем в форме уравнения Ляпунова,

227

состоит в том, чтобы матрица MS RM размерности n n была

симметричной и положительно определенной, то есть имела n

положительных собственных значений.

Условие 12.2 (УС12.2). Условие разрешимости матричного уравнения (12.13), используемого в задачах исследования динамических свойств непрерывных систем при стохастических экзогенных воздействиях стационарных в широком смысле в форме уравнения Ляпунова, состоит в том, чтобы матрица

MS RM размерности n n была симметричной, положительно полуопределенной и имела ранг, равный рангу матрицы входа системы.

Применительно к матричному уравнению (12.14), которое в задачах синтеза формирователей сигналов управления и устройств динамического наблюдения принимает вид матричного уравнения Сильвестра, эти условия принимают вид.

Условие 12.3 (УС12.3). Условие разрешимости матричного уравнения (12.14) в форме (12.42), используемого в задачах синтеза формирователей сигналов управления в форме уравнения Сильвестра,

состоит в том, чтобы столбцы MS RM i матрицы MS RM размерности n n принадлежали образу матрицы V

MS RM i Jm V ;i

 

.

 

1,n

(12.44)

В случае выполнения условия 12.3 в форме (12.44) решение

(12.42) уравнения (12.14) относительно матрицы

W представимо в

форме

 

1

(12.45)

W V TV V T (MS RM ) .

Условие 12.4 (УС12.4). Условие разрешимости матричного уравнения (12.14) в форме (12.43), используемого в задачах синтеза устройств динамического наблюдения в форме уравнения Сильвестра,

состоит в том, чтобы столбцы матрицы M T RT ST M T размерности

n n принадлежали образу матрицы W T

 

 

M T RT ST M T Jm W T ;i 1,n

.

 

(12.46)

i

 

 

 

 

 

В случае выполнения условия 12.4 в форме (12.46) решение

(12.43) уравнения (12.14) относительно матрицы V представимо в

форме

 

 

 

 

 

1

.

 

 

 

(12.47)

V (MS RM )W T WW T

 

 

 

Примеры и задачи

 

 

 

 

 

12.1.*Определить, разрешимо

ли

матричное

уравнение

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

M,

если матрицы

S, R и U

принимают значения

228

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]