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

МОТС ИТМО

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

11. ДИНАМИЧЕСКИЕ СИСТЕМЫ (ОБЪЕКТЫ) ПРИ ЭКЗОГЕННОМ ВОЗДЕЙСТВИИ

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

Таким образом, рассматриваются проблемы, связанные с изучением реакции динамической системы (ДС) непрерывной и дискретной по времени версий модельного представления, фиксируемой как на ее выходе, так и по вектору состояния в случае подачи на ее вход экзогенного воздействия. Полагается, что исследуемая система линейная (или локально линейная) так, что ее векторно-матричное описание (ВМО) имеет вид

 

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

 

 

 

 

 

 

k 1

 

 

 

 

 

 

 

 

 

 

x(k) F k x(0)

F (k 1 i) Gg(i),

 

 

 

 

 

 

i 0

 

 

 

 

 

 

 

 

 

 

 

k 1

(k 1 i)

 

 

y(k)

 

 

 

k x(0)

 

 

 

 

 

C

F

C

F

Gg(i) .

i 0

(11.5)

(11.6)

Из приведенных выражений (11.3) – (11.6) видно, что движения непрерывной и дискретной систем как по вектору состояния, так и по вектору выхода содержат свободные составляющие, порожденные ненулевым начальным состоянием x(0)

x

 

 

(t) x(t, x(0) 0) eFt x(0),

(11.7)

св

 

 

 

 

 

 

 

 

 

 

 

 

 

y

св

(t) Cx

св

(t) Ce Ft x(0),

(11.8)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(k) x(k, x(0) 0)

 

k x(0),

(11.9)

x

 

 

F

св

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(k)

 

 

 

 

 

k x(0),

 

 

y

св

Cx (k)

С

F

(11.10)

 

 

 

 

 

 

св

 

а также вынужденные составляющие, порожденные внешними

воздействиями g(t) и g(k)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

 

xв (t) x(t, g(t), x(0) 0) eF (t )Gg( )d ,

(11.11)

 

 

 

 

 

 

 

 

0

 

yв

xв

yв

 

 

t

 

 

 

 

 

 

 

 

 

 

 

(t) Cxв (t) Ce F (t )Gg( )d ,

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k 1

 

 

 

 

(k) x(k, g(k), x(0)

0) F

(k 1 i) Gg(i) ,

 

 

 

 

 

 

 

 

 

i 0

 

 

 

 

 

 

k 1

 

 

 

 

 

 

 

 

 

 

(k) Cxв (k)

C F

(k 1 i) Gg(i) .

 

 

 

i 0

 

 

 

 

 

 

 

(11.12)

(11.13)

(11.14)

Нетрудно видеть, что аналитическое вычисление вынужденных составляющих движений в непрерывной и дискретной системах с помощью выражений (11.12) и (11.14) представляет заметные трудности, так как требует необходимости взятия интеграла и формирования суммы. При этом вычисление установившихся

компонентов yу (t) и

yу (k) реакций непрерывной и дискретной ДС,

определяемых в силу соотношений

 

yу (t) lim yв (t), yу (k) lim yв (k) ,

(11.15)

t

k

 

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

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

190

11.1. Динамические непрерывные системы (объекты) при конечномерном экзогенном воздействии. Уравнения Сильвестра

Рассматривается непрерывная линейная система (11.1), в которой экзогенное воздействие g(t) формируется с помощью автономной

конечномерной ДС вида (10.13), имеющей ВМО вида

 

z(t) Ez(t) , z(0) z(t)

 

t 0 ,

g(t) Pg z(t) .

 

(11.16)

 

 

 

 

 

Объединение систем

(11.1)

и

(11.16)

образует агрегированную

систему

с вектором

состояния

~

x

T

, z

T T

задаваемую

x

 

col x, z ,

системой векторно-матричных описаний

 

 

 

 

~

~~ ~

 

T

(0), z

T

 

T

,

 

 

 

(11.17)

x (t) Fx (t), x (0) x

 

 

(0)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

~ ~

~ ~

 

 

~ ~

 

 

~ ~

(11.18)

x(t) Cx x (t), y(t)

Cy x (t), z(t) Cz x (t), g(t) Cg x (t) .

Утверждение 11.1 (У11.1). Матричные компоненты ВМО (11.17) и (11.18) имеют представления

~

F

GPg

~

I n n

~

F

 

, Cx

0n l , C y C 0m l ,

~

0

E

~

 

Pg .

Cz 0l n

Il n , Cg 0m n

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

Сформируем производную

(11.19)

□ (11.20)

~ вектора x (t)

