Теория автоматического управления. Часть II. Дискретные системы. Учебное пособие
.pdf
8) Найти оригинал y[i].
Соединения звеньев 1) В отличие от непрерывных систем для случая последова-
тельного соединения звеньев с общей передаточной функцией
k |
k |
W (S) Wi (S) , W (z) Wi (z) , |
|
i 1 |
i 1 |
где передаточная функция W (z) должна определяться по общей весовой функции WП (S).
2) Параллельное соединение звеньев
k |
k |
W (S) Wi (S) , |
W (z) Wi (z) . |
i 1 |
i 1 |
2.8 Характеристики замкнутой импульсной системы
Для определения д.п.ф. замкнутой импульсной системы можно использовать правила структурных преобразований типовых соединений, сформулированные для непрерывных систем. Но при этом следует помнить, что:
1)обычные правила структурных преобразований справедливы для импульсных систем, лишь если каждая ветвь типового соединения представляет собой типовую импульсную цепь, состоящую из идеального квантователя (на входе цепи) и непрерывной части;
2)при иной структуре цепи и всего типового соединения эквивалентная д.п.ф. определяется более сложными правилами. Для основной схемы одноконтурной импульсной системы д.п.ф.
по каналу д-х.
Рис. 2.3. Система с единичной обратной связью
Пусть для системы с единичной обратной связью (рис. 2.3) определена (для общего случая 0) передаточная функция ра-
41
зомкнутой системы W (z, ). Тогда изображение выходной вели-
чины |
|
Y (z, ) W (z, )X (z,0), |
(2.14) |
где X (z,0) - изображение ошибки, так как ИЭ реагирует на |
|
значения Х в дискретные моменты времени t iT ( 0). При
0 |
имеем X (z,0) (z,0) Y(z,0) , подставляя его в (2.13), |
|||||||||||||||||||||||
получим: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Y (z,0) |
|
W (z,0) |
|
|
|
|
|
(z,0) Ф(z,0) (z,0), |
|||||||||||||||
|
1 W (z,0) |
|||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
||||||||||||||||
|
X (z,0) |
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
(z,0) ФX (z,0) (z,0), |
|||||||
|
1 W (z,0) |
|||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
||||||||||||||||
|
Y (z) |
|
W (z) |
|
|
|
(z) Ф(z) (z) , |
|
||||||||||||||||
|
1 W (z) |
|
||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
Ф(z) Y (z) |
|
|
|
|
|
W (z) |
; |
(2.15) |
|||||||||||||||
|
1 |
W (z) |
||||||||||||||||||||||
|
|
|
|
|
G(z) |
|
|
|
||||||||||||||||
|
X (z) |
|
|
1 |
|
|
|
|
|
|
|
|
(z) ФX (z) (z) |
, |
||||||||||
|
1 W (z) |
|||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
Ф (z) |
X (z) |
|
|
|
|
|
|
|
|
1 |
|
|
(2.16) |
||||||||||
|
|
1 W (z) |
||||||||||||||||||||||
|
X |
|
|
|
|
|
G(z) |
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
где Ф(z) - передаточная функция замкнутой системы,
ФX (z) - передаточная функция замкнутой системы по ошибке.
W (z) K(z) - д.п.ф. разомкнутого контура, представляю-
D(z)
щего собой (в данной схеме) типовую импульсную цепь. Характеристическое уравнение импульсной системы
1 W (z) 0
или в развернутых формах
D(z) K(z) F(z) 0
Характеристическое уравнение: |
|
|
|
a zn a zn 1 |
.... a |
0. |
|
0 |
1 |
n |
|
42
Условием применимости полученных формул является требо-
вание |
равенства 0 приведенной весовой |
функции в момент |
t 0, |
П (0) 0. Для этого в системах с бесконечно короткими |
|
импульсами в виде -функций требуется, |
чтобы степень числи- |
|
теля передаточной функции WH (S) по крайней мере на два была
меньше степени знаменателя.
В системах с конечными по длительности импульсами достаточно чтобы разность была бы не меньше, чем 1.
Передаточные функции W (z), Ф(z), ФX (z) - могут быть ис-
пользованы для оценки устойчивости и качества импульсных систем.
Если 0, то учитывая, что X (z,0) - изображение ошибки
Y (z, ) W (z, ) (z,0) . Это выражение практически не
1 W (z,0)
используется.
Кроме того, W (z,0) 1zW (z,1) .
Для случая неединичной обратной связи, рис. 2.4.
Рис. 2.4. Система с неединичной обратной связью
Wgy(z) |
|
W1(z) |
|
|
1 Z{W (S)W (S)} . |
||||
|
1 |
2 |
|
|
Пример. Определим характеристики замкнутой импульсной системы, разомкнутый контур которого соответствует цепи, содержащей «ключ», фиксатор и идеальный интегратор.
Подставляя точную д.п.ф. (2.13) в формулы (2.15) и (2.16), получим соответствующие д.п.ф. замкнутой системы:
43
Ф |
(z) Y (z) |
|
kT |
|
; |
|
|||
(z kT 1) |
|
||||||||
д y |
|
G(z) |
|
|
|||||
Ф |
(z) |
Х (z) |
|
|
(z 1) |
. |
|||
|
(z kT 1) |
||||||||
|
д х |
|
G(z) |
|
|||||
(2.17)
(2.18)
Характеристическое уравнение системы:
z kT 1= 0.
Найдем операторное уравнение динамики системы по каналу д-х. Разделив предварительно числитель и знаменатель д.п.ф.
(2.18) на z, получим |
|
X (z)[1 (kT 1)z 1] G(z)[1 z 1]. |
(2.19) |
Уравнению (2.19) соответствует разностное уравнение в ре-
куррентной форме (при Т=1):
x[i] [1 kT]x[i 1] д[i] д[i 1].
3.УСТОЙЧИВОСТЬ ИМПУЛЬСНЫХ СИСТЕМ
3.1Основное условие устойчивости
Динамические свойства импульсных систем с амплитудной модуляцией во многом аналогичны динамическим свойствам непрерывных систем. Поэтому методы анализа таких импульсных систем являются аналогами соответствующих методов исследования непрерывных систем.
Устойчивость импульсной системы управления, как и устойчивость непрерывной системы, определяется характером еѐ свободного движения. Импульсная система устойчива, если свободная составляющая переходного процесса с течением вре-
мени затухает, т.е. если
ilim yc[iT ] 0
Свободная составляющая yc[iT ] является решением однородного разностного уравнения
(a zn a zn 1 |
|
a )Y (z) 0 |
(3. 1) |
|
0 |
1 |
|
n |
|
Решением уравнения (3.1) при отсутствии у него одинаковых корней z представляет собой сумму:
44
|
y [iT ] |
n |
|
|
|
|
|
|
C zi |
|
(3.2) |
||||
|
c |
|
k 1 |
k |
k |
|
|
|
|
|
|
|
|
|
|
Ck – постоянные интегрирования, |
зависящие от начальных |
||||||
условий; zk – корни характеристического уравнения. |
|
||||||
a zn a zn 1 |
|
a |
|
= 0 |
(3.3) |
||
0 |
1 |
|
|
n |
|
|
|
Из выражения (3.2) |
видно, |
что при i→∞ решение |
yc(iT ) |
||||
стремится к нулю лишь в том случае, если все корни zk |
по моду- |
|
лю меньше единицы, т.е. если |
|
|
| zk | 1, |
k 1;2; ;n |
(3.4) |
Запись (3.4) выражает общее условие устойчивости: для устойчивости импульсной системы необходимо и достаточно, чтобы все корни характеристического уравнения системы находились внутри круга единичного радиуса с центром в начале координат
Если хотя бы один корень zk располагается на окружности единичного радиуса, то система находится на границе устойчивости. При | zk |>1 система неустойчива.
Таким образом, единичная окружность в плоскости корней zk
является границей устойчивости и, следовательно, играет такую же роль, как мнимая ось в плоскости корней pk (рис. 3.1, б). Этот
вывод вытекает также из основной подстановки метода z- преобразования. Действительно, пусть pk k j k . Тогда
z e pkT e kT e j kT |
|
k |
|
и требование | zk |<1 сводится к неравенству: |
|
| e kT | 1, |
(3.5) |
откуда следует известное в теории непрерывных систем условие устойчивости: k <0.
Для суждения об устойчивости импульсных систем можно использовать обычные критерии устойчивости линейных систем, но при этом приходится учитывать лишь некоторые особенности импульсных систем. Так, для того чтобы применить критерий
45
Гурвица, необходимо предварительно в уравнении (3.3) произвести замену переменной z на переменную w путем подстановки
Z (w 1)/(w-1) |
(3.6) |
И получить преобразованное характеристическое уравнение
n |
n-1 |
|
(3.7) |
a0w |
a1w |
an 0 |
Корням уравнения (3.3), расположенным в плоскости корней внутри единичного круга (см. рис. 3.1, а), теперь будут соответствовать корни преобразованного уравнения (3.7), находящиеся в плоскости корней wk слева от мнимой оси. Действительно, если
| zk |<1, то модуль числителя в выражении (3.6) должен быть меньше модуля знаменателя, т.е.| wk 1| | w-1|. А это возможно лишь в том случае, если вектор wk расположен в левой полуплоскости (см. рис. 3.1, в).
3.2 Частотный критерий устойчивости Михайлова
При использовании критерия Михайлова в характеристический полином F (z) подставляют z e j T , изменяют w от 0 до π/T и в комплексной плоскости строят годограф вектора F( e j T ). Импульсная система устойчива, если при возрастании w от 0 до π/T характеристический вектор F( e j T ) повернется против часовой стрелки на угол nπ. Если годограф характеристического вектора проходит через начало координат, то система находится на границе устойчивости.
а |
б |
|
с |
|
|
|
|
|
|
|
|
Im Z |
Im p |
|
-1 |
+1 Im w |
|
|
|
-1 |
+1 |
|
|
|
|
|
|
|
|
|
|||||
|
|
|
w |
|
|
|
|
|
|
|
|
|
|
|
|
- |
|
|
|
|
|
|
|
|
|
|
|
1 w |
w+1 |
|
1 |
|
1 |
|
|
|
|
|
|
|
+ |
|
w |
- |
|
|
|
|
|
|
|
|
|
w |
|
|
|||
|
|
|
|
|
|
w |
|
|
|
||
1 |
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
Re w |
||
Re Z |
|
|
Re p |
|
|
|
|
|
|
|
|
0 |
|
0 |
-3 |
-2 |
-1 |
0 |
|
1 |
2 |
3 |
|
|
|
|
|
||||||||
Рис. 3.1. Области устойчивости в плоскости корней
46
Годографы вектора F( e j T ) для устойчивой и неустойчивой системы второго порядка показаны на рис. 3.2, а.
Имеем характеристический полином:
F(z) a |
zn a zn 1 |
... a |
|
0 |
1 |
|
n |
Подставим z e j |
при Т=1 и учитывая, что |
||
e j cos jsin , в результате получим
a0(cosn jsinn ) ... an U ( ) jV ( ) F( j ) ,
Изменяя частоту от 0 до получим кривую – годограф вектора Михайлова.
Необходимое условие устойчивости:
Характеристический полином
F( j ) 0 0;
F( j ) 0;при четном n
F( j ) 0; при нечетном n
Отметим, что импульсные системы второго и даже первого порядка, в отличие от непрерывных систем такого же порядка, могут быть неустойчивыми при положительных коэффициентах характеристического уравнения. Это объясняется тем, что фиксатор, содержащийся обычно в контуре импульсной системы, вносит дополнительное отставание по фазе.
3.3 Частотный критерий устойчивости Найквиста
Критерий Найквиста для импульсных систем формулируется так же, как и для непрерывных систем: система устойчива, если
АФХ W( e j T ) устойчивого разомкнутого контура не охватывает точку (-1; j0).
Характеристики устойчивой импульсной системы и неустойчивой системы, находящейся на границе устойчивости, показаны штриховыми линиями на рис. 3.2, б.
47
Устойчивость разомкнутого контура импульсной системы определяется устойчивостью еѐ непрерывной части: если последняя устойчива, то и весь контур (включая импульсный элемент) устойчив.
Следует иметь в виду, что, хотя импульсный элемент не влияет на устойчивость разомкнутого контура, он существенно влияет на устойчивость и качество замкнутой системы. При малых периодах повторения частотная характеристика разомкнутого контура совпадает с частотной характеристикой непрерывной части, и устойчивость импульсной системы полностью определяется свойствами непрерывной части. По мере увеличения периода повторения у большинства систем уменьшается предельный передаточный коэффициент и ухудшаются динамические свойства. Однако на некоторые структурно-неустойчивые непрерывные системы и на системы с запаздыванием, АФХ которых заходит в правую полуплоскость, импульсный элемент оказывает стабилизирующее действие. Для таких систем рекомендуется период повторения выбирать из условия
0 – частота, при которой АФХ непрерывной части пересекает положительную мнимую ось Q( ).
jQ( ) |
F e j T |
jQ( ) |
|
|
|
|
|
|
|
|
|
||
|
0 |
|
|
0 |
|
|
|
1 |
|
|
|
||
|
|
|
|
|
||
0 |
P( ) |
|
1 |
P( ) |
||
T |
||||||
|
|
|
|
|
||
|
|
|
|
|
|
|
|
T |
2 |
|
W |
e j T |
|
|
|
|
|
|||
а) |
|
б) |
|
|
|
Рис. 3.2. Критерий Михайлова (а) и Найквиста (б) для импульсной системы
48
Пример. Найти модуль АФХ импульсной системы, состоящей из ИЭ, генерирующего короткие прямоугольные импульсы относительной длительности j 1 и непрерывной части
|
|
|
|
|
|
|
|
|
|
|
|
WH (S) |
|
|
|
k |
|
|
|
. |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
TS 1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
Передаточная функция формирователя WФ(S) |
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
WПН (S) |
|
k |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
TS 1 |
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
Для того, чтобы получить W (z) воспользуемся z- |
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||
преобразованием. |
|
k |
|
|
|
z |
; |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
W (z) k z |
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
T |
|
S 1/T |
T |
|
z d |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
d e 1/T |
|
|
|
|
|
~ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
~ |
|
|
|
|
|
|
|
|
~ |
|
|
|
|
|
|
|
|
|
|
|
|
|||||
подставим вместо z e j cos jsin |
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||
W ( j ) k |
|
|
|
cos j sin |
|
|
|
|
k |
(cos j sin )((cos d) j sin ) |
|
||||||||||||||||||||||||||||||
cos j sin d |
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||
|
|
|
|
|
T |
|
|
T |
|
|
|
|
|
|
|
|
|
|
(cos d)2 |
sin2 |
|
|
|
|
|||||||||||||||||
|
k |
|
cos (cos d ) j sin (cos d ) j cos sin j2 sin2 |
|
|
|
|||||||||||||||||||||||||||||||||||
T |
|
|
|
|
|
|
|
cos2 2 cos d d 2 |
sin2 |
|
|
|
|
||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||
|
k |
|
1 d cos jd sin |
|
k |
|
|
|
|
1 d cos |
|
|
j |
|
d sin |
|
|
|
|||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
1 d 2 2cos d |
|
|
|
|
2d cos d 2 |
|
2d cos d 2 |
||||||||||||||||||||||||||||||||
|
T |
|
|
T |
|
1 |
|
1 |
|
||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
k |
|
|
1 d cos |
|
|
|
2 |
|
|
|
|
|
|
|
|
|
d sin |
2 |
|
|
|
|
||||||||||||||
A( ) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
T |
|
|
|
|
|
|
|
|
2 |
|
|
2d cos d 2 |
|
|
|
|
|||||||||||||||||||||||||
|
|
|
|
|
|
1 2d cos d |
|
|
|
|
1 |
|
|
|
|
|
|||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||
|
1 2d cos d 2 cos2 d 2 sin2 |
|
k |
|
|
|
|
|
1 |
|
|
|
|
|
|||||||||||||||||||||||||||
|
|
|
|
|
(1 2d cos d 2 )2 |
|
|
|
|
|
|
|
|
T |
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
1 2d cos d 2 |
|
|
||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||
Пример. Определим с помощью критерия Гурвица предельное значение передаточного коэффициента импульсной системы, рассмотренной в разделе 2.5.
Характеристическое уравнение системы
Z kT -1 0
49
После подстановки (3.5) уравнение принимает вид
kTw 2- kT 0
и для него условие устойчивости Гурвица заключается, как известно, в положительности коэффициентов, т.е.
|
kT 0 |
|
2 - kT 0 |
a0 |
и a1 |
отсюда допустимые пределы измерения передаточного коэффициента 0<k<2T.
4.КАЧЕСТВО ИМПУЛЬСНЫХ СИСТЕМ
4.1Оценка качества импульсных систем
Качество импульсных систем управления характеризуется такими же показателями, как и качество непрерывных систем:
-точностью в установившихся режимах;
-длительностью и перерегулированием переходного процес-
са.
Длительность и перерегулирование оценивают непосредственно по переходной характеристике. Переходная характеристика импульсной системы строится гораздо проще, чем для непрерывной системы. Для этого записывают z-изображение выходной величины при единичном ступенчатом воздействии
Y(z) [z /(z -1)]Ф(z)
Азатем по изображению находят оригинал – решетчатую функцию y[iT ].
В простых случаях функцию можно найти при помощи
таблиц обратного z-преобразования, разложив предварительно изображение Y (z) на простые дроби.
В тех случаях, когда разложение на дроби связано с трудностями, целесообразно разложить функцию Y (z) в степенной ряд
по отрицательным степеням z (делением числителя на знаменатель):
Y (z) c z i c |
c z 1 |
c |
z 2 ... c z i ... |
(4.1) |
|
i |
0 |
1 |
2 |
i |
|
50
