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

vostrikov

.pdf
Скачиваний:
133
Добавлен:
23.01.2019
Размер:
12.99 Mб
Скачать

7.2. Динамические характеристики линейных импульсных систем

231

Запишем систему уравнений (7.2) для одноканального объекта в скалярной форме

x1(k

1)

a11x1(k)

a12 x2 (k) ...

a1n xn (k)

b1u(k),

x2 (k

1)

a21x1(k)

a22 x2 (k ) ...

a2n xn (k )

b2u(k ),

 

 

 

 

 

(7.15)

xn (k

1)

an1x1(k)

an2 x2 (k) ...

ann xn (k)

bnu(k),

y(k)

c1x1(k) c2 x2 (k ) ... cn xn (k ).

 

Системе уравнений (7.15) соответствует структурная схема, представленная на рис. 7.17.

b1

u(k) ...

bn

x1(k 1)

z 1

 

a11

 

...

 

a1n

 

...

...

 

xn (k 1)

z 1

an1

...

ann

x1(k)

...

xn (k)

c1

... y(k)

cn

Рис. 7.17. Структурная схема, соответствующая уравнениям (7.15)

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

232

Глава 7. ЛИНЕЙНЫЕ ИМПУЛЬСНЫЕ СИСТЕМЫ

7.2.9.ПЕРЕХОД ОТ ПЕРЕДАТОЧНЫХ ФУНКЦИЙ

КСТРУКТУРНЫМ СХЕМАМ

Рассмотрим первый вариант перехода к структурной схеме от передаточной функции

 

y(z)

 

b

 

zn 1

... b z

b

 

 

W (z)

 

 

n 1

 

1

0

.

(7.16)

u(z)

 

zn

an 1zn 1

... a1z a0

 

 

 

 

Передаточную функцию (7.16) представим в виде произведения двух передаточных функций

1

W (z) A(z) B(z)

 

 

 

 

1

 

 

 

 

 

 

b

zn 1

... b z

b .

(7.17)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

zn a

zn 1

 

 

... a z

a

 

 

n 1

 

 

1

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n 1

 

1

0

 

 

 

 

 

 

 

 

 

Введем новую переменную (z) , как показано на рис. 7.18:

 

 

 

u(z)

 

 

 

 

 

(z)

 

 

 

 

 

y(z)

 

 

 

 

 

1

 

 

 

 

B(z)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A(z)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рис. 7.18. Структурное представление системы (7.17)

Для каждого из звеньев запишем операторное уравнение

zn

a

zn 1

...

a z

a

(z) u(z),

 

n 1

 

 

 

1

 

0

 

(7.18)

 

 

 

zn 1

 

 

zn

2

 

y(z)

b

 

b

2

... b z

b (z).

 

n 1

 

n

 

 

1

0

От операторной формы уравнений (7.18) перейдем к их записи в дискретном времени:

(k

n) u(k )

a0

(k )

a1

(k

1)

...

an 1

(k

n 1),

(7.19)

y(k)

bn 1 (k

n

1)

bn 2

(k

n

2)

...

b1 (k

1) b0 (k ).

 

По уравнениям (7.19) построена структурная схема рис. 7.19. Полученная структурная схема позволяет перейти к модели системы

7.2. Динамические характеристики линейных импульсных систем

233

в переменных состояния. Выход каждого звена задержки обозначим как переменную состояния

x1(k) (k), x2 (k) (k 1), xn (k) (k n 1) ,

это позволяет представить разностные уравнения (7.19) в виде

 

 

 

x1(k

1)

 

x2 (k),

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x2 (k 1)

x3 (k),

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

xn (k 1)

a0 x1(k) a1x2 (k) an 1xn (k ) u(k ),

 

 

 

 

 

 

y(k) b0 x1(k) b1x2 (k) bn 2 xn 1(k ) bn 1xn (k ).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

bn 1

 

...

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

b1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ε(k

 

n)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ε(k 1)

 

 

 

ε(k)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u(k)

 

 

 

 

 

 

 

ε(k

n 1)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z 1

 

 

 

 

 

z 1

 

 

 

 

 

b0

 

 

 

 

 

 

 

 

 

 

z 1

 

...

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

an 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

...

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(7.20)

y(k)

Рис. 7.19. Структурное представление передаточной функции в соответствии с системой уравнений (7.19)

Систему разностных уравнений (7.20) можно представить в вектор- но-матричной форме (7.2) с матрицами

0

1

0

 

0

0

 

0

0

1

 

0

0

 

A

, B

, C

b0 b1 bn 2 bn 1 .

0

0

0

 

1

0

 

a0

a1

a2

an 1

1

 

234

Глава 7. ЛИНЕЙНЫЕ ИМПУЛЬСНЫЕ СИСТЕМЫ

