
книги из ГПНТБ / Юсупбеков Н.Р. Автоматизация технологических процессов производства растительных масел
.pdf
|
Уъ — fb (.x w -M.0; |
|
|
(224) |
|
|
Уо~/в(-^12« -Из)! |
|
|
(225) |
|
Уравнения (220 —225) |
представлены |
в виде ли |
|||
нейных регрессионных зависимостей |
типа |
|
|||
У = а0 |
П |
( г =1, 2......... я). |
(226) |
||
|
|||||
Заметим, что линейные зависимости типа (226) мо |
|||||
гут не удовлетворять требованиям адекватности |
мате |
||||
матической модели реальному объекту. |
Если же ис |
||||
пользовать идеи |
кусочно-линейной |
аппроксимации |
|||
характеристик процесса [37] |
и процедуру |
корректиров |
|||
ки модели с целью приспособления ее |
к текущему |
||||
состоянию объекта, |
то можно достичь требуемой |
адек |
ватности и при линейном представлении уравнений
связи |
(219). |
|
|
|
|
|
|
|
|
|
|
||
В уравнении (226) неизвестные ап, а1— коэффициен |
|||||||||||||
ты регрессии. При известных значениях у, |
х |
(имеется |
|||||||||||
в виду, |
что |
статистика |
достаточная) |
величины ап, |
at |
||||||||
определяют |
методами |
классического |
регрессионного |
||||||||||
или корреляционного |
анализов с использованием |
ЭВМ. |
|||||||||||
Достаточную |
статистику |
о поведении |
параметров |
||||||||||
у, л: набирают путем |
организации непосредственно |
на |
|||||||||||
объекте |
регистрационного эксперимента [38]. |
|
|
|
|||||||||
Результаты эксперимента |
заносятся |
в табл. 20. |
|
||||||||||
|
|
|
|
|
|
|
|
|
Т а б л и ц а |
20 |
|||
Н о м е р э к с |
|
У |
|
X, |
|
|
|
|
|
|
|
|
|
п е р и м е н т а |
|
|
|
|
х к |
|
|
|
хп |
|
|||
|
|
|
|
|
|
|
|
|
|||||
1 |
|
|
У 1 |
|
Х 11 |
|
|
Х К, |
|
|
|
ч |
|
|
|
|
|
|
|
|
|
|
|
||||
2 |
|
|
Уз |
|
* 1 2 |
|
|
ч |
|
|
|
ч |
|
|
|
|
|
|
|
|
|
|
|
||||
N |
|
|
У н |
|
Х 2 Ц |
|
|
x k N |
|
|
|
x n N |
|
|
|
|
|
|
|
|
|
|
|
||||
В этой |
таблице у — значения выходного |
параметра |
|||||||||||
объекта при /-ом эксперименте. |
объекта |
при |
|||||||||||
у — значения |
выходного |
параметра |
|||||||||||
у- ом |
эксперименте |
(или |
в у- ый момент |
времени) |
|||||||||
у = 1, |
N.. |
х и — значение i-ой входной |
величины при |
у'-ом эксперименте.
311 |
181 |
Стало быть, табл. 20 представляет собой запись входных и выходных переменных эксперимента.
Для достаточности статистики число экспериментов {количество строк в табл. 20) должно отвечать следу ющему условию [39].
|
|
|
> 20, |
(227) |
|
где |
« — число |
входных переменных в |
представлении |
||
(226). |
|
регистрационного |
эксперимента |
||
При организации |
|||||
[38] |
следует обратить |
внимание на тот факт, что у |
и |
||
|
1,«) на |
одной |
строке таблицы |
фиксируются |
с |
учетом инерционных свойств объекта по соответству ющим каналам вход — выход, ибо практически всем технологическим объектам присуще явление запазды вания.
Рассматриваемому процессу также свойственно чис тое запаздывание. Определение времени запаздывания и его учет необходимы для разделения моментов съема данных о входных и выходных параметрах объекта.
При организации регистрационного эксперимента не обходимо правильно определить интервалы времени между замерами значений параметров. Чрезмерное увеличение этого интервала уменьшает точность ап проксимации характеристик объекта уравнениями, по лученными на основе собранной информации; умень шение приводит к частым опросам датчиков, так как очередные замеры могут быть сильно коррелирован ными друг с другом, информативность полученной та ким образом статистики окажется низкой. Поэтому значение интервала для конкретного технологического процесса рассчитывается и определяется, исходя из специфических особенностей процесса: время запазды вания по каналам регулирования, скорость изменения характеристик оборудования, скорость изменения спе цифических параметров самого моделируемого процесса, возможность работы информационных каналов (ско рость отбора информации, время лабораторных анали зов и т. д.).
Технологической процесс производства хлопкового масла и шрота в каждый момент времени представлен нами в виде машинных значений 19 параметров, среди которых имеются данные о расходах, давлениях, тем
182
пературах и, наконец, лабораторные анализы. Длитель ность анализов остается узким звеном в окончатель ной фиксации текущей статистики.
Последние исследования времени запаздывания по каналам регулирования и учет длительности лабора торных анализов для данного объекта позволили опре делить интервал между замерами очередных значений параметров, составляющий 2 часа.
Таким образом, с учетом соответствующего време ни запаздывания собрана статистика с объекта по ме тоду пассивного эксперимента [39] в течение 15 дней (по 12 групп замеров в сутки). За это время на Янгиюльском масложиркомбинате провели эксперименты при различных режимах функционирования технологи ческого процесса, обеспечивая весь диапазон измене ния расходов исходного сырья и других продуктов.
Суточный режим объекта эквивалентировался 12 группами значений упомянутых 19 параметров, т. е. фиксировался материал, состоящий из 19X12=228 слов. Следовательно, 15-дневная статистика содержала 3420 слов. Итак, исследуемый объект при моделировании можно заменить устройством, на котором записаны числа, состоящие из значений выходных и входных па раметров объекта. На целесообразность такого подхода уже указывалось в ряде работ [40].
В качестве указанного устройства наиболее разум но использовать запоминающее устройство электрон ной машины. Мы использовали память ЭВМ М-20.
Собранный с помощью регистрационного экспери мента и зафиксированный в памяти ЭВМ М-20 стати стический материал обработали с помощью стандарт ных программ, получили в явном виде регрессионные уравнения (220) —(225) и провели их подробный ста тистический анализ.
Уравнения сведены в табл. 21. Эти уравнения мож но принять за приближенную модель объекта. Оче видно, что такая модель может не соответствовать, не быть адекватной процессу в каждый конкретный мо мент времени. Со временем она будет расходиться с реальным объектом, так как характеристики изменя ются под воздействием ряда случайных неконтроли руемых возмущающих факторов. Поэтому модель не обходимо периодически корректировать. В конкретной
183
Т а б л и ц а 21
Номер |
|
|
|
|
а1 |
|
|
|
|
|
|
урав |
«0 |
|
|
|
Ох |
|
«г |
|
г |
|
|
нения |
|
«а |
|
|
«7 |
|
|
|
|||
220 |
1 7 ,8 4 |
2 , 6 9 |
— 0 ,9 7 7 |
0 , 1 6 |
- 0 , 0 6 |
0 , 1 5 |
0 ,0 8 9 |
0 , 8 9 |
5 , 5 8 |
1,0 |
0 ,9 1 |
221 |
1 9 ,7 6 |
— 1 , 3 3 |
0 ,4 1 |
— 0 ,0 1 |
0 ,0 1 6 |
0 , 0 6 |
— 0 ,0 1 |
0 , 0 7 8 |
9 , 2 0 |
1 ,0 |
0 , 9 5 |
|
«0 |
«8 |
«, |
|
*11 |
- |
- |
- |
|
|
|
222 |
— 0 , 9 6 |
0 , 8 2 |
- 0 , 0 1 4 |
— 0 , 2 9 |
0 ,0 0 3 |
|
|
|
9 , 8 8 |
1 ,3 5 |
0 , 9 5 |
223 |
— 1 ,8 2 |
0 , 3 9 |
0 ,0 0 3 |
— 0 , 0 2 ! |
0 ,0 0 0 4 |
|
|
|
1 ,9 2 5 |
1 ,3 5 |
0 ,7 0 5 |
|
«о |
«.а |
«■а |
- |
- |
- |
- |
- |
|
|
|
224 |
0 , 4 4 |
— 0 ,0 2 4 |
—0,022 |
|
|
|
|
|
1 ,6 0 9 |
1 ,4 5 |
8 ,6 2 5 |
225 |
9 2 , 0 |
3 , 5 4 3 |
3 5 ,9 6 |
|
|
|
|
|
3 , 4 9 |
1 , 5 |
0 ,8 4 9 |
постановке корректируются параметры модели, образу ющие множество коэффициентов регрессионных урав нений [а'\. Корректировку проводят по схеме метода обучающейся модели с применением алгоритма ста тистического поиска с пересчетом, используя текущую информацию о состоянии объекта [41; 42; 43]. В про цессе приспосабливания модели минимизируется функ ция качества (Q), зависящая от степени рассогласова ния выходов модели и объекта.
|
|
Рис. |
41. Схема |
метода обучающейся модели. |
|||
на |
Схема |
метода |
обучающейся модели [42] изображена |
||||
рис. |
41, |
|
|
|
|
||
где: X |
= |
(хх, х г, |
. . . , |
х п) — общий вход объекта и его |
|||
|
|
|
|
|
|
модели; |
|
|
У = |
(У\, |
У2, |
• •• . Ут) — ВЫХОД объекта; |
|||
|
У' = |
(Уь У2, |
• • • . Ут) — выход модели; |
||||
|
|
|
|
|
|
Q— функция |
качества; |
|
А = |
(аъ а 2, |
. . . , |
а'2) — вектор параметров модели. |
|||
ся |
Смысл метода сводится к следующему. Выдвигает |
||||||
исходная |
гипотеза |
об объекте, |
сформированная в |
виде определенной математической структуры, прини
маемой за первоначальную модел_ь объекта. |
При |
об |
|
щем выходе X выходы объекта У и модели |
У |
со |
|
поставляются между собой, а |
вычислительное устрой |
||
ство ВУ формирует величину |
Q = Q ( y , У'), |
характе |
ризующую степень несоответствия выходов объекта и его математической модели. При этом модель выби рают так, чтобы величина Q была минимальной, а в идеальном случае равнялась бы нулю.
185
Таким образом, задача построения модели объекта сводится к минимизации функции Q, определяющей степень точности модели. Минимизация производится
оптимизатором, так подбирающим вектор А' |
коэффи |
|
циентов модели, чтобы величина Q была минимальной. |
||
Начальное условие работы по этой схеме — наличие |
||
исходной гипотезы, |
сформулированной в виде некото |
|
рой математической |
структуры. Для нашего |
объекта, |
включающего в себя процессы форпрессования и экс тракции,—это совокупность регрессионных уравнений, приведенных в табл. 21.
При минимизации функции качества Q применяют различные алгоритмы поиска.
Существующие методы поиска разделяют на два класса: детерминированные (регулярные) и статисти ческие (случайные) методы поиска.
В первом случае алгоритм поиска представляет си стему действия, строго предопределяющуюся сложив шейся ситуацией. Это значит, что в одинаковых ситуа циях система действий будет также одинаковой, в про тивоположность случайным алгоритмам, допускающим систему действий в тождественных случаях.
Имеется много методов регулярного поиска: по очередного изменения параметров (метод Гаусса-Зей- деля), градиента, наискорейшего спуска и другие. Пе рейдем сразу же к характеристике статистического (случайного) поиска [43; 47].
Допустим, в пространстве оптимизируемых парамет ров из какого-то исходного состояния система делает случайный шаг. Если он не уменьшает функцию ка чества, то происходит возврат в исходное состояние. При уменьшении функции качества новое состояние принимается за исходное и из него делается следую щий случайный шаг. Процедура повторяется, пока не получаем требуемое минимальное значение функции качества.
Каждое состояние системы при случайном блуждении порождает определенное значение функции ка чества и, естественно, множество таких состояний дает на выходе системы некоторое множество функций ка чества (Q).
Очевидно, задача многопараметрической оптимиза ции в этом случае — выделить подмножества {Q*} с
186
{Qj, соответствующие некоторым требованиям, предъяв ляемым к оптимизируемой системе. Например, к этому подмножеству относят все значения функции качества меньше заданного или же равные нулю.
Основное достоинство метода — наличие элемента случайности, вводимого при поиске. Это позволяет вести систему в лучшем направлении, минуя скользя щие режимы, к которым часто приводят регулярные методы поиска. Случайный поиск не попадет в ловуш ку из-за элемента случайности, выводящего систему из затруднительных положений. Примером распростра ненной ловушки является „хребет“ для метода ГауссаЗейделя и „седло" для градиентного метода [43].
Рассмотрим шаговый алгоритм случайного поиска, использованный при корректировке модели исследуемо го объекта.
Поиск с пересчетом. Блок-схема алгоритма приве дена на рис. 42. При возврате системы после неудач ного случайного шага в какую-то ранее пройденную точку на ней не производится повторного вычисления функции качества и его удобно применять для опти мизации систем с мало изменяющейся функцией ка чества.
Рекуррентная формула смещения на г -f-1 -ом шаге
имеет вид а] , —щ + Ааг+1;
* _ *
а;, если Q (а[) < Qi-t
Д а ;+1 =
— Аа\ + а£, если Q(a,) > Q*_х
Для корректировки модели исследуемого объекта остановимся на алгоритме случайного поиска с перес четом. Это связано с тем, что из-за сильной инерцион ности объекта однажды рассчитанная в точке прост
ранства оптимизируемых параметров [а'\ функция качества Q практически остается постоянной за время работы алгоритма в контуре настройки.
Для работы контура корректировки необходима те кущая информация с объекта. При решении задачи непрерывной идентификации с помощью цифровой си стемы, подключенной непосредственно к действующему объекту [44], получение новой информации связано с опросом датчиков и преобразующих устройств, т. е. с необходимостью учитывать скорость работы информа
187
нионной части цифровых систем. Поэтому, хотя быст родействие УВМ само по себе большое, скорость на стройки модели в конечном счете зависит от скорости получения и преобразования информации с объекта. Это в какой-то степени нейтрализует высокоскорост ные качества вычислительных блоков управляющих машин.
Представление объекта эквивалентной статистикой, записанной заранее в памяти машин, необязательно управляющей, позволяет при идентификации перейти от реального масштаба времени к шкале машинного времени. Так, если цифровая система идентификации, связанная с объектом, вынуждена дожидаться очеред ного замера параметров через определенные интервалы времени, то при представлении объекта эквивалентной статистикой режимы выбирают из памяти машин, т. е.
“ I
Рис. 42. Блок-схема алгоритма с пересчетом.
188
скорость моделирования определяется лишь быстро действием ЭВМ или УВМ.
Поэтому в общем случае в последнем варианте на идентификацию затрачивают намного меньше времени. Эквивалентирование объекта статистикой в памяти ЭВМ. выгодно и тем, что пока не доказана необходимость установки на объекте средств вычислительной техники на производственных предприятиях, можно не устанав ливать сложных дорогостоящих технических средств.
Принятый в работе метод эквивалентирования объ екта позволил исключить все затруднения, связанные, например, с длительностью лабораторных анализов или необходимостью учета инерционных характеристик ка налов. Кроме того, в этом случае удобно по шкале машинного времени исследовать дрейф коэффициентовполученных регрессионных уравнений. Этот дрейф свя зан с изменением характеристик и режимов работы объекта. Его исследование необходимо, так как коли чественные оценки, получаемые при этом, нужны для эффективной организации настройки модели объекта с использованием алгоритма случайного поиска.
Разумеется, такое эквивалентирование целесообраз но, когда необходимо получить наиболее полные дан ные об объекте, математически описать его и выявить необходимость создания цифровой системы оптималь ного управления. Если все этапы уже пройдены, нд объекте функционирует система, модель корректируеется на основе текущей информации, получаемой че рез систему. Иными словами, работы по оптимальному управлению с использованием периодически коррек тируемой модели ведут через цифровую систему не посредственно на объекте.
С использованием эквивалентной статистики, накоп ленной в памяти машины, установлена практическая, возможность применения процедуры настройки для выбранного реального объекта и работоспособность алгоритма случайного поиска в контуре такой наст ройки.
Из совокупности уравнений (220) —(225) корректи ровке подвергли уравнение (222). Корректировку про водили по методике, изложенной в работе [41], с при менением алгоритма с пересчетом при величине шага; поиска а = 0,1. Результаты корректировки приведены
189
на рис. 43 для уравнения (222). На рис. 44 показаны изменения функции качества для уравнения (222). По оси абсцисс отмечено число шагов поиска, а по оси ординат — значения корректируемых коэффициентов регрессионных уравнений. Значения коэффициентов, помеченных индексом (к), соответствуют ко-нцу процес са настройки, когда обеспечивается такое множество
коэффициентов \а'\ , что уравнение по определенному критерию оказывается адекватным реальному объекту.
Организация настройки моделей поиском в прост ранстве параметров требует дальнейшего детального изучения возникающих при этом специфических воп росов: масштабирования шагов поиска; целесообразную организацию функции качества; выбор того или иного алгоритма из множества алгоритмов; обеспечение тре буемой надежности и достоверности статистического ма териала; создания быстродействующих систем получе ния и первоначальной переработки информации и др.
Решение этих вопросов, а также использование бо лее быстродействующих средств вычислительной тех ники позволяет резко повысить эффективность предла гаемых здесь методов и процедур настройки моделей технологических объектов управления.
УПРАВЛЕНИЕ ПРОЦЕССОМ ПРИ ПОМОЩИ ЛОГИЧЕСКОГО АВТОМАТА-СОВЕТЧИКА (ЛАС)
Для управления рассматриваемым процессом не тре буется высокой точности регулировки управляющих органов. Вполне достаточно той точности, которую соблюдает опытный оператор.
В этом случае возможно непрерывные величины за менить дискретными, разбив диапазон изменения пара метров на несколько участков. Тогда для реализации алгоритма управления вместо системы управлений ре шают логическую задачу оптимального выбора страте гии управления в соответствии с определенной ситуа цией в каждом конкретном случае. Во многих сложных ситуациях при одновременном отклонении нескольких параметров задача управления усложняется. Успешное ее решение зависит от опыта и интуиции оператора. Можно намного упростить задачу, формализовав и при менив логический автомат-советчик (ЛАС), реализую щий формализованный алгоритм управления.
190