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

Цифровые системы автоматического регулирования

..pdf
Скачиваний:
8
Добавлен:
05.02.2023
Размер:
1.38 Mб
Скачать

Соответствующие матрицы в привычных обозначениях будут иметь вид

0,1

0

,

0,25

,

C 0,1

0,5 ,

D 1 .

A

0

 

B

 

 

0,5

 

0,75

 

 

 

 

Если у передаточной функции присутствуют кратные полюсы, например полюс zi имеет кратность p, то в разложении

W(z) на простые дроби будут слагаемые

p 1

Ap j

 

 

,

(z z

) p j

j 0

i

 

 

которые представлены на диаграмме состояний (рис. 4.8).

В этом более общем случае матрица А представляет собой матрицу Жордана, а уравнения состояния имеют жорданову каноническую форму.

Последовательную декомпозицию лучше всего использо-

вать для действительных полюсов и нулей импульсной передаточной функции. В этом случае представим W(z) в виде

W (z)

K (z c1)(z c2 ) ... (z cm )

,

(4.56)

 

 

(z z )(z z

2

) ... (z z

n

)

 

 

 

1

 

 

 

 

 

где ci (i = 1,2,…,m) — нули W(z);

zj

(j = 1,2,…,n) — полюсы

W(z), а для физической реализуемости требуется, чтобы выполнялось условие m n.

Выражение (4.56) представим в виде

W z KW1 z W2 z Wn z ,

где

 

 

 

z c

 

 

 

 

1 c

z 1

 

 

 

W (z)

 

k

 

 

 

k

 

при k m;

(4.57)

 

 

 

 

 

 

 

k

 

z zk

 

1 zk z 1

 

 

 

 

 

 

 

 

W (z)

1

 

 

 

 

 

z 1

при m 1 k n.

(4.58)

 

 

 

 

 

 

 

 

k

 

z zk

 

1

zk z 1

 

 

 

 

 

 

 

 

 

91–

Диаграммы состояний для выражений (4.57) и (4.58) приведены на рис. 4.9,а и б соответственно.

A1

A2

A3

Ap

z 1

zi

z –1

 

 

1

 

 

 

 

zi x1(k)

1 z 1

x2(k)

zi 1

z –1

x3(k)

zi

1

xp (k)

Рис. 4.8. Параллельная декомпозиция при кратных полюсах

Диаграмма состояния всей системы, описываемой передаточной функцией (4.56), будет являться последовательным соединением диаграмм, представленных на рис. 4.9.

92–

1

1

 

z –1

ck

 

 

 

 

 

zk

а

1

 

z–1

1

 

 

 

 

 

zk

б

Рис. 4.9. Последовательная декомпозиция

Пример 4.8. Опять возьмем систему из примера 4.6. Передаточную функцию системы представим в форме (4.56):

 

W (z)

 

 

 

z2 0,2z

 

 

z(z 0,2)

 

W (z)W (z),

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z2 0,6z 0,05

 

(z 0,1)(z 0,5)

 

 

1

2

 

 

 

 

 

 

 

 

 

 

где

W (z)

 

 

z

 

1

 

 

; W (z)

z 0,2

 

1 0,2z 1

.

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

z

0,1

 

1 0,1z 1

2

z

0,5

 

1 0,5z 1

 

 

 

 

 

 

 

 

 

Строим диаграмму состояния в виде последовательно соединенных цепочек:

 

 

1

 

 

 

1

 

r(kT) 1

 

–1

 

1

 

z

–1

x1(kT)

 

 

 

z

x2(kT)

 

 

 

 

 

 

 

 

 

 

 

 

 

–0,2

 

 

 

 

 

 

 

 

 

 

 

 

0,1

 

 

 

0,5

 

93–

По диаграмме состояния записываем уравнения состояния и уравнение выхода:

x1 (k 1)T 0,5x1(kT ) 0,1x2 (kT ) r(kT ); x2 (k 1)T 0,1x2 (kT ) r(kT );

y(kT ) (0,5 0,2)x1(kT ) 0,1x2 (kT ) r(kT ).

Окончательно матрицы в уравнениях системы будут иметь вид

0,5

0,1

,

1

,

C 0,3

0,1 ,

D 1 .

A

0

 

B

 

0,1

 

1

 

 

 

 

4.3.4.Диаграмма состояний цифровых систем

снепрерывной частью

Вобщем случае цифровые системы могут содержать как цифровые, так и аналоговые элементы. Интерфейсом между ними является устройство выборки и хранения.

