
MS_bak_220400_220700
.pdf
Рис. 3.13. Структурная схема каскадной системы управления:
WP1 ,WP2 ,WO1 ,WO2 - передаточные функции регуляторов и каналов управления внутреннего и внешнего контуров; y 3 - задающее воздейст-
вие; e1 ,e2 - величины рассогласования; u1 , u2 - управляющие воздей-
ствия; y1 , y2 - выходы объекта
Расчет цифровых регуляторов каскадной системы может быть выполнен двумя методами. Первый - метод декомпозиции, когда с помощью численных методов оптимизации проводится расчет регулятора внутреннего контура в составе одноконтурной схемы (см. пункт 3.5), а затем расчет регулятора внешнего контура по системе уравнений, описывающей динамику каждого элемента каскадной схемы.
Второй - метод свертки, при котором после расчета настроек регулятора внутреннего контура проводится свертка каскадной схемы в одноконтурную с эквивалентным объектом с помощью оператора сдвига z для расчета регулятора внешнего контура.
Рассмотрим метод декомпозиции. Для оптимизации настроек q02 ,q12 регулятора внешнего контура по критерию (3.44)
необходимо составить уравнения (3.70) для расчета динамических характеристик системы и квазианалитические выражения (3.71) для расчета производных критерия в диапазоне тактов
i mc, N :
u2 |
u2 |
q2 |
y3 |
y |
2 |
q2 |
y |
3 |
|
y2 |
|
, |
(3.70) |
||
i |
i 1 |
0 |
|
i |
|
i |
|
1 |
|
i 1 |
i 1 |
|
|
||
u1 |
u1 |
q1 |
u2 |
y1 |
q1 u2 |
|
y1 |
, |
|
||||||
i |
i 1 |
0 |
|
i |
i |
|
|
1 |
i 1 |
|
i 1 |
|
|
|
101

