МОТС ИТМО
.pdfимеет вид
(9.44)
Как и в непрерывном случае рассмотрим три возможные структурыA алгебраического спектра собственных значений матрицы A .
1.Простая вещественная структура собственных значений матрицы: |
|||||||||||||||||||||||||||||||||||||
|
|
|
|
|
i : det |
|
|
I |
|
0; |
|
i |
|
j ;i j; Jm |
|
|
i, j 0;i, j 1, n . В этом |
||||||||||||||||||||
A |
|
|
A |
|
|
|
|||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
случае матрица A |
|
|
приводима к диагональному представлению |
||||||||||||||||||||||||||||||||||
|
|
diag |
|
ι ;i |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
1, n |
с помощью матрицы M преобразования подобия в |
|||||||||||||||||||||||||||||||||
|
|
|
|
|
|
1 |
|
diag |
|
|
;i |
|
|
|
1 . Воспользуемся свойством |
||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
1, n |
||||||||||||||||||||||||||||
форме |
A |
M |
M |
M |
|
ι |
M |
||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
сохранения матричного условия подобия для матричных функций от |
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
в форме f |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1.Тогда становится справедливым |
|||||||||||||||||||||||||||||||||||||||||||||||||||||
матриц |
|
|
A |
|
Mf |
|
M |
|||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
представление решения (9.44) в виде |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||||||||||
|
|
k |
|
|
|
|
|
|
|
|
|
|
|
|
k |
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
k |
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||
x k A |
|
|
|
|
|
|
|
|
|
|
|
|
M |
x 0 Mdiag |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||||||||||
|
|
x 0 M |
|
|
|
|
i ;i 1, n M |
|
|
x 0 , |
|
|
|
|
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
k |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||
y k C x k C Mdiag i ;i 1, n M |
|
|
|
|
x 0 . |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
x k |
|
|||||||||||||||||||||||||||||||||||||||||||||||||||||||
Переход к векторным и матричным нормам для нормы |
|
|
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
вектора состояния x k |
c учетом свойств норм позволяет записать |
|||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
x k |
|
|
|
|
|
diag |
|
ki ;i |
|
|
|
1 x 0 |
|
|
|
|
|
diag |
|
ki ;i |
|
|
|
|
|
1 |
|
|
|
x 0 |
|
|
|
|
|||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
M |
|
1, n |
M |
|
|
M |
|
|
1, n |
|
|
M |
|
|
|
|
|
x 0 |
|||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
C |
|
|
|
|
|
|
kM |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
M |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(9.45) |
2. Простая комплексно–сопряженная структура собственных значений матрицы:
|
|
|
|
|
: det |
|
|
|
|
|
0; |
|
|
|
|
|
|
|
|
|
|
|
|
|
l |
|
cos l |
j sin l ;l 1, n 2 . В |
|||||||||||||||
A |
|
l |
|
I |
A |
|
l |
l |
j |
l |
|
|
|
||||||||||||||||||||||||||||||
этом |
случае |
|
матрица |
|
|
|
|
|
|
|
|
|
приводима |
к |
блочно–диагональному |
||||||||||||||||||||||||||||
|
|
|
|
A |
|
||||||||||||||||||||||||||||||||||||||
представлению |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
cos l |
sin l |
|
|
|
|
|
||||||||||||||||||
|
|
|
|
|
|
|
|
|
l |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||
~ |
|
|
|
~ |
|
|
l |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||
diag |
|
|
l |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
l |
|
|
|
|
|
|
|
|
|
|
|
|
;l 1, n 2 |
|
с |
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
l |
|
|
|
|
|
|
|
|
|
|
|
|
sin |
|
cos |
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||
|
|
|
|
|
|
|
|
l |
|
~ |
|
|
|
|
|
|
|
|
|
|
|
l |
|
|
|
|
|
l |
|
|
|
|
|||||||||||
помощью |
|
|
матрицы |
|
|
|
|
|
|
|
|
|
преобразования |
|
подобия |
в форме |
|||||||||||||||||||||||||||
|
|
|
|
M |
|
|
|
|
|
|
~ |
~ |
~ |
~ |
|
|
|||||
A M M 1 |
Mdiag |
||||
|
|
|
|
|
|
~ l
l l
|
|
|
|
|
|
|
|
|
~ |
|
Тогда |
||||
|
|
|
|||||
l |
;l 1, n 2 |
M 1. |
|||||
l |
|
|
|
|
становится справедливым представление решения (9.44) в виде
|
|
|
|
|
~ |
|
|
~ |
|
|
|
|
|
|
k cos k l |
|
|
sin k l |
|
|
|
|
|
|
~ |
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||
|
|
|
|
|
l |
|
|
|
||||||||||||||||||||||||||||
x k A k x 0 Mdiag |
kl |
|
sin k |
|
|
|
cos k |
|
|
;l 1, n 2 |
M |
1x 0 , |
||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
l |
l |
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
~ |
|
|
|
|
|
|
|
k cos k |
|
|
|
|
sin k |
|
|
|
|
|
|
~ |
|
|
||||||||
|
|
|
|
|
~ |
|
|
|
|
|
|
|
|
|
|
|
|
1x 0 . |
||||||||||||||||||
y k C x k |
CMdiag |
kl |
l |
|
|
|
l |
|
|
|
|
|
l |
;l 1, n 2 |
|
M |
|
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
sin k |
l |
|
cos k |
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
l |
|
|
|
|
|
Если теперь перейти к векторным и матричным нормам, то для нормы x k вектора состояния x k c учетом свойств норм в том
149
числе и того, что |
|
cos k |
l |
|
sin k |
l |
|
|
|
|
1 при |
|
|
x 0 |
|
1 получим |
||||||||||||||||||||||||
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
sin k |
l |
|
cos k |
l |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
~ |
|
~ |
|
|
|
|
k cos k l |
sin k l |
|
|
|
|
|
~ |
|
|
~ |
|
|
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||
|
x k |
|
|
|
|
|
l |
|
|
|
|
|
C M kM . |
|||||||||||||||||||||||||||
|
Mdiag |
kl |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
;l |
1, n 2 |
M 1 |
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
sin k |
|
cos k |
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
l |
|
|
|
|
|
|
|
|
l |
|
|
|
|
|
|
|
|
|
(9.46)
3. Кратная вещественная структура собственных значений (простоты
ради будем как и в непрерывном случае полагать кратность |
n ) |
||||||||||||||
|
|
|
|
|
i : det |
|
I |
|
0; Jm 0;i 1, n . В |
|
|||||
матрицы: |
A |
|
|
A |
этом |
||||||||||
случае |
матрица |
|
|
|
|
приводима |
к |
жордановой |
форме |
||||||
|
A |
||||||||||||||
J col 1 0 0, 0 |
1 0, , 0 |
0 |
1, 0 0 0 |
. |
Тогда становится справедливым представление решения (9.44) в виде
|
|
|
|
|
|
|
1 |
k k k 1 2! 2 |
k k 1 k n 2 n 1 ! n 1 |
||||
|
|
|
|
|
|
|
|
|
|
k |
k k 1 k n 3 n 2 ! n 2 |
||
|
|
|
|
|
|
|
0 |
1 |
|||||
|
|
|
|
|
|
|
|
|
|
|
k k 1 k n 4 n 3 ! |
n 3 |
|
|
|
|
|
|
|
|
|
0 |
0 |
1 |
|
||
|
|
|
|
|
|
||||||||
x k TJ k T 1x 0 T k |
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
0 |
0 |
|
k |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
0 |
0 |
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
k 1
y k Cx k CTJ T x 0
|
|
|
|
|
|
|
|
T |
1x 0 . |
|
|
|
|
|
|
|
|
|
|
|
|
|
Если теперь в полученном выражении перейти к векторным и |
||||||||||||||||||||||||||
матричным нормам, то для нормы |
|
|
x k |
|
|
|
вектора состояния |
x k c |
|||||||||||||||||||
|
|
|
|
||||||||||||||||||||||||
учетом свойств норм |
при |
|
|
|
x 0 |
|
|
|
|
1 по аналогии с непрерывным |
|||||||||||||||||
|
|
|
|
||||||||||||||||||||||||
случаем при кратных вещественных собственных числах получим |
|||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
n 1 |
|
|
k! |
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
x k |
|
TJ k T 1x 0 |
C T k |
|
|
|
|
|
|
|
i! i . |
(9.47) |
||||||||||||||
|
|
|
|
|
|
|
|
||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
i 0 |
|
k 1 i ! |
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Анализ решений автономной системы (9.43) для трех возможных структур алгебраического спектра собственных значений матрицы состояния позволяет сделать вывод о том, что необходимым и достаточным условием выполнения предельных переходов для
зависимых переменных в виде векторов состояния и выхода |
||
дискретного объекта lim x x 0 , k 0 и lim y x 0 , k 0 |
является |
|
k |
k |
|
удовлетворение значений модулей всех собственных чисел неравенству i 1; i 1, n . Это же неравенство является условием устойчивости исходной модели ЛДДО (9.18)–(9.19).
150
9.4. Структурные свойства объектов управления: управляемость и наблюдаемость, управляемость и наблюдаемость собственных значений матрицы состояния над бесконечными полями. Каноническое структурное представление Р.Калмана. Полнота моделей ВСВ и ВВ динамических объектов
Остановимся на важных структурных свойствах динамических объектов как непрерывных, так и дискретных, которые удалось установить специалистам по теории систем лишь при изучении моделей ВСВ этих объектов. Ниже ограничимся рассмотрением двух базовых структурных свойств объектов управления (ОУ):
управляемости и наблюдаемости.
Определение 9.7 (О9.7). Непрерывный/дискретный динамический |
|
объект с парой матриц (A, B) / (A, B) |
называется полностью |
управляемым, если его можно из произвольного начального состояния
x0 x(t) |t t0 / x0 |
x(k) |k k0 |
перевести за конечное время в |
|
произвольное |
конечное |
состояние |
xk x(t) |t tk / xk x(k) |k kk |
применив подходящим образом выбранное управляющее воздействие
(возможно, даже неограниченное). |
|
□ |
|||||||
Определение 9.8 (О.9.8). Непрерывный/дискретный динамический |
|||||||||
объект |
с |
матрицами |
( A, C) / |
|
, |
|
|
называется |
полностью |
A |
C |
||||||||
наблюдаемым |
на |
|
интервале |
наблюдения |
|||||
T {t : t0 |
t tk }/T {k : k0 |
k kk }, если |
его состояние x(t) / x(k) |
может быть определено на основе наблюдений за выходом y(t) / y(k) (а возможно, и входом u(t) / u(k) ) в течение интервала наблюдения. □
Сформулируем критерии управляемости и наблюдаемости,
которые инвариантны относительно специфики динамических объектов.
Утверждение 9.1 (У9.1). (Критерий управляемости 1 (КУ1))
Объект с парой матриц ( A, B) / ( A, B ) является полностью управляемым
тогда и только тогда, когда матрица управляемости объекта,
построенная в силу матричного соотношения
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
n 1 |
|
|
(9.48) |
|
|
|
|
|
|
|
|
|||||||||
W |
у |
B |
A |
B |
A |
B |
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
имеет ранг, равный n dim x , т.е. |
|
||||||||||||||
rang Wу dim x n . |
□ (9.49) |
Доказательство сформулированного утверждения проведем на примере дискретного объекта управления с использованием его рекуррентного модельного представления (9.19).
Поставим задачу перевода дискретного ОУ (9.19) из произвольного ненулевого начального состояния x(0) за n dim x –
интервалов дискретности (тактов управления) в желаемое конечное
151
x(n). Вычислим последовательность управляющих воздействий,
образующих «стратегию управления», осуществляющих этот перевод. Для этой цели воспользуемся суммарной моделью дискретного объекта (9.26)
x k Ak x 0 Ak 1Bu 0 Ak 2Bu 1 ABu k 2 Bu k 1
k1
Ak x 0 Ak 1 i Bu i .
i0
Запишем последнее выражение для момента k n в форме, поменяв
при этом порядок суммирования компонентов, |
|
x(n) A n x(0) B AB ... |
A n 1B uT (n 1)uT (n 2)...uT (0) T . (9.50) |
Разрешим полученное выражение относительно вектора «стратегии управления», осуществляющего перевод за n тактов ОУ из начального состояния x(0) в желаемое конечное x(n) , в результате получим
T |
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
uT (n 1)uT (n 2)...uT (0) |
B AB ... |
An 1B |
x(n) An x(0) . |
||||||||||||
Нетрудно видеть, что |
ключевым |
моментом, гарантирующим |
существование искомого управления является обратимость матрицы,
которая имеет место только при выполнении условия (9.48). |
■ |
|||||
Критерий управляемости 2 (КУ2). Пара матриц ( A, B) является |
||||||
полностью управляемой тогда и только тогда, когда матрица |
||||||
Q W W T |
|
|
|
|
|
(9.51) |
у у |
|
|
|
|
|
|
является положительно определенной, т.е. имеет все |
строго |
|||||
положительные |
собственные |
|
значения |
|||
уi 0, уi {Q}, : det у I Q 0; |
|
|
|
|
|
|
i 1, n. |
|
|
||||
Для оценки положительной определенности матрицы |
Q полезно |
|||||
воспользоваться |
критерием |
Сильвестра |
положительной |
определенности матриц. Суть критерия Сильвестра состоит в следующем.
Пусть произвольная квадратная матрица E имеет вид
E11 |
E12 |
E1n |
|
E |
E |
E |
|
E 21 |
21 |
|
2n . |
|
|
|
|
|
En2 |
|
|
En1 |
Enn |
Рассмотрим определители подматриц (главные диагональные миноры) матрицы Е
1 E11; |
2 |
E |
E |
|
|
det 11 |
12 |
|
; |
||
|
|
E21 |
E22 |
|
|