Модель системы в переменных состояния (7.20) будем называть «прямой» формой (в литературе встречается другое название – «первая каноническая форма»).

Рассмотрим второй вариант перехода к структурной схеме от передаточной функции (7.16), для чего запишем соответствующее ей разностное уравнение в операторной форме

zn a

zn 1 ... a z

a y(z)

b

zn 1 ... b z

b u(z) . (7.21)

n 1

1

0

n 1

1

0

Перейдем к записи уравнения (7.21) в дискретном времени:

y(k n) a0 y(k) b0u(k) a1y(k 1) b1u(k 1)

an 1y(k n 1)

bn 1u(k

n 1) .

(7.22)

В уравнении (7.22) сдвинем аргумент на один шаг назад:

y(k n 1) a0 y(k 1) b0u(k 1) a1y(k) b1u(k)

an 1y(k n 2) bn 1u(k n 2)

 

и введем новую переменную x1(k)

a0 y(k 1)

b0u(k

1) , тогда

y(k n 1) x1(k) a1 y(k) b1u(k)an 1 y(k n 2) bn 1u(k n 2).

Процедуру сдвига аргумента и замены переменной выполним n раз. В новом базисе разностное уравнение (7.22) можно представить в век- торно-матричной форме

x1(k 1)

a0 xn (k) b0u(k),

 

 

 

(7.23)

xn (k 1)

xn 1(k) an 1xn (k) bn 1u(k ),

 

y(k) xn (k)

7.2. Динамические характеристики линейных импульсных систем

 

 

235

с соответствующими матрицами

 

 

 

 

 

0

0

0

 

a0

 

b0

 

 

 

 

1

0

0

 

a1

 

b1

 

 

 

 

A 0 1 0 a2

, B , C 0

0 0 1 .

 

 

 

 

bn 2

 

 

 

 

0

0

0

 

an 1

 

bn 1

 

 

 

 

Модель системы в переменных состояния (7.23) будем называть

«транспонированной» формой.

 

 

 

 

 

 

По уравнениям (7.23) построена структурная схема (рис. 7.20).

 

 

 

 

 

 

...

 

bn 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

b1

 

 

 

 

 

 

 

u(k)

 

x1(k 1)

x1(k)

x2 (k 1)

x2 (k)

xn 1(k)

xn (k 1)

xn (k)

y(k)

 

 

z 1

 

z 1

...

z 1

 

 

b0

 

 

 

 

 

 

 

 

 

 

 

 

...

an 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a1

 

 

 

 

 

 

 

 

a0

 

 

 

 

 

 

 

Рис. 7.20. Структурное представление передаточной функции в соответствии

 

 

 

с системой уравнений (7.23)

 

 

 

ПРИМЕР 7.6

Задано линейное разностное уравнение, необходимо перейти к системе разностных уравнений

0,5y(k 2) 1, 2 y(k 1) 1,5y(k) 0, 4u(k 1) u(k) .

236 Глава 7. ЛИНЕЙНЫЕ ИМПУЛЬСНЫЕ СИСТЕМЫ

Исходное разностное уравнение приведем к нормированному виду, для этого коэффициент перед y(k + 2) должен быть равен единице:

y(k 2) 2, 4 y(k 1) 3y(k) 0,8u(k 1) 2u(k) ,

далее разрешим это уравнение относительно y(k + 2):

y(k 2)

2, 4 y(k 1) 3y(k) 0,8u(k 1) 2u(k) .

Выполним сдвиг аргумента в последнем уравнении назад на один шаг:

y(k 1)

3y(k 1)

2u(k 1)

2, 4 y(k) 0,8u(k),

введем новую переменную x1(k)

 

3y(k 1)

2u(k 1) , тогда

y(k

1) x1(k)

2, 4 y(k)

0,8u(k) .

Процедуру сдвига аргумента и замены переменной выполним для последнего уравнения еще раз:

y(k) x1(k 1) 2, 4 y(k 1) 0,8u(k 1) , y(k) x2 (k) .

Запишем систему разностных уравнений для новых переменных х1 и х2 :

x1(k 1)

3x2 (k ) 2u(k ),

x2 (k

1)

x1(k ) 2, 4x2 (k ) 0,8u(k ),

y(k)

x2 (k).

В итоге получаем систему разностных уравнений в «транспонированной» форме.

7.3. УСТОЙЧИВОСТЬ ЛИНЕЙНЫХ ИМПУЛЬСНЫХ СИСТЕМ

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

Устойчивость свойство объекта или системы с течением времени приходить в равновесное состояние.

Равновесное состояние такое состояние, в котором все переменные состояния неизменны.

7.3. Устойчивость линейных импульсных систем

237

Запишем исходную систему разностных уравнений для объекта управления:

x(k 1) Ax(k) Bu(k),

(7.24)

