книги из ГПНТБ / Козин В.З. Методы исследований в обогащении учеб. пособие
.pdfПроверка |
адекватности модели |
|
|
|
|
||||||||
|
Sk |
~ ? |
(ÿi ~ ÿi-pacz) |
|
|
|
|
|
|||||
Число степеней |
свободы |
|
|
|
|
|
|
||||||
|
/ _ а/ |
( l + z H t + i ) |
|
= op ( * + г ) ( з и } - ю |
|||||||||
|
J* |
|
|
|
|
|
X |
|
|
|
|
2 |
|
По ранее |
полученным данным ошибка воспроизводлмости |
||||||||||||
i l = |
2,25. |
|
|
|
|
|
|
|
|
|
|
|
|
|
V |
|
_ |
|
|
|
|
= 5 1 8 |
= 51,8 |
- |
|
||
|
|
|
- |
|
J* |
|
|
|
|||||
|
|
|
|
|
|
|
TÖ“ |
|
|
|
|
||
.Расчетные значения критерия Фишера |
|
|
|||||||||||
|
Г |
|
|
о аа. |
|
_ |
51,8 |
= |
2з . |
|
|
||
|
Грpact |
|
'ПГ- |
|
|
"2 Ж |
|
|
|
|
|||
|
|
|
|
|
й |
|
|
|
|
|
|
||
Табличное значение критерия Фишера равно |
2,54, т.е, |
||||||||||||
модель неадекватна. |
|
|
|
|
|
|
|
|
|||||
Запишем, в итоге, уравнение в натуральных координатах, |
|||||||||||||
используя преобразования. |
|
|
|
|
|
|
|||||||
|
- |
|
_ Я'і - іо |
|
|
|
ѵ _ 9 л - а |
|
|||||
|
|
|
|
|
|
г |
_ |
ъ |
- |
и |
|
|
|
|
|
|
|
|
|
|
3 |
|
8 |
. |
' |
■ |
|
ÿ * |
15i>7* |
|
|
|
|
|
|
|
ъ м с р ь |
+ |
|||
T 0 ,1 2 & ^ - 0 ,1 5 7 '^ + 0 , 0 1 7 7 ^ 5 < ' |
|
||||||||||||
+ |
|
|
ÿ ity s |
+ |
°>i5 e |
< ^ 3 |
+ |
Q;° |
|
• |
|||
У к а э + с е м |
е щ е |
р а ь ? |
ч т о |
|
к о э с р с р і - щ ц е н т ы |
||||||||
при |
Оуг 5 |
c j , 1!L |
и |
<^, £ г |
н |
е з н |
а ч и |
м ы . |
|
|
ІА. Регрессионный анализ при ротатабепыюм |
|
||||||||
|
планировании |
второго порядка |
|
|
||||||
После вычисления коэффициентов регрессии находят |
||||||||||
связанную с ними сумму квадратов |
|
|
|
|
||||||
S = a 0-(0lj)+ |
^ a r f U j ) * |
■k |
|
|
|
|||||
Z |
OLu’(Cjÿ). /7.49/ |
|||||||||
с числом степеней |
свободы |
|
HJ • |
J |
|
|
||||
|
г _ |
|
|
|
, r |
j |
|
|
|
|
|
(&+z)(-k+i: |
|
|
|
/7.50/ |
|||||
|
J |
~ |
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
Остаточная сумма |
квадратов |
|
|
|
|
|
|
|||
|
S R |
= |
ЗоЛц |
“ |
S .= |
(у# ) - S . |
|
/7.51/ |
||
С числом степеней свободы |
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
/7.52/ |
Если |
нулевые опыты |
используются для |
вычисления |
5 ß |
||||||
с числом степеней |
свободы |
П.0 —і |
, то |
|
|
|||||
|
( |
_ |
./ |
(è+2)(Â +l) |
( м |
' |
/7,53/ |
|||
|
J R |
-J1 - |
s------2"---------- ('г° Ü |
|||||||
Sjj |
может быть найдена и таким |
образом, как |
это |
сделано |
впримере 7.2.
7.15.Разбиение на ортогональные блоки
Так как редко удается поставить все необходимые
эксперименты в одинаковых условиях, целесообразно разби ение плана па ортогональные блоки.
Условия могут нарушаться, например, потому что при ходится готовить новые порции реагентов, брать несколь ко отличные пробы руды и т.п.
Разбиение на ортогональные блоки позволяет устранить влияние низкочастотного дрейфа.
Например, матрица ротатабепьного плана для двух пере менных может быть представлена в виде двух ортогональ ных блоков.
Таблица 7 ' Матрица, разбитая на ортогональные блоки
; |
|
Л |
Х і |
X z |
|
|
|
|
|
Блок |
1 |
1 |
+ 1 |
+1 |
|
|
2 |
- 1 |
+1 |
|
|
3 |
+1 |
- 1 |
|
|
4 |
- 1 |
- 1 |
|
|
5 |
0 |
0 |
|
|
6 |
0 |
0 |
Блок |
2 |
7 |
-1,41 |
0 |
|
|
8 |
+1,41 |
0 |
|
|
9 |
0 |
-1,41 |
|
|
ю - |
0 |
+1,41 |
|
|
И |
0 |
0 |
|
|
1 2 |
0 |
0 |
7.16.Исследование почти стационарной области, •представленной полиномами второро порядка
Известно, что любое уравнение второго порядка общего вида путем преобразования координат может быть приведе
но |
к канонической форме, |
например |
• |
|
|
|||
|
|
ÿ = |
ß d x j |
+ ß 2 |
|
|
||
|
Еозможнс. несколько вариантов сечений гиперповерхно |
|||||||
сти |
описываемой такими |
каноническими уравнениями |
в |
|||||
зависимости от |
знака и величины коэффициента |
ß^ и |
В^ |
|||||
|
а / Эллипсы, |
ßji |
и |
ß 2 |
имеют |
одинаковые |
знаки. |
|
Гиперповерхность имеет явный максимум или минимум. |
|
|||||||
|
б/ Гиперболы, |
и |
ß 2 |
имеют |
разные знаки. Гипер |
|||
поверхность имеет |
седло. |
|
бесконечности/. Один |
|||||
|
в/ Параболы |
/центр |
находится в |
из коэффициентов равен нулю. Локально гиперповерхность отклика представляет собою возрастающее повышение "гребень".
г/ Параллельные прямые. Один из коэффициентов равен нулю /условно центр находится в тобой точке оси/. Гиперповерхность стационарное возвышение /невозраста-
юідий гребень/.
Так как переход к канонической форме требует значи тельных преобразований, рекомендуется по полученному уравнению нанести несколько сечений на миллиметровой бумаге, по возможности графически оценив фокальные точки, после чего отнести гиперповерхность к одному из четырех типов и принимать решения, руководствуясь свой ствами этих типов гиперповерхностей.
7,17. Планирование отсеивающих экспериментов
Если заранее.неизвестно, какие факторы могут оказать
ся существенно влияющими на критерий |
оптимизации .це |
|
лесообразно рассмотреть все возможные |
факторы и по |
|
пытаться отсеять |
малозначимые. |
|
Это значит, что |
на первых этапах исследования может |
рассматриваться большое число факторов.
Для отсеивания малозначимых факторов используется три метода :
а/ Насыщенные планы.
Насыщенными называются планы, у которых все эффек ты воздействия заменены новыми факторами. Следователь но, предполагается, что доминирующими являются линей ные эффекты. Например, для числа факторов '3 насыщен ным будет план
Можно видеть, что насыщенный план можно построить
ч
не для любого числа факторов, а лишь для 3,7,15,31 и
т.д. т.е. для |
чисел Я |
- і . |
"к = 11,19,23,35. |
Имеются |
разработки |
и для |
б/ Сверхнасыщенные планы Если ввести дополнительное условие - число сущест
венно влияющих факторов сравнительно мало по сравнению с первоначально исследуемым количеством, можно перей ти к сверхнасыщенным планам /число степеней свободы для них отрицательно/.
матрицу уровней факторов по столбцам /пользуясь, напри мер, таблицей случайных чисел/. Ранг' матрицы /число наблюдений / / / выбирается существенно меньше чисщ ■эффектов, взятых под подозрение .
После выполнения экспериментов результатъ! наносят на диаграммы рассеяния.
Визуально выделив значимые Факторы, вычисляем ли
нейные эффекты по формулам |
5 - |
|
U J |
м |
’ |
после чего корректируем результаты, |
вычитая из результа |
|
тов полученных на уровне + 1 для |
j -го |
фактора значе |
ния эффекта, после чего строим вновь диаграмму рассея ния и т.д. до тех пор, пока не будут выделены все значи
мые |
факторы, т.е. разброс оставшихся будет соизмерим |
|
с допустимой ошибкой в оценке |
^ |
|
в/ |
Последовательное отсеивание. |
|
Если заранее известно, что |
число значимых факторов |
мало, то можно использовать последовательное отсеива ние но группам факторов. Для этого факторы группируют ся по несколько .л эти группы рассматриваются как новые факторы. Среди них выделяются значимые, например с помощью насыщенных планов, а незначимые группы сразу отбрасываются и т.д.
/Этот метод можно сравнивать с методом поиска од ной фальшивой монеты среди, например, ста других/.
Р а з д е л УШ. РЕГРЕССИОННЫЙ И КОРРЕЛЯЦИ ОННЫЙ АНАЛИЗЫ
8.1. Уравнения для определения коэффициентов
' Обычно, при постановке лабораторных экспериментов, либо при обработке результатов наблюдений за производ ственным процессом отыскиваемые связи представляют в графическом виде, так как при этом анализ становится наглядным. Для удобства анализа получаемых графиков их представляют в виде лилейных, либо релинеііпых уравнений,
Такое представление становится почти необходимым при изучении влияния многих переменных, так как в этом случае графики могут передать лишь некоторое сечение изучаемых гиперповерхностей, но никогда не могут дать представления о всей поверхности в целом. Наконец, по требность в отыскании связи в виде некоторых уравне ний возникает при получении большого числа эксперимен тальных данных, графический анализ которых неточен и затруднителен.
Дополнительным преимуществом математических мето дов является однозначность результатов обоснованной оценки точности получаемых уравнений.
Уравнения связи используются для расчетов, предска зывающих с некоторой точностью появление определенных значений выходных показателей при заданных входных.
Если имеется набор экспериментальных данных X и ÿ , то с помощью математического аппарата можно полу чить уравнения почти любого вида. Например, в виде по линома любой степени, экспоненты, гиперболы, логарифми ческой функции либо в виде функций,обратных указанным и т,д. Однако, как правипо, математически выбор формы связи не решается. Можно лишь проверить, насколько хо рошо удовлетворяет выбранная форма связи эксперимен тальным точкам. Таким образом, выбор формы связи или,
иными словами, аналитического вида подучаемой закономер ности зависит от профессиональной подготовки исследова теля, Для этого выбора он должен располагать некоторы ми теоретическими предпосылками. Если же исследова тель не имеет возможности заранее выбрать вид связи, приходится руководствоваться формой расположения точек на графике. Можно лишь утверждать, что в практически наблюдающихся ■диапазонах изменения переменных промыш
ленных процессов уравнение |
вида |
а + ё х применимо |
в подавляющем большинстве |
случаев /примерно в 80-90 |
|
случаях из 100/, Использование |
полинома второго поряд |
|
ка ^ = а + ё х + с х г доводит эту |
цифру почти до 1 0 0 , |
Естественно, что научение влияния изменения входной
величины X в широком диапазоне, что и наблюдается в Лабораторной практике, приводит чаще всего к нели нейным уравнениям, принципиально любой формы. Так как одному и тому же значению аргумента в наборе экспериментальных данных соответствует несколько зна чений искомой функции, то получаемая связь описывает
лишь изменение среднего значения |
|
у [ |
при изменении X . |
|||||||||||
т.е. в процессе обработки найдем связь |
|
|
|
|
|
|||||||||
|
|
|
|
ÿ L = OL + |
ê x c . |
|
|
|
|
|
|
|
||
Для определения коэффициентов |
CL |
и |
& |
использу |
||||||||||
ем метод |
наименьших квадратов. Формулы для опреде |
|||||||||||||
ления |
ОС. |
а |
£ |
находятся |
из |
условия |
|
|
|
|
||||
|
|
|
_ |
^ = Z ( Ч і - У і ) |
|
|
m i n , |
/ 8 .1 / |
||||||
Тогда, |
если |
|
с=і |
* |
° |
|
|
|
|
|
||||
|
—CL + оХ.^ , |
условием минимума |
будет |
|||||||||||
равенство нулю частных производных по |
X |
и |
£ |
вы |
||||||||||
ражения / 8 .1 / |
І |
|
|
|
Л |
|
|
|
|
= |
||||
|
jj = |
|
|
|
|
|
|
|
||||||
|
JÜr |
= |
ê . |
Z ( u i |
- a |
- ê |
x |
i ) ( - X i ) |
= 0 . |
|
,ѣ’і ' |
|||
|
dg |
|
|
o 1. |
|
|
|
|
|
|
|
|
|
|
Отсюда получаем систему уравнений, позволяющих |
||||||||||||||
найти |
CL |
|
и |
ê , |
V |
|
|
V |
|
; |
|
|
|
|
|
|
|
|
|
= я Ѵ + й г |
X i |
|
|
|
|||||
Если же ищется связь |
у.^ = а |
+ é>JC£ + С at£ |
, то |
|
/аз/ |
|||||||||
систе |
||||||||||||||
ма уравнений для |
нахождения |
CLt ё |
и |
С. |
будет |
|
||||||||
|
' £ y L = |
ct'JÏ + ê £ x L+ C ^ x l |
; |
|
|
|
||||||||
|
< |
|
t Ly L - a ^ x |
L + g % . x l + C l . x f |
■ |
|
/8t4/ |
f X L^L = « f * ? + |
X? + C £ X * . |
Аналогично можно получить уравнения и для полино
мов более высокого порядка.
Таким образом, для определения интересующих нас коэффициентов необходимо вычисление соответствующих сумм и решения системы уравнений / 8 ,3 / и / 8 .4 /.
Можно пользоваться приведенными выражениями, но при ручной обработке данных используется ряд приемов, упрощающих и облегчающих расчеты. Покажем это на примере.
П р и м е р 8.1. Пусть имеется ряд наблюдений над двумя величинами: гранулометрическим составом пульпы,
поступающей на флотацию Г |
в % - 2 0 0 |
меш |
и содержа |
|
нием меди |
s хвостах основной Флотации |
i f |
в %. Для |
|
нахождения |
уравнения связи |
нанесем на пасте милли |
метровки точки, соответствующие отдельным наблюдениям Г а 1Г и соответственно выбранных осях так, чтобы . точки располагались в пределах всего графика, рис.8 .1 ;
При этом мы получим поле корреляции. Затем оси |
Г и |
|||
Т/' |
разбиваем |
на интервалы, число которых К |
ориен |
|
тировочно можно определить |
по формуле /8,5/. в зависи |
|||
мости |
от числа |
наблюдений у |
. |
|
|
|
К = 1 + 1 , 3 |
|
/ 8 .8 / |
К округляется до ближайшего целого значения. В нашем
примере |
V |
= 1 0 0 и |
К = 8 . По оси ординат принимаем |
|||
7 интервалов, а по оси абсцисс - |
9 /для удобства/. |
|||||
Заменяем натуральные координаты точек упрощенными. |
||||||
Для этого интервалы |
по оси |
Г |
нумеруем от -4 |
до +4, . |
||
а по оси |
гг |
от -2 |
до +4. |
При |
этом координату |
0 при |
даем столбцу и строке с наибольшим числом наблюдений, либо близким к нему. Все дальнейшие расчеты, таким образом, проводятся в упрощенных осях- Г 1 и ТГ"1 . Под
считываем число точек в каждой |
подученной клетке и |
записываем в ней посредине соответствующую цифру. |
|
Умножаем эту цифру на координату |
и записываем в |
правом верхнем углу клетки. |
|
Теперь можно перейти к подсчету сумм. Для этого
необходимо заполнить соответствующие строки под табли |
|
||||||||||
цей. . |
|
М |
получается |
суммированием |
в столбце |
|
|
||||
Строка |
|
|
|||||||||
всех цифр, |
записанных посредине клеток. Например для |
|
|||||||||
столбца |
с координатой |
Г = / —1 / |
П = 1 + 1 + 1 |
+ 5 |
+ |
|
|||||
+ 2 = 10. Суммируя цифры в строке |
Ш получим V =100. |
|
|||||||||
Строку |
пГ-,_ |
получим, умножая цифры в строке |
Л |
на |
|
||||||
значения упрощенных координат |
Гс‘ |
. Например , для |
|
||||||||
того же столбца П Гі |
= 10 х /-1 / |
= -10. |
Сумма |
цифр |
|
||||||
в этой строке |
даст |
= 43. Строку П(f/7'j^получим, |
|
|
|||||||
умножая цифры в строке ПГі |
на значения |
Гс |
, т.е. |
|
|||||||
для нашего примера |
|
-10 х /-1 / = 10. Сумма |
|
|
|||||||
цифр в строке даст Е(Гс)&= 507. |
|
|
|
|
|
|
|||||
Таким же образом можно получить суммы £(/Ѵ)3^ |
|
’ |
|||||||||
требующиеся для системы /8.4/. |
|
|
|
|
|
||||||
Суммируя по столбцам, цифры в правых верхних углах |
|
||||||||||
клеток, |
получим цифры для |
строки |
. Например, |
для |
|
||||||
столбца |
с |
координатой |
Г ' |
= /-1 / Ѵ}' =4+2+1+0-2 = 5. |
; |
||||||
Умножая цифры в строке |
|
на значение координат Г |
|
||||||||
получим цифры строки |
Гі 'Ѵ'і |
. Суммируя строчки, |
полу |
|
|||||||
чим /Lût |
|
= 14 яЕГі'Ѵі' ~ |
|
|
|
|
|
|
|
||
Если |
разделить цифры VI ' |
на цифры строки |
Л |
,то |
|
получим данные для|построения эмпирической линии регре
ссии VÎ |
.Н а нашем графике она представляет собой ло |
|
маную линию. |
_ |
|
Найдем |
значения CL' и ê |
в уравнении t?l'= OL +ê f |
для чего в систему /8.3/ подставим соответствующие сум-
мы- |
|
Г i k |
= m a ' + J +ъё* ; |
|
||
|
I |
[ - 7 |
8 = |
Jf5a' + 507ê’ . |
т ж |
|
Отсюда |
^ |
|
О! |
= —0,174. Сиедаватеаьио в |
||
CL |
= 0,21 |
|
О |
|||
упрощенных координатах |
|
|
|
|||
|
|
|
- |
0,21 - в , і 7k f i . |
/а. 7/ |
п |
|
ft |
$ В 1ІО ß 2И5 9 9 |
100 = /У |
|||||
па' %Гі! |
-2Н-І8-М-Ш 0п |
SO27 i& |
кЪ = %Гі! |
||||||
пГі‘*ШІг |
9В5 И 0 |
0 |
ОП ВО М М |
507 = № ' і а |
|||||
|
|
т іrï |
і м |
в |
ь |
іг |
і |
н -h -» |
М |
|
|
J.7îîrc |
-ko-н -м |
I B-il-W |
- 7 i * z W . |
||||
|
|
jfL _ |
ij <1,5 |
fl |
OßtßOßtß-l\m,7 |
|
|||
|
|
ï ï i - n1 i ß m ^ m ß ü ß H ß m -ü .5 m t |
|
||||||
|
|
tvî-гг|аДОaß ДОДОо ,01,010.5307і |
|
||||||
|
|
ІѴі-Ѵ‘А |
Piß<ДО0.21,3Ö0,20J5 2,376,392tßi |
||||||
T~ эмпирическая |
линия |
peepeeeuuj |
|
||||||
й |
- |
т е о р е т и ч е с к а я |
лини$ь р е в р е с с и и ; |
|
|||||
Й |
- |
границы |
З обер и тел ь н ы хинтйрёалаё |
для |
|||||
|
|
средник |
а н а ч е н ш і; |
|
и н т е р в а л о в Зля |
||||
іу - |
гр ан и ц ы |
доверительных |
|||||||
|
|
и н З и ё и З у а л ь н ы х |
|
н а й л ію З е н и і. |
|
Рис. 81. Попе корреляции и вычисление статистических характеристик.