
МОТС ИТМО
.pdf11. ДИНАМИЧЕСКИЕ СИСТЕМЫ (ОБЪЕКТЫ) ПРИ ЭКЗОГЕННОМ ВОЗДЕЙСТВИИ
Прежде чем приступать к изложению материала, вынесенного в заголовок раздела объясним, почему в нем наряду с понятием объект появилось понятие система. Напомним, что системой называется функциональное объединение объекта управления и устройства управления, которое наделяет образованную таким образом систему необходимыми динамическими показателями качества процессов в переходном и установившемся режимах движения. Исходный объект управления, взятый в отдельности, может обладать плохими динамическими показателями и даже быть неустойчивым, поэтому ставить задачу исследования его показателей при экзогенном (внешнем) воздействии возможное, но бессмысленное занятие. Поэтому приводимый ниже материал сориентирован в основном на исследование динамических систем.
Таким образом, рассматриваются проблемы, связанные с изучением реакции динамической системы (ДС) непрерывной и дискретной по времени версий модельного представления, фиксируемой как на ее выходе, так и по вектору состояния в случае подачи на ее вход экзогенного воздействия. Полагается, что исследуемая система линейная (или локально линейная) так, что ее векторно-матричное описание (ВМО) имеет вид
|
x(t) Fx(t) Gg(t), x(0) x(t) |t 0 , y(t) Cx(t) |
(11.1) |
|||||||||||||||||||||
для непрерывного случая и |
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
x(k 1) Fx(k) Gg(k), x(0) x(k) |k 0 , y(k) Cx(k) |
(11.2) |
|||||||||||||||||||||
для дискретного случая. В (11.1) и (11.2) |
x – вектор состояния ДС, |
||||||||||||||||||||||
x Rn ; g |
– вектор экзогенного (внешнего) |
воздействия, y |
– вектор |
||||||||||||||||||||
|
|
|
|
|
|
|
|
|
Rn n ;G, |
|
|
||||||||||||
выхода; |
g, y Rm ; F, |
|
F |
– матрицы |
состояния; |
F, |
F |
G |
– |
||||||||||||||
матрицы |
входа; G, |
G |
Rn m ; C, |
C |
|
– |
матрицы |
выхода; C, |
C |
Rm n . |
|||||||||||||
Непрерывное время t и дискретное k |
связаны соотношением t t k , |
где t - интервал дискретности.
В разделе 9 показано, что решение системы (11.1) относительно
векторов состояния |
x(t) и выхода y(t) |
может быть представлено в |
форме |
|
|
|
t |
|
x(t) eFt x(0) eF (t )Gg( )d , |
(11.3) |
|
|
0 |
|
|
t |
|
y(t) Ce Ft x(0) |
Ce F (t )Gg( )d . |
(11.4) |
|
0 |
|
В свою очередь решение системы (11.2) относительно векторов состояния x(k) и выхода y(k) принимает вид
189
~
Матрица F представляет собой блочную матрицу (11.19), матричные компоненты которой можно связать матричным уравнением Сильвестра. Для этой цели докажем утверждение.
Утверждение 11.2 (У11.2). Пусть выбором начальных состояний x(0) и z(0) обеспечивается выполнение векторного равенства
x(t) Tz(t) t , |
|
|
|
|
|
|
|
|
(11.27) |
||||||||
тогда T (n l) |
|
– матрица подобия (в общем случае особого) |
|||||||||||||||
удовлетворяет матричному уравнению Сильвестра |
|
||||||||||||||||
TE FT GPg . |
|
|
|
|
|
|
|
|
□(11.28) |
||||||||
Доказательство утверждения строится на дифференцировании |
|||||||||||||||||
(11.27) по времени, которое дает равенство |
|
|
|||||||||||||||
x(t) Tz(t) t , |
|
|
|
|
|
|
|
|
(11.29) |
||||||||
и на последующей подстановке (11.27) и |
(11.29) в соотношения |
||||||||||||||||
(11.1) |
для |
x(t) и (11.16) для |
z(t) , |
в результате чего возникают две |
|||||||||||||
цепочки векторно-матричных соотношений |
|
|
|||||||||||||||
x(t) Fx(t) Gg(t) FTz(t) GPg z(t) (FT GPg )z(t), |
(11.30) |
||||||||||||||||
x(t) Tz(t) TEz(t) . |
|
|
|
|
|
|
(11.31) |
||||||||||
Векторно-матричные соотношения (11.30) и (11.31) делают |
|||||||||||||||||
справедливыми матричные равенства |
|
|
|
||||||||||||||
TE FT GPg TE FT GPg . |
|
|
■(11.32) |
||||||||||||||
Матричное уравнение Сильвестра (11.28) позволяет представить |
|||||||||||||||||
матрицу |
~ |
(11.19) в форме |
|
|
|
|
|
|
|||||||||
F |
|
|
|
|
|
|
|||||||||||
~ |
|
F GPg |
F TE FT |
|
|
|
(11.33) |
||||||||||
F |
|
|
|
|
|
|
|
|
|
|
. |
|
|
|
|||
|
|
|
|
0 |
|
|
E |
|
|
0 |
|
E |
|
~ |
|
|
|
Представление |
|
матрицы |
вида |
(12.33) |
позволяет |
||||||||||||
|
F |
||||||||||||||||
сформулировать и доказать следующее утверждение. |
|
||||||||||||||||
Утверждение 11.3 (У11.3). Показательная (степенная) матричная |
|||||||||||||||||
функция |
|
|
~ |
~k |
|
от |
матрицы |
~ |
вида (11.33), где |
k – целое |
|||||||
f (F ) F |
|
F |
|||||||||||||||
положительное число представима в форме |
|
|
|||||||||||||||
~ k |
|
F k |
TE k |
F k T |
|
|
|
|
|
|
|||||||
F |
|
|
|
|
|
0 |
|
|
E k |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
. |
|
|
|
|
(11.34) |
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
Доказательство |
утверждения |
строится |
на непосредственном |
||||||||||||||
вычислении целых степеней матрицы |
~ |
|
|
||||||||||||||
F путем их перемножения. В |
|||||||||||||||||
результате получим |
|
|
|
|
|
|
|
|
|
||||||||
~1 |
|
~ |
|
F |
TE FT |
|
|
|
|
|
|
||||||
F |
|
|
F |
|
|
|
E |
; |
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
~ |
2 |
|
|
~~ |
|
F |
TE FT F |
TE FT |
|
|
|
||||||
F |
|
|
FF |
|
|
E |
|
|
E |
|
|
|
|||||
|
|
|
|
|
|
|
0 |
|
0 |
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
192 |
|
|
|
|
F 2 |
|
FTE F 2T TE2 FTE |
F 2 |
TE2 F |
2T |
|
|
|
|||||||||||
= |
|
0 |
|
|
|
E |
2 |
|
|
|
|
E |
2 |
|
|
; |
|
|
||
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|||||
~3 |
|
~2 |
~ |
F 2 |
TE2 F 2T F TE FT |
|
|
|
|
|
|
|
|
|||||||
F |
|
F F |
|
|
|
E |
2 |
|
E |
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
0 |
|
|
|
0 |
|
|
|
|
|
|
|
|
|||
|
F 3 |
|
F 2TE F 3T TE3 F 2TE |
F 3 |
TE3 F 3T |
|
|
|
|
|||||||||||
|
|
0 |
|
|
|
E |
3 |
|
|
E |
3 |
|
; |
|
|
|
||||
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
||||
База индукции сформирована, поэтому для искомой матричной |
||||||||||||||||||||
функции |
|
~ |
~k |
можно записать |
|
|
|
|
|
|
|
|
|
|||||||
f (F ) F |
|
|
|
|
|
|
|
|
|
|||||||||||
~ ~ ~ |
F k 1 |
TE k 1 F k 1T F TE FT |
|
|
|
|
|
|
||||||||||||
F k F k 1F |
|
0 |
|
|
E k 1 |
|
0 |
E |
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
■(11.35) |
||||||
F k |
F k 1TE F kT TE k F k 1TE |
F k |
TE k F kT |
|
||||||||||||||||
|
||||||||||||||||||||
|
|
|
|
|
|
|
|
E k |
|
|
|
|
|
E k |
|
. |
|
|||
|
0 |
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
Примечание 11.1 (ПР11.1). Нетрудно видеть, что положения
утверждения 11.3 позволяют сделать следующее обобщение. Если
~
матричная функция от матрицы f (F ) представляет собой степенной
конечный или бесконечный ряд по целым положительным степеням
матрицы |
~ |
F GPg |
|
F TE FT |
|
|
|
||||||||
F |
|
|
|
|
|
|
|
|
|
, то оказывается справедливым |
|||||
|
|
0 |
|
E |
|
0 |
|
E |
|
|
|
|
|||
представление |
|
|
|
|
|
|
|
|
|
|
|
|
|
||
~ |
f (F ) |
Tf (E) f (F )T |
|
|
|
(11.36) |
|||||||||
f (F ) |
0 |
|
|
|
f (E) |
. |
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
||||||
Представление (11.36) матричной степенной функции от матрицы |
|||||||||||||||
делает справедливыми положения следующего утверждения. |
~ |
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Утверждение |
11.4 |
(У11.4). Матричная экспонента |
eFt , где |
||||||||||||
матрица |
|
~ |
|
F GPg |
|
|
F TE FT |
имеет вид |
|
||||||
|
F |
|
0 |
|
E |
|
|
0 |
E |
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
||||
~ |
eFt |
|
TeEt eFtT |
|
|
|
|
|
|
□■(11.37) |
|||||
eFt |
|
|
|
eEt |
|
|
. |
|
|
|
|
||||
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Если в (11.22) – (11.24) подставить (11.26), в котором |
|||||||||||||||
использовать |
(11.37), |
|
то |
|
получим |
для |
|
переменных x(t), y(t), z(t) |
|||||||
представление |
|
|
|
|
|
|
|
|
|
|
|
|
|
||
x(t) eFt x(0) (TeEt eFtT )z(0) TeEt z(0) eFt (x(0) Tz(0)); |
(11.38) |
||||||||||||||
y(t) CeFt x(0) C(TeEt |
eFtT )z(0) CTeEt |
CeFt (x(0) Tz(0)); |
(11.39) |
||||||||||||
z(t) eEt z(0). |
|
|
|
|
|
|
|
|
|
|
|
|
|
(11.40) |
193
Анализ выражений (11.38) – (11.40) позволяет получить представления всех компонентов движения системы (11.1) при конечномерном экзогенном воздействии на ее входе.
Вынужденные составляющие движения по состоянию, выходу и ошибке (t) g(t) y(t)
x (t) x(t, g(t), x(0) 0) x(t, z(0) 0, x(0) 0) |
TeEt z(0) eFtTz(0); |
(11.41) |
||||
в |
|
|
|
|
|
|
y |
в |
(t) Cx |
в |
(t) CTe Et z(0) Ce FtTz(0); |
|
(11.42) |
|
|
|
|
|||
|
|
(t) (P |
|
CT )eEt z(0) Ce FtTz(0). |
|
|
в |
|
|
|
|||
|
g |
|
|
|
Установившиеся составляющие движения по состоянию, выходу и ошибке
xу |
(t) lim xв (t) TeEt z(0); |
|
|
|
|
|
|
|
|
|
(11.43) |
|||||||
|
|
|
t |
|
|
|
|
|
|
|
|
|
|
|
|
|
||
y |
у |
(t) Cx |
у |
(t) CTeEt z(0); |
у |
(t) g(t) y |
у |
(t) (P CT )eEt z(0). (11.44) |
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
g |
|
|
||||
|
|
Переходные составляющие движения по состоянию, выходу и |
||||||||||||||||
ошибке |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
x |
пер |
(t) x |
у |
(t) x (t) eFtTz(0), |
y |
пер |
(t) Cx |
пер |
(t) CeFtTz(0) |
пер |
(t). (11.45) |
|||||
|
|
|
|
|
|
в |
|
|
|
|
|
|
|
Из сравнения представлений (11.40) и (11.43) становится понятным математическое «содержание» матрицы T , состоящее в том, что она представляет собой матрицу подобия между процессами по вектору состояния источника конечномерного экзогенного воздействия (ИКЭВ) и установившейся составляющей вектора состояния исследуемой динамической системы так, что устанавливаются соотношение
x |
у |
(t) TeEt z(0) Tz(t). |
(11.46) |
|
|
|
|
Если обратиться к общему решению (11.38) для |
x(t) , то из него |
можно видеть, что соотношение (11.46) устанавливаются с начального момента t 0, если начальные соотношения подчинить условию x(0) Tz(0).
Таким образом, для вычисления реакции непрерывной линейной системы (11.1) на произвольное конечномерное экзогенное воздействие можно воспользоваться приводимым ниже алгоритмом.
Алгоритм 11.1 (А11.1)
1.Построить векторно-матричное описание непрерывной системы (11.1) с матрицами (F,G,C);
2.Построить векторно-матричное описание источника экзогенного конечномерного воздействия (11.16) с матрицами (Е,Pg);
3.Решить аналитически или численно матричное уравнение Сильвестра (11.28) относительно матрицы подобия T ;
194
4. С помощью соотношений (11.41) – (11.45) решить задачу анализа реакции непрерывной системы на конечномерное экзогенное воздействие.
Остановимся теперь на ситуации, когда ИКЭВ формирует
типовые экзогенные воздействия, ограничившись случаями ступенчатого и гармонического экзогенных воздействий. Первое типовое воздействие используется для оценки динамических свойств
исследуемой ДС по кривой переходного процесса, второе – по амплитудным частотным характеристикам.
Случай 11.1 (Сл11.1) ступенчатого экзогенного воздействия
характеризуется матрицами E и Pg модельного представления ИКЭВ (11.16) в силу (11.19), имеющими вид
E 0m m , |
Pg Im m . |
|
|
|
|
(11.47) |
||
Подстановка матриц вида(11.47) в уравнение Сильвестра (11.28) |
||||||||
дает |
|
|
|
|
|
|
|
|
FT G , что позволяет для матрицы T записать |
|
|||||||
|
T F 1G. |
|
|
|
|
|
(11.48) |
|
Матрицы E вида (11.47) и T |
вида (11.48) дают для матричной |
|||||||
|
~ |
|
|
|
|
|
|
|
экспоненты eFt , записанной в форме (11.37) представление |
|
|||||||
~ |
eFt |
TeEt eFtT |
|
eFt |
eFt I |
F 1G |
(11.49) |
|
eFt |
|
|
|
|
|
. |
||
|
eEt |
|
|
|
|
|||
|
0 |
|
|
0 |
I |
|
|
Подстановка матричной экспоненты вида (11.49) и выделение из
полного движения вынужденной составляющей по вектору состояния и выхода исследуемой системы дает
|
|
x (t) eFt |
I F 1Gz(0) H |
x |
(t)z(0); z(0) argmin |
|
|
|
z(0) |
|
|
|
|
1 , |
(11.50) |
||||||||||||
|
|
|
|
|
|
||||||||||||||||||||||
|
|
в |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
y |
в |
(t) Cx |
в |
(t) C eFt |
I F 1Gz(0) CH |
x |
(t)z(0) H |
y |
(t)z(0). |
(11.51) |
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
В выражениях (11.50) и (11.51) матрицы |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
H |
x |
(t) eFt I F 1G и H |
y |
(t) C eFt I F 1G CH |
x |
(t), |
|
|
|
|
|
|
|
(11.52) |
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
представляют |
собой |
матричные |
|
переходные |
функции |
соответственно по состоянию и выходу. Если возникает необходимость скаляризации задачи, состоящей в исследовании динамических свойствi, j го сепаратного канала ДС, связывающего i й выход yi (t) с j м
входом |
g j (t) , то в этом случае в переходной матрице по выходу |
||||||||
H y (t) H (t) (или просто в |
|
переходной матрице) выделить |
i, j й |
||||||
элемент |
(t) Ci eFt I F 1G |
|
i, j |
|
, |
|
|||
|
|
|
|
|
|||||
H |
ij |
j |
1,m |
(11.53) |
|||||
|
|
|
|
|
|
|
|
||
где |
Ci ,G |
j |
соответственно i я строка матрицы выхода C |
и j й |
|||||
|
|
|
|
|
|
|
|
|
столбец матрицы входа G.
Скаляризация векторных процессов в системах типа «многомерный вход – многомерный выход» (МВМВ) строится с
195

