- •Н. Б. Борковский Математические методы в экологии
- •Часть I. Численные методы
- •Тема 1 9
- •Тема 2 14
- •Тема 3 30
- •3.0 Постановка задачи 30
- •1. Постановка проблемы
- •2. Выбор или построение математической модели
- •3. Постановка вычислительной задачи
- •4. Предварительный анализ свойств вычислительной задачи
- •5. Выбор или построение численного метода
- •6. Алгоритмизация и программирование
- •7. Отладка программы
- •8. Счет по программе
- •9. Обработка и интерпретация результатов
- •10. Использование результатов и коррекция математической модели
- •1.0 Источники приближенных чисел
- •1.1 Источники и классификация погрешностей
- •1.2 Абсолютная и относительная погрешности вычисления
- •1.3 Погрешности арифметических операций
- •1.4 Обратная задача теории погрешностей
- •Постановка задачи
- •Метод половинного деления (метод «вилки», метод дихотомии)
- •Метод Ньютона (метод касательных)
- •Метод хорд (метод пропорциональных частей)
- •Метод итерации (метод последовательных приближений)
- •Постановка задачи
- •Решение систем линейных уравнений методом итерации
- •Решение систем нелинейных уравнений методом итерации
- •Решение систем линейных уравнений – «подводные камни»
- •4.0 Интерполирование – постановка задачи
- •4.1 Интерполирование полиномами
- •4.2 Аппроксимация функций. Метод наименьших квадратов
- •4.3 Интерполирование сплайнами
- •4.4 Интерполирование – примеры и «подводные камни»
- •4.5 Сглаживание данных (скользящее среднее)
- •5.0 Постановка задачи
- •5.1 Формула трапеций
- •5.2 Формула Симпсона
- •5.3 Формула Ньютона–Котеса
- •5.4 Формула Гаусса
- •5.5 Связь формул интегрирования
- •5.6 Численное интегрирование. Сравнение методов
- •5.7 Численное интегрирование – «подводные камни»
- •6.0 Введение
- •6.1 Обыкновенные ду. Метод Эйлера
- •6.2 Обыкновенные ду. Основные методы решения
- •6.3 Обыкновенные ду. Методы прогноза и коррекции
- •6.5 Обыкновенные ду. Сравнение методов
- •6.6 Оду. Влияние неточности исходной информации
- •6.7 Оду. Представление о конечных разностях
- •6.8 Оду. Краевые задачи – введение
- •8.0 Ду в частных производных. Постановка задачи
- •8.1 Ду в частных производных. Эллиптические уравнения
- •8.2 Ду в частных производных. Гиперболические уравнения
- •8.3 Ду в частных производных. Параболические уравнения
- •9.0 Методы безусловной оптимизации. Введение
- •9.1 Методы безусловной оптимизации. Классификация методов
- •9.2 Методы безусловной оптимизации нулевого порядка
- •9.3 Методы безусловной оптимизации первого порядка
- •9.4 Методы оптимизации второго порядка
- •Тема 10
- •10.0 Метод Монте–Карло. Некоторые задачи
- •10.1 Проблема получения случайных чисел
- •10.2 Общая схема метода
- •10.3 Метод Монте–Карло. Пример
- •10.4 Вычисление кратных интегралов
- •Число испытаний n не зависит от размерности интеграла i0
9.0 Методы безусловной оптимизации. Введение
Решение многих теоретических и практических задач сводится к отысканию экстремума скалярной функции f(x) n-мерного векторного аргумента x.
Под x будем понимать вектор-столбец:
Вектор-строка получается путем транспонирования:
Оптимизируемую функцию f(x) называют целевой функцией, или критерием оптимальности.
В дальнейшем без ограничения общности будем говорить об отыскании минимума функции f(x)
f(x) → min.
Вектор x*, доставляющий минимум целевой функции, называют оптимальным.
Задачу максимизации можно заменить эквивалентной задачей минимизации и наоборот.
Рассмотрим это на примере функции одной переменной. Если x* – точка минимума функции y = f(x), то для функции y = –f(x) она – точка максимума. Т. е. min f(x) = –max (–f(x)) (рис. 9.1).
|
Рис. 9.1. Экстремумы |
Сказанное справедливо и для функции многих переменных:
min (f(x1, …, xn )) = –max (–f(x1, …, xn )).
Так что далее речь будет идти только о минимизации.
Точка x* доставляет глобальный минимум функции одной переменной f(x), заданной на числовой прямой X, если x*X и f(x*) f(x) для всех x X.
Точка x* называется точкой строгого глобального минимума, если неравенство выполняется как строгое.
Если же в выражении f(x*) f(x) равенство возможно при x x*, то реализуется нестрогий минимум, а под решением понимают множество x* = {x X: f(x) = f(x*)} (рис. 9.2).
|
Рис. 9.2. Глобальный минимум: а – строгий; б – нестрогий |
Точка x* X доставляет локальный минимум функции f(x) на множестве X, если при некотором, достаточно малом > 0 для всех x x*, x* X, удовлетворяющих условию |x – x*| , выполняется неравенство f(x*) f(x).
Если неравенство строгое, x* является точкой строгого локального минимума.
Все определения для максимумов функции получаются заменой знаков в неравенствах на обратные.
На рис. 9.3. показаны экстремумы функции одной переменной f(x) на отрезке [a, b].
|
Рис. 9.3. Экстремумы функции. Точки x1, x3, x5 – локальные максимумы; точки x4, x6 – локальные минимумы; точка x2 – глобальный минимум; точка x7 – глобальный максимум |
Проиллюстрируем трудности оптимизации на примере очень простых функций.
Минимумы любого типа отсутствуют, когда функция не ограничена снизу
|
f(x) = x
|
|
Даже если она ограничена снизу, минимум может отсутствовать
|
f(x) = e –x
|
|
Возможно бесконечное число локальных и глобальных минимумов
|
f(x) = sin x
|
|
Возможно бесконечное число локальных и ни одного глобального минимума |
f(x) = x + 2sinx
|
|
Если функция или ее производная разрывна, ситуации могут быть очень своеобразные: глобальный минимум есть, а локального нет! В точке минимума производная |
|
|
В точке x = 0 производная
|
f(x) = x3
|
|
Седловая точка функции двух переменных – по одной переменной достигнут максимум, по другой – минимум |
z(x, y) = x2 – y2
|
|

!
,
а минимума нет