
- •Министерство общего и профессионального образования
- •Российской Федерации
- •Алтайский Государственный Университет
- •Математический факультет
- •Кафедра ткпм
- •Барнаул,1999 год п.1. Дихотомия
- •Варианты заданий.
- •П. 2. Метод золотого сечения
- •Если же число вычислений значений функций f(X) заранее жестко
- •П.3.Метод фибоначчи.
- •Числа Фибоначчи определяются соотношениями
- •П.4. Метод сканирования.
- •Упражнения
- •П.5. Метод градиентного спуска.
- •Варианты заданий:
- •П.6.Метод сопряженных градиентов.
- •П.7.Упражнения.
- •Упражнение10. Доказать, что выпуклая на отрезке [a,b] функция f(х), отличная от константы, не может достигать своей верхней грани внутри отрезка [a,b].
- •Упражнение17. Применить метод штрафных функций к задачам
- •Литература
Если же число вычислений значений функций f(X) заранее жестко
задано и равно n,
то процесс на этом заканчивается.В
качестве решения задачи можно принять
пару x*n,f(x*n),
при этом (x*n,xn)
max{bn-x*n
,x*n-an}=(
-1)(b
n
-a
n)/2=( (
-1)/2)n(b-a )=A
n.
С помощью дихотомии за n=2k вычислений значений функций мы получили точку x*n с погрешностью
(x* n ,X* ) ≤ 2–n/2(b-a- ) < 2-n/2 (b-a)=B n.
Отсюда имеем A
n /B
n=2
(
+1))n=(0.87)
n.
Очевидно,что уже при небольших n (числе итераций) преимущество золотого сечения перед дихотомией становится ощутимым.
Порядок выполнения лабораторной работы и варианты задания такие же, как в методе дихотомии.
П.3.Метод фибоначчи.
Можно показать,что для решения задачи одномерной минимизации оптимальным является метод Фибоначчи основанный на использованиизнаменитых чисел Фибоначчи. При достаточно большом количестве итераций окончательный интервал [a n ,b n ] интервал неопределенности в мет оде золотого сечения лишь на 17% больше, чем в методе Фибоначчи, однако организация вычислительного процесса значительно проще.
Числа Фибоначчи определяются соотношениями
F 1=1; F2=2; Fn+2=Fn+1 + F n,n=1,2,3,....
Метод Фибоначчи относится к классу симметричных методов и определяется заданием на отрезке [a,b] точки x1=a+(b-a)F n /Fn+2 или симметричной ей точки x2=a+b-x1=a+(b-a)Fn+1/Fn+2. Описанные выше методы допускают естественное обобщение.Начальное разбиение определяется точками
x1=a + (1- )(b-a), x2=a + (b-a).
Для метода дихотомии =1/2 + /(b-a)2.
Для метода золотого сечения =( -1)/2=0.618033989....
Для метода Фибоначчи =Fn+1/Fn+2,1- =F n / Fn+2.
Утверждение1.Количество необходимых при решении задачи минимизации вычислений значений унимодальной функции,гарантирующих достижение точности, равно числу n, удовлетворяющему неравенствам (b-a)/Fn+2 <=ε <=(b-a)/Fn+1.
Порядок выполнения лабораторной работы и варианты задания такие же,как в методе дихотомии.
П.4. Метод сканирования.
Пусть функция y=f(x) является унимодальной на некотором промежутке. Предположим,что произвольная точка x0 этого промежутка является исходной для поиска точки x* локального минимума и число -заданная точность нахождения х*. Обозначим через h произвольное приращение аргумента х и, сделав один шаг от точки x0,получим новое значение аргумента х1=x0+h.
Сравним значения функции у0=f(x0) и у1=f(х1)=f(x0+h). Возможны три различных продолжения в приближении к точке х*.
у1< у0– произошло уменьшение значения функции. Тогда примем в качестве нового стартового значения x0(1)=х1 и сделаем шаг h от этой точки x0(1) к точке х1(1),т.е. х1(1)=x0(1) +h. Если окажется у1 (1) <у0(1),то снова сделаем шаг h от новой стартовой точки х0(2)=x1(1) и т.д. На некотором к-ом шаге произойдет увеличение значения функции, т.е. у1 (k) <у0(k), и если при этом h<, то принимаем х*x0(k) с погрешностью h. В противном случае полагаем,что точка x0**==x(k) является исходной для продолжения вычислений по следующей схеме 2.
у1> у0 – значение функции возрастло. В этом случае полагаем, что начальной точкой вычислений является точка x0**=x1,меньшим шагом в продолжении счета – величина
h*=-h/к, где к- некоторое целое положительное число, к 2. Далее производим вычисления по схеме 1. или 2., вплоть до достижения заданной точности.
у1=у0. В этом практически маловероятном случае естественно либо принять x*=(x0+х1)/2 при достижении заданной точности h либо следовать схеме 2.
Поиск минимума функции одной переменной указанным методом представляет собой колебательный процесс, совершающийся около точки х* локального минимума функции f(х) с непрерывно уменьшающейся амплитудой.
Описанный метод поиска локального минимума называют методом сканирования. Порядок выполнения лабораторной работы и варианты задания такие же, как в методе дихотомии.
Пример. Найти точку х* локального минимума функции f(x)=2x2-lnx методом сканирования с точностью =0.1, полагая х0=0.25; h=0.2 и K=4.
Решение.
Последовательность вычислений представим
в виде пар координат
,
где i=0;1,
а верхний индекс k
является номером очередного шага
сканирования. Имеем:
h=0.2>=0.1;
х0=0.25; х1=х0+h=0.45;
;
;
;
Здесь при вычислениях
с шагом
от точки
=0.65
до точки
=0.5
происходило уменьшение значений целевой
функции, а в точки
=0.45
значение функции возросло. Следовательно,
можно принять минимальное значение
функции достигается в точке x*=
=0.5
c
погрешностью 0.05.