
- •Численные методы (Вычислительная математика) Санкт-Петербург
- •Оглавление
- •Глава 1. Элементарная теория погрешностей........................10
- •Глава 2.Численные методы решения уравнений……………24
- •Глава 3. Приближенное представление функций.................50
- •Предисловие
- •Введение
- •Глава 1. Элементарная теория погрешностей
- •1.1. Приближенные числа и их погрешности
- •1.2. Правила приближенных вычислений и оценка погрешностей при вычислениях
- •1.3. Неустранимая погрешность вычисления значения функции
- •1.4. Обеспечение заданной точности результата
- •Глава 2.Численные методы решения уравнений
- •2.1. Итерационные методы решения нелинейных уравнений
- •2.2. Метод половинного деления
- •2.3. Метод хорд
- •2.4. Метод касательных (метод Ньютона).
- •2.5. Комбинированный метод хорд и касательных
- •2.6.Метод итераций
- •2.7. Решение систем линейных уравнений методом итераций.
- •2.9. Метод итераций для нелинейных систем уравнений
- •Глава 3. Приближенное представление функций (аппроксимация, интерполяция)
- •3.1. Интерполирование функции алгебраическими многочленами
- •3.2. Единственность интерполяционного многочлена
- •3.3. Интерполяционная формула Лагранжа
- •3.4. Погрешность интерполяционной формулы Лагранжа
- •3.5. Погрешность линейной интерполяции
- •3.6.Конечные разности и их свойства
- •3.7.Интерполяционная формула Ньютона (для равностоящих узлов)
- •3.10. Эмпирические формулы и метод наименьших квадратов
- •Глава 4. Численное дифференцирование и интегрирование
- •4.1. Дифференцирование на основе интерполяционной формулы Ньютона
- •4.2. Дифференцирование на основе формулы Стирлинга
- •4.3. Дифференцирование на основе интерполяционной формулы Лагранжа
- •Глава 5. Численное интегрирование функций
- •5.1. Квадратурные формулы Ньютона-Котеса
- •5.2. Формула трапеций
- •5.3. Формула Симпсона (формула парабол)
- •5.4. Формулы Ньютона – Котеса высших порядков
- •Глава 6. Численное интегрирование обыкновенных дифференциальных уравнений
- •6.1. Метод Эйлера
- •6.2. Метод Рунге – Кутта
- •6.3. Метод Адамса
- •Глава 7. Численные методы оптимизации
- •7.1. Градиентные методы
- •7.2. Метод Франка – Вулфа
- •7.3. Метод наискорейшего спуска
- •7.4. Метод штрафных функций
- •Литература
7.3. Метод наискорейшего спуска
Метод предназначен
для приближенной минимизации функции
,
определенной и выпуклой во всем
пространстве, т.е. для приближенного
решения задачи без ограничений:
.
Описание алгоритма метода наискорейшего спуска.
Выбор начального приближения. Начальное приближение
выбирается произвольно, но выгоднее выбирать начальное приближение так, чтобы значение функции
было по возможности меньшим.
Условие оптимальности (критерий остановки алгоритма). Алгоритм заканчивает работу, если приближение
оптимально, т.е.
. При практическом применении метода задается малое число
, определяющее точность вычислений и тогда условие остановки алгоритма
или аналогичное ему условие
. При выполнении этого условия вычисления прекращаются и минимальное значение
и
.
Выбор направления спуска. За направление спуска выбирается направление наискорейшего убывания функции в точке , это направление, противоположной направлению градиента в этой точке. Следовательно, координаты вектора
, задающего направление спуска, определяются как частные производные в точке
:
.
Выбор шага спуска
. Для определения шага введем функцию шага
, значения которой совпадают со значениями вдоль направления, задаваемого вектором
. В таком случае она имеет вид
, где
. Шаг определяется из условия минимума функции , т.е.
. Приближенно эту точку можно найти как точку экстремума из условия
Построение очередного приближения. Очередное приближение вычисляется следующим образом:
. Вычисления продолжаются до тех пор, пока не выполнится критерий оптимальности.
Описываемый алгоритм спуска сходится, если исследуемая функция удовлетворяет следующим условиям:
Функция определена, выпукла и непрерывна во всем пространстве;
;
Функция имеет всюду непрерывные частные производные первого и второго порядков.
При выполнении
перечисленных условий (1)-(3) функция
имеет точку абсолютного минимума
.
Последовательные приближения
образуют минимизирующую последовательность,
т.е.
.
Если строго выпуклая функция
имеет единственную точку минимума
,
то приближения
сходятся к этой точке:
.
7.4. Метод штрафных функций
Рассмотрим задачу
минимизации функции
,
которая определена во всем пространстве
и является выпуклой во всей области
определения, при наличии ограничений.
Метод штрафных функций сводит решение
этой задачи к решению некоторой
последовательности задач без ограничений.
Для этого вводится функция
,
где M
– неотрицательный
параметр (
).
Функция
обладает следующими свойствами:
1.
,
для всех значений x,
входящих в область допустимых планов
;
2.
вне области планов
;
3.
если
:
4.
,
если
.
Введенная таким
образом функция
называется штрафом, сумма
является штрафной функцией задачи. В
области допустимых планов штрафная
функция совпадает с целевой функцией
,
вне ее
При этом исходная
задача заменяется задачей минимизации
штрафной функции без ограничений. Если
оптимальный план
,
полученный при минимизации штрафной
функции, лежит в области планов
исходной задачи, то он является ее
оптимальным планом, поскольку области
планов штрафной и исходной функций
совпадают. Если оптимальный план
задачи минимизации штрафной функции
не входит в область планов исходной
функции, то оптимальное значение штрафной
функции не больше оптимального значения
целевой функции исходной задачи
.
В основе алгоритма решения задачи минимизации методом штрафных функций лежит следующая теорема.
Теорема. Если
функция определена, выпукла, непрерывна
во всем пространстве и
,
а функция
непрерывна как функция переменной x,
то задачи минимизации исходной целевой
и штрафной функций имеют оптимальные
планы,
и
.