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

книги / Теория автоматического управления. Дискретные системы

.pdf
Скачиваний:
3
Добавлен:
19.11.2023
Размер:
72.74 Mб
Скачать

Различают три типа желаемых ЛЧХ дискретных систем:

1) cистемы с астатизмом 2-го порядка типа 2-1-0 (рис. 7.1, а)

 

k

ε

 

 

(1jλT / 2)

 

T

 

 

W ( jλ) =

 

 

(1+ jλτ1 )

(1+ jλT / 2)

1+ jλ

 

TΣ

;

( jλ)

2

 

 

 

 

 

2

 

 

а

б

б

Рис. 7.1. Типовые желаемые ЛЧХ дискретных систем

2) системы 1-го порядка (следящие системы) типа 1-2-1-0

 

k

v

 

(1+ jλτ )

 

(1jλT / 2)

 

T

 

 

W ( jλ) =

 

 

1

 

(1+ jλT / 2)

1+ jλ

 

TΣ

;

jλ

(1+ jλT1 )

 

 

 

 

 

2

 

 

121

3) статические системы регулирования 0-1-2-1-0

 

k(1+ jλτ1 )

 

(1jλT / 2)

 

T

 

W ( jλ) =

 

 

 

1+ jλ

 

TΣ .

(1+ jλT1 )(1+ jλT2 )

(1+ jλT / 2)

2

 

 

 

 

 

7.1.1. Построение желаемых ЛЧХ по заданной точности

При малых значениях периода квантования в области средних и малых частот характеристики непрерывных и дискретных систем практически не отличаются (ω ≈ λ), поэтому построение желаемых частотных характеристик дискретных систем основано на методике для непрерывных систем. При построении исходят из требования ограничения максимальной ошибки системы управления при воспроизведении задающего воздействия, как наиболее общего показателя качества.

Пусть задающее воздействие на систему v(t)=vmaxsin(ωkt). Из передаточной функции замкнутой системы по ошибке имеем

emax =

 

vmax

 

 

 

 

.

 

1+W ( jωk)

Поскольку emax << ymax и

 

W ( jωk)

 

>>1,

 

 

из (7.1) будет

 

 

 

 

 

v

emax (max ).

W jωk

(7.1)

максимальная ошибка

(7.2)

Из этого условия вытекает требование к низкочастотной асимптоте амплитудно-частотной характеристики статических систем: левее точки ωk=λk должно выполняться условие, чтобы при λ < λk

Lk) ≤ 20 lg

vmax

.

(7.3)

 

 

e

 

 

max

 

Для среднечастотного участка частотных характеристик ис-

пользуют условия ограничения скорости и ускорения изменений

122

входного сигнала. Пусть на входе системы известны максимальные по модулю значения первой (скорости) и второй (ускорения) производных входного сигнала v(t). Для эквивалентного гармонического этому сигналу режима vэ(t)=vэmaxsin(ωэt), который характери-

 

 

 

 

 

 

′′

2

 

зуется производными vmax =ωэvэmax ,

vmax =ωэvэmax ,

можно опреде-

лить частоту и амплитуду:

 

 

 

 

 

 

 

 

 

 

 

v′′

 

 

 

(v

)2

 

 

ω

 

=

max

,

v

=

max

.

(7.4)

 

 

 

э

 

 

эmax

 

′′

 

 

 

 

 

 

vmax

 

 

 

vmax

 

Если ωэ<<T1 ээ), то на плоскости ЛЧХ можно отметить точку Аk, определяющую «запретную зону» (рис. 7.2).

а

б

Рис. 7.2. Запретные зоны для астатических (а) и статических (б) систем

Если амплитуда первой производной равна vmaxи остается не-

изменной, то уменьшение амплитуды ее производной второго порядка в соответствии с (7.4) приведет к снижению ωээ и увеличению vэmax . При этом точка Ак будет перемещаться вдоль асимптоты

с наклоном 1 влево, причем, если vmax′′ устремить к нулю, то λэ 0. Это соответствует режиму с постоянной скоростью v′=vmax′′ , и тогда имеем

v′′

emax = max , (7.5)

Kv

где Kv добротность по скорости.

123

Если рассмотреть режим, когда амплитуда второй производной гармоники равна максимальному значению vmax′′ и амплитуда первой производной ниже заданного значения vmax, то аналогич-

ными предыдущему случаю рассуждениями можно установить, что точка Аk будет перемещаться в область более высоких частот, двигаясь по асимптоте с наклоном 2. Тогда этот режим соответствует режиму с постоянным ускорением, для которого имеем