полагаем при этом, что управляющее воздействие неизменно и ограниченно:

u(0) u(1) ... u(k n) u const .

Поскольку в состоянии равновесия все переменные состояния в выражении (7.24) неизменны, уравнение равновесия принимает вид

x0

Ax0

Bu ,

(7.25)

где x0 – положение равновесия.

 

 

 

Введем новые координаты

(k)

x(k)

x0 – отклонение от поло-

жения равновесия. Пример движения

системы в отклонениях к поло-

жению равновесия приведен на рис. 7.21.

 

 

n (k)

 

 

 

 

 

(0)

1(k)

Рис. 7.21. Пример движения системы в отклонениях к положению равновесия

Преобразуем исходное уравнение (7.24) в уравнение в отклонениях от положения равновесия:

x(k)

 

(k) x0

x(k 1) (k 1) x0 ,

(k 1) x0 A (k) x0

BU 0 ,

(k

1)

A (k)

Ax0 BU 0

x0 .

238

Глава 7. ЛИНЕЙНЫЕ ИМПУЛЬСНЫЕ СИСТЕМЫ

Последняя скобка равна нулю в состоянии равновесия в соответствии с (7.25), поэтому

(k 1) A (k).

(7.26)

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

7.3.1. ОБЩЕЕ УСЛОВИЕ УСТОЙЧИВОСТИ

Линейная импульсная система (7.26) является устойчивой, если

для ее процессов выполняется условие

lim (k) 0 .

(7.27)

k

Утверждение: для того чтобы линейная импульсная система была устойчивой, необходимо и достаточно, чтобы модули собственных чисел матрицы А были меньше единицы:

| zi | 1; i 1, n ,

где zi собственные числа матрицы A (корни характеристического уравнения).

Докажем это утверждение, для чего запишем характеристическое уравнение объекта

A(z) det(zI A) zn a

zn 1

a z

a 0 .

n 1

 

1

0

Будем искать частные решения разностного уравнения в отклоне-

ниях от положения равновесия (7.26) в виде

 

i

(k)

zk

i

,

(7.28)

 

i

 

 

 

 

 

 

 

i1

где zi – корень характеристического уравнения; i

... – соответ-

in

ствующий этому корню собственный вектор.

7.3. Устойчивость линейных импульсных систем

239

Общее решение разностного уравнения (7.26) представляет собой линейную комбинацию частных решений (7.28) и имеет вид

 

n

 

 

 

(k )

C zk

i

.

(7.29)

 

i i

 

 

i

1

 

 

 

Как следует из анализа выражения (7.29), для выполнения предельного соотношения (7.27) необходимо, чтобы каждая компонента общего решения разностного уравнения стремилась к нулю, тогда и вся сумма будет стремиться к нулю.

Если все корни характеристического уравнения по модулю меньше единицы, то выполняется условие

lim zk

 

 

0, i 1, n

k

i

 

 

 

 

 

и, следовательно, выполняется соотношение (7.27), что и требовалось доказать.

Для нахождения собственных векторов i частное решение (7.28) подставим в однородное разностное уравнение (7.26):

zk

1

Azk

i

,

 

i

i

i

 

 

из которого следует

 

 

 

 

 

(zi I

A) i

0 .

(7.30)

Поскольку в выражении (7.30) матрица в круглых скобках вырожденная (так как zi – собственное число матрицы A), для вектора i существует бесконечное множество решений. Для нахождения любого из них произвольно задают одну из компонент этого вектора. Из системы исходных уравнений исключают зависимое уравнение, остальные компоненты вектора i находят по сформированной системе уравнений (n 1)-го порядка.

Для нахождения неизвестных констант Сi используют начальные условия:

 

n

 

(0)

C z0

.

 

i i i

 

i

1

 

Это система из n уравнений для нахождения констант Ci , i 1, n .

240 Глава 7. ЛИНЕЙНЫЕ ИМПУЛЬСНЫЕ СИСТЕМЫ

7.3.2. ГЕОМЕТРИЧЕСКАЯ ИНТЕРПРЕТАЦИЯ ОБЩЕГО УСЛОВИЯ УСТОЙЧИВОСТИ

Изобразим плоскость корней линейной импульсной системы

(рис. 7.22).

Очевидно, линейная импульсная система (ЛИС) устойчива, если все ее корни лежат в круге единичного радиуса.

 

Im z

 

j

1

1 Re z

 

j

Рис 7.22. Область устойчивости в пространстве корней

Процедура анализа устойчивости линейной импульсной системы:

1)записать характеристическое уравнение det(zI – A) = 0;

2)найти корни zi .

3) проанализировать zi по критерию | zi |< 1, i 1, n.

7.3.3. БИЛИНЕЙНОЕ ПРЕОБРАЗОВАНИЕ

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

Соседние файлы в предмете Основы Теории Управления