Математические задачи исследования операций, которые основаны на нелинейном программировании
Пример 4: Алгоритмом штрафных функций решить задачу
Этот пример показан на рис. 1.9. Поиск проводится из начальной точки (–2;–7) с параметрами a0=0,1 и b=2.
Здесь, ка
к и в предыдущем примере, на первой итерации алгоритма движение направлено в сторону безусловного минимума целевой функции. Это объясняется небольшим начальным значением параметра штрафа, при котором основное влияние на направление оказывает целевая функция. С возрастанием a направление поиска ориентируется на условный экстремум.
Пример 5: Алгоритмом штрафных функций решить задачу
Этот пример приведен на рис. 1.10. Поиск производился при a0=1 и b=10. Такие параметры обусловили другой характер движения к условному минимуму: первая итерация уже не приводит в окрестность безусловного экстремума и траектория не имеет резких изменений направления поиска.
Таким образом, выбор параметров поиска имеет существенное влияние на эффективность алгоритма.
5. Метод барьерных функций
В отличие от метода штрафных функций, данный метод применим к задачам с ограничениями только в виде неравенств.
|
Исходная задача на условный экстремум задается в виде
f(x) → min; (1.25)
ji(x) £ 0, . (1.26)
Она преобразуется в задачу безусловной минимизации вспомогательной функции
Q(x) = f(x) + mB(x),
где B(x) – барьерная функция, m – параметр барьера.
Обязательное условие: внутренность области не должна быть пустой (имеются точки, в которых "ji (x) < 0).
Барьерная функция строится так, чтобы она была неотрицательной и непрерывной на допустимом множестве и стремилась к бесконечности при приближении изнутри к границе:
Как и в случае штрафной функции, существует несколько конструкций B(x), удовлетворяющих этим условиям. Но в основном используется барьерная функция в виде
(1.27)
Понятно, что решение вспомогательной задачи зависит от значения параметра барьера. Покажем влияние m на результат минимизации Q.
Примеры
Пример 1: f(x) = x → min; j(x)=3 – x £ 0.
Барьерную функцию строим согласно (1.27). Тогда вспомогательная функция имеет вид Находим точку минимума Q: Отсюда
Следовательно, с уменьшением m точки минимума вспомогательной функции приближаются к минимуму исходной задачи. Геометрическая иллюстрация решения приведена на рис. 1.11.
Как хорошо видно на рис. 1.11, при оптимуме исходной задачи на границе допустимого множества последовательность точек минимума Q с уменьшением m приближается к оптимальному решению изнутри допустимой области. По этой причине метод барьеров называют методом внутренних штрафов.
В связи с возможными трудностями поиска при малых значениях m решается не одна, а последовательность вспомогательных задач с уменьшающимися значениями параметра барьера.
Алгоритм.
1. Выбрать начальную точку x0 так, чтобы "ji(x0)<0; задать точность e, начальное значение m0 и число b Î (0, 1).
2. Минимизировать Q(x) одним из методов безусловной оптимизации, в результате чего определяется .
3. Проверить: если , то остановиться, приняв за оптимальное решение задачи.
4. Положить , за начальную точку принять и вернуться к шагу 2.
Значение m0 можно брать из интервала [2, 10]. Важное замечание касается пункта 2 алгоритма: в процессе поиска минимума вблизи границы из-за дискретности шагов возможен выход за допустимую область, где барьерная функция становится отрицательной, что повлечет расхождение поиска. Поэтому необходима явная проверка на допустимость точек на каждом шаге при минимизации Q.
Пример 2: f(x) = (x1–2)4+( x1–2x2)2® min; j(x)=x2£ 0.
Решение находим, используя вспомогательную функцию
Q=(x1–2)4+(x1–2 x2)2 –
За начальную точку возьмем допустимую точку (0;1), значения m и b принимаем равными 10. Результаты поиска алгоритмом барьерных функций представлены в табл. 1.2 и на рис. 1.12.
Таблица 1.2
Результаты поиска алгоритмом барьерных функций
№ итерации |
m |
x1 |
x2 |
f |
Q |
mB |
1 |
10 |
0.7079 |
1.5315 |
8.3338 |
18.0388 |
9.705 |
2 |
1.0 |
0.8282 |
1.1098 |
3.8214 |
6.1805 |
2.3591 |
3 |
0.1 |
0.8989 |
0.9638 |
2.5282 |
3.1701 |
0.6419 |
4 |
0.01 |
0.9294 |
0.9162 |
2.1291 |
2.3199 |
0.1908 |
5 |
0.001 |
0.9403 |
0.9011 |
2.0039 |
2.0629 |
0.0590 |
6 |
0.0001 |
0.94389 |
0.89635 |
1.9645 |
1.9829 |
0.0184 |
Другие рефераты на тему «Экономико-математическое моделирование»:
- Статистические методы обработки данных
- Проектирование модели для составления оптимального рациона кормления скота
- Оптимизация сетевой модели комплекса производственных работ
- Математические задачи исследования операций, которые основаны на нелинейном программировании
- Прогнозирование и регулирование развития производственной инфраструктуры
Поиск рефератов
Последние рефераты раздела
- Выборочные исследования в эконометрике
- Временные характеристики и функция времени. Графическое представление частотных характеристик
- Автоматизированный априорный анализ статистической совокупности в среде MS Excel
- Биматричные игры. Поиск равновесных ситуаций
- Анализ рядов распределения
- Анализ состояния финансовых рынков на основе методов нелинейной динамики
- Безработица - основные определения и измерение. Потоки, запасы, утечки, инъекции в модели