книги из ГПНТБ / Козин В.З. Методы исследований в обогащении учеб. пособие
.pdfРис. 9.d. Экспоненциальная /1 / автокорреляционная функция Rj.^ (С) = {•<è~z/v и экспоненциальнокосинусная /2 / Roi*, (t) = 1 'Й- ^ т>Соіо)‘Сд ля содер жания металла в руде; вааимокоррепяционная
Функция /3/ между содержанием металла в руде и в хвостах
регрессивного анализа. Это, естественно, снижает точ ность описания, и, кроме того, требует учета методоло
гических особенностей. |
|
|
В ы б о р |
ч а с т о т ы |
з а м е р о в . |
Частота замеров при получении корреляционных зависи мостей процесса должна быть такой,, чтобы время между
замерами |
~Ь^ |
было равно или больше времени |
корре |
||
ляции ‘Ск |
: |
|
|
|
|
|
|
^ à |
■ |
■ |
/9 .2 1 / |
Приэтом выполняется одно из основных требований
применения статистических методов - независимость ис пытаний., /выборки/.
Требордние независимости выборки основывается на том, что выборочная совокупность должна наиболее полно и объективно, ртрщ.лт^ генеральную. Ооьзктивног.ть - это результат такого способа отбора, когда каждая
величина имеет такую же возможность быть отобранной, как и любая другая;.или, иначе, вероятность для того или иного варианта попасть в выборку равна вероятности этого варианта в составе генеральной совокупности. При
этом функция распределения вероятностей для выборочных данных при достаточном объеме выборки приближается к Функции распределения вероятностей генеральной совокуп ности.
Для стационарных случайных процессов требование /9,21/ явпяется излишним. Действительно, непрерывная реализация достаточной длины в этом случае дозволяет получить функцию распределения вероятностей, как угодно близко приближающуюся к действительной функции распре деления, т.е, для стационарного случайного процесса ко нечная реализация является объективной "выборкой" иссле дуемой переменной.
Таким образом, для стационарного процесса даже не возникает вопроса об ограничении интервала временя меж ду замерами в смысле независимости. Вопрос о частоте замеров возникает лишь в связи с необходимостью наи более точного представления функции по данным фискретных замеров или проб, так как для важнейших технологи ческих параметров используются приборы дискретного контроля.
Требование /9.21/, будучи справедливым при получении статистической модели любого процесса, в случае ста ционарного случайного процесса и инерционного объекта не только не обязательно, но и уменьшает возможность полного и точного исспедоваішя процесса.
При практическом применении корреляционного анализа задача ставится к получению на первом этапе коэффици ентов корреляции между всеми параметрами., измеряемыми в процессе исследования, и параметрами и показателем процесса, для которого находится уравнение парной или множественной корреляции. Коэффициент корреляции в обычном стати сгическом смысле - это значение норми рованной корреляционной функции при фиксированиям С - временном сдвиге между замерами. Применяя корреляцион-
иый |
анализ для получения |
математической |
ліидели ста- |
тики |
процесса, Iвыбираем и |
"замораживаем' |
ъ , после |
чего фактически вычисляем значение нормированной корре ляционной функция.
Таким образом, при снятии ограничений /9.21/ интер вал времени между соседними замерами может быть вы бран достаточно малым.
Коэффициент корреляции в случае непрерывного про цесса, являясь фиксированным значением взаимнокорреля ционной функции, зависит от частоты входного сигнала и свойств самого объекта. На коэффициент корреляции ока зывает влияние инерционность объекта, переменное время запаздывания и изменение характеристик объекта по време ни, Поэтому, должно быть поставлено условие проведения опыта в наккратчайший срок. При этом будет подучено некоторое "мгновенное" частное уравнение связи, но этс уравнение раскроет истинное влияние параметров на ход процесса.
Заметим, что цель достигается в том случае, если
|
|
|
T < T UJM, |
/9.22/ |
|
где |
Т |
- |
время |
эксперимента; |
|
|
Тіумвремя, |
за которое характеристики объекта |
|||
|
|
|
изменяются существенным |
образом. |
|
Если |
же |
Т * TUJM , то статистическим путем полу |
|||
чать точное |
уравнение нельзя. |
|
|||
Следует отметить, что при постановке эксперимента |
|||||
с частыми |
замерами, когда существенно сокращается |
время эксперимента, есть возможность считать многие неучитываемые возмущающие факторы постоянными, что заметным образом может повысить точность модели, а Также имеется возможность наглядного изучения измере ния исследуемых фактов по данным дискретных проб.
Определение длины реализации, необходимой дпя полученияматематической модели. Целью наблюдения является получение как можно более точной статичес кой модели. Применительно к корреляционному анали
зу это соответствует как можно более точному опре делению коэффициента корреляции.
Ошибка в определении коэффициента корреляции оце нивается в зависимости от среднеквадратичного отклоне ния коэффициента корреляции 3*, , вычисляемого по формуле
|
|
|
5 ъ |
1 - г 1 |
|
/9.28/ |
|
|
|
- |
|
||
где |
|
|
|
V Ï Ï H T ' |
|
|
П - число данных. |
|
|
||||
Для фиксированного значения |
1 величина |
будет |
||||
зависеть от |
И |
. Рассмотрим график функции |
ѴТѵТ |
|||
в зависимости |
от |
П |
. Можно видеть, что при |
|||
д > |
ІОСЬ200 замеров |
падение функции незначительно и |
||||
не может компенсировать затраты |
труда на получение |
данных. Следовательно, практически, число данных может
находиться именнф интервале |
П = 100+200. |
Выбор интервала времени между аамерами в соседних |
|
точках технологической схемы. |
|
При дискретной регистрации |
/измерения/ непрерывных |
случайных переменных на входе и выходе инерционного объекта, интервалы времени между замерами должны быть равны эквивалентным запаздываниям соответствую щих каналов объекта.
Эквивалентное запаздывание - это временный сдвиг, при котором взаимокоррепяционная функция имеет макси мум, т.е. теснора связи между входным и выходным сигналом наибольшая. При известных параметрах объек та и свойствах входного сигнала величина эквивалентно го запаздывания может быть определена расчетным путем.
Полное эквивалентное запаздывание для инерционного объекта с чистым запаздыванием складывается из двух
величин: |
( |
|
кВ |
+ ^ , |
/9.24/ |
где |
er1 |
|
- |
эквивалентное запаздывание инерционной |
|||||||
L эк| |
|||||||||||
|
|
|
|
|
части |
объекта; |
|
|
|
||
|
|
|
|
- |
чистое |
запаздывание объекта. |
|
|
|||
|
Эквивалентное запаздывание для инерционного звена |
||||||||||
первого |
порядка |
|
|
р |
2 То |
|
|
||||
|
|
|
|
|
|
І'вЛ |
’ |
|
|||
|
|
|
|
|
|
= То‘ Т |
Т о + т |
|
|||
где |
То |
- |
|
|
|
|
Т о - Т |
|
|
|
|
постоянная времени объекта; |
|
|
|||||||||
|
X |
- |
параметр |
автокорреляционной функции |
вход |
||||||
|
|
|
|
ного сигнала, записанной в виде |
R(T) = T6 L^T |
||||||
|
При |
То » |
Т |
величина |
почти не зависит от |
||||||
параметра |
|
Т |
н стремится к значению |
|
|
||||||
|
|
|
|
|
|
|
|
= т„іп г |
= о.ѳэ т„. |
/925/ |
Для инерционного звена второго и более высокого порядка эквивалентное запаздывание определяется по бо лее сложным формулам. Поэтому целесообразнее инер ционные участки объекта высокого порядка представ лять в виде инерционного звена первого порядка с чис тым запаздыванием.
Влияние нестационарности и инерционности объекта
на статистические связи. Неоташюпарность динамичес
ких параметров процесса приводят к тому, что экспери
ментально получаемый коэффициент корреляции оказыва ется заниженным. Сравнительно просто это можно пока зать для объектов чистого запаздывания и синусоидаль ных входных сигналов.
Пусть |
X |
( f ) |
= а |
iLn |
( u ) t + *f) |
; |
|
ÿ |
(+■) |
= a |
ôùn |
І |
+ %) > |
т.е. фактическая связь между переменными со сдвигом на время запаздывания будет ÿ = л ■
♦
При получении уравнения связи выбирают для постро |
|
||||
ения корреляционного поля некоторый постоянный сдвиг |
|
||||
Ÿ . Обозначим |
Д І'і = Ÿx “ Ÿ |
. Найдем среднее |
|
||
значешіе |
дня |
tl |
замеров, |
соответствующих одному |
|
и тому же |
значению |
X |
|
|
|
~ ~гГ ' <z |
|
+ У)‘ |
(“>k +'f)' b** & |
> |
При симметричном /в частности, нормальном/ законе рас пределения
Z ü n A Ÿ c ^ O ; |
U |
X C o s ù Y i = i < 1 , |
L |
|
*• |
Следовательно |
|
|
|
|
/B2ä/ |
Таким образом, снимая экспериментально зависимость
ÿ - X. в условиях переменного времени запаздывания, мы получили бы вместо существующей у~&Х при &<І Коэффициент корреляции при этом получился бы также меньше 1. Этот вывод можно распространить и на стацио нарную случайную функцию.
Поэтому для получения более точных уравнений связи необходимо вести эксперимент при достаточно малых коле баниях объемного расхода пульпы.
Инерционность объекта также приводит к снижению коэффициента корреляции, т.ѳ. нарушается равенство
* = ъ » ■ приведенное в рааделе УШ. Особенно сильно это проявляет
ся в случае наличия у входных факторов резких всплесков или провалов.
Р а з д е л X. ЦЕЛЕВЫЕ ФУНКЦИИ
10.1Виды целей и ограничения
Целевая функция - это строгая формулировка цепи, которой нужно достичь в результате исследований. Жела
тельно, чтобы |
эта формулировка была |
математической. |
||
Если |
- некоторые выходные показатели, то |
|||
целью может, быть стабилизация у |
при произвольном % |
|||
|
conit, |
|
/ 10. 1/ |
|
Экстремизация |
/минимизация либо максимизация/ ÿ |
при |
||
произвольном |
£ |
е д х г , |
|
|
|
и —^ |
|
/ 10.2/ |
|
|
|
|
|
|
и достижение условного экстремума |
|
|
||
|
и, — - |
е><Ьг ; |
|
|
|
- |
c o m t . |
|
/10.3/ |
Общей целью управления промышленным технологическим |
процессом является производство материальных благ соот ветствующих данной технологии при наименьших затратах общественного труда. Эта бо льшая цель применительно к конкретному технологическому процессу нуждается в более узкой формулировке, так как вычислить затраты обществен ного труда практически невозможно. Кроме того, любой технологический процесс связан со многими ограничениями, которые обязательно должны учитываться при формулировке дели. Наконец, в большинстве случаев изменение и упро щение математической формулировки цели управления свя зано с необходимостью учета реальных возможностей рас четов целевой функции в приемлемое для цепей управления время. Желательно чтобы любая упрощенная запись цели управления достаточно строго вытекала из общей записи.
Рассмотрим наиболее характерные ограничения, накла
дываемые |
на управление процессом. |
|
Важнейшим ограничением является ограничение, накла |
||
дываемое |
на содержание |
основного металла в кон |
центрате. |
|
|
Обычно это ограничение типа неравенства |
|
||||||
|
|
|
ß |
* |
ß n A , |
|
/10.4/ |
здесь |
ßnA |
- плановое содержание |
основного металла, |
||||
|
|
|
в концентрате. |
|
|
||
Часто накладываются ограничения на содержание jön |
|||||||
каких-либо |
примесей в концентратах |
|
|
||||
|
|
|
ß i -П |
^ JàLnnpeg, |
|
/10.6/ |
|
Здесь ß i n |
пред ~ предельное значение |
содержания |
С-ой |
||||
|
|
|
примеси |
в концентрате. |
|
||
Другим, играющим большую роль, является ограничение, |
|||||||
накладываемое |
на производительность по руде |
|
|||||
|
|
|
0-р |
^ |
О. р пл |
|
/10.Ѳ/ |
Здесь |
|
па |
- плановая |
производительность по руде. |
Плановые и предельные значения параметров могут быть переменными и связанными либо с условиями обогащения, либо с изменением условий работы поставщиков и потреби телей сырья.
Целый ряд параметров процесса имеет технические ограничения, например, расход воды Gê , реагента ÿ,p , циркуляционная нагрузка Q. ц ;
|
^ |
max |
|
/10.7; |
|
|
|
||
9'Р |
^ |
9'і: тах |
> |
/ 10.8/ |
|
||||
0-ц |
^ |
Q.4 ma x |
|
/10,9/ |
Наконец, ограниченными являются возможности техно |
||||
логического персонала |
по управлению, |
|
10.2.Общий подход к выбору целевой функции
Сматематической точки зрения конечной целью 'про
изводства является достижение экстремума некоторйй функции /функционала / f
|
F - F ( X , Y ,2 |
, и ,t |
) |
t |
/ і о д о / |
|
где |
Ъ ~ время^ X - |
возмущающие, |
Y |
—управляющие |
||
|
воздействия, LL |
- помехи, 2 |
- |
выходные показатели-. |
||
|
, |
137 |
|
|
|
|
На |
все |
составляющие функции |
F |
наложены ограни |
|||||
чения типа |
равенств или неравенств. В общем виде |
||||||||
|
|
X і min |
< |
X |
« |
X i |
max |
; |
/ 10.11/ |
|
|
Y tn Ln |
« |
Y |
« |
Y |
max |
|
/ 10. 12/ |
|
|
-2 min |
4 |
г |
« |
Ï |
щ а х . |
/ 10, 13/ |
|
|
|
|
|
|
|
|
|
|
|
Если процесс является непрерывным, все входные и вы |
|||||||||
ходные |
величины, как |
л критерий эффективности |
/10.1 0/ |
являются функциями времени. Нели вычисляется некоторый сводный показатель F по измеряемым входным и вы ходным величинам, то необходимо компенсировать погреш ности, связанные с динамическими свойствами объекта. Одним из способов компенсация является сдвиг реализаций на величину транспортного запаздывания и усреднение из меряемых величин в течение промежутка времени.
Детальный анализ целевых функций дан в работах Л.А.Барского. В его работах сформулированы и некоторы. требования, которые могут быть предъявлены к целевой функции. Это: эффективность с точки зрения конечной цели, количественность, универсальность, статистическая э<рфективность, наличие физического смысла, простота и легкость вычисления и т.п.
В связи с наличием большого числа предложенных фор мул для целевых функций, выбор нужной формулы сложная задача. Дело в том, что после того, как предложена новая формула, ее часто начинаю’! рассматривать в отрыве от тех условий, для которых она сформулирована. Условия же могут быть различными. Различными могут быть и возмож ности получения для вычисления целевой Функции необхо димых данных.
Действительно, пусть, например, в лабораторных усло виях изучается флотируемость какой-либо руды. Иссдецователь знает, что его рекомендации послужат основой для строительства обогатительной фабрики. Он, естественно, жепал бы наиболее точно оценить свое решение в соотве'і ствни с дальнейшей целью переработки руды на фабрике,
где, как известно, учитывается не только содержание полезного ископаемого в продуктах обогащения, но я затраты на реагенты, электроэнергию, руду, а также многочисленные ограничения. Однако, при постановке лабораторных опытов исследователь не располагает дан ными ни о конкретных свойствах рабочих аппаратов, іш о качестве реагентов, которые будут использоваться на фабрике, ни о фактических затратах на добычу руды и т.п. Использование ориентировочных оценок позволяет отбросить некоторые резко различающиеся между собой
решения, но не дает возможности выполнить точпые коли чественные расчеты. В таком случае исследователь вынуж ден обращаться к формулам, которые объединяли бы име ющиеся в его распоряжении данные и позволяли бы срав нивать результаты отдельных опытов между собой. Приме ром такого критерия является единый критерий /метод/ оценки эффективности обогатительных процессов, предло женный Н.Г.Тюпенковым.
/10.14/
Несмотря на четкость физического представления оцен киэффективности процесса по Н.Г.Тюренкову, её нельзя применять к собственно процессам обогащения, т.к. про цесс будет одинаково эффективен, если извпечёте породы в концентратах будет уменьшаться, а потери металла с хвостами - увеличиваться.
Подобные целевые функции принято называть критери ями оценки технологической эффективности. Обычно это комбинации нескольких выходных показателей обогати
тельных процессов. |
Если сС - содержание |
металла |
в |
||
руде; |
о/3 - то |
же в концентрате; |
^ |
- то же |
в |
хвостах; |
% —выход концентрата, |
то достаточно широ |
ко используемыми фгормулами в качестве целевых функций являются1: извлечение
/10.15/