- •Вычислительная математика
- •210100 "Управление и информатика в технических системах"
- •680035, Хабаровск, ул. Тихоокеанская, 136.
- •Прямые методы решения систем линейных алгебраических уравнений Метод Гаусса
- •Метод Гаусса с выбором главного элемента по столбцам
- •Итерационные методы решения систем линейных алгебраических уравнений
- •Нормы матриц
- •Итерационные методы решения систем линейных алгебраических уравнений Метод простых итераций
- •Задания
- •Описание программы, входных и выходных параметров
- •Метод Зейделя
- •Контрольные вопросы
- •Лабораторная работа 2 решение нелинейных уравнений
- •Решение алгебраических уравнений методом половинного деления
- •Метод отделения корней
- •Метод итераций
- •Метод поразрядного приближения
- •Задания
- •Описание программы, входных и выходных параметров
- •Контрольные вопросы
- •Точечное среднеквадратичное приближение функций многочленами Чебышева
- •Задания
- •Контрольные вопросы
- •Лабораторная работа 4 интерполирование функций
- •Интерполяционная формула Лагранжа
- •Интерполяция кубическими сплайнами
- •Задание
- •Контрольные вопросы
- •Лабораторная работа 5
- •Описание программы, входных и выходных параметров
- •Контрольные вопросы
- •Лабораторная работа 6 Численное интегрирование
- •Квадратурные формулы интегрирования
- •Метод трапеций
- •Метод прямоугольников
- •Метод Симпсона (парабол)
- •Квадратурная формула Гаусса
- •Задание
- •Лабораторная работа 7 численное решение дифференциальных уравнений
- •Методы Рунге-Кутта
- •Метод Эйлера
- •Метод Эйлера-Коши
- •Метод Рунге-Кутта 4-го порядка
- •Задание
- •Контрольные вопросы
- •Лабораторная работа 8 поиск экстремумов функции
- •Сужение промежутка унимодальности
- •Численные методы поиска минимума функции одной переменной Метод половинного деления
- •Метод золотого сечения
- •Метод сканирования
- •Задание
- •Контрольные вопросы
- •Библиографический список
- •Оглавление
Численные методы поиска минимума функции одной переменной Метод половинного деления
Пусть
определен отрезок
,
которому принадлежит точка локального
минимума
,
и функция
является унимодальной на этом отрезке.
Используем
точки, расположенные симметрично
относительно середины отрезка
:
,
где
.
Тогда точки
и
принадлежат отрезку
,
и, следуя рассмотренной выше схеме
сужения промежутков унимодальности,
получим новый суженный отрезок
и оценим его длину
в каждом из трех возможных случаев:
.
.
.
Таким
образом, после первого шага преобразований
найден новый отрезок унимодальности,
длина которого уменьшилась. Название
метода (половинного
деления)
мотивировано тем, что если величина
очень мала, то отрезок уменьшается почти
вдвое (в случаях 1 и 2).
Вычисления
продолжаются до тех пор, пока не будет
достигнута требуемая точность
.
Метод золотого сечения
Точка
являетсязолотым
сечением
отрезка
,
если отношение длины
всего отрезка к длине
большей части равно отношению большей
части к длине меньшей части
,
т. е.
.
Аналогично, точка
,
симметричная точке
относительно середины отрезка, является
вторым золотым сечением этого отрезка.
Так
как точки
и
симметричны относительно середины
отрезка
,
то можно записать
,
(8.1)
где
.
Свойство золотого сечения: точка
одновременно является золотым сечением
отрезка
,
а другая точка
золотым сечением отрезка
.
При поиске минимума функции используется следующий алгоритм:
На исходном отрезке
по формуле (8.1) при
найдем точки
и
,
а затем разность
.Вычисляются значения функции
и
,
и по схеме сужения промежутка
унимодальности образуется суженный
отрезок
.На полученном отрезке
находятся два сечения
и
.
При этом возможны три случая:
.
.
.
По приведенной схеме находятся отрезки
,
и т. д., с учетом того, что в случаях 1) и
2) значение
или
целевой функции уже получено на
предыдущем шаге
.Точность приближенного равенства
на
-м
шаге вычислений можно оценить неравенством
,
где
.
Метод сканирования
Пусть
функция
является унимодальной на некотором
промежутке. Предположим, что произвольная
точка
этого промежутка является исходной для
поиска точки локального минимума и
число
- заданная точность нахождения
.
Обозначим через
произвольное приращение аргумента
и, сделав один шаг от точки
,
получим новое значение аргумента
.
Сравним
возможные значения функции
и
.
Возможны три варианта продолжения
приближения к точке минимума.
произошло
уменьшение значения функции. В качестве
нового стартового значения принимается
.
Вычисления по этой схеме продолжаются
до тех пор, пока не произойдет увеличение
значения функции, т. е.
,
и если при этом
,
то принимаем
с погрешностью
.
В противном случае полагаем, что точка
является исходной для продолжения
вычислений по следующей схеме 2.
значение
функции возросло. В этом случае полагаем,
что начальной точкой вычислений является
точка
,
а меньшим шагом для продолжения
вычислений – величина
,
где
- некоторое целое число,
.
Далее производим вычисления по схеме
1 или 2 до достижения требуемой точности.
(маловероятно).
Принимается либо
при достижении требуемой точности
,
либо следовать схеме 2.
Поиск минимума функции одной переменной данным методом представляет собой колебательный процесс около точки локального минимума функции с непрерывно уменьшающейся амплитудой.
Для
определения точки локального минимума
метод сканирования применим без
предварительного нахождения промежутков
унимодальности функции. С помощью метода
сканирования можно найти различные
точки локального минимума, если целевая
функция
имеет не один минимум.
