Файл: Лебедкин, В. Ф. Проектирование систем управления обогатительными производствами.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