Тема 6. Приближенные методы поиска экстремальных значений функций одной переменной.

Функция y = f(x) имеет в некоторой окрестности точки x^* локальный минимум, если для любой точки x, отстоящей от нее не более чем на некоторую величину  (т.е. | x - x^*| < ), справедливо неравенство f(x) < f(x^*). Функция y = f(x) имеет в некоторой окрестности точки x^** глобальный минимум на некотором интервале [a, b], если для любой точки x из этого интервала (т.е. x  [a, b]), справедливо неравенство f(x) < f(x^**). Точку достижения функцией глобального минимума можно определить посредством выбора из всех точек локальных минимумов или одной из границ интервала той, в которой значение функции будет минимальным. Аналогично можно опредилить локальный и глобальный максимум функции одной независимой переменной. Точки минимума и максимума функции имеют общее название - точки экстремума функции. Ограничимся рассмотрением только дифференцируемых функций.

Необходимым условием наличия в некоторой точке локального экстремума функции является равенство нулю ее производной. Это условие не является достаточным, поэтому такие точки являются только подозрительными на наличие локального экстремума. Достаточным условием наличия в них экстремума является тот факт, что этой точке вторая производная функции будет отличной от нуля (или при переходе через эту точку ее первая производная меняет знак). При этом

• функция имеет минимум, если ее вторая производная в подозрительной на экстремум точке положительна (или в этой точке ее первая производная меняет знак с “+” на “-“)

• функция имеет максимум, если ее вторая производная в подозрительной на экстремум точке отрицательна (или в этой точке ее первая производная меняет знак с “-” на “+“).

Если вторая производная в подозрительной на экстремум точке равна нулю (или в этой точке ее первая производная не меняет знак) - для решения вопроса о наличии в ней экстемума необходимы более сложные исследования.

Поскольку искать корни призводной функции и определять в них знаки вторых ее производных не совсем просто - имеются численные методы приближенного нахождения как глобальных, так и локальных экстремумов функций (с заданной точностью). Для выполнения поиска локальных экстремумов они требуют предварительного выбора интервала такого поиска, т.е. определения границ интервала, внутри которого этот экстремум содержится. А именно,

• для поска локального минимума необходимо выбирать интервал, на котором вторая производная функции знакопостоянная и положительная

• для поска локального максимума необходимо выбирать интервал, на котором вторая производная функции знакопостоянная и отрицательная.

В противных случаях - достаточно проверять конечные точки интервала.

Если вторая производная внутри интервала не является знакопостоянной - применение численных методов затруднено (точне полученный результат может оказаться неверным). Следует также отметить, что задаваемое пользователем значение точности для поиска экстремума тем или иным численным методом должно быть на порядок большим точности вычисления значения функции.

Метод сетки предназначен для поиска глобального экстремума функции на заданном интервале. Суть метода состоит в том, что весь интервал [a, b] покрывается “сеткой” - множеством точек, удаленных друг от друге на более чем на величину заданной точности (>0), а затем в каждой из этих точек производится расчет значений функции. Далее из них выбирается минимальное (или максимальное), которое (вместе с соответствующей ей абсциссой (т.е. значением переменной x) и будет найденным значением глобального экстремума заданной функции на заданном интервале. Выбор шага разбивки интервала производится так, чтобы выполнялось неравенство а это возможно при .

Метод равномерного поиска предназначен для поиска локального экстремума функции на заданном интервале. Суть метода состоит в том, что весь интервал [a, b] аналогичо предыдущему методу покрывается “сеткой” - множеством точек, удаленных друг от друге на более чем на величину заданной точности (>0), а затем начиная с левого конца заданного интервала (точки x=a) последовательно в каждой из этих точек производится расчет значений функции и сравнение его с расчитанным на предыдущем шаге. Если значения не убывают (при поиске максимума) или не возрастают (при поиске минимума) процесс вычислений повторяют до достижения конца интерваая. Если же тенденция изменения значений нарушается - процесс перебора точек прекращают и предыдущую точку принимают в качестве решения задачи.

Метод последовательного приближения предназначен для поиска локального экстремума функции на заданном интервале. Способ дости-жения цели - в некотором смысле аналогичен ранее описанному методу равномерного поиска. Но имеются два существенных отличия

• шаг поиска (приближения) не является постоянным, а изменяется от некоторого начального значения до значения, меньшего заданной точности (по абсолютной величине) по мере приближения к искомому значению x. Начальное значение шага выбирается равным , где значение k выбирается таким, чтобы было выполнено двойное неравенство 5  k  10. В процессе приближения к точке экстремума как знак шага, так и его величина будут изменяться

• начиная с левого конца заданного интервала (точки x=a) последовательно в каждой из этих точек производится расчет значений функции и сравнение его с расчитанным на предыдущем шаге. Если значения не убывают (при поиске максимума) или не возрастают (при поиске минимума) процесс вычислений повторяют до достижения конца интервала. Если же тенденция изменения значений нарушается - производят изменение величины и знака шага, полагая

,

и из последней точки перебор продолжают. Если перед очередной сменой значения шага его абсолютная величина окажеться меньше заданной точности, процесс перебора точек прекращают и предпоследнюю точку принимают в качестве решения задачи.

Процес поиска значения экстремума этим методом удобно проводить с занесением результатов расчетов в таблицу следующего вида

i	x	f(x)	h
0	x0=a		h0
1	x1=x0+h
…

Подробный алгоритм поиска экстремума по этому методу приведен в на рис. 7.1.