В современной интерпретации этот метод был предложен Кифером, хотя корни его восходят к математику XIII века Фибоначчи и даже к более ранним исследованиям Евклида.
Будем считать, что для исследования исходного интервала неопределенности единичной длины было выделено N экспериментов, N - 1 из которых выполнены. В последнем эксперименте требуется исследовать интервал неопределенности LN-1опт держащий точку-эксперимент Е+N-1 для которой получено самое большое значение y при первых N - 1 испытаниях (рис. 15-10). Положение точки E+N-1 целиком зависит от стратегии поиска. Последнюю точку EN можно выбрать произвольно, но не ближе, чем на расстоянии ε от точки EN-1. Ситуация на этом последнем шаге полностью совпадает с обычным поиском при наличии двух экспериментов. Минимальное значение последнего интервала LNопт получается, если точки ЕN-1+ и EN расположить симметрично на расстоянии ε/2 относительно середины интервала LN-1опт. Поэтому значения LN-1опт и LNопт связаны соотношением
LN-1опт=2LNопт-ε. (15-6)
Рис. 15-10. Пояснение к методу Фибоначчи
Продолжая процедуру движения с конца, аналогичную динамическому программированию, введем в рассмотрение интервал неопределенности LN-2опт, который получается в результате N-2 экспериментов из общего числа N. Внутри этого интервала содержится эксперимент Е+N-1 дающий наибольшее значение y при первых N-2 экспериментах. В этом интервале следует произвести следующую пару экспериментов. Лучшее значение из двух обозначим через ЕN-1+. Другая точка из этой пары, которую обозначим через ЕN-1-, будет граничной между частью LоптN-2 интервала LоптN, исключаемой и частью LN-2опт, оставляемой для рассмотрения. Но в начале эксперимента неизвестно, какое из двух значений EN-1 будет наибольшим и обозначится как ЕN-1+ а какое наименьшим и обозначится как Эти два значения могут меняться местами. Поэтому обе точки ЕN-1+ и ЕN-1- располагаются на одинаковом расстоянии LN-1опт от концов интервала. Так как эти точки расположены симметрично относительно середины интервала LоптN-2 и одна из них становится EN-1+ (причем это может быть любая из двух точек EN-1), то каждая из них должна находиться на расстоянии LNопт от одного конца интервала и на расстоянии LN-1опт от другого конца интервала. При этом учитывается, что ранее была доказана формула (15-6). Поэтому
LN-2опт=LN-1опт+LNопт
С учетом формулы (15-6) напишем:
LN-2опт=(2LNопт-ε)+LNопт=3LNопт-ε
Аналогично можно доказать для любого 1<k<N, что
Lk-1опт=Lkопт+Lk+1опт
Отсюда при k=N-2
Для удобства записи формул введем последовательность чисел Фибоначчи:
F0=F1=1
Fk=Fk-1+Fk-2, k=2,3...,
в соответствии с которой имеем:
F2=2, F3=3, F4=5, F5= 8 и т. д.
С помощью этих чисел можно записать следующую формулу для оптимальных интервалов неопределенности:
LN-kопт=Fk+1LNопт-Fk-1ε (15-7)
Полагая длину исходного интервала равной единице, получаем:
Отсюда находим оптимальный интервал после N опытов:
Как видно из приведенной позднее табл. 15-1, для уменьшения интервала неопределенности более чем в 100 раз по методу Фибоначчи необходимо провести 11 опытов, по методу дихотомии- 14, а при пассивном методе - 198 [см. формулу (15-5)]. При методе Фибоначчи каждая последующая точка выбирается симметрично по отношению к точке, которая осталась от предыдущего эксперимента и попала в оставшийся интервал. Первый эксперимент следует выбирать на расстоянии L2опт от одного из концов начального интервала, безразлично от левого или правого, в силу симметрии. Чтобы получить точное выражение для L2опт воспользуемся формулой (15-7), учитывая, что N-k=2 и k=N-2:
L2опт=FN-1LNопт-FN-3ε.
Исключив из этого соотношения LNопт, с помощью ^формулы (15-8) получим:
Используя тождества FN-1=FN-2+FN-3, FN=FN-1+FN-2 которые следуют из самого определения чисел Фибоначчи, напишем:
Можно доказать [Л. 92], что правая часть этого равенства равна (-1)n, поэтому формулу (15-9) можно переписать в виде
Формула (15-11) позволяет начать процедуру поиска, если известны ε и N. В этом состоит недостаток метода: требуется заранее знать необходимое число опытов. Очень часто задается только величина интервала неопределенности, которую надо знать к концу N-го опыта. Эта особенность придает последовательному методу Фибоначчи свойства пассивного поиска. Она напоминает взаимоотношения классической процедуры Неймана - Пирсона и последовательного метода Вальда в теории проверки статистических гипотез (см. гл. 3). Иногда поэтому обращаются к методу дихотомии, который хотя и обладает меньшей эффективностью по сравнению с методом Фибоначчи, но не требует предварительного знания числа опытов.
Однако существует другой метод поиска, который почти также эффективен, как метод Фибоначчи, но не требует априорного задания числа опытов N. Имеется в виду известный издавна метод золотого сечения.