- •Дискретные системы автоматического управления. Основные положения.
- •Математические методы описания дискретных систем Решетчатые функции
- •В соответствии с этим значение рассматриваются на полуинтервале
- •Функции f[nT, t] являются функциями двух переменных (аргументов) n и , поэтому целесообразно обозначать эти функции как
- •Дискретное преобразование Лапласа.
- •Свойства z-преобразования.
- •Передаточные функции дискретных систем
- •Определение процессов в импульсных системах при типовых воздействиях.
- •Далее определим вычеты в полюсах передаточной функции ф(z,). Для простых полюсов получим
- •Учитывая, что ф(z,) является дробно-рациональной по отношению к переменной z и обозначая
- •Анализ устойчивости дискретных систем.
- •Критерий Рауса-Гурвица.
- •Частотные характеристики
- •Построение лчх дискретных систем Построение низкочастотной части лчх
- •Построение высокочастотной части цифровых систем с экстраполятором нулевого порядка.
- •Синтез цифровых систем автоматического регулирования. Обеспечение заданной точности.
- •Расчет дискретных корректирующих средств.
Свойства z-преобразования.
Свойство линейности. Для краткости записи приводим при =0. Равенства справедливы и при 0.
Если решетчатые функции f1[n], f2[n]… fk[n] являются оригиналами и их изображения соответственно F1(z), F2(z)…Fk(z), то справедливо равенство
Z{f[n]}= F(z), где -произвольная постоянная.
Пример. Найти изображения тригонометрических решетчатых функций sin n и cos n (n0)
Z{ sin n}= Z{ljn- l-jn}= [Z{ ljn }-Z{ l-jn }]
Используя результаты предыдущего примера, получим
Z{ sin n }= [ - ]=[]=
(T=1)
Аналогично для косинуса
Z{ cos n }=
2. Теорема смещения. Для функции времени f(t-), где - положительное число, причем 0< и f(t-)0 при t<<
Z{f(t-)= Z-(1+m)F(z, 1+-)},
Здесь m- целая, - дробная часть числа , а
F(z, )= Z{f(t)} (*)
Если <1, то изображение равно
Z{f(t-)}=Z-mF(z, -).
В частном случае запаздывание может составлять целое число периодов дискретности т.е.
=mT. Тогда
Z{f(t-mT)=z-mF(z, )}
3. Теорема об умножении оригинала f(t) на экспоненту.
Если для оригинала f(t) изображение определяется формулой
F(z,)=Z{f(t)}(*) определяемой в дискретные моменты времени t=nT+, то
Z{eT f(t)}= dF(,), где - комплексное число, а d=eT
4. Теорема об умножении оригинала на смежную функцию.
Пусть для оригинала f(t) изображение определяется формулой (*), тогда для =0 имеем
Z{tmf(t)}=(-1)m eqT=Z
5. Изображение конечных разностей.
Z{kf[n,]}=(z-1)kF(z,)-
Если дискретной последовательности f[n, ] соответствуют нули в первых k точках оси времени, т.е.
f[0, ]=f[1, ]=…f[(k-1), )]= 0, то f[0,]=0 =1,k-1
Тогда Z{kf[n,]}=(z-1)kF(z,).
Для обратной разности справедливо выражение
Z{k f[n, ]}= ()k F(z, ),
6. Изображение конечной суммы.
Если для оригинала f(t) изображение определяется формулой (*), то
Z{ f[,]}=
7. Сверка решетчатых функций.
Если Z{f1[n,]}=F1(z,);
Z{f2[n,]}=F2(z,), то
Z{ f1[,]f2[(n-)],}= Z{f1[(n-),]f2[,]}=F1(z, )F2(z, )
8. Теорема о конечном значении оригинала.
Если f[n, ] - оригинал, а F(z, )- изображение, то
lim f[n, ]=lim F(z,)
n z1
9. Теорема о начальном значении оригинала. При тех же условиях f[0,]=lim F(z,).
z
Передаточные функции дискретных систем
Чтобы найти передаточные функции дискретных систем, как и в случае непрерывных систем, необходимо первоначально определить передаточную функцию разомкнутой системы. Структурная схема на рисунке.
g x
T T
Передаточная функция
W(Z)=,
Где X(z), G(z ) - изображение функций времени x(t), g(t). Рассмотрим первоначально простейший случай, когда D(z)=1. В этом случае последовательно соединены идеальный импульсный элемент, экстраполятор и непрерывная часть системы. Экстраполятор и непрерывная часть в совокупности образуют приведенную непрерывную часть, передаточная функция которой определяется выражением
W(s) =Wэ(s) W0(s)
Импульсной переходной функцией или функцией веса приведенной непрерывной части системы является оригинал изображения W(s)
k (t)=L-1[W(s)]
Тогда можно записать
W (Z) =Z{k(t)}=Z{W(s)}.
Таким образом, чтобы найти передаточную функцию W(z) необходимо определить Z-преобразование передаточной функции приведенной непрерывной части.
Пусть в системе используется экстраполятор нулевого порядка. Тогда
Wэ(s)=L{kэ(t)}= .
kэ
1
0 T T t
kэ(t)=1(t)-1(t-T)
L[1(t)]=
L[1(t-T)]=e-TS(), т.е. Wэ(S)=,
(если импульсный элемент формирует импульсы с последовательностью , то Wиэ(S)=)
В этом случае
W(s) =(1-e-Ts) =W1(e-Ts) .
Для отыскания Z-преобразования функции времени, заданной преобразованием Лапласа в такой форме можно воспользоваться выражением
Z{W(s)}=W1(z)Z{ }
Поэтому
W(z)=Z{ }=
Z-преобразование Z { } находится путем разложения рациональной дроби на простые с дальнейшим использованием таблиц Z-преобразований.
Пример.
Пусть =
Тогда пользуясь таблицей, будем иметь
W(z)=K Z{} =K [] ,
где d=e-()
Для большинства цифровых систем управления D(z)1. В этом случае
W (z)=D(z)Z{W(s)}),
т.к. ЭВМ или цифровое управляющее устройство, осуществляя модуляцию последовательности входных -функций, не изменяет дискретной природы сигналов.
В любой момент времени nT ЭВМ в соответствии с алгоритмом работы определяет и выдает на выход числовую величину, получаемую в общем случае по значениям входного и выходного сигналов в данный и предшествующие моменты времени.
x1 x2
x2(n)= ix1(n-i)- ix2(n-i) (*)
где i- целое положительное число
Учитывая, что
Z{f[n-i]}=Z-IF(z), где F(z)=Z{f[n]}
Можно записать, что
D(Z)= и применяя к (*) Z-преобразование
X2(z)= iz-iX1(z)- iz-iX2(z) или
X2(z)= Iz-i=X1(Z)iz-i,
где 0=1
Откуда D(Z)=
Умножив числитель и знаменатель на zl , получим
D (Z)=
Пример
Пусть ЭВМ реализует функцию корректирующего устройства с алгоритмом
X2[n]=d1x1[n]+d2x1[n]=d1x1[n]+d2{x1[n]-x[n-1]}=
=(d1+d2)x[n]-d2x[n-1]=0x[n]-1x[n-1]
X2[Z]= 0X1[z]- 1z-1X1[z];
Тогда
D(Z)= 0-1z-1=
Зная D(z) всегда можно определить W(z).
Рассматривая полученные выражения, следует учитывать особый вид дискретной передаточной функции, которую не следует путать с обычной передаточной функцией непрерывного звена, т.к. W(z)W(s) при замене в последней символа s на z.
Для систем с единичной обратной связью, если W(z) представляет собой передаточную функцию разомкнутой системы
Ф(z)=
и передаточную функцию по ошибке
Ф(z)=
и по возмущению
Фf(z)=
Знаменатель рассмотренных передаточных функций замкнутой системы называется ее характеристическим полиномом. Обычно изображения входных сигналов и передаточные функции представляют собой дробно-рациональные функции z. Они позволяют использовать различные оценки качества систем регулирования.