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

Мехатроника и динамика мини-роботов

.pdf
Скачиваний:
2
Добавлен:
24.11.2025
Размер:
4.17 Mб
Скачать

Учитывая соотношение

 

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