|
|
|
|
y1 |
|
|
|
a1 y1 |
a1 y1 |
|
|
|
|
b1u1 , |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||
|
|
|
|
|
i d1 1 |
|
|
|
1 i d1 |
|
|
|
2 i d1 1 |
|
|
|
|
|
i |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||
|
|
|
|
y2 |
|
|
|
|
|
|
a |
2 y |
2 |
|
|
|
|
|
a2 y2 |
|
|
|
|
|
|
b2 y1 |
|
|
, |
|
|
|
|
||||||||||||||||||||||
|
|
|
|
|
i d1 d 2 2 |
|
|
|
|
1 i d1 d 2 1 |
|
|
|
2 i d1 d 2 |
|
|
|
|
|
|
|
i d1 1 |
|
|
|
|
|
||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
i mc, N . |
||||||
|
|
Расчет динамики проводится при следующих начальных |
|||||||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
y3 |
0 |
|
при |
|
|
i mc |
|
|
|
y 3 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||
условиях: |
|
|
|
3 |
|
|
|
|
|
|
|
i mc |
, |
|
|
- величина задающего |
|||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
i |
|
|
|
y |
|
|
при |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||
воздействия; |
|
|
|
|
|
u2 |
|
0, |
i 1,mc 1; |
|
|
|
u1 0, |
i 1,mc 1; |
|||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
i |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
i |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
y1 |
0, |
|
i 1,mc d1 |
; y2 |
0, |
|
i 1,mc d1 d2 1 . |
||||||||||||||||||||||||||||||||||||||||||||||||
i |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
i |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Дифференцируя уравнения (2.22), получим: |
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||||||||
|
|
ui2 |
|
|
ui2 1 |
|
|
y3 |
y |
2 |
q 2 |
yi2 |
q |
2 |
|
|
yi2 1 |
|
, |
|
|
|
|
|
|
(3.71) |
|
|
|||||||||||||||||||||||||||
|
|
q 2 |
|
q 2 |
|
|
|
|
|
|
|
|
|
q 2 |
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
i |
|
|
i |
|
|
|
0 q 2 |
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||
|
0 |
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
ui2 |
|
|
ui2 1 |
|
q 2 yi2 |
|
y |
3 |
y 2 |
q 2 |
yi2 1 |
|
, |
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||
|
|
q 2 |
|
q 2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
0 q 2 |
|
|
|
|
i 1 |
|
|
|
i 1 |
|
|
|
|
|
1 q 2 |
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||
|
1 |
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
||||||
|
1 |
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
1 |
|
|
||||||||||||
|
|
ui |
|
|
|
|
ui 1 |
q |
1 |
|
ui |
|
|
yi |
|
q1 |
|
|
ui 1 |
|
|
|
yi 1 |
, |
|||||||||||||||||||||||||||||||
|
|
q 2 |
|
|
|
|
|
q 2 |
q 2 |
|
|
q 2 |
|
|
q 2 |
||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
q 2 |
|
|
|
|
|
0 |
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||
|
0 |
|
|
0 |
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|||||||||
|
1 |
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
1 |
|
|
||||||||||||
|
|
ui |
|
|
|
|
ui 1 |
q |
1 |
|
ui |
|
|
yi |
|
q1 |
|
|
ui 1 |
|
|
|
yi 1 |
, |
|||||||||||||||||||||||||||||||
|
|
q 2 |
|
|
|
|
|
|
|
q 2 |
q 2 |
|
|
q 2 |
|
|
q 2 |
||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
q 2 |
|
|
|
|
|
0 |
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||
|
1 |
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|||||||
|
|
y1 |
|
|
|
|
|
|
|
|
|
|
|
y |
1 |
|
|
|
|
|
|
|
y1 |
|
|
|
|
|
|
|
|
u1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
i d 1 1 |
a |
1 |
|
|
|
i d 1 |
a |
1 |
|
|
i d 1 1 |
b1 |
|
|
|
|
i |
, |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||
|
|
|
1 |
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||
|
|
q 2 |
|
|
|
|
|
|
|
q 2 |
|
|
|
|
|
q |
2 |
|
|
|
|
|
|
q 2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
y1 |
|
|
|
|
|
|
|
|
|
|
|
y |
1 |
|
|
|
|
|
|
|
y1 |
|
|
|
|
|
|
|
|
u1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
i d 1 1 |
a |
1 |
|
|
|
i d 1 |
a |
1 |
|
|
i d 1 1 |
b1 |
|
|
|
|
i |
|
, |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||
|
|
q 2 |
|
|
|
|
|
1 |
|
|
q 2 |
|
|
|
|
2 |
|
q |
2 |
|
|
|
|
|
|
q 2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
y 2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
y 2 |
|
|
|
|
|
|
|
|
y |
2 |
|
|
|
|
|
|
|
|
|
|
|
|
y |
1 |
|
|
|
|
|
|
||||||||
|
|
i d1 d 2 2 |
a 2 |
|
|
|
i d 1 d 2 1 |
a 2 |
|
|
|
i d 1 d 2 |
b2 |
|
|
|
i d 1 1 |
, |
|||||||||||||||||||||||||||||||||||||
|
|
q 2 |
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
q 2 |
|
|
|
|
|
2 |
|
|
|
q 2 |
|
|
|
|
|
|
|
|
|
|
|
|
q 2 |
|
|
|
|
|
||||||||||
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|||
|
|
y 2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
y 2 |
|
|
|
|
|
|
|
|
y |
2 |
|
|
|
|
|
|
|
|
|
|
|
|
y |
1 |
|
|
|
|
|
|
||||||||
|
|
i d1 d 2 2 |
a 2 |
|
|
|
i d 1 d 2 1 |
a 2 |
|
|
|
i d 1 d 2 |
b2 |
|
|
|
i d 1 1 |
, |
|||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||||||
|
|
q 2 |
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
q 2 |
|
|
|
|
|
2 |
|
|
|
q 2 |
|
|
|
|
|
|
|
|
|
|
|
|
q 2 |
|
|
|
|
|
||||||||||
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
i mc, N .
102
Начальные условия:
|
|
u |
2 |
0, |
|
|
|
|
i |
|
|||
|
|
q |
2 |
|
|
|
|
|
|
0 |
|
|
|
|
|
u1 |
0, |
|
||
|
|
|
i |
|
||
|
|
q |
2 |
|
||
|
|
|
|
|||
|
|
|
0 |
|
|
|
|
y1 |
0, |
|
|||
|
|
i |
|
|
|
|
|
q 2 |
|
||||
|
|
|
||||
|
0 |
|
|
|
|
|
y2 |
0, |
y |
2 |
|||
i |
i |
|||||
q 2 |
|
|
|
|
q |
2 |
0 |
|
|
|
|
1 |
|
u2 |
|
|
|
|
|
|
|
||
0, |
i 1,mc 1 ; |
|||||||||
|
i |
|||||||||
|
q 2 |
|
|
|
|
|
|
|
||
1 |
|
|
|
|
|
|
|
|||
|
u1 |
|
|
|
|
|
|
|
||
|
0, |
i 1,mc 1 ; |
||||||||
|
i |
|||||||||
|
q 2 |
|||||||||
|
|
|
|
|
|
|
|
|||
1 |
|
|
|
|
|
|
|
|||
y1 |
|
|
|
|
|
|
|
|||
0, |
i 1,mc d1; |
|||||||||
|
i |
|
||||||||
q 2 |
||||||||||
|
|
|
|
|
|
|
||||
1 |
|
|
|
|
|
|
|
|
||
|
|
|
|
|||||||
0, |
i 1,mc d1 d 2 1 . |
Выполнение второго этапа оптимизации в каскадной системе полностью совпадает с расчетом регулятора одноконтурной ЦСУ (см. пункт 3.5.).
После определения оптимальных настроек (табл. 3.8) проводится исследование динамических характеристик каскадной ЦСУ (рис. 3.14).
|
|
|
|
|
Таблица 3.8 |
|
|
Параметры оптимальной каскадной ЦСУ |
|||||
|
|
|
|
|||
Параметры моделей |
Оптимальные настройки |
Такт |
|
|||
объекта управления |
цифровых регуляторов |
квантования |
|
|||
|
1 |
|
1 |
|
сигналов |
|
|
WO |
|
WP |
|
|
|
a1 |
1,501454 |
q1 |
0,661289 (т/ч)/ С |
T0 1,0 мин |
|
|
1 |
|
0 |
|
|
|
|
a1 |
0,571428 |
q1 |
0,581424 |
(т/ч)/ С |
|
|
2 |
|
1 |
|
|
|
|
b1 0,0907 С/(т/ч) |
|
|
|
|
|
|
d1 4 |
|
|
|
|
|
|
|
W 2 |
|
W 2 |
|
|
|
|
O |
|
P |
|
|
|
a2 |
1,856459 |
q2 |
5,028216 |
С/% |
|
|
1 |
|
0 |
|
|
|
|
a2 |
0,861244 |
q2 |
4,909592 С/% |
|
|
|
2 |
|
1 |
|
|
|
|
b2 |
0,000598 %/ С |
|
|
|
|
|
d2 30 |
|
|
|
|
|
|
|
|
|
|
|
103 |

