Файл: Лебедкин, В. Ф. Проектирование систем управления обогатительными производствами.pdf
ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 16.10.2024
Просмотров: 124
Скачиваний: 0
в работе [134], средняя квадратичная погрешность определения ве личины при ступенчатой экстраполяции принимает максимальное значение
^х=2[Кх{0)-Кх{х)} |
+ ^ |
х , |
(Ш.205) |
где /Ск(0)/Сж (т) — значение корреляционной |
функции в |
соответст |
|
вующие моменты сдвига. |
|
|
|
Если погрешность, рассчитанная |
по формуле (Ш.205), превос |
ходит заданную, то можно применить более сложные методы, ос нованные на вычислении полинома, (параболическая [85], триго нометрическая [149], статистическая [48] экстраполяции).
Например, нами при разработке алгоритма оптимизации про цесса обогащения в тяжелых суспензиях обогатительной фабрики Зыряновского свинцового комбината использовалась статистиче ская экстраполяция. Сущность этой экстраполяции по п точкам заключается в определении функции вида
? т с (0= |
2 Рі I х ('/) - |
тЛ+т*< |
(ІІІ.206) |
|
; = 1 |
|
|
где РІ — коэффициенты |
многочлена; |
x(ti)—значение |
случайной |
функции в момент tf, тх |
— математическое ожидание. |
|
Коэффициенты рі определяются из условия минимизации по грешности экстраполяции решением п уравнений:
2 |
РіКх (ts — ti)=Kx |
{t — ts) при 5 = 1 , 2 . . . Л, |
* |
=1 |
|
где Кх(т.) —корреляционная функция.
Были определены аппроксимирующие функции при статисти
ческой |
экстраполяции по одной, двум и трем точкам |
и получены |
|||||||||||||
OS |
|
1 |
|
|
ошибки |
предсказания. |
На |
рис. |
111.28 |
||||||
средние |
квадратические |
||||||||||||||
0,5 |
|
у |
и |
! |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|||
0,3 |
|
1 |
|
|
|
|
1 |
- я |
ft &1 |
> |
|
||||
|
|
|
|
|
|
г* |
|
||||||||
0,2 |
п |
Г |
1йп |
|
|
|
|
|
h |
|
|||||
0,1 |
|
|
|
|
\V г* г |
|
|
|
|
||||||
|
Z |
4 S |
в W 12 |
П 16 |
IS |
20 22 |
24 26 28 30 32 34 36 3d 40 |
42 44 46 |
±,емен |
||||||
|
|
Рис. |
Ш.28. Действительная и |
предсказанная |
функции ссг п р е д : |
|
|||||||||
|
р е а л и з а ц и я |
с л у ч а й н о й |
ф у н к ц и и |
|
v(t)\ |
2 |
— п р е д с к а з а н н а я |
ф у н к ц и я |
( с т а т и с т и ч е с к а я |
э к с т р а п о л я ц и я по д в у м т о ч к а м )
167
показано изменение во времени аппроксимирующей функции при экстраполяции по двум точкам.
Сравнение средних квадратических ошибок показывает, что если переход от ступенчатой экстраполяции к статистической по одной точке не увеличивает точности, то использование статисти ческой экстраполяции по двум и трем точкам существенно сни жают ошибки предсказания.
Кроме того, методы статистической экстраполяции хороши для стационарных эргодических процессов. В этом случае достаточно просто определяются ошибки экстраполяции. Если же процессы не могут быть сведены к классу стационарных, то модели прогноза определяются адаптивными методами [239], причем предпочти тельно применять их для стационарных процессов, так как объем вычислительной работы в этом случае значительно меньше.
Рассматривая задачи в классе адаптивных, обычно имеем дело
со структурой, |
в которой на вход поданы значения на |
n, |
п—1, |
||||||||
п — 2 ... |
шагах |
(глубина, очевидно, должна подбираться в |
про |
||||||||
цессе исследования), |
а |
на выходе получаем |
значение |
п+1 |
или |
||||||
в общем |
случае |
на |
n + k |
шагов |
вперед |
(см. рис. III.19). |
При |
этом |
|||
задача |
сводится |
к |
идентификации |
безынерционных |
объектов |
||||||
[239]. |
|
|
|
|
|
|
|
|
|
|
|
|
Поиск |
экстремума |
целевой |
функции |
|
|
|||||
Рассмотрим |
частный |
случай, |
когда |
|
оценка |
тесноты |
связи |
ме |
жду содержанием металлов в легкой фракции очень высока. При этом, применяя уравнение (III.173), допускаются несущественные ошибки с точки зрения потери эффективности. Этот случай имеет практическую ценность для монометаллических и некоторых по лиметаллических руд. Естественно, задача поиска экстремума имеет смысл только тогда, когда процесс находится в области по ложительной эффективности, т. е. когда выполняются ограничения г / 2 ^ г / * и а 2 < а | п р е д .
Возможны также случаи, когда выполняется первое требование и не выполняется второе. Тогда следят только за изменением г/*
{у* и агпред — функции времени, определяемые с помощью прогно
зирующих моделей).
Возможен и другой частный случай, когда при достаточно ап риорной информации о коэффициенте fo и незначительных его из менениях относительно среднего, очевидно, оптимизация состоит в наблюдении за величиной <z*opt , определяемой по формуле
4 о р < : = а 2 п р е д + / 0 _ |
( I I L 2 0 7 ) |
Цель управления в этом случае заключается в минимизации функционала
• / ( 0 = Л 1 а 2 о р Л ( а 2 о Р , - а 2 ) 2 ) - |
(III.208) |
168
Рассмотрим вариант, когда изменениями fo пренебречь нельзя. При этом известен вид функции а2 = /(г/г), но неизвестно значение параметров этой связи в каждый момент времени.
Для вывода технологического процесса в окрестность точки
экстремума достаточно решить |
задачу идентификации |
«2 = !{уг) |
и затем, пользуясь выражением |
для определения a*opt, |
вывести |
процесс в окрестность экстремума. Таким образом, недостаточность априорной информации приводит нас к необходимости совмещать изучение объекта с управлением им.
Связь между переменными представляется в следующем виде:
* 2 = / о + / У 2 . |
(Ш.209) |
Пусть в начальный момент состояние процесса соответствует точке 0 (рис. III.29). При этом известны переменные а,%, г/2, но
Рис. III.29. Итерационная процедура поиска оптимума (детерминированный вариант) :
/ , / / — ф у н к ц и и э ф ф е к т и в н о с т и , п р и м е н я е |
|
|
|
|
|
|||||||
м ы е |
с о о т в е т с т в е н н о |
при I и |
I I |
и т е р а ц и я х |
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
0 |
|
|
Уг |
неизвестны |
значения |
коэффициентов |
fa и fi . Предположим, |
что |
||||||||
fo = 0. По уравнению |
(Ш.209) |
на этом |
шаге найдем коэффициент f |
|||||||||
и по (III.207) определим a*2opt |
• Сделаем |
шаг в этом направлении |
||||||||||
(если « 2 > а * о р ( , уменьшим выход легкой |
фракции, соблюдая |
усло |
||||||||||
вие у2^У* |
|
, если a,2<a*opt |
увеличим выход легкой фракции). |
|
||||||||
|
В новом |
состоянии (точка |
/ ) , получив |
информацию |
о величине |
|||||||
«2 и уг, уточним коэффициенты fo и f при помощи |
следующего |
|||||||||||
итерационного |
процесса: |
|
|
|
|
|
|
|||||
/ о |
[ л ] = / о |
\п- |
1] +ъ |
[«] («2 Щ - / о [п- |
1] - / [ * - 1] у2 |
[л]); |
(III.210) |
|||||
/ і |
[ Л ] = / І |
[ л - 1 ] + 7 2 [ я ] ( а 2 [ л ] — / о [ я — 1 ] — / [ * — 1 ] У 2 [ л ] ) |
Уг[л], |
|||||||||
|
|
|
|
|
|
|
|
|
|
(Ш.211) |
||
где Yi [п], Уг [п] — скаляры, зависящие от п |
и удовлетворяющие |
схо |
||||||||||
димости итерационного |
процесса. |
|
|
|
|
|
||||||
|
Итерационный процесс строится по критерию минимума мате |
|||||||||||
матического ожидания квадратического |
отклонения |
|
|
|
||||||||
|
|
|
|
|
J(f)=M |
{ ( а 2 - / 0 - / у 2 ) 2 ) . |
(Ш.212) |
169
Остановимся на особенности приведенного алгоритма оптими зации. Как уже указывалось, функция а 2 = /(г/2) представляется кусочно-линейной формой, причем известно (из физических предпо сылок), что коэффициенты углов наклона и свободные члены ку сков рассматриваемой характери стики подчиняются неравенству
|
Функция |
критерия |
эффектив |
|||
|
ности в этом случае также будет |
|||||
|
иметь разрыв первого рода. При |
|||||
|
этом возможны случаи, когда оп |
|||||
Рис. III.30. К пояснению единствен |
тимальное |
значение критерия эф |
||||
фективности |
будет |
достигнуто |
||||
ности оптимума при f ô ^ / ô ' и f1 С f [ I |
||||||
именно в |
точке разрыва |
(рис. |
||||
|
III.30). При использовании |
выше |
приведенного алгоритма оптимизации в случае кусочно-линейной аппроксимации необходимо, очевидно, соблюдать еще одно огра-
\Руда
Легкая |
Разделение» |
|
Хвосты, |
|
|
фракция |
Флота |
концентраты |
|
||
|
в |
тяжелых\ |
ция |
|
|
|
суспензиях |
|
|
||
|
|
|
|
||
\Идентифика- |
а, Ъ |
|
Вычисление |
*2npl |
|
ция_ _ |
- |
|
|
В,*С„ |
|
|
|
"гпрЬ-Ч |
'2пр M ] |
||
|
|
33, |
33. |
||
|
|
|
|
||
|
|
|
|
Прогноз |
|
Идентификация |
|
|
|
||
|
|
|
|
*2преЭ |
|
|
|
f, |
Вычисление |
|
|
|
|
|
<?r{Kv+fo) |
|
|
|
|
Уг |
|
c2opt |
|
|
|
Управление |
с учетом |
|
|
|
|
ограничения |
|
|
|
Рис. III.31. Структура алгоритма оптимизации процесса |
|||||
обогащения руд в тяжелых |
суспензиях: |
|
Э З ь Э З г , . . ., Э 3 , — э л е м е н т ы з а д е р ж к и
170
ниченне
где у**— координата точки разрыва.
Если значение у** известно заранее, то задача оптимизации
решается просто, однако, когда координаты точки разрыва неиз вестны, возникает необходимость применения более сложных алго ритмов управления, в которых координаты разрыва должны опре деляться на основе обработки текущей информации. В структуру алгоритма управления необходимо включить алгоритмы опознава
ния разрыва. |
|
Структура алгоритма оптимизации |
процесса обогащения руд |
в тяжелых суспензиях показана на рис. |
III.31. |
II 1.6. ОПТИМИЗАЦИЯ ПРОЦЕССОВ ИЗМЕЛЬЧЕНИЯ И КЛАССИФИКАЦИИ
Выше была рассмотрена постановка задачи управления циклом измельчения и классификации. При ее формализации возможны различные варианты построения алгоритмов управления. При на личии статической модели процесса в виде уравнения множествен ной регрессии [179, 180] и рассматривая данный передел в виде самостоятельного контура или подсистемы управления общей си стемы управления фабрикой, задачу можно свести к максимизации частного критерия эффективности [34, 35], например количество готового класса на сливе классификатора или гидроциклона Qu- Отбор информации об основных параметрах процесса возможен с использованием серийных приборов; для измерения гранулометри ческого состава могут быть использованы гранулометры [77, 83].
В формализованном виде задачу можно сформулировать так [178]:
|
Q 7 4 = Q ß 7 4 — m a x ; |
|
|
||
|
ai<^=f(Xlf |
Xj, Vh |
Yl X.yù<a2; |
|
|
|
* i < 7 < £ 2 ; |
Q<cu |
|
|
|
где $n = f(Xi, |
X2., Yи Y2., XiYi) |
—уравнение множественной |
регрес |
||
сии, связывающее возмущающие ХІ |
И управляющие Yi |
воздействия |
|||
с процентным |
содержанием |
готового класса; ai, ai, |
Ь\, |
Ъі, ci — |
числовые значения параметров (ограничения).
Используя, например, алгоритм управления, методом перебора можно при данном значении величин возмущающих воздействий найти оптимальные управления.
Близкая к этой формализация задачи возникает и тогда, когда измельчительный агрегат работает при достаточно большой общей нагрузке. В этом случае необходимо наложить дополнительное ог раничение на допустимую производительность, под которой следует
171