Мехатроника и динамика мини-роботов
.pdfУчитывая соотношение
|
H T z 1 |
n |
I z 1H T 1, |
|
|
||||
n 0 |
|
|
|
|
выражение z-преобразования для вектора состояния дискретной системы запишем в виде
V z I z 1H T 1 |
0 . |
(2.34) |
|
|
|
|
|
Применяя обратное z-преобразование к уравнению (2.34), для вектора состояния получим
nT Z 1 I z 1H T 1 |
0 , |
(2.35) |
|
|
|
|
|
что является общим решением уравнения состояния, из которого могут быть определены переменные состояния системы в последовательные моменты квантования, если только известны расширенная матрица перехода и матрица В и заданы начальные условия. Уравнение (2.35) позволяет систематизировать и унифицировать анализ динамического поведения дискретных систем.
Анализ устойчивости системы
Пусть матрица H T разделена на блоки согласно размерностям векторов m и x:
T |
0 |
|
, |
(2.36) |
|
H T |
T |
|
|
||
|
T |
|
|
||
где T – квадратная матрица.
70
Тогда
1 |
0 |
z 1 |
T |
0 |
|
|
|
I z 1H T |
|
|
T |
|
|
||
0 |
1 |
|
|
T |
|
||
(2.37)
I z 1 T |
||
|
z 1 T |
|
|
||
|
|
|
Так как по определению |
|
|
V z |
M z |
, |
|
|
|
|
X z |
|
.
I z 1 T
0 m 0 |
|
, |
|
|
|
x 0 |
|
|
то уравнение (2.34) можно записать как |
|
|
|
|
|
|
|
|
||||||||||
M z |
|
|
|
|
|
I z |
1 |
T |
|
1 |
|
|
|
|
||||
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
I z 1 T |
|
1 T |
|
|
z 1 T |
|
I z 1 T |
|
|
|
||||
X z |
z 1 |
|
|
I |
1 |
1 |
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(2.38) |
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
m 0 |
. |
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
x 0 |
|
|
|
|
|
|
|
|
|
|||
При выводе этого уравнения было использовано равенство (2.37). Отсюда находим z-преобразования для векторов состояния входа и процесса соответственно:
M z |
|
I z 1 T |
|
1 m 0 |
(2.39) |
|
|
и
X z z 1 I z 1 T 1 T M z I z 1 T 1 x 0 . (2.40)
71
Первое слагаемое в правой части z-преобразования (2.39) присутствует благодаря вектору M, а второе слагаемое – благодаря ненулевым начальным условиям. Второе слагаемое равно нулю, если в начальный момент времени система находится в состоянии покоя. Последовательность дискретных значений переменных состояния и переменных состояния процесса можно получить по формулам (2.39) и (2.40), применяя к ним обратное z-преобразование.
Если система в начальный момент времени находится в состоянии покоя, то z-преобразование для вектора состояния процесса принимает вид
X z z 1 I z 1 T 1 T M z .
Характеристическое уравнение системы имеет вид |
|
det I z 1 T 0. |
(2.41) |
Для устойчивости системы необходимо, чтобы корни уравнения (2.40) находились в пределах единичной окружности наплоскости z.
Учитывая равенство (2.36), выражение (2.30) для вектора состояния t можно записать в виде
m t |
t nT |
0 |
m nT |
|
||
|
|
|
|
|
. |
|
x t |
t nT |
t nT x nT |
|
|||
Следовательно, в интервале |
|
|
|
|
||
|
|
nT t n 1 T |
|
|
|
|
вектор состояния процесса определяется выражением |
|
|||||
x t t nT m nT t nT x nT , |
(2.42) |
|||||
где значения векторов m nT и x nT могут быть легко найдены при использовании выражений (2.39) и (2.40). Полагая n 0,1, 2,
72
в уравнении (2.42), находим функции времени для переменных состояния процесса.
Пр и м е р 2.2
Сцелью иллюстрации описанного выше метода z-преобразования для анализа дискретных систем рассмотрим разомкнутую систему.
Схема простой импульсной системы показана на рис. 2.12. Предполагается, что входное воздействие имеет вид единичной ступенчатой функции и период прерывания равен T, с. Требуется определить реакцию на выходе системы.
Рис. 2.12. Простая импульсная система
Из рассмотрения схемы системы в переменных состояния, показанной на рис. 2.13, следует, что дифференциальные уравнения состояния системы имеют вид
m1 0; m2 0; x1 m2 x1.
Рис. 2.13. Схема в переменных состояния для простой импульсной системы
Уравнения переходных состояний системы находим в виде m1 nT m1 nT ;
m2 nT m1 nT ; x1 nT x1 nT .
73
Матрицы А и В записываются следующим образом:
0 |
0 |
0 |
|
|
|
|
0 |
|
; |
A 0 |
0 |
|||
|
|
1 |
|
|
0 |
1 |
|
||
|
1 |
0 |
0 |
|
B |
|
0 |
|
|
1 |
0 |
|
||
|
|
0 |
|
|
|
0 |
1 |
|
|
Матрицу перехода находим в виде
1 |
|
0 |
|
|
|
0 |
|
|
|
|
||
|
|
1 |
|
|
|
0 |
|
|
|
|
||
Ф t 0 |
|
|
|
|
. |
|
|
|||||
|
1 e |
t |
|
e |
t |
|
|
|
||||
0 |
|
|
|
|
|
|
|
|||||
Отсюда находим матрицу H T : |
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
1 |
|
|
|
|
0 |
|
0 |
|
H T Ф T B |
|
|
1 |
|
|
|
|
0 0 |
|
|||
|
|
|
|
|
|
|
||||||
|
|
|
e |
T |
|
|
0 |
e |
T |
|||
и |
|
1 |
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
T |
1 |
|
0 |
|
|
|
|
|
|
|||
|
|
|
; |
|
|
|
|
|
|
|||
|
|
1 |
|
0 |
|
|
|
|
|
|
||
t e T ; |
|
|
|
|
|
|
|
|||||
t 1 |
e T |
0 |
; |
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
1 z 1 T |
|
|
|
|
1 |
|
|
|
; |
|
|
|
1 z |
|
|
|
0 |
|
|
||||||
|
|
|
z 1 |
|
|
1 |
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
1 z 1 T 1 z 1e T .
74
Так как
X z z 1 1 z 1 T 1 T 1 z 1 T 1 m 0 ,
то окончательно имеем
z 1 e T
X z z 1 z e T .
Отметим, что решение этой простой задачи можно получить классическим методом. В данном случае применение обобщенного подхода не дает каких-либо преимуществ. Заметим также, что данный пример был приведен с целью проиллюстрировать применение рассмотренного внастоящемпараграфеметода анализа дискретных систем.
П р и м е р 2.3
Замкнутая система
Требуется найти выход импульсной системы, структурная схема которой приведена на рис. 2.14. Входное воздействие имеет вид единичной ступенчатой функции.
Рис. 2.14. Простая импульсная система с обратной связью
Схема системы в переменных состояния изображена на рис. 2.15. Из рассмотрения этой схемы находим, что дифференциальные уравнения состояния системы имеют вид
m1 0; m2 0; x1 m2 x1.
75
Рис. 2.15. Схема в переменных состояния для импульсной системы с обратной связью
Уравнения переходных состояний находим в виде
m1 nT m1 nT ;
m2 nT m1 nT ;
x1 nT x1 nT .
Отсюда находим матрицу А и В:
0 |
0 |
0 |
|
|
|
0 |
0 |
|
; |
A 0 |
|
|||
|
1 |
|
|
|
0 |
1 |
|
||
1 |
0 |
0 |
|
|
|
0 |
|
|
|
B 1 |
1 . |
|||
|
0 |
1 |
|
|
0 |
|
|
||
В качестве иллюстрации рассмотрим определение расширенной матрицы перехода другим методом. Обращаясь к схеме системы в переменных состояния, запишем
76
m1 1 m1 0 ;
m2 m1 x1 m2 0 ; x1 m2 0 x1 .
Решение последнего дифференциального уравнения имеет вид
x1 1 e m2 0 e x1 0 .
Отсюда находим расширенную матрицу перехода:
1 |
0 |
|
|
0 |
|
|
|
|
1 |
|
|
0 |
|
|
|
Ф 0 |
|
|
. |
|
|||
|
1 e |
|
e |
|
|
||
0 |
|
|
|
|
|||
Матрицу H T получаем в виде |
|
|
|
|
|
|
|
|
1 |
|
|
0 |
|
0 |
|
|
1 |
|
|
0 |
|
1 |
|
H T Ф T B |
|
|
|
. |
|||
|
e |
T |
0 |
1 2e |
T |
||
1 |
|
|
|
||||
Матрицу H T разобьем на блоки таким образом, что
T 1; |
|
|
|
0 |
1 |
|
; |
T |
|
|
|
0 |
1 2e T |
|
|
|
|
|
|
|
1 |
T |
. |
1 |
e T |
|
|
77
Следовательно, |
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
1 |
|
1 |
|
|
|
1 |
|
|
1 |
z 1 |
2e T 1 |
1 z 1 |
|
|
1 |
z |
T |
|
|
|
|
|
|
|
|
|
|
|
|
; |
|
|
|
1 |
2e |
T |
|
|
|
|
||||||||
|
|
|
|
|
1 z |
|
|
|
|
|
0 |
1 |
|
|
||
|
|
|
|
|
|
|
1 |
|
|
|
||||||
1 z 1 T 1 1 1z 1 .
Если x1 0 , то из формул (2.37) и (2.38) следует, что
|
|
|
z |
|
z |
|
|
|
|
|
z |
1 |
|
|
|
|
T |
|
|
|
|
1 |
|
T |
|
|
|||||||
M2 |
|
|
|
|
|
1 |
|
|
2e |
|
|
1 1 z |
|
1 e |
|
|
, |
||||||||||||||||
|
X |
1 |
z |
z 1 z 2e |
T |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 e |
T |
|
|
|
|
||||
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
откуда |
|
|
|
|
|
|
|
|
|
z 1 e T |
|
|
|
|
|
|
|
|
|
||||||||||||||
|
|
|
|
|
X1 z |
|
|
|
|
|
|
|
|
. |
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
z |
1 z 2e T |
1 |
|
|
|
|
|
|
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||
|
Используя классический метод, получаем |
|
|
|
|
|
|
|
|
||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
e |
Ts |
|
|
1 |
e |
T |
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
; |
|
|
|
|
|
|
||||||||||
|
|
|
|
|
G z Z |
|
|
|
|
|
|
|
|
z |
e T |
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
s |
s 1 |
|
|
|
|
|
|
|
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
R z |
|
|
|
|
|
z |
|
|
; |
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
z |
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
C z |
|
|
|
G z |
|
; |
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
1 G z |
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
z 1 e T
R z z 1 z 2e T 1 ,
что совпадает с результатом, полученным выше.
78
Если начальное условие задано в виде
x1 0 c0 ,
то при соответствующей подстановке в формулы (2.39) и (2.40) будем иметь
1 z 1
откуда
|
|
M2 z |
|
|
|
z |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
z 1 z 2e T |
1 |
|
|
|
|
|
|
|
||||||||||||
|
|
|
X1 z |
|
|
|
|
|
|
|
|
|
|||||||||||||
2e |
T |
|
|
|
|
z |
1 |
|
T |
|
|
|
|
|
|
z |
|
|
|
|
|
1 |
|
|
|
|
1 1 |
|
1 e |
|
|
|
|
|
|
|
|
c0 |
1 z |
|
|
, |
|||||||||
|
|
1 e |
T |
|
|
|
|
|
|
z 2e |
T |
1 |
|
|
c0 |
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
X1 |
z |
|
|
|
z 1 e T |
|
|
|
|
|
c0 z |
|
|
. |
|
|
|
|
|||||||
z |
1 z 2e T |
1 |
z 2e T |
1 |
|
|
|
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
Из изложенного видно, что по удобству и эффективности рассматриваемый общий подход превосходит классический в анализе сложных систем, а также в оптимальном синтезе.
79