Рис. 3.14. Динамические характеристики оптимальной каскадной ЦСУ температуры на верхней контрольной тарелке колонны ректификации с коррекцией по концентрации изоамиленов в дистилляте:
104
а) y1 - температура; u2 - корректирующее воздействие по температуре;
u1 - расход флегмы; б) y 3 - задание по концентрации; y 2 - концентрация
Рассмотрим метод свертки. После идентификации дискретных моделей объекта по каналам управления и настройки регулятора внутреннего контура проводится свертка структурной схемы управления и определяется передаточная функция эквивалентного объекта внешнего контура, в которой отсутствует внутренняя связь между выходом регулятора и входом объекта внутреннего контура:
экв |
|
y |
2 |
|
|
|
1 |
1 |
|
|
2 |
|
z |
|
z |
|
|
WP |
z WO |
z |
|
z . (3.72) |
|||
WO |
|
|
|
|
WO |
|||||||
u |
2 |
z |
1 WP1 z WO1 z |
Для этого проводится переход прямым преобразованием от разностных уравнений (3.70) к соответствующим передаточным функциям. Переход осуществляется путем описания входных и выходных величин системы оператором сдвига z :
|
1 |
|
|
u1 z |
q01 q11z 1 |
|
|
|
Q1 z |
|
||||||||||
WP z |
|
|
|
|
|
|
|
|
|
|
, |
|
|
(3.73) |
||||||
e1 z |
1 z 1 |
|
P1 z |
|||||||||||||||||
WO1 z |
|
y1 z |
|
|
|
b1z d 1 1 |
|
|
|
|
|
B1 z |
|
|||||||
|
|
|
|
|
|
, |
|
|||||||||||||
u1 z |
1 a11z 1 a21 z 2 |
|
A1 z |
|
||||||||||||||||
WO2 z |
|
|
y 2 z |
|
|
|
|
b2 z d 2 1 |
|
|
|
|
|
B2 z |
|
|||||
|
|
|
|
|
|
, |
|
|||||||||||||
|
y1 z |
1 a12 z 1 a22 z 2 |
A2 z |
|
где Q1 z , P1 z , B1 z , A1 z , B2 z , A2 z - полиномы.
Подставляем передаточные функции объекта внутреннего, внешнего контура и регулятора внутреннего контура в эквивалентную передаточную функцию (3.72):
105