v′′

emax = max , (7.6)

Kε

где Kε – добротность системы по ускорению.

Определим «запретную зону» для статических систем аналогично построению ее для астатических. Координаты точки Аk определим также аналогично при эквивалентном гармоническом сигнале на входе, при этом эквивалентная круговая частота должна быть выбрана из условия ωэvmax =vmax. Тогда координаты точки Аk будут,

как показано на рис. 7.2, б:

 

 

 

vmax

 

 

 

ωэ1 =

vmax

,

L(ωэ1 )=20lg

.

(7.7)

 

 

 

 

 

v

 

e

 

 

 

max

 

max

 

Если у эквивалентного гармонического сигнала vэ max=vmax и

, то при уменьшении ωэ точка Аk будет перемещаться

ωэvmax <vmax

влево по асимптоте с 0 наклоном. При vэ max<vmax и

ωэvmax =vmax ,

наоборот, точка Аk будет двигаться по асимптоте с наклоном 1 до тех пор, пока vэ′′<vmax′′ , т.е. пока выполняется условие ωэvэmax vmax′′ .

v′′

Отсюда находим абсциссу точки Аk на асимптоте 1 ω = max . Да-

э2 vmax

лее наступает режим с постоянным ускорением в соответствии с (7.4), и точка Аk движется до оси частот.

124

7.1.2. Построение желаемых ЛЧХ с заданным запасом устойчивости

Для синтеза дискретных систем удобно применение частотных методов оценки запасов устойчивости, при этом пользуются показателем колебательности систем в замкнутом состоянии, значения которого могут быть установленными в пределах 1,3–1,7. Связь запасов устойчивости с колебательностью может оцениваться с использованием любых частотных характеристик, например, W(ejωТ), W(jλ). Условием соблюдения требуемых запасов устойчивости, как известно, является запрет на вхождение частотных характеристик внутрь М-окружностей (линий равных значений колебательности), как показано на рис. 7.3.

Рис. 7.3. Желаемые ЛПЧХ с требуемыми запасами устойчивости

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

При построении желаемых логарифмических псевдочастотных амплитудных характеристик используются расчетные формулы их параметров, обозначенных в соответствии с рис. 7.1–7.3 и связанных с показателем колебательности для «симметричных» типовых характеристик:

при ν = 2 (см. рис. 7.1, а, 7.2, а, 7.3)

125

ω =

1

, τ =

 

1

 

М

,

τ=

 

1

 

М

;

τ

ω М 1

 

ω

 

М +1

э

1

 

 

 

 

 

 

1

 

 

ср

 

 

 

 

 

ср

 

 

 

h =

τ1

, ω =

 

K

;

 

 

(7.8)

 

 

 

 

τ

 

э

 

ε

 

 

 

 

 

при ν = 1 (см. рис. 7.1, б, 7.2, б)

 

 

 

 

 

 

ω =

 

Kν

 

,

K T =

M 2 + M M 2 1

 

, τ =

1

 

 

М

;

(7.9)

 

τ

2

 

ω

М 1

э

 

 

 

 

ν 1

 

 

1

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

э

 

 

 

 

 

при ν = 0 (см. рис. 7.1, в, 7.2, б)

 

 

 

 

 

 

 

 

 

ω =

 

K

 

,

 

KT

=

 

M 2 + M M 2 1

, τ =

 

1 М

 

.

(7.10)

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

TT

 

T

2

 

 

ω

 

М 1

э

 

 

 

 

 

1

 

 

 

 

 

1 2

 

 

1

 

 

 

 

 

 

э

 

 

 

 

 

7.2. Синтез регуляторов дискретных систем

При проектировании дискретных систем автоматического управления имеется множество структурных решений задачи, в том числе традиционных.

На рис. 7.4 приведены типовые структурные схемы включения регуляторов в дискретных системах.

а

б

Рис. 7.4. Структурные схемы скорректированных дискретных систем (см. также с. 127)

126

в

г

Рис. 7.4. Окончание

Вдискретных системах применяются как аналоговые, так и цифровые регуляторы, в виде последовательного соединения в прямой цепи (рис. 7.4, а, б) или в виде внутренних обратных связей

(рис. 7.4, в, г).

7.2.1.Синтез аналогового регулятора в прямой цепи

Всоответствии с (4.16) запишем импульсную передаточную функцию системы как комплексную частотную функцию разомкнутой системы:

 

 

 

1

 

 

 

 

 

 

 

 

Wр* ( jω)=

 