Прежде чем получать диаграмму состояния всей системы, представим диаграмму состояния экстраполятора нулевого

порядка. Обозначая входной и выходной сигналы экстраполятора через u*(t) и h(t) соответственно, для интервала времени

kT t k 1 T имеем

h t u kT .

Применяя преобразование Лапласа к последнему выражению, получаем следующее соотношение:

H (s) u(kT ) .

(4.59)

s

 

Диаграмма состояния для уравнения (4.59) состоит из единственного ребра с весом s–1, соединяющего вершины H(s)

и u(kT) (рис. 4.10).

94–

Далее получаем диаграмму состояния всей системы, которая будет являться соединением диаграммы состояния цифровой части с диаграммой состояния непрерывной части. Это соединение осуществляется через диаграмму состояния экстраполятора нулевого порядка. Момент времени t0 для задания начальных условий в диаграмме состояния непрерывной части полагаем равным kT, как и при выводе уравнения (4.13) на основе выражения (4.12).

u(kT)

s –1

H(s)

Рис. 4.10. Диаграмма состояний экстраполятора нулевого порядка

Пример 4.9. Рассмотрим получение диаграммы состояния и уравнений состояния для системы, содержащей как непрерывные, так и дискретные элементы (рис. 4.11).

r(t) e(t)

e*(t)

 

u*(t)

 

h(t)

 

y(t)

Цифровая

ЭНП

Непрерывная

 

 

 

 

 

 

часть

 

 

 

часть

 

 

 

 

 

 

 

 

 

Рис. 4.11. Структурная схема системы для примера 4.9

Цифровая часть описывается импульсной передаточной функцией

D(z) U (z)

 

1

,

(4.60)

z2 5z 6

E(z)

 

 

 

а непрерывная часть — передаточной функцией

W (s)

Y (s)

 

5

 

.

(4.61)

H (s)

s

2

 

 

 

 

95–

Диаграмма состояния для выражения (4.60) ранее уже составлялась (см. рис. 4.4), а диаграмма состояния для непрерывной части приведена на рис. 4.12.

x1(t0)

 

 

 

 

 

s –1

 

H(s)

5

 

s –1

 

1

Y(s)

X1(s)

–2

Рис. 4.12. Диаграмма состояний непрерывной части системы к примеру 4.9

Соединяя диаграммы состояния, представленные на рис. 4.4, 4.10 и 4.12, получим диаграмму состояния всей сис-

темы (рис. 4.13).

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x1(kT)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

s 1

 

 

 

 

 

 

 

 

 

 

r(kT) 1

z –1

 

z –1 h(kT) s –1

5

s –1

 

1 Y(s)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

e(kT)

 

 

x3(kT)

x2(kT) H(s)

 

 

 

X1(s)

 

 

5

 

 

 

 

 

 

 

 

2

 

 

 

6

Рис. 4.13. Диаграмма состояний системы к примеру 4.9

На полученной диаграмме состояний учтено, что t0 kT и

x1 kT y kT .

Применяя формулу Мейсона для выходных вершин X1(s), x2(kT), x3(kT), получим

96–

X

 

 

(s)

1 x (kT )

5

x (kT ),

 

 

 

 

s(s 2)

 

 

 

1

 

s 1

 

 

 

2

 

 

 

x

 

 

k 1 T

x

kT ,

 

 

 

(4.62)

2

 

 

 

 

3

 

 

 

 

 

 

x

 

k 1 T x

kT 6x

kT

5x

kT r kT .

3

 

 

 

 

1

 

2

 

3

 

Переходя в первом из уравнений (4.62) к функции времени, запишем

x

(t) x

(kT )

5

1 e 2(t t0 ) x

(kT ).

1

1

 

2

 

2

 

 

Полагая в последнем уравнении

t0 kT , а

t k 1 T и

учитывая оставшиеся уравнения (4.62), получим уравнения состояния системы:

 

(k 1)T

 

 

1

5

1 e 2T

0

 

 

 

 

 

x

 

 

2

x (kT )

0

1

(k 1)T

 

 

 

 

 

 

1

(kT )

0

r(kT ).

x

 

0

 

0

1

x

2

 

 

 

1

 

6

 

 

2

 

 

 

x

(k 1)T

 

 

 

5

x

(kT )

1

 

3

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Уравнение выхода будет иметь вид y kT x1 kT .

4.4. Управляемость и наблюдаемость цифровых систем