|
|
|
|
Q1 z B1 z |
|
|
|
|
||||||||
|
y 2 z |
|
|
|
|
|
|
|
|
|
|
|
B2 z |
|||
WOэкв z |
|
|
P1 z |
A1 z |
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
||||||
u2 z |
1 |
Q1 |
z B1 |
z |
|
A2 z |
||||||||||
|
|
|
|
|
|
|
|
(3.74) |
||||||||
|
|
|
P1 |
z |
A1 |
z |
P1 z A1 z Q1 z B1 z A2 z .
После подстановки полиномов из формул (3.73) и проведения аналитических преобразований получим:
W экв z |
y 2 z |
|
|
|
|
|
|
|
1z d 1 d 2 2 |
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
O |
|
|
u2 z |
1 z 1 |
z 2 |
z 3 |
z 4 |
|
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
2 |
|
3 |
|
|
4 |
|
(3.75) |
||
|
|
|
|
|
|
|
|
|
|
|
|
|
z d 1 d 2 3 |
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
, |
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
z 5 |
z d 1 1 |
z d 1 2 |
|
z d 1 3 z d 1 4 |
|
|
|
||||||||||||||||||||
|
|
|
|
5 |
|
6 |
|
|
|
|
|
|
7 |
|
|
|
|
8 |
|
|
9 |
|
|
|
|
|
|
|
где |
1 |
q1b1b2 ; |
|
2 |
q1b1b2 ; |
|
|
a1 |
a2 1; |
|
|
|
|
|
|
|||||||||||||
|
|
|
0 |
|
|
|
|
|
1 |
|
|
1 |
|
1 |
|
1 |
|
|
|
|
|
|
|
|||||
a1 |
a1 |
a2 a1 |
a1 |
; a1 |
a2 a1 |
a2 a1 |
a2 a1 |
a2 |
; |
|||||||||||||||||||
2 |
|
|
2 |
1 |
|
1 |
1 |
|
|
2 |
|
|
3 |
2 |
|
1 |
2 |
1 |
1 |
|
2 |
1 |
2 |
|
||||
a2 a1 a2 a1 |
a2 a1 |
; a |
2 a1 |
; q1b1 ; |
|
|
|
|
|
|
||||||||||||||||||
4 |
|
|
1 |
2 |
2 |
2 |
|
|
2 |
1 |
|
|
5 |
|
2 |
2 |
|
6 |
0 |
|
|
|
|
|
|
|
||
q1 |
a2 q1b1 ; a2 q1 a2 q1b1 ; |
a2 q1 . |
|
|
|
|
||||||||||||||||||||||
7 |
|
|
1 |
1 |
0 |
|
|
|
8 |
|
|
|
|
1 |
1 |
|
|
2 |
0 |
|
9 |
|
2 |
1 |
|
|
|
|
Используя обратное преобразование z , находим разностное уравнение эквивалентного объекта:
y 2 |
y 2 |
|
y 2 |
|
y 2 |
|
||||||
i d 1 d 2 2 |
|
|
1 i d 1 d 2 1 |
2 |
|
i d 1 d 2 |
|
|
3 i d 1 d 2 1 |
|
||
y 2 |
|
|
y 2 |
|
|
y 2 |
|
|
y 2 |
(3.76) |
||
4 i d 1 d 2 2 |
5 i d 1 d 2 3 |
|
|
6 i d 2 1 |
7 i d 2 |
|
||||||
y 2 |
|
, y 2 |
u2 |
u2 |
. |
|
|
|||||
8 i d 2 1 |
|
9 i d 2 2 |
|
1 i |
|
2 i 1 |
|
|
|
106

