Vedyakov_articles / Операторный метод анализа и синтеза ЛСУ
.pdfТогда допустима передаточная функция (2.41), если исходные требования удовлетворяют дополнительному условию
|
δ f ≥ |
|
10ωy2 |
. |
|
|
|
(2.44) |
||||||||
|
|
δ |
y |
ω2 |
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
f |
|
|
|
|
|
|||
В противном случае, но при условии |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
1 |
|
, δ f > |
1 |
|
|
(2.45) |
||||||||
ωf > µ |
|
|
|
|
|
|||||||||||
T |
|
µ3 |
|
|
||||||||||||
|
|
|
|
3 |
|
|
|
|
|
|
|
|
|
|
||
можно использовать функцию вида |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
Wp (s) = |
|
|
|
|
|
|
kp (T2s +1) |
, |
(2.46) |
|||||||
s(T1s |
+1)(T3s +1)(T4s +1) |
|||||||||||||||
|
|
|
||||||||||||||
с теме же параметрами, но при |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
T |
|
|
= |
|
1 |
T . |
|
|
|
(2.47) |
|||||
|
|
|
µ |
|
|
|
||||||||||
|
|
4 |
|
|
|
3 |
|
|
|
|
|
|||||
Все записанные выше формулы вытекают из вида асимптотических ЛАХ, представленных на рис.2.5.
Когда найдена желаемая передаточная функция разомкнутой системы Wp (s) , то
для того, чтобы обеспечить выполнение тождества (2.16) необходимо положить
W2 (s) =Wy−1 (s)Wp (s) . |
(2.48) |
На этом этапе проверяется условие строгой реализуемости передаточной функции регулятора. Последнее означает, что степень числителя не должна превышать степени знаменателя. Если степень числителя передаточной функции W2 (s) выше сте-
пени знаменателя, то можно воспользоваться законом
W (s) =W −1 |
(s)W (s) |
1 |
, |
|
1 |
>> ω . |
(2.49) |
|
∏(Tis +1) |
Ti |
|||||||
ε |
p |
|
c |
|
||||
i
После определения передаточной функции регулятора следует найти характеристический полином замкнутой системы (2.18) и по его корням оценить устойчивость системы. Рассчитать передаточную функцию разомкнутой системы Wp (s) =Wy (s)W2 (s) и определить запас устойчивости по амплитуде и фазе. Запасом
21
Lp (ω) |
|
|
|
-20 дБ/дек |
|
|
|
−20 lg δy |
|
−lg δ1 |
|
|
|
|
|
20 lg δf |
lg ωy |
lg ωf |
lg ω |
|
|
|
|
а) ЛАХ передаточной функции (2.40) |
|
||
Lp (ω) |
|
|
|
-20 дБ/дек |
|
|
|
−20 lg δy |
ωc |
1/T3 |
|
|
|
||
1/T1 |
lg ωy 1/T2 |
lg ωf |
lg ω |
20lg δf |
|
|
|
b) ЛАХ передаточной функции (2.41) |
|
||
Lp (ω) |
|
|
|
-20 дБ/дек |
|
|
|
−20 lg δy |
|
1/T3 1/T4 |
|
|
ωc |
|
|
1/T1 |
lg ωy 1/T2 |
lg ωf |
lg ω |
20lg δf |
|
|
|
c) ЛАХ передаточной функции (2.46) |
|
||
Рис. 2.5. ЛАХ допустимых передаточных функций
22
устойчивости по амплитуде называется величина ∆L = −Lp (ω) , |
где ω таково, что |
||||||
ϕp (ω) = −π. Величина ∆ϕ = π+ϕp (ωc ) называется запасом устойчивости по фазе. |
|||||||
Для удовлетворительной работы системы необходимо выполнение условий |
|||||||
|
∆L ≥ 6 дБ, ∆ϕ ≥ 300 . |
|
|
|
|
(2.50) |
|
Если последние требования не выполнены, то следует видоизменить желаемую ЛАХ |
|||||||
разомкнутой системы и заново определить передаточную функцию W2 ( p) . |
|||||||
В заключении раздела строится ЛАХ и ЛФХ передаточной функции разомкнутой |
|||||||
системы Wp (s) . |
|
|
|
|
|
|
|
2.4. Построение электронной модели регулятора |
|
|
|
||||
Построение электронной модели регулятора включает три этапа. На первом этапе |
|||||||
осуществляется преобразование модели вход – выход регулятора |
|
||||||
|
|
|
|
~ |
|
|
|
|
u =W0 ( p)(−W1( p) y +W2 ( p)ε) |
|
|
|
|||
к модели вход-состояние-выход |
|
|
|
|
|
|
|
& |
|
~ |
|
|
|
~ |
(2.51) |
x = Ax + B1 y + B2 ε , u = Cx + D1 y + D2 ε , |
|||||||
где матрицы A, B1 , B2 , C, D1 , D2 |
такие, что |
|
|
|
|
|
|
|
−W ( p)W ( p) = C(Ip − A)−1 B + D |
, |
|
||||
|
0 |
1 |
|
1 |
1 |
|
|
|
W ( p)W ( p) = C(Ip − A)−1 B + D |
, |
|
||||
|
0 |
2 |
|
2 |
2 |
|
|
а также построение по уравнениям (2.51) структурной схемы на элементарных |
|||||||
звеньях: сумматор, интегратор, усилитель. На z1 |
R1 |
|
|
C |
|||
следующем этапе в построенной структурной |
z2 |
R2 |
|
|
|||
схеме сумматоры, интеграторы и усилители |
|
|
R |
||||
|
|
|
|
||||
заменяются на блоки состоящие из операци- |
|
|
|
|
|||
zk |
Rk |
|
|
z |
|||
онного усилителя, резисторов и конденсато- |
|
|
|
|
|||
ров, электрическая схема которого представ- |
|
|
|
|
|
||
лена на рис. 6. Данный |
блок |
осуществляет |
|
|
|
|
|
преобразование сигналов |
z1 , z2 ,K, zk в сиг- |
|
Рис. 2.6. Электронный блок |
||||
|
|
|
|
|
|||
|
|
|
|
|
|
|
23 |
нал z по формуле
1 |
|
k |
1 |
|
|
|
|
|
+ Cp z = −∑ |
|
zi . |
(2.52) |
|
|
Ri |
|||||
R |
|
i=1 |
|
|
||
На последнем этапе, составляются уравнения электронной модели регулятора и находятся значения сопротивлений резисторов и емкости конденсаторов включенных в схему.
После построения электронной модели необходимо скорректировать значения коэффициентов передаточных функций регуляторов в соответствии с существующими номинальными значениями параметров электронных элементов, использованных при построении электронной модели.
2.5. Исследование замкнутой системы управления
Исследование замкнутой системы управления производится по структурной схеме на рис. 2.3 и включает проведение следующих вычислительных экспериментов:
1)построение процесса y(t) при y (t) =1(t), f (t) = 0 ;
2)построение процессов ε(t), y (t) при y (t) = a1t, f (t) = 0 ;
3)построение процессов ε(t), y (t) при y (t) = ay cos(ωy t), f (t) = 0 ;
4) построение процессов |
ε(t), f (t) при y (t) = 0, f (t) = |
a |
f cos(ω f t) ; |
5) построение процессов |
y(t), y (t), f (t) при y (t) = a y cos(ωy t), f (t) = a f cos(ω f t) . |
||
По результатом экспериментов требуется определить время переходного процесса, перерегулирование, а также максимальные по модулю значения установившейся ошибки в каждом эксперименте.
24
3. ПРИМЕР РАСЧЕТА СИСТЕМЫ УПРАВЛЕНИЯ
3.1. Постановка задачи управления
Рассмотрим процедуру синтеза регулятора для объекта управления, структурная схема которого представлена на рис. 3.1. Будем считать, что следящая система должна удовлетворять следующим требованиям:
1)запас устойчивости по амплитуде не менее 6 дБ, а запас устойчивости по фазе не менее 300;
2)установившаяся ошибка ε∞ (t) отработки сигнала
y* = a0 + a1t, | a1 |≤ a1 = 0.5 |
(3.1) |
|||
должна удовлетворять условию |
|
|||
|
ε∞ (t) |
|
≤ ε1 = 0.05 ; |
(3.2) |
|
|
|||
3) установившаяся ошибка ε∞(t) отработки сигнала |
|
|||
y* = ay cos(ωyt), | ay |≤ ay = 0.25, | ωy |≤ ωy = 1 |
(3.3) |
|||
должна удовлетворять условию |
|
|||
| ε∞ (t) |≤ ε y = 0.05 ; |
(3.4) |
|||
Рис. 3.1. Структурная схема объекта управления
25
4) установившаяся ошибка ε∞(t) , вызванная наличием помехи |
|
f = a f cos(ωf t), | a f |≤ a f = 0.1, | ωf |≥ ωf =500 |
(3.5) |
должна удовлетворять условию |
|
| ε∞(t) |≤ εf = 0.01. |
(3.6) |
3.2. Анализ объекта управления
По структурной схеме определим передаточную функцию. Введем дополнительные переменные как показано на рис. 3.2 и представим переменную y(t) относи-
тельно входа u(t) :
y = x |
+ x |
|
= |
1 |
x |
|
+ x |
|
= |
|
1 + p |
x |
|
= |
p +1 1 |
(u + y) . |
(3.7) |
||
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
p p + 2 |
||||||||||||
1 |
|
2 |
|
p |
2 |
|
2 |
|
|
p |
2 |
|
|
|
|||||
Рис. 3.2. Структурная схема объекта управления |
|
||
Откуда следует |
|
|
|
( p2 + 2 p − p −1) y = ( p +1)u |
(3.8) |
||
и передаточная функция имеет вид |
|
|
|
W (s) = |
s +1 |
. |
(3.9) |
|
|||
|
s2 + s −1 |
|
|
Характеристический полином системы s2 + s −1 = 0 имеет корни с положительной вещественной частью и, следовательно, объект управления неустойчив. Корень чис26
лителя является нулем передаточной функции и равен -1. Таким образом, данная передаточная функция является минимально-фазовой и решение задачи управления существует.
Логарифмическая амплитудная частотная характеристики (ЛАХ) и логарифмическая фазовая частотная характеристики (ЛФХ) представлены на рис. 3.3.
3.3. Решение задачи стабилизации
Представим сигнал управления в виде
|
|
|
u =W0 ( p)(u1 +u2 ) , |
(3.10) |
где W ( p) = b0 ( p) |
- дробно – рациональная функция. |
|
||
0 |
a0 |
( p) |
|
|
|
|
|
||
Рис. 3.3. ЛАХ и ЛФХ объекта управления
Определим алгоритм формирования переменной u1 для того, чтобы решить задачу стабилизации. Для этого выберем алгоритм формирования сигнала u1 в виде
u1 = −W1( p) y , |
(3.11) |
27
где W1( p) = b1( p) . Подставляя (3.9), (3.11) в уравнение (2.2) и разрешая его относи- a1( p)
тельно выходной переменной, найдем уравнение замкнутой системы
|
|
|
|
|
|
|
y =Wy ( p)u2 , |
|
|
|
|
|
|
|
|
|
|
где |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
Wy ( p) = |
W ( p)W0 ( p) |
|
= |
|
( p +1)b0 ( p)a1( p) |
|
|
. |
(3.12) |
||||||||
1 +W ( p)W ( p)W ( p) |
( p2 + p −1)a ( p)a ( p) +b( p)b ( p)b ( p) |
||||||||||||||||
0 |
1 |
|
|
|
0 |
1 |
|
|
|
|
|
0 |
1 |
|
|
||
Рассчитаем полиномы |
a0 ( p),a1( p),b0 ( p),b1( p) |
из условия |
|
|
|
|
|||||||||||
|
|
|
|
b( p)b0 ( p)a1( p) |
|
= |
|
|
1 |
|
, |
|
|
(3.13) |
|||
|
|
a( p)a ( p)a ( p) +b( p)b ( p)b ( p) |
a |
y |
( p) |
|
|
|
|||||||||
|
|
|
0 |
1 |
|
0 |
1 |
|
|
|
|
|
|
|
|
||
где ay ( p) - произвольный устойчивый полином степени n − m = 2 −1 =1. |
|
||||||||||||||||
Пусть полином ay ( p) = p +1, тогда |
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
a |
( p) =b |
( p) = ( p +1)n−1 = p +1 |
|
|
|
|
|
(3.14) |
||||||
|
|
|
1 |
0 |
|
|
|
|
|
|
|
|
|
|
|||
и, используя следующее тождество |
|
|
|
|
|
|
|
|
|
|
|
||||||
|
a( p)a ( p) +b( p)b ( p) =b( p)( p +1)2n−m−1 , |
|
|
|
(3.15) |
||||||||||||
|
|
|
0 |
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
найдем коэффициенты многочленов a0 ( p),b1( p) .
Тогда, подставляя соответствующие числовые значения в уравнение (3.15), получаем
( p2 + p −1)(a |
0,1 |
p + a |
0,0 |
) + ( p +1)(b |
p + b |
) = ( p +1)2 ( p +1) . |
(3.16) |
|
|
1,1 |
1,0 |
|
|
Откуда следует
a |
0,1 |
p3 + (a |
0,0 |
+ a |
0,1 |
+b |
|
) p2 |
+ (a |
0,0 |
− a |
0,1 |
+b |
+b |
) p + (b |
− a |
0,0 |
) = |
|
|
|
1,1 |
|
|
|
1,0 |
1,1 |
1,0 |
|
|
|||||||
|
|
|
|
|
|
|
= p3 +3 p2 +3 p +1. |
|
|
|
|
(3.17) |
||||||
Приравнивая члены при соответствующих степенях, получаем
a0,1 =1,
a0,0 + a0,1 +b1,1 =3 ,
a0,0 − a0,1 +b1,0 +b1,1 =3,
b1,0 − a0,0 =1,
28
откуда следует
a0,0 +b1,1 = 2 , b1,0 +b1,1 + a0,0 = 4 ,
b1,0 =1 + a0,0 .
Производя простые преобразования, находим коэффициенты полиномов a0 ( p) и b1( p) . В нашем случае получилось: a0,1 = a0,0 =1 и b1,1 =1, b1,0 = 2 . Таким образом передаточные функции регуляторов имеют вид:
W |
= |
s +1 |
и W (s) = |
s + 2 |
, |
||||
|
|
||||||||
0 |
|
s +1 |
1 |
|
|
s +1 |
|
||
|
|
|
|
|
|
|
|||
а передаточная функция замкнутой системы: |
|
|
|
||||||
|
|
Wy (s) = |
|
1 |
|
. |
|
|
|
|
|
|
s +1 |
|
|
||||
|
|
|
|
|
|
|
|
||
ЛАХ и ЛФХ замкнутой системы представлены на рис. 3.4. Структурная схема системы управления приведена на рис. 3.5. Результаты компьютерного моделирования представлены на рис. 3.6 иллюстрируют асимптотическую устойчивость системы. При постановке эксперимента были выбраны ненулевые начальные условия на интеграторе с выходом x1 .
3.4. Расчет передаточной функции регулятора |
|
Управляющее воздействие объектом (3.13) строится в форме: |
|
~ |
(3.18) |
u =W2 ( p)ε , |
где W2 ( p) = b2 ( p) . a2 ( p)
Для расчета передаточной функции регулятора воспользуемся методом, предполагающим построение желаемой передаточной функции замкнутой системы, удовлетворяющей выше перечисленным требованиям и синтез регулятора на соответствии между логарифмическими частотными характеристиками.
29
Рис. 3.4. ЛАХ и ЛФХ замкнутой системы
Рис. 3.5. Структурная схема системы управления
30