4.4.1. Понятия управляемости и наблюдаемости

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

97–

переменные состояния системы и на ее выход. Дуальным является понятие наблюдаемости, под которой понимают возможность восстановления некоторого начального состояния по наблюдениям (измерениям) выхода системы при известном входном сигнале.
Формализуем понятия наблюдаемости и управляемости для цифровой системы, описываемой, например, уравнениями
(4.7), (4.8):
x k j 1 A k j x k j B k j r k j ;
y k j C k j x k j D k j r k j .
Система называется полностью управляемой по состоя-

нию, если для произвольного начального момента времени k0 существует последовательность неограниченных входных воздействий r ki , i 1,2,..., N 1, переводящая каждое на-

чальное состояние x k0 в некоторое конечное состояние x kN за ограниченное время kN k0 .

Система называется полностью управляемой по выходу,

если для произвольного начального момента времени k0 существует последовательность неограниченных входных воздействий r ki , i 1,2,..., N 1, с помощью которой может

быть достигнуто некоторое конечное значение выхода y kN из произвольного начального значения y k0 за ограниченное

время kN k0 .

Система называется абсолютно управляемой (по состоянию или по выходу), если она полностью управляема (по со-

стоянию или по выходу) для всех k0 и всех kN k0 .

Полная управляемость, как видно из определений, требует управляемости для всех начальных условий (но необязательно для всех конечных состояний). Абсолютная же управляемость подразумевает также и все конечные состояния, то есть является более сильным понятием. Если элементы матриц А(k), B(k), C(k), D(k) являются аналитическими функциями дис-

98–

кретного времени k, то понятия полной и абсолютной управляемости совпадают. Тем более они совпадают для стационарных систем.

Система называется полностью наблюдаемой, если для некоторого момента времени k0 состояние x k0 может быть

определено по известным выходным y(k) и входным r(k) сигналам для k0 k kN , где kN — ограниченное время.

Система называется глобально (абсолютно) наблюдаемой,

если она полностью наблюдаема для всех k0 и всех kN k0 .

Для стационарных систем понятия полной и глобальной наблюдаемости совпадают.

4.4.2. Определение управляемости по уравнениям динамики

Анализ уравнений состояния позволяет ответить на вопрос об управляемости системы, о чем утверждает следующая теорема.

Теорема 4.1. Для того чтобы система, описываемая уравнениями состояния (4.7), была полностью управляемой по состоянию, необходимо и достаточно, чтобы ранг матрицы управляемости Q был равен n. Матрица управляемости Q в блочной записи имеет вид

Q Q0 Q1 QN 1 , (n N p), (4.63)

где Qi Φ kN ,ki 1 B ki , n p , i 0,1,2,..., N 1.

Иногда установление ранга матрицы довольно утомительное занятие, поэтому равенство ранга матрицы Q величине n можно определить по невырожденности соответствующей матрицы Грама:

N 1

 

 

G Q QT Qi QTi

 

 

i 1

 

 

N 1

 

 

Φ(kN , ki 1)B(ki )BT (ki )ΦT (kN ,ki 1).

(4.64)

i 1

99–

Теорема 4.2. Для полной управляемости по выходу системы, описываемой уравнением (4.7), необходимо и достаточно, чтобы матрица Т управляемости по выходу имела ранг, равный m. Матрица Т задается выражением

 

T T0

T1

TN , (m (N 1) p),

(4.65)

 

C(k

N

)Φ(k

N

,k

)B(k ), 0 i N 1,

 

где T

 

 

i 1

i

 

 

 

),

 

 

i N.

 

i

D(k

N

 

 

 

 

 

 

 

 

 

 

Управляемость по выходу можно установить и по невырожденности соответствующей матрицы Грама:

N 1

TTT TiTiT .

i 0

Для стационарных систем определение управляемости по состоянию и по выходу облегчается. Кроме того, для стационарного случая существует ряд теорем, которые часто более удобны в практическом применении.

Теорема 4.3. Линейная стационарная система

x(ki 1) Ax(ki ) Br(ki );

(4.66)

y(ki ) Cx(ki ) Dr(ki ),

 

полностью управляема по состоянию, если

и только если

(n N p)-мерная матрица

 

Q B AB A2B AN 1B

 

 

имеет ранг n или матрица Грама QQT является невырожденной.

Теорема 4.4. Линейная стационарная цифровая система (4.66) полностью управляема по состоянию, если и только если строки матрицы

zE A 1 B

линейно независимы.

100–