После свертки структурной схемы описание системы сводится к двум уравнениям (3.77) (рис. 3.15) - уравнению регулятора внешнего контура (первое уравнение системы (3.70)) и уравнению эквивалентного объекта (3.76):
u2 u2 |
q2 |
y3 |
y2 |
q2 |
y3 |
y2 |
, |
(3.77) |
||||||
i |
i 1 |
|
0 |
i |
i |
|
|
1 |
i 1 |
|
i 1 |
|
|
|
y 2 |
|
y 2 |
|
|
y 2 |
|
y 2 |
|
||||||
i d 1 d 2 2 |
|
1 i d 1 d 2 1 |
2 |
i d 1 d 2 |
|
|
3 i d 1 d 2 1 |
|
||||||
y 2 |
|
|
y 2 |
|
|
y 2 |
|
|
y 2 |
|
||||
4 |
i d 1 d 2 2 |
|
5 i d 1 d 2 3 |
|
6 |
i d 2 1 |
7 i d 2 |
|
||||||
y 2 |
|
, y 2 |
u |
2 u2 |
. |
|
|
|||||||
8 |
i d 2 1 |
|
9 i d 2 2 |
|
1 i |
2 i 1 |
|
|
|
Рис. 3.15. Структурная схема системы с эквивалентным объектом
Оптимизация настроек q02 ,q12 регулятора внешнего конту-
ра по критерию минимум интегрально-квадратичной ошибки (3.44) осуществляется с использованием уравнений (3.77) также как и в одноконтурной схеме (см. пункт 3.5).
3.8. Синтез комбинированной системы цифрового управления по критерию минимум интегральноквадратичной ошибки и из условия инвариантности
Комбинированные системы (рис. 3.16, 3.17) применяются в тех случаях, когда управляемый параметр технологического объекта существенно зависит не только от управляющего воздействия, но и от контролируемого возмущения.
107

