- •Методы оптимизации
- •Методы оптимизации
- •1. Основные понятия и определения.
- •1.1 Место и роль методов оптимизации при моделировании и решении прикладных задач.
- •1.2 Общая постановка прикладной задачи.
- •1.3 Классификация оптимизационных методов и моделей.
- •1.4 Основные этапы построения оптимизационных моделей.
- •2. Линейное программирование.
- •2.1 Пример постановки задачи линейного программирования.
- •2.2 Общая постановка задачи линейного программирования.
- •2.3 Методика составления экономико-математических моделей.
- •2.4 Геометрическая интерпретация решения злп.
- •Выполнить самостоятельно. Задание №1.
- •2 .6 Симплекс-метод (см) решения злп
- •2.7 Отыскание опорного и оптимального решения злп с использованием табличного алгоритма с заменой базисных переменных.
- •Выполнить самостоятельно.
- •2.8 Двойственные задачи линейного программирования (дзлп).
- •2.9 Алгоритм составления двойственных злп:
- •3. Теорема двойственности
- •3. Целочисленное программирование.
- •3.1 Постановка задачи целочисленного программирования (зцп)
- •3.2 Метод отсечения (метод Гомори).
- •3.3 Алгоритм решения злцп
- •4. Нелинейное программирование.
- •4.1Классические методы оптимизации.
- •4.2 Метод множителей Лагранжа.
- •5. Методы направленного поиска экстремумов.
- •5.1 Методы определения экстремума унимодальной функции.
- •5.2 Методы определения локального экстремума функции нескольких переменных
2.7 Отыскание опорного и оптимального решения злп с использованием табличного алгоритма с заменой базисных переменных.
Алгоритм составления симплексных таблиц (СТ), рассмотрим на примере решения задачи отыскания max.
Пример 2.6
Линейная функция:
F=2x1+3x2→max
Система ограничений
x 1+3x2≤18
2x1+x2≤16
x2≤5
3x1≤21
x1, x2≥0
Приведем эту систему к каноническому виду, введя дополнительные переменные х3, х4, х5, х6:
х1+3х2+х3=18
2х1+х2+х 4=16
х2+х5=5
3х1+х6=21
Целевую функцию представим в виде:
F-2x1-3x2=0;
Заполняем первую симплексную таблицу: в ней х3, х4, х5, х6 – основные переменные (базис). Последняя строка называется оценочной.
Таблица №1.
Базис
Свободный член
Переменные
Оценочные отношения
х1
х2
х3
х4
х5
х6
х3
18
1
3
1
0
0
0
6
х4
16
2
1
0
1
0
0
16
х5
5
0
1
0
0
1
0
5
х6
21
3
0
0
0
0
1
F
0
-2
-3
0
0
0
0
←разрешающая строка
↑
разрешающий столбец
Проверяем выполнение критерия функции на max – первый опорный план не оптимальный, так как в F коэффициенты при x1 и x2 < 0
Выбираем наибольший по модулю отрицательный коэффициент F, который определяет разрешающий столбец.(второй столбец)
Делим свободные члены на коэффициенты разрешающего столбца, определяем оценочные отношения. И выбираем строку в качестве разрешающей, где это отношение минимальное min {6,16,5,∞}=5. На пересечении разрешающего столбца и разрешающей строки находиться разрешающий элемент a23 = 1.
Для построения таблицы №2 в качестве основной переменной мы выбираем х2, так как она образует разрешающей столбец таблице 1.
Переход к новому плану осуществляется пересчетом симплексной таблице (СТ) методом Жордана Гаусса.
Таблица №2.
Базис
Свободный член
Переменные
Оценочные отношения
х1
х2
х3
х4
х5
х6
х3
3
1
0
1
0
-3
0
3
х4
11
2
0
0
0
1
0
11/2
х2
5
0
1
0
0
1
0
х6
21
3
0
0
0
0
1
7
F
15
-2
0
0
0
3
0
← разрешающая строка
↑
разрешающий столбец
Построение 2ой таблицы:
1)Заменим переменные в базисе с х5 на х2.
2)Делим элементы разрешающей строки х5 (табл.1) на разрешающий элемент, результаты занесем в строку х2, но в таблицу №2.
3)В остальных клетках разрешающегося столбца (табл.1) записываем 0.
4)Остальные клетки заполняем по правилу прямоугольника:
НЭ=СТЭ-(АВ)/РЭ
НЭ- новый элемент.
СТЭ- старый элемент.
РЭ- разрешающий элемент.
А,В- эл-ты старого плана, образующие прямоугольник со старым эл-том и разрешающим элементом.
СТЭ А
В РЭ
b1=18-(3х5)/1=3
а11=1-(3х0)/1=1
b2=16-(1х5)/1=11 и т.д.
Критерии оптимальности опять не выполнен, так как F имеет коэффициент -2<0
- наибольший отрицательный по модулю коэффициент |-2| определяет разрешающий столбец x1
- min (3;11/2;∞;7)=3. Следовательно, 1ая строка разрешающая а11 - разрешающий элемент.
Таблица №3
Базис
Свободный член
Переменные
Оценочные отношения
х1
х2
х3
х4
х5
х6
х1
3
1
0
1
0
-3
0
3/-3=-1
х4
5
0
0
-2
1
5
0
5/5
х2
5
0
1
0
0
1
0
5/1
х6
12
0
0
-3
0
9
1
12/9
F
21
0
0
2
0
-3
0
← разрешающая строка
↑
разрешающий столбец
В таблице 3 критерий оптимальности вновь не выполнен. Разрешающий столбец x5, разрешающая строка x4, разрешающий элемент 5.
1)х1 вместо х3.
2)В строке х1 делим все на 1.
Таблица №4
Базис |
Свободный член |
Переменные |
Оценочные отношения |
|||||
х1 |
х2
|
х3
|
х4
|
х5
|
х6
|
|||
х1 |
6 |
1 |
0 |
-1/5 |
3/5 |
0 |
0 |
|
Х5 |
1 |
0 |
0 |
-2/5 |
1/5 |
1 |
0 |
|
х2 |
4 |
0 |
1 |
2/5 |
-1/5 |
0 |
0 |
|
х6 |
3 |
0 |
0 |
3/5 |
-9/5 |
0 |
1 |
|
F |
24 |
0 |
0 |
4/5 |
3/5 |
0 |
0 |
|
Новая СТ№4 – критерий оптимальности выполнен – оптимальное базисное решение X(6,4,0,0,1,3)
F=24 =max
Вспоминая экономический смысл всех переменных, логично сделать следующие выводы.
Прибыль принимает максимальное значение Fmax=24 при реализации 6 единиц продукции P1 (x1=6) и 4 единиц продукции P2(x4=4). Дополнительные переменные x3, x4, x5, x6 показывают остатки ресурсов каждого вида. При оптимальном плане производства x3=x4=0, то есть остатки ресурсов S3 и S4 равны 0, а остатки ресурсов S5 и S6 равны соответственно 1 и 3 единицам.