Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
6 семестр.doc
Скачиваний:
0
Добавлен:
01.04.2025
Размер:
331.26 Кб
Скачать

Если же число вычислений значений функций 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% больше, чем в методе Фибоначчи, однако организация вычислительного процесса значительно проще.

Числа Фибоначчи определяются соотношениями

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. Для метода дихотомии =1/2 + /(b-a)2.

  2. Для метода золотого сечения =(  -1)/2=0.618033989....

  3. Для метода Фибоначчи =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. у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.

  2. у1> у0 – значение функции возрастло. В этом случае полагаем, что начальной точкой вычислений является точка x0**=x1,меньшим шагом в продолжении счета – величина

h*=-h/к, где к- некоторое целое положительное число, к 2. Далее производим вычисления по схеме 1. или 2., вплоть до достижения заданной точности.

  1. у10. В этом практически маловероятном случае естественно либо принять x*=(x01)/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; х10+h=0.45; ; ;

;

Здесь при вычислениях с шагом от точки =0.65 до точки =0.5 происходило уменьшение значений целевой функции, а в точки =0.45 значение функции возросло. Следовательно, можно принять минимальное значение функции достигается в точке x*= =0.5 c погрешностью 0.05.