Рис. 3.16. Пример комбинированной системы управления температурой на контрольной тарелке (к.т.) колонны ректификации при контролируемом возмущении по расходу экстрагента (растворителя): ТЕ, FE – датчики температуры, расхода; ТС – регулятор температуры (контур управления); FC – компенсатор расхода (контур возмущения)
Рис. 3.17. Структурная схема комбинированной системы управления:
WP ,WK ,WO ,WOx - передаточные функции компенсатора, регулятора и каналов управления, возмущения; y 3 - задающее воздействие; x - контролируемое возмущение; e - величина рассогласования; u,uk - управляющие воздействия от регулятора и компенсатора; u - суммарное управляющее воздействие; y - измеряемый выход объекта
108
Синтез цифровой комбинированной системы управления (рис. 3.17) заключается в определении структур и настроек основного регулятора и компенсатора возмущения. Расчет цифрового регулятора проводиться по основному каналу в составе одноконтурной схемы (см. пункт 3.5). Синтез цифрового компенсатора выполняется одним из численных методов оптимизации в составе комбинированной схемы или аналитически из условия инвариантности.
Для численной оптимизации настроек компенсатора сначала выбирается его структура и составляется модель динамики
комбинированной ЦСУ: |
|
|
|
|
|
|
|
|
|
|
|||||||||||
u |
|
u |
q |
0 |
y3 y |
|
q |
y3 |
y |
i 1 |
, |
(3.78) |
|||||||||
i |
|
|
i 1 |
|
|
|
|
|
i |
|
i |
|
|
1 |
i 1 |
|
|
|
|||
u |
k i |
|
q x |
x |
i |
q x |
x |
i 1 |
, |
|
|
|
|
|
|||||||
|
|
|
0 |
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|||
u i |
|
ui |
uk i , |
|
|
|
|
|
|
|
|
|
|
|
|||||||
yi d 1 a1 yi d a2 yi d 1 bu i , |
|
|
|
|
|||||||||||||||||
y x |
|
|
a x y x |
|
bx x |
i |
, |
|
|
|
|
|
|
||||||||
|
i dx 1 |
|
|
1 |
|
|
i dx |
|
|
|
|
|
|
|
|
|
|
||||
y |
|
|
y |
|
|
|
|
y x |
|
|
, |
|
|
|
|
|
|
|
|||
|
|
i d 1 |
|
|
|
|
|
|
|
|
|
||||||||||
|
i d 1 |
|
|
|
i d 1 |
|
|
|
|
|
|
|
|
где q0 ,q1 - оптимальные настройки основного регулятора, определяемые по модели основного канала (см. пункт 3.5); q0x ,q1x - искомые настройки цифрового компенсатора 1-го порядка без интегральной составляющей (табл. 2.1); a1 ,a2 ,b,d и a1x ,b x ,dx -
параметры моделей объекта.
Расчет проводится при нулевых начальных условиях:
|
|
0 |
при i mc |
|
|
0 |
при |
i mc |
y3 |
|
|
|
при i mc , |
||||
|
|
при i mc |
; x |
|
x |
|||
i |
y3 |
|
i |
|
|
i mc |
||
|
|
|
|
|
|
0 |
при |
где y3 , x - величины входных воздействий (ступенчатого и им-
пульсного); mc - переменная, принимающая наибольшее значение из порядков n и nx ( mc =2);
109

|
|
|
|
|
|
|
|
|
|
|
|
|
|
ui 0, uk |
0, u |
0, |
i 1,mc 1; |
||||||
|
|
|
|
i |
i |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
y |
i |
0, |
i 1,mc d ; |
y x 0, |
i 1,mc dx ; |
||||||
|
|
|
|
|
i |
|
|
|
|
|
yi 0, i 1,mc d ,
d - наименьшее число тактов чистого запаздывания из d и dx . Оптимизация настроек q0x ,q1x выполняется по критерию
минимум интегрально-квадратичной ошибки (3.79) между заданием y 3 и суммарным выходом объекта y методом покоординатного спуска (см. пункт 3.5):
N |
|
S yi3 yi 2 . |
(3.79) |
i mc |
|
После оптимизации настроек регулятора и компенсатора (табл. 3.9) проводится расчет динамических характеристик комбинированной ЦСУ (рис. 3.18).
|
|
|
|
Таблица 3.9 |
|
|
Параметры оптимальной комбинированной ЦСУ |
||||
|
|
|
|
||
Параметры моделей объ- |
Оптимальные настройки |
Такт |
|
||
|
екта управления |
регулятора и компенсатора |
квантова- |
|
|
|
|
|
|
ния |
|
|
WO |
|
WP |
сигналов |
|
|
|
|
|
||
a1 |
1,501454 |
q0 0,661289 (т/ч)/ С |
T0 1,0 |
|
|
a2 0,571428 |
q1 |
0,581424 (т/ч)/ С |
мин |
|
|
b 0,0907 С/(т/ч) |
|
|
|
|
|
d 4 |
|
|
|
|
|
|
W x |
|
W |
|
|
|
O |
|
K |
|
|
a x |
0,665132 |
q x |
0,207585 (т/ч)/ С |
|
|
1 |
|
0 |
|
|
|
bx |
0,36092 С/(т/ч) |
q x |
0,648909 (т/ч)/ С |
|
|
|
|
1 |
|
|
|
dx 2 |
|
|
|
|
110