WЭ (jω+ jrωи )Wp (jω+ jrωи )WНЧ (jω+ jrωи )=

 

 

 

 

 

 

 

T r=−∞

 

 

 

 

 

 

 

 

1

sin

ω+ rωи

T

 

j(ω+rωи )

T

(7.11)

 

 

2

 

 

 

 

 

 

 

2

 

=

 

T

 

 

 

 

e

 

 

Wp (jω+ jrωи )WНЧ (jω+ jrωи ).

 

 

ω+ rω

 

 

 

 

 

T r=−∞

 

 

 

 

и

 

 

 

 

 

 

 

 

 

 

 

2

 

T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

127

Поскольку в большинстве систем управления регулятор имеет

 

sin

ω+ rωи

T

свойства фильтра нижних частот, а отношение

 

2

 

на

ω+ rω

 

 

 

и T

 

 

 

 

2

 

 

низких частотах приблизительно равно единице, то выражение (7.11) можно аппроксимировать только первым членом ряда (r=0). Это означает, что экстраполятор нулевого порядка заменяется звеном чистого запаздывания, равным половине периода квантования, как было показано на рис.1.9. Таким образом, получаем эквивалентную непрерывную систему с передаточной функцией в разомкнутом состоянии:

Wрc (р)=WР (р)WНЧ (р)ер

Т

 

2

.

(7.12)

По желаемой передаточной функции, обеспечивающей дискретной системе требуемые показатели качества, с учетом (7.12) определяется структура и параметры непрерывного регулятора.

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

7.2.2. Синтез цифрового регулятора в прямой цепи

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

Пример 7.1. Определить структуру цифрового регулятора систе-

мы с передаточной функцией объекта

WНЧ (р)=

 

K

 

, если

p T p +1

 

 

(

1

)

 

K=100 c–1, T1=0,1 c, период квантования сигнала T=0,01 c. Запаздывание в регуляторе и ошибку от квантования по уровню не учитывать. Показатель колебательности в системе не должен превышать 1,5.

128

Передаточная функция разомкнутой исходной системы в фор-

 

 

 

 

 

 

 

 

 

 

 

W (z)

 

z 1

WНЧ (р)

 

ме

z-преобразования

 

будет

=

 

 

Z

 

 

=

 

z

p

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T

 

T 1d

)

 

 

T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= K

1 (

, где

d = e T1 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z 1

 

z d

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

После подстановки параметров системы получим

 

 

 

 

 

 

 

 

 

W (z)

=

0,04837z + 0,04679

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(z 1)(z 0,9048)

 

 

 

 

 

 

 

На рис. 7.5 приведен переходный процесс исходной замкнутой системы, показатели качества которой оказываются неудовлетворительными.

Рис. 7.5. Переходная характеристика исходной системы

129

Выполним синтез цифрового фильтра с помощью билинейно-

го преобразования. Используя подстановку z = 11+ww в передаточ-

ной функции, получим W (w)=

(0,952 + 0,00158w)(1w).

 

 

 

 

 

 

2w(0,952 +1,9048w)

 

Сделаем замену w = j

T

λ

и получим комплексную псевдо-

2

 

 

 

 

 

 

 

частотную

функцию

разомкнутой

исходной

системы

W (jλ)=

100

(1+ 0,000083 jλ))10,005 jλ.

 

 

jλ

 

 

 

(1+ 0,1 jλ)

 

 

 

Построение желаемой ЛАПЧХ проводим из условия, что передаточный коэффициент разомкнутой системы (добротность по скорости) и постоянная времени объекта управления неизменны. Это означает, что ордината Ak «запретной зоны» (см. рис. 7.2, а) известна и равна L1), с контрольной частотой, равной сопрягающей

частоте ω1.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Базовая частота λ0 =

Kν

=

 

100

= 31,6 c1.

 

 

 

 

 

 

 

 

 

0,1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T1

 

 

 

 

 

 

 

 

 

 

Постоянная времени, соответствующая левой границе средне-

частотной асимптоты ЛАПЧХ,

определится как

τ =

1

 

M

 

=

λ0

 

M 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

=

1

 

 

1,5

 

 

= 0,055 c.

 

 

 

 

 

 

 

 

 

 

 

 

31,6

 

 

1,5 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Частота среза желаемой ЛАПЧХ определится из соотношения

τ

=

1

 

 

M

.

При известных

значениях M и

τ

частота среза

λср M 1

1

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

λср = 54,5 c1.

Верхняя граница среднечастотного диапазона определяется из условия «симметричности» желаемой псевдочастотной ампли-

130