Файл: Козин В.З. Методы исследований в обогащении учеб. пособие.pdf
ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 28.06.2024
Просмотров: 151
Скачиваний: 1
где ЗСс& |
и 3-сН |
- |
верхнее и га»жнее значение |
I -го |
||
|
|
|
|
Фактора; |
|
|
|
|
ftli |
- |
чиспо |
шагов. |
|
Число |
узлов |
сетки |
равно |
+ |
|
|
(1 - |
число |
факторов. |
'•'* |
|
||
Этот метод поиска весьма прост, но и весьма трудо |
||||||
емок. Действительно, |
если имеете« всего 3 фактора и |
каждый из них нужно испытать в 10 узлах /такова требу емая точность определения необходимого значения фак тора/, то необходимо выполнить 103 = 1000 опытов, что в .'подавляющем большинстве случаев нереально.
Однако, этот случай следовало бы рассмотреть, так как часто все же, пренебрегая точностью, для небольшою
числа факторов такой поиск планируют, ибо он исключитель но прост. Мы не рекомендуем этот метод дли использова ния.
Более эффективным является локальный поиск на сетке. От предыдущего он отличается тем, что опыты осуществ ляются не все, а лишь те, которые способствуют движению к дели, другими словами выбор следующего варианта опыта доплееч зависеть от результатов предыдущих /обычно со седних/ опытов.
Отметим, что число необходимых опытов по принципу локального поиска на сетке зависит от выбора начальной тички. При удачное выборе это число может быть весьма малым /в пределах - 5/, в менее удачных случаях это чиспо увеличивается, но все равно остается меньшим по сравнению с перебором на сетке.
Прибііи чнтепьно число шагов в локальном поиске на сетке /в сам JM неблагоприятном случае/ может быть оце нено так.
' Пусть в одномерном опыте придется пройти всю ось, тогда чиспо опытов
вдвумерном /пусть зигзагообразно/ мы пройдем
т± + т Р + 1
ер.т.чмальных на дачном іпаге опытов и по 4 = 2 к 2 соседних опытов дпн каждой оптимальной точки.В трех
мерном |
|
|
|
|
|
|
tn i + m t + m 5 + |
3 и |
по 6 = 2 ?: 3 для каждой |
из |
этих |
||
точек н т.д. |
|
|
|
|
|
|
|
Всего надо сделать |
|
|
|
|
|
|
2. И |
( г п + і ) |
опытов |
|
|
|
|
Следовательно, |
для |
приведенного ранее примера |
В = 3 |
||
и |
ГН =1 —1C, |
необходимо |
выполнить всего |
18Û опы |
тов /по сравнению с 1000/ и то лишь в-самом неблагопри ятном случае.
Правда, это достигнуто следующей ценой: предполага ем, что найденный локальный экстремум совпадает с гло бальным, т.е. функция одноэкстремальна, В противном слу чае і.іы рискуем ошибиться.
Но выигрыш в числе опытов столь велик, что очень части идут на такое предположение. Однако, число опытов мож но сокращать и далее.
11.2Пропорциспальный поиск
Вени г процедуре покалыюго поиска на сетке принять не
постоянный шаг, а выбирать его также в зависимости от результатов предыдущих опытов, то выигрыш бѵдет боль шим.
Пусть шаг сетки будет пропорциональным производной функции на исследуем ом участке. Практически
|
+ ( х + & # -і)-4 -(х) |
|
|
- К |
'11.3/ |
Здесь мробь |
соответствует производной, |
вычисленной |
по приращениям функции н аргумента, ь |
К -выбранный |
|
нами коз|н«шиент |
дропоршюнапьно тк. |
|
Иногда этот метод называют методом градиента. Выбор К - сло/к ный процесс. Мочено вначале сказать,
что чем больше К - тем быстрее движение к экстремуму, но в случае неточных данных или неблагоприятных свойств функции /склоны с резко переменной крутизной/ большое К
может привести к "качке", причем процесс "качки" может стать далее расходящимся. Поэтому еще Ньютон предложил депать К переменным. Так, чем меньше ско рость изменения функции /вторая производная/, тем боль ше К можно взять /допустимы широкие шагп/, чем быстрее меняется характер функции, тем осторожнее необходимо двигаться.
Так, если произвести два измерения функции с шагом Л и 2Ù , то приблизительно молено вычислить вторую производную
t ± |
~ |
1 |
r»(X+U)~X(x+à) |
2(Х+*)~У(*П / 11.4/ |
|||
д х і |
" |
л |
L |
д |
|
д |
J |
Тогда^ взяв |
|
|
|
|
|
||
|
|
|
К ~ |
± |
Hit |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
дх*- |
менщу скоростью движе |
|
обеспечим разумный компромисс |
ния и сходимостью. /Знак минус берут при достижении максимума, плюс - минимума, чтобы обеспечить положи
тельность |
К |
/. |
|
|
выбиратьd |
Однако, в данной формуле необходимо снова |
|||||
Обычно рекомендуют |
= 1/2, |
но его можно |
увеличить |
||
почти до |
1. |
|
|
|
|
Практически |
экспериментатор |
выбирает вначале доста |
точно большое К с целью быстрейшего движения и при наступлении "качки" начинает его снижать, руководству ясь интуицией. Правда, многое зависит от стоимости и процолж ительности опытов. И всегда можно подсчитать, что дает выигрыш в затратах и времени - некоторый пере
расход в числе опытов за счет незнания второй производном, либо более точный*расчет правда также при перерасходе в числе опытов.
Наконец, заметим, что если бы мы могли эксперимен
тальным путем |
установить вид и коэффициенты функции |
||
У , то |
для определения экстремума достаточно было бы |
||
решить |
уравнение |
—Q |
|
Это |
такліе |
один из |
возможных вариантов поиска. |
11.3.Метод стохастической аппроксимации
Двигаясь к экстремуму в условиях высокого уровня помех /больших ошибок игмерения/ вполне возможны ошибки в определении градиента , в особенности ' вблизи экстремума, когда частные производные преближаются к нулю.
Например, пусть ошибка измерения рагна Z .тогда в исходной точке будет получено значение
= ÿ ( x o ) + £ o .
После осуществления шага, будет измерено
у , * ( х 0 + л ) = у ( х о + й ) + £ і .
Тогда вместо истинного значения производной
Ц- = т ^ (зс»+й)- ^ х”)]
получим
Ц - = |
|
|
|
^ + - 8 ^ . |
В зависимости |
от величины |
можно двинуться |
||
даже в противоположном |
направления. |
|
|
|
Конечно, пока |
Ж - |
велика по сравнению с ошибкой |
||
|
Эх |
|
|
|
в среднем движение будет к экстремуму /хотя и более |
||||
медленным, чем при отсутствии ошибок/, |
но как только |
|||
станет выполняться условие |
|
|
||
|
|
i t ~ |
д |
° Ѵ |
движение прекрати гоя и экстремума мы |
не достигнем. |
Процесс поиска экстремума превратится в колебательные движения со случайной амплитудой вокруг экстремума.
Рассмотрим |
формулу, иллюстрирующую это явление J-8J, |
Пусть ^ = |
—0,5^ X z |
Максимум функции лежит в точке
ij* ——L х ~ 0 при х = 0 .
Правило |
пропорционального |
поиска |
таково |
|
||||||||
X t+1 = |
|
+ Ѳ[-Lxt + |
= X± (-І- SL) + QS,± . |
|
||||||||
Здесь |
|
|
- |
ошибка |
в определении производной. |
|||||||
Для |
ряда последовательных шагоЬ' получим |
|
||||||||||
|
|
X L = Х0 ( і - B L ) + Q £ , o } |
|
|
||||||||
|
|
X z = Х 0 ( i - 8 L ) l + Ѳ - ( І - 0 U £ o + Ö £ i ; |
|
|||||||||
JCW = Хл (i - QL)'V+ [ Ѳ ( i - |
|
|
+ ... + B ( i S L ) C Z+ Ö ^ J |
|||||||||
Так как экстремум в начале координат, то значение Х^ И |
||||||||||||
есть и отклонение от экстремума. Если |
І1~Ѳ^І< І |
то |
||||||||||
при J / — |
|
|
( i - Q L ) " — |
0 |
. Но второе слагаемое |
|||||||
/в квадратных скобках/ ведет себя по законам случайных |
||||||||||||
чисел» Правда, |
благодаря тому, |
что |
£, |
= Ü, сумма |
в |
|||||||
квадратных скобках также стремится к нулю, однако, под |
||||||||||||
считаем |
дисперсию. |
|
|
|
|
|
|
|
||||
Дисперсия суммы случайных чисел равна сумме диспер |
||||||||||||
сий, тогда |
|
|
|
|
|
|
|
|
|
|
||
|
D |
M |
|
= |
e * [ ( i - e J L ) * (/'~i)'Di: + . . . + ( i - B L ) i Df: + D j . . |
|||||||
Считая, что |
(i~BL)>0 |
и применяя формулу суммы |
||||||||||
бесконечно убывающей геометрической прогрессии получим |
||||||||||||
|
D |
M |
|
= |
в * Р & |
_ _ Ѳ П ь |
|
|
||||
|
|
- (i |
e |
|
Z L t i - * £ ) |
|
||||||
|
|
|
|
|
1 - U - |
ôi)t)* |
|
|
||||
Но дисперсия производной на каждом шаге равна |
|
|||||||||||
0,'с'° " |
TW |
« |
|
а д » |
°Л) |
• |
|
|
||||
|
|
|
и |
ІХ М) - |
|
|
|
|
||||
Следовательно, ошибка |
в определении экстремума |
|
||||||||||
тем больше, |
чем |
больше |
/что очевидно \-(\ тем |
|
||||||||
больше, |
чем |
больше |
Ѳ |
, т.е, |
чем |
быстрее мы; желаем |
||||||
прийти |
к .-экстремуму. |
Тем |
больше |
чем |
меньше Л |
,т.е. |