МОТС ИТМО
.pdfУтверждение 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 ( ) ГR[ГR (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