помощью мажорант и минорант, являющихся для t максимальным и минимальным сингулярными числами переходной матрицы H (t).
Случай 11.2 (Сл11.2) гармонического экзогенного воздействия характеризуется матрицами E и Pg модельного представления ИКЭВ
(11.16) в силу (10.26), имеющими вид
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
E diag |
|
|
|
|
|
, P |
diag P |
1 |
0 ; j 1, m (11.54) |
||||||||
E |
j |
|
|
|
j |
|
; j 1, m |
||||||||||
|
|
|
j |
0 |
|
|
|
g |
gj |
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Прежде, чем формировать аналитическое представление реакции динамической системы на гармоническое воздействие следует заметить, что частотные методы применяются для исследования установившихся составляющих движения системы, которые для векторов состояния и выхода имеют соответственно вид (11.43) и (11.44), в которых требуется вычислить матрицу T из решения
матричного уравнения Сильвестра(11.28) и матричную экспоненту eEt . Для решения уравнения Сильвестра (11.28) с матричными
компонентами вида (11.54) запишем его в столбцовой форме |
|
||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
TE |
FT |
GPg ; 1,2m. |
|
|
|
|
|
|
|
|
|
|
|
|
(11.55) |
||||||||
Выделим |
в |
(11.55) |
по два соседних |
столбца |
с индексами |
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
||||||||||||||
2 j 1, 2 j ( j 1,m) |
в составе каждой из правых матриц E,T , Pg . |
||||||||||||||||||||||
Тогда с |
учетом |
структуры (11.54) матриц E и Pg |
для |
смежных |
|||||||||||||||||||
столбцов T2 j 1,T2 j матрицы T получим два уравнения Сильвестра |
|||||||||||||||||||||||
T2 j j |
FT2 j 1 G j , T2 j 1 j |
FT2 j 0. |
|
|
|
|
|
|
|
|
|
(11.56) |
|||||||||||
Решение системы (11.56) относительно фрагмента матрицы T , |
|||||||||||||||||||||||
составленного из двух смежных столбцов принимает вид |
|
|
|
||||||||||||||||||||
T |
|
( 2I F 2 ) 1 F |
|
I G |
|
|
|
|
|
|
|
|
|||||||||||
T |
|
; j 1,m . |
|
|
|
|
(11.57) |
||||||||||||||||
2 j 1 |
2 j |
|
|
|
j |
|
|
j |
j |
|
|
|
|
|
|
|
|
|
|
|
|||
Матрица T , как решение матричного уравнения Сильвестра |
|||||||||||||||||||||||
(12.28) в итоге принимает вид |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
T row T |
|
|
; j |
|
row ( 2I F 2 ) 1 F |
|
|
I G |
; j |
|
. |
(11.58) |
|||||||||||
|
T |
1, m |
j |
1, m |
|||||||||||||||||||
2 j 1 |
|
2 j |
|
|
|
|
|
|
|
j |
|
|
|
|
|
j |
|
|
|
|
|
Отметим, что результат (11.58) получен для разночастотного случая гармонических экзогенных воздействий, подаваемых на входы динамической системы МВМВ - типа, если же режим воздействия
одночастотный так, что j ,( j 1, m) , в этом случае матрица T
принимает вид
T ( 2 I F 2 ) 1 F I G. (11.59)
Матричная экспонента eEt для источника разночастотного гармонического воздействия в силу (10.28) записывается в форме
196


Алгоритм 11.2(А11.2)
1. Формирование векторно–матричного описания (11.2) исследуемой дискретной системы с матрицами F, G, C ;
2. Формирование векторно-матричного описания (10.62) источника дискретного конечномерного экзогенного воздействия с
матрицами E, Pg ;
3.Формирование агрегированной автономной дискретной
системы |
|
~ |
|
|
|
|
|
|
|
|
|
|
|
|
|
~ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
~ |
|
|
|
|
~ |
|
|
|
|
|
||||||
~ |
|
|
|
~ |
|
|
|
|
|
|
|
|
|
|
~ |
|
|
|
|
|
|
|
|
|
|
|
~ |
|
|
~ |
k ; |
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||
x k 1 Fx k ; x(k) Cx x k ; y(k) Cy x k ; z(k) Cz x |
(11.63) |
||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
~ |
~ |
|
|
|
|
|
|
|
|
|
|
|
~ |
~ |
k , |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
g(k) Cg x k ; (k) |
C x |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||
где |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
x (k) x |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
T |
|
|
~ |
|
|
|
|
F |
|
GP |
|
|
~ |
|
|
~ |
|
|
|
|
|
||||||||||||||
|
(k) z |
|
|
|
|
|
|
|
|
;Cx I |
0 ; C y C |
0 ; |
|
||||||||||||||||||||||||||||||||||
T |
T |
k |
; F |
|
|
|
|
|
|
|
g |
|
|
||||||||||||||||||||||||||||||||||
~ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
~ |
|
|
|
|
|
|
~ |
|
|
|
|
|
|
|
|
|
|
|
|
|
~ |
|
0 |
|
|
E |
|
|
|
|
|
|
|
|
|
(11.64) |
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
Cz 0 |
I ;Cg |
0 |
Pg ;C |
C P g |
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||||
|
|
позволяющей |
для |
|
агрегированного вектора |
состояния |
~ |
||||||||||||||||||||||||||||||||||||||||
|
|
|
x k |
||||||||||||||||||||||||||||||||||||||||||||
записать |
|
|
|
~ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
~ |
|
|
|
|
k |
~ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(11.65) |
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
x k F |
x 0 ; |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
4. Введение априорной «гипотезы» о соблюдении векторного |
|||||||||||||||||||||||||||||||||||||||||||||
подобия |
|
x k |
|
|
z k , k 0, |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(11.66) |
|||||||||||||||||||||
|
|
|
|
|
T |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
где T – n l - матрица подобия;
5. На основе (11.2), (11.62) и (11.66) установление того, что
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
матрица T |
удовлетворяет |
|
матричному |
|
уравнению |
Сильвестра вида |
|||||||||||||||||||||||||||
(11.28), записываемого для дискретного случая в форме |
|||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
T E F T G Pg |
; |
|
|
|
|
|
|
|
|
(11.67) |
||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
~ |
|
|
|
|
|
|
|
|
|
6. Представление матрицы |
F вида (11.64) на основе уравнения |
||||||||||||||||||||||||||||||||
Сильвестра (11.67) в форме |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
~ |
|
|
|
|
|
|
|
|
|
|
|
|
|
T E FT |
|
|
|
|||||||||||||||
|
F |
|
F |
G |
Pg |
F |
|
; |
(11.68) |
||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||
|
|
|
|
|
0 |
|
|
E |
0 |
|
|
|
|
|
|
E |
|
|
|
7. С использованием положений утверждения 11.3 представление
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
k |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
~ |
|
|
показательной (степенной) функции F |
в блочной матричной форме |
||||||||||||||||||
|
~ |
|
|
|
k |
|
|
|
|
k |
|
k |
|
|
|
|
|||
|
k |
F |
T |
E |
F |
T |
|
(11.69) |
|||||||||||
|
|
|
|||||||||||||||||
F |
|
|
|
|
|
|
|
k |
|
|
|
; |
|
||||||
|
|
|
0 |
|
|
|
|
|
E |
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
8. Формирование аналитических представлений |
векторных |
переменных исследуемой дискретной системы (11.2) и источника дискретного конечномерного экзогенного воздействия (11.62) с
198