MS_bak_220400_220700
.pdf37.Провести сравнительную оценку схем несвязанного и связанного регулирования. От чего зависит выбор той или иной схемы?
38.Назвать основные показатели качества регулирования. Для каких целей они определяются?
39.Как выполняется расчет показателей качества при моделировании с использованием дискретных моделей?
40.Как выполняется структурный синтез регуляторов и компенсаторов по показателям качества (в каскадной, комбинированной, несвязанной и связанной системах)?
41.Какие факторы вызывают дрейф динамических характеристик технологических объектов (приведите примеры)?
42.Для каких целей исследуется чувствительность систем регулирования к изменению параметров модели объекта?
43.Назвать основные отличия цифровых и непрерывных систем управления.
44.Пояснить формирование комплексного критерия (среднеквадратичной суммы ошибки регулирования и отклонения управляющей величины от установившегося значения) для выбора оптимального такта квантования сигналов.
45.Как скажется изменение длительности такта квантования на параметрах дискретных моделей объектов и регуляторов?
46.Как выполнить пересчет параметров дискретных моделей объектов и регуляторов при изменении такта квантования?
151
3.15. Варианты заданий к лабораторной работе “Идентификация дискретной динамической модели по имитационной модели объекта”
Цель работы:
изучение методов построения дискретных динамических моделей, используемых при синтезе цифрового управления, и идентификация параметров моделей объектов регулирования, описываемых конечно-разностными уравнениями.
Постановка задачи.
Для указанного варианта исходных данных (табл. 3.16) построить дискретную динамическую модель объекта регулирования, соответствующую апериодическому звену 1-го, 2-го или 3-го порядка. Используя МНК составить систему уравнений и решить ее относительно определяемых параметров дискретной модели для двух заданных структур (графа 11). При отсутствии экспериментальных данных реального объекта регулирования используется имитационная модель объекта, заданная в виде уравнения (3.2) (графы 4 – 8). Дополнительно заданы коэффициент помехи (графа 9) и длительность такта квантования (графа 10). Составить схему алгоритма и программу расчета параметров дискретной модели для двух заданных структур по данным имитационной модели (при наличии помехи и без нее), оценить адекватность модели по критерию Фишера (приложение А) и выбрать структуру модели с наибольшим критерием.
Содержание отчета:
1.Название и цель работы, номер варианта.
2.Математическая формулировка задачи, содержащая все формулы, выводы и пояснения для решения поставленной задачи.
3.Схема алгоритма расчета с комментариями.
4.Листинг программы.
5.Анализ результатов (приводятся численные значения параметров модели объекта для двух структур уравнений при идентификации по имитационной модели (при наличии помехи и без
152
нее), критерия Фишера, графики переходных процессов имитационной модели объекта и восстановленной модели (при наличии помехи и без нее)) с конкретными выводами.
Поставленную задачу можно представить следующей структурной схемой (рис.3.39).
Рис.3.39. Структурная схема задачи идентификации модели объекта по имитационной модели
Например, примем в качестве объекта регулирования имитационную модель в виде конечно-разностного уравнения второго порядка с чистым запаздыванием с наложением случайной помехи i на выходную переменную (см. вариант 2):
|
|
yi a1 yi 1 a2 yi 2 |
bui d 1 , |
(3.103) |
|||
|
|
y u y |
i |
p , |
|
(3.104) |
|
|
|
i |
|
i |
|
|
|
где y |
i |
- выход объекта без помехи; |
y u |
- измеряемое значение |
|||
|
|
|
|
i |
|
|
выходного сигнала с помехой i ; p - коэффициент помехи, оп-
ределяющий уровень помехи на каждом такте T0 квантования;
i 2random 1001 / 1000 1.
153
Варианты заданий
|
Технологический |
Канал объекта регулирования |
Пар. мо- |
№ |
объект (ТО) |
|
дели ТО |
|
|
|
|
вар |
|
|
Первая |
|
|
|
пост. |
|
|
|
врем. |
|
|
|
|
T1 |
|
|
|
|
|
1 |
|
2 |
3 |
4 |
|
|
|
|
|
1 |
Теплообменник |
Расход греющего пара - температура смеси |
4,21 |
|
|
|
|
на выходе теплообменника |
мин |
|
|
|
|
|
2 |
Колонна ректификации |
Расход теплоносителя в кипятильник - |
81 мин |
|
|
|
|
содержание компонента смеси на к.т. |
|
|
|
|
|
|
3 |
Полимеризатор |
Расход рассола в рубашку полимеризатора - |
6 мин |
|
|
|
|
температура смеси |
|
|
|
|
|
|
4 |
Колонна |
экстрактив- |
Расход флегмы - содержание компонента |
41,8 |
|
ной ректификации |
смеси в дистилляте |
мин |
|
|
|
|
|
|
5 |
Теплообменник |
Расход греющего пара - температура смеси |
8 мин |
|
|
|
|
на выходе теплообменника |
|
|
|
|
|
|
6 |
Колонна ректификации |
Расход теплоносителя в кипятильник - |
12 мин |
|
|
|
|
температура смеси в кубе колонны |
|
|
|
|
|
|
7 |
Паровой котел |
Расход топлива - температура в котле |
153,5 с |
|
|
|
|
|
|
8 |
Полимеризатор |
Расход рассола в рубашку полимеризатора - |
78 мин |
|
|
|
|
конверсия мономера |
|
|
|
|
|
|
9 |
Теплообменник |
Расход греющего пара - температура смеси |
4,21 |
|
|
|
|
на выходе теплообменника |
мин |
|
|
|
|
|
10 |
Колонна ректификации |
Расход теплоносителя в кипятильник - |
3,59 |
|
|
|
|
температура смеси в кубе колонны |
мин |
|
|
|
|
|
11 |
Теплообменник |
Расход греющего пара - температура смеси |
7,86 |
|
|
|
|
на выходе теплообменника |
мин |
|
|
|
|
|
12 |
Колонна |
экстрактив- |
Расход флегмы - содержание компонента |
11,7 |
|
ной ректификации |
смеси на в.к.т. |
мин |
|
|
|
|
|
|
13 |
Теплообменник |
Расход греющего пара - температура смеси |
8 мин |
|
|
|
|
на выходе теплообменника |
|
154
Таблица 3.16
к выполнению работы
|
|
Параметры модели ТО |
|
Такт |
Порядки |
Метод |
|
|
|
|
|
|
кван- |
уравне- |
решения |
|
|
|
|
|
|
|
|
Вторая |
Третья |
Коэффи- |
Время |
Коэф. |
това- |
ний вос- |
системы |
пост. |
пост. |
циент |
чистого |
поме- |
ния |
станав- |
лин. |
врем. |
врем. |
усиления |
запазд. |
хи |
|
ливаемых |
уравне- |
|
|
|
|
|
|
|
|
T2 |
T3 |
k |
|
p |
T0 |
моделей |
ний |
|
|
|
|
||||
5 |
6 |
7 |
8 |
9 |
10 |
11 |
12 |
|
|
|
|
|
|
|
|
- |
- |
12,1 |
2,3 мин |
0,61 |
0,1 |
1;2 |
Правило |
|
|
ºС/(т/ч) |
|
|
мин |
|
Крамера |
|
|
|
|
|
|
|
|
27 |
- |
-5 %/(т/ч) |
27 мин |
0,25 |
1 мин |
1;3 |
Метод |
мин |
|
|
|
|
|
|
Гаусса |
|
|
|
|
|
|
|
|
2 мин |
- |
4 ºС/(м3/ч) |
2 мин |
0,2 |
0,1 |
2;3 |
Правило |
|
|
|
|
|
мин |
|
Крамера |
|
|
|
|
|
|
|
|
15,3 |
10,2 |
-5,81 |
16 мин |
0,29 |
2 мин |
1;2 |
Метод |
мин |
мин |
%/(т/ч) |
|
|
|
|
Гаусса |
|
|
|
|
|
|
|
|
2 мин |
- |
1 ºС/(т/ч) |
3,8 мин |
0,05 |
0,2 |
1;3 |
Правило |
|
|
|
|
|
мин |
|
Крамера |
|
|
|
|
|
|
|
|
2,44 |
- |
4,5 |
3 мин |
0,225 |
0,25 |
2;3 |
Метод |
мин |
|
ºС/(т/ч) |
|
|
мин |
|
Гаусса |
|
|
|
|
|
|
|
|
24 с |
15 с |
1,77 |
0 с |
0,089 |
10 с |
1;2 |
Правило |
|
|
ºС/( м3/ч) |
|
|
|
|
Крамера |
77 |
- |
1,5 |
6 мин |
0,075 |
2 мин |
1;3 |
Метод |
мин |
|
%/(м3/ч) |
|
|
|
|
Гаусса |
3,19 |
- |
12,1 |
0,5 мин |
0,61 |
0,1 |
2;3 |
Правило |
мин |
|
ºС/(м3/ч) |
|
|
мин |
|
Крамера |
2,8 |
- |
7,3 |
4 мин |
0,37 |
0,5 |
1;2 |
Метод |
мин |
|
ºС/(т/ч) |
|
|
мин |
|
Гаусса |
|
|
|
|
|
|
|
|
- |
- |
3,05 |
2 мин |
0,15 |
0,25 |
1;3 |
Правило |
|
|
ºС/(т/ч) |
|
|
мин |
|
Крамера |
|
|
|
|
|
|
|
|
5,48 |
4,9 |
-4,33 |
20 мин |
0,22 |
2 мин |
2;3 |
Метод |
мин |
мин |
%/(т/ч) |
|
|
|
|
Гаусса |
|
|
|
|
|
|
|
|
2 мин |
- |
1 ºС/(м3/ч) |
2,1 мин |
0,05 |
0,1 |
1;2 |
Правило |
|
|
|
|
|
мин |
|
Крамера |
|
|
|
|
|
|
|
|
155
|
Технологический |
Канал объекта регулирования |
Пар. мо- |
№ |
объект (ТО) |
|
дели ТО |
|
|
|
|
вар |
|
|
Первая |
|
|
|
пост. |
|
|
|
врем. |
|
|
|
T1 |
|
|
|
|
1 |
2 |
3 |
4 |
|
|
|
|
14 |
Колонна ректификации |
Расход теплоносителя в кипятильник - |
73,84 |
|
|
температура смеси в кубе колонны |
мин |
|
|
|
|
15 |
Полимеризатор |
Расход катализатора - температура в реак- |
41,72 |
|
|
торе |
мин |
|
|
|
|
16 |
Теплообменник |
Расход греющего пара - температура смеси |
4,21 |
|
|
на выходе теплообменника |
мин |
|
|
|
|
17 |
Колонна ректификации |
Расход теплоносителя в кипятильник - |
73,84 |
|
|
температура смеси в кубе колонны |
мин |
|
|
|
|
18 |
Колонна ректификации |
Расход теплоносителя в кипятильник - |
35,18 |
|
|
температура смеси в кубе колонны |
мин |
|
|
|
|
19 |
Полимеризатор |
Расход катализатора - вязкость выходного |
45,45 |
|
|
продукта |
мин |
|
|
|
|
20 |
Колонна ректификации |
Расход теплоносителя в кипятильник - |
65 мин |
|
|
содержание компонента смеси в кубе кол. |
|
|
|
|
|
21 |
Теплообменник |
Расход греющего пара - температура смеси |
4,21 |
|
|
на выходе теплообменника |
мин |
|
|
|
|
22 |
Колонна ректификации |
Расход теплоносителя в кипятильник - |
73,84 |
|
|
температура смеси в кубе колонны |
мин |
|
|
|
|
23 |
Колонна ректификации |
Расход теплоносителя в кипятильник - |
10,12 |
|
|
температура смеси в кубе колонны |
мин |
|
|
|
|
24 |
Колонна экстрактив- |
Расход флегмы - содержание компонента |
41 мин |
|
ной ректификации |
смеси в дистилляте |
|
|
|
|
|
25 |
Колонна ректификации |
Расход теплоносителя в кипятильник - |
12 мин |
|
|
температура смеси в кубе колонны |
|
156
Окончание табл. 3.16
|
|
Параметры модели ТО |
|
Такт |
Порядки |
Метод |
|
|
|
|
|
|
кван- |
уравне- |
решения |
|
|
|
|
|
|
|
|
Вторая |
Третья |
Коэффи- |
Время |
Коэф. |
това- |
ний вос- |
системы |
пост. |
пост. |
циент |
чистого |
поме- |
ния |
станав- |
лин. |
врем. |
врем. |
усиления |
запазд. |
хи |
|
ливаемых |
уравне- |
|
|
|
|
|
|
|
|
T2 |
T3 |
k |
|
p |
T0 |
моделей |
ний |
|
|
|
|
||||
5 |
6 |
7 |
8 |
9 |
10 |
11 |
12 |
|
|
|
|
|
|
|
|
15,13 |
- |
0,08 |
10 мин |
0,04 |
1 мин |
1;3 |
Метод |
мин |
|
ºС/(т/ч) |
|
|
|
|
Гаусса |
|
|
|
|
|
|
|
|
17,28 |
- |
-1,03 |
8 мин |
0,052 |
1 мин |
2;3 |
Правило |
мин |
|
ºС/(м3/ч) |
|
|
|
|
Крамера |
3,19 |
- |
12,1 |
5 мин |
0,61 |
0,5 |
1;2 |
Метод |
мин |
|
ºС/(м3/ч) |
|
|
мин |
|
Гаусса |
15,13 |
- |
0,08 |
15 мин |
0,04 |
1 мин |
1;3 |
Правило |
мин |
|
ºС/(т/ч) |
|
|
|
|
Крамера |
|
|
|
|
|
|
|
|
2,88 |
- |
1,04 |
9 мин |
0,052 |
1 мин |
2;3 |
Метод |
мин |
|
ºС/(т/ч) |
|
|
|
|
Гаусса |
|
|
|
|
|
|
|
|
20,35 |
- |
0,4 |
12 мин |
0,02 |
1 мин |
1;2 |
Правило |
мин |
|
г/см3/(т/ч) |
|
|
|
|
Крамера |
40 |
15 мин |
-0,9 %/(т/ч) |
10 мин |
0,045 |
2 мин |
1;3 |
Метод |
мин |
|
|
|
|
|
|
Гаусса |
|
|
|
|
|
|
|
|
- |
- |
12,1 |
7 мин |
0,61 |
0,5 |
2;3 |
Правило |
|
|
ºС/(м3/ч) |
|
|
мин |
|
Крамера |
15,13 |
- |
0,08 |
8 мин |
0,004 |
1 мин |
1;2 |
Метод |
мин |
|
ºС/(т/ч) |
|
|
|
|
Гаусса |
|
|
|
|
|
|
|
|
2,77 |
- |
3,12 |
8 мин |
0,16 |
0,5 |
1;3 |
Правило |
мин |
|
ºС/(т/ч) |
|
|
мин |
|
Крамера |
|
|
|
|
|
|
|
|
19,17 |
0,698 |
-0,13 |
20 мин |
0,007 |
2 мин |
2;3 |
Метод |
мин |
мин |
%/(т/ч) |
|
|
|
|
Гаусса |
|
|
|
|
|
|
|
|
- |
- |
2 ºС /(т/ч) |
2 мин |
0,2 |
0,2 |
1;2 |
Правило |
|
|
|
|
|
мин |
|
Крамера |
|
|
|
|
|
|
|
|
157
Расчет значения случайной помехи i реализуется стан-
дартным оператором – функцией int random int num в язы-
ке Cu (генератор случайных чисел). Помеха принимает зна-
чения в пределах от минус 1 до плюс 1 случайным образом.
С помощью имитационной модели рассчитываются значения точек переходного процесса объекта регулирования с тактом квантования T0 . Для расчета кривой разгона по имитационной модели необходимо знать параметры a1, a2 , b, d (определяются пересчетом по исходным данным (графы 4 – 8), такт квантования T0 , коэффициент помехи p и задаться нулевыми начальными условиями по входу и выходу объекта, а также величиной входного ступенчатого воздействия u .
Для расчета параметров восстанавливаемых моделей двух структур (графа 11) составляются системы уравнений и решаются указанным методом (графа 12) (см. пункт 3.2).
Для оценки адекватности моделей по критерию Фишера используется формула (3.28).
158
3.16. Варианты заданий к лабораторной работе “Оптимизация настроек цифровых регуляторов в составе одноконтурной схемы”
Цель работы:
овладение методикой синтеза и исследования одноконтурной цифровой системы управления (ЦСУ).
Постановка задачи.
Для указанного варианта исходных данных (табл. 3.17) провести оптимизацию градиентным методом настроек цифровых регуляторов, соответствующих заданным непрерывным законам регулирования (графа 10). В качестве исходных данных заданы параметры дискретной динамической модели объекта (графы 4 – 8) и длительность такта квантования (графа 9). Составить схему алгоритма и программу расчета настроек цифровых регуляторов для двух заданных структур, построить динамические характеристики синтезированных оптимальных ЦСУ (по заданию и возмущению) и оценить эффективность их работы по показателям качества (3.89) – (3.93).
Содержание отчета:
1.Название и цель работы, номер варианта.
2.Математическая формулировка задачи, содержащая все формулы, выводы и пояснения для решения поставленной задачи.
3.Схема алгоритма расчета с комментариями.
4.Листинг программы.
5.Анализ результатов (приводятся численные значения оптимальных настроек цифровых регуляторов для двух законов регулирования, графики динамических характеристик синтезированных оптимальных ЦСУ, значения показателей качества) с конкретными выводами.
Для составления математического описания задачи оптимизации настроек регуляторов см. пункт 3.5. Формулы расчета показателей качества приведены в пункте 3.11.
159
Таблица 3.17
Варианты заданий для выполнения лабораторной работы
№ |
Объект |
Канал ОР |
|
Параметры модели объекта |
|
|
Законы |
||
вар. |
регулиро- |
|
|
|
|
|
|
|
регули- |
|
вания (ОР) |
|
a1 |
a2 |
a3 |
b |
d |
T0 |
рования |
1 |
2 |
3 |
4 |
5 |
6 |
7 |
8 |
9 |
10 |
1 |
Тепло- |
Расход греющего пара – |
0,976247 |
– |
– |
0,287411 |
6 |
0,1 |
ПИ, |
|
обменник |
температура смеси |
|
|
|
С/(т/ч) |
|
мин |
ПД |
2 |
Колонна |
Расход гр. пара в кипят. – |
1,952505 |
-0,952941 |
– |
-0,002179 |
10 |
1,0 |
П, |
|
ректиф. |
конц. комп. на к.т. колонны |
|
|
|
%/(т/ч) |
|
мин |
ПИД |
3 |
Полиме- |
Расход рассола в рубашку – |
1,936719 |
-0,9375 |
– |
0,003125 |
5 |
0,1 |
И, |
|
ризатор |
температура смеси |
|
|
|
С/(м3/ч) |
|
мин |
ПИ |
4 |
Колонна |
Расход флегмы – конц. |
2,676247 |
-2,383372 |
0,706259 |
-0,005032 |
4 |
2,0 |
ПИ, |
|
экстр.рект. |
комп. в дистилляте |
|
|
|
%/(т/ч) |
|
мин |
ПД |
5 |
Тепло- |
Расход греющего пара – |
1,886667 |
-0,888889 |
– |
0,002222 |
3 |
0,2 |
П, |
|
обменник |
температура смеси |
|
|
|
С/(т/ч) |
|
мин |
ПД |
6 |
Колонна |
Расход гр. пара в кипят. – |
1,888339 |
-0,890241 |
– |
0,008551 |
6 |
0,25 |
И, |
|
ректиф. |
конц. комп. в кубе колонны |
|
|
|
С/(т/ч) |
|
мин |
ПИД |
7 |
Паровой |
Расход топлива – темпера- |
2,253742 |
-1,661497 |
0,400507 |
0,012828 |
0 |
10,0 с |
ПИ, |
|
котел |
тура в котле |
|
|
|
С/(м3/ч) |
|
|
ПИД |
|
160 |
|
|
|
|
|
|
|
|