Цифровые системы автоматического регулирования
..pdfСоответствующие матрицы в привычных обозначениях будут иметь вид
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–
нию, если для произвольного начального момента времени 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–