состояния агрегированной системы из компонентов x(t), z(t) , тогда на

основании (11.1) и (11.16) можно построить цепочку матричных равенств

~

x(t)

 

Fx(t) Gg(t)

 

Fx(t) GPg z(t)

 

x (t)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0x(t) Ez(t)

 

 

z(t)

 

 

0x(t) Ez(t)

 

 

 

 

 

F GPg

x(t) ~

 

 

 

 

 

 

 

 

 

 

 

 

~

 

 

 

 

 

 

0

E

 

 

 

Fx (t).

 

 

 

 

 

 

 

z(t)

 

 

 

 

 

 

векторно-

(11.21)

Из сравнения (11.17) и (11.21) следует справедливость первого равенства (11.19). Запишем для переменных x(t), y(t), z(t), g(t)

следующие представления

x(t) Ix(t) 0z(t) I

 

0 x

T

(t), z

T

 

 

T

~ ~

 

 

 

 

 

(t)

Cx x (t) ,

y(t) Cx(t) 0z(t) C

 

0 x

T

(t), z

T

 

 

 

 

 

T

 

~ ~

 

 

 

(t)

Cy x (t) ,

z(t) 0x(t) Iz(t) 0

I x

T

(t), z

T

 

 

 

 

T

 

~ ~

 

 

(t)

 

Cz x (t) ,

g(t) 0x(t) Pg z(t) 0

 

Pg x

T

(t), z

T

 

 

 

T

~ ~

 

 

 

 

 

(t)

Cg x (t) .

Из (11.22) – (11.25) следуют остальные соотношения в

(11.20).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Представления (11.17),

 

(11.18)

 

 

для

 

векторов x(t)

(11.22)

(11.23)

(11.24)

(11.25)

(11.19) и

иy(t)

исследуемой динамической системы (11.1) позволяют записать

~

~

~

~

 

Ft ~

Ft ~

(11.26)

x(t) Cxe

x (0),

y(t) Cy e

x (0) .

191

~

Матрица 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

 

Et

 

E jt

cos(

t)

sin(

t)

 

 

 

 

 

 

 

e

 

 

j

 

j

 

 

 

 

, (11.60)

 

exp(Et) diag e

 

sin( jt)

cos( jt)

; j 1, m

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

для источника одночастотного гармонического воздействия принимает вид (11.60), в котором следует положить j , ( j 1, m) .

В

итоге для векторов

состояния x(t) , выхода y(t) и ошибки

(t) для

разночастотного

и одночастотного гармонического

воздействий соответственно можно записать

x(t) Tz(t) TeEt z(0)

row ( 2j I F 2 ) 1 F

x(t) ( 2 I F 2 ) 1 F

cos(

t)

 

j

 

j I G j diag

sin( j t)

 

 

 

 

I G diag

eEt j

cos( t)

 

 

 

 

 

sin( t)

y(t) Cx(t) CTz(t) CTe Et z(0), (t) (P CT )z(t).

sin(

t)

 

 

 

 

 

j

 

 

 

 

 

cos( j t)

; j 1, m z(0);

 

 

 

 

 

 

 

 

 

sin( t)

 

 

 

; j 1, m ;

cos( t)

 

 

(11.61)

Завершая рассмотрение проблемы конструирования аналитического представления реакции динамической непрерывной системы (11.1) в общем случае МВМВ – типа следует заметить, что в силу (11.46) матрицы T ( ) , CT ( ) и P CT ( ) выполняют функции

передаточных матриц (функций) для случая вещественного гармонического экзогенного воздействия.

Скаляризация векторных процессов в системах МВМВ-типа строится с помощью мажорант и минорант, являющихся для максимальным и минимальным сингулярными числами передаточных матриц T ( ) и CT ( ) .

11.2. Динамические дискретные системы (объекты) при конечномерном экзогенном воздействии. Уравнение Сильвестра

Рассматривается линейная дискретная система (11.2), в которой

экзогенное воздействие g(k) формируется с помощью

автономной

конечномерной ДС вида (10.37), имеющей ВМО вида

 

z(k 1)

 

 

 

 

; g(k)

 

 

 

Ez(k); z(0) z(k)

k

0

Pg z(k)

(11.62)

 

 

 

 

 

 

 

Построение аналитического представление реакции дискретной системы (11.2) опирается на возможности матричного уравнения Сильвестра и строится по той же схеме, что и в случае непрерывных ДС. В связи со сказанным процедуру формирования аналитического представления реакции дискретных динамических систем на конечномерные экзогенные воздействия представим конспективно в виде алгоритма 11.2.

197

Алгоритм 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