Файл: Козин В.З. Методы исследований в обогащении учеб. пособие.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/

Здесь мробь

соответствует производной,

вычисленной

по приращениям функции н аргумента, ь

К -выбранный

нами коз|н«шиент

дропоршюнапьно тк.

 

Иногда этот метод называют методом градиента. Выбор К - сло/к ный процесс. Мочено вначале сказать,

что чем больше К - тем быстрее движение к экстремуму, но в случае неточных данных или неблагоприятных свойств функции /склоны с резко переменной крутизной/ большое К


может привести к "качке", причем процесс "качки" может стать далее расходящимся. Поэтому еще Ньютон предложил депать К переменным. Так, чем меньше ско­ рость изменения функции /вторая производная/, тем боль­ ше К можно взять /допустимы широкие шагп/, чем быстрее меняется характер функции, тем осторожнее необходимо двигаться.

Так, если произвести два измерения функции с шагом Л и , то приблизительно молено вычислить вторую производную

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 +

= (-І- 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

«

 

а д »

°Л)

 

 

 

 

 

и

ІХ М) -

 

 

 

 

Следовательно, ошибка

в определении экстремума

 

тем больше,

чем

больше

/что очевидно \-(\ тем

 

больше,

чем

больше

Ѳ

, т.е,

чем

быстрее мы; желаем

прийти

к .-экстремуму.

Тем

больше

чем

меньше Л

,т.е.