3.3. Критерий оптимальности Критерий оптимальности полученного опорного решения:
если в последней строке симплексной таблицы – строке целевой функции - все оценки неотрицательны
p0,
то полученное опорное решение оптимально.
Оптимальное решение не существует, если в симплексной таблице есть хотя бы одна отрицательная оценка, столбец которой не содержит положительных элементов, в этом случае целевая функция f не ограничена в области допустимых решений:
Пример 9.
Проверить критерий оптимальности по заданной симплекс – таблице:
Базис
| х1
| х2
| х3
| х4
| х5
| bi0
| bi0 / bip
| x1
| 1
| 0
| 4/5
| -1/5
| 0
| 140
|
| x2
| 0
| 1
| -3/5
| 2/5
| 0
| 120
|
| x5
| 0
| 0
| 1
| -1
| 1
| 100
|
| f
| 0
| 0
| -2
| 3
| 0
| 13200
|
|
Решение.
Найденное в предыдущем примере опорное решение
Х=(140, 120, 0, 0, 100)
не будет оптимальным, так как в строке целевой функции есть отрицательная оценка 4= -2 0.
3.4. Невырожденные ЗЛП, алгоритм их решения Опорное решение ЗЛП называется невырожденным, если все базисные переменные строго больше нуля.
Задача называется невырожденной, если все ее опорные решения невырожденные.
Алгоритм решения невырожденной ЗЛП.
Приводим систему ограничений к единичному базису ( с помощью метода Жордана- Гаусса )
Находим приведенное выражение для целевой функции f и заполняем исходную симплексную таблицу.
Если все оценки p0, то делаем вывод, что полученное решение оптимально и процесс решения задачи окончен.
Хопт=(b10,b20, ... , br0,0,0, ... ,0),
fmax = b00.
Если есть отрицательная оценка p<0 и все элементы столбца над этой оценкой отрицательны или равны нулю, то это означает, что f неограниченна, max f и задача решения не имеет.
Если есть отрицательная оценка p<0 и хотя бы один из элементов столбца над ней неотрицателен bip>0, то переходим к лучшему опорному решению следующим образом:
выбираем ведущий p- столбец из условия: наибольшая по модулю оценка p<0 и хотя бы один из элементов bip>0;
выбираем ведущую q- строку по наименьшему симплексному отношению .
в симплексной таблице выбираем ведущий элемент bqp (элемент, стоящий на пересечении q- ой строки и p- того столбца) и выполняем итерацию метода Жордана- Гаусса. Таким образом переходим к новому базису с базисной переменной xp.
находим новое опорное решение и проверяем выполнения
критерия оптимальности.
Пример 10.
Найти оптимальное решение по заданной симплекс – таблице:
Базис
| х1
| х2
| х3
| х4
| х5
| bi0
| bi0 / bip
| x1
| 1
| 0
| 4/5
| -1/5
| 0
| 140
|
| x2
| 0
| 1
| -3/5
| 2/5
| 0
| 120
|
| x5
| 0
| 0
| 1
| -1
| 1
| 100
|
| f
| 0
| 0
| -2
| 3
| 0
| 13200
|
|
Решение.
Найденное в предыдущем примере опорное решение
Х1=(140, 120, 0, 0, 100)
не будет оптимальным, так как в строке целевой функции есть отрицательная оценка 4= -2 0.
Перейдем от исходного опорного решения Х1 к лучшему Х2.
В качестве ведущего столбца выберем третий столбец по отрицательной оценке 3=-2.
Для положительных элементов третьего столбца составим симплексные отношения и запишем их в последнем столбце симплексной таблицы.
В качестве ведущей строки возьмем строку с минимальным симплексным отношением- это третья строка, так как
==100,
таким образом, ведущий элемент b33=1.
С этим ведущим элементом проведем одну итерацию метода Жордана- Гаусса, а именно.:
все элементы ведущей строки разделим на ведущий элемент
(на 1);
в ведущем столбце все элементы кроме ведущего заменим нулями (запишем единичный столбец);
все остальные элементы симплекс-таблицы преобразуем по правилу « определителя второго порядка », получим
новую симплекс – таблицу:
Базис
| х1
| х2
| х3
| х4
| х5
| bi0
| bi0 / bip
| x1
| 1
| 0
| 0
| 3/5
| -4/5
| 220
|
| x2
| 0
| 1
| 0
| -1/5
| 3/5
| 180
|
| x5
| 0
| 0
| 1
| -1
| 1
| 100
|
| f
| 0
| 0
| 0
| 1
| 2
| 13400
|
| Новое опорное решение
Х2=(220, 180, 100, 0, 0)
будет оптимальным, так как в строке целевой функции все оценки положительные p0, причем значение целевой функции fmax= f(X2)=13400. В следующей задаче составим математическую модель и найдем оптимальное решение симплексным методом.
Пример 11.
Составить оптимальный план производства изделий двух видов А и В, обеспечивающий максимальную стоимость их реализации, если на
изготовление единицы изделия А требуется затратить а1=2 кг сырья первого типа, а2=3 кг сырья второго типа и а3=1 кг сырья третьего типа. Для единицы изделия В требуется b1=1 кг сырья первого типа, b2=4 кг сырья второго типа и b3=3 кг сырья третьего типа. Производство обеспечено сырьем каждого типа в количестве 400 кг, 900 кг, 600 кг соответственно. Стоимость единицы изделия А составляет 60 руб., а единицы изделия В- 40 руб.
Решение.
Построим математическую модель ЗЛП.
Пусть х1- количество изделий вида А,
х2- количество изделий вида В.
Тогда потребуется затратить сырья первого типа в количестве а1х1=2х1 кг на изделия вида А и b1х2=х2 кг на изделия вида В. Значит, общие затраты сырья первого типа по плану составят а1х1+b1х2=(2х1+х2) кг. Аналогично, сырья второго типа в килограммах потребуется а2х1+b2х2=3х1+4х2 и сырья третьего типа а3х1+b3х2=х1+3х2. При реализации изделий будет получено f=(60х1+40х2) рублей.
Нужно найти такие х1, х2, чтобы выполнялись условия данной задачи по запасам сырья трех типов:
и при этом функция f=60х1+40х2
должна достигать максимума.
Решим задачу симплексным методом
Запишем эту задачу в канонической форме, введя дополнительные балансовые переменные х3, х4, х5, которые имеют смысл остатков сырья соответственно первого, второго и третьего типов.
Тогда система ограничений имеет вид:
и приведенное выражение для целевой функции f f-60x1-40x2=0
Очевидно, что система ограничений приведена к единичному базису:
базисные переменные- х3, х4, х5,
свободные переменные- х1, х2.
Так как свободные члены (400, 900, 600) строго больше нуля, то исходное опорное решение получим при х1=х2=0 равное
Х0=(0, 0, 400, 900, 600)
Составим исходную симплекс- таблицу.
Базис
| х1
| х2
| х3
| х4
| х5
| bi0
| bi0 / bip
| x3
| 2
| 1
| 1
| 0
| 0
| 400
| 400:2=200
| x4
| 3
| 4
| 0
| 1
| 0
| 900
| 900:3=300
| x5
| 1
| 3
| 0
| 0
| 1
| 600
| 600:1=600
| f
| -60
| -40
| 0
| 0
| 0
| 0
|
|
В строке целевой функции обе оценки свободных переменных(-60;-40) отрицательны. Следовательно, решение Х0=(0, 0, 400, 900, 600) не является оптимальным.
В каждом столбце с отрицательной оценкой все элемент- положительные, следовательно, задача имеет оптимальное решение.
Перейдем от исходного опорного решения Х0 к лучшему Х1.
В качестве ведущего столбца выберем столбец, в котором наибольшая по абсолютной величине отрицательная оценка 1=-60, это первый столбец.
Для положительных элементов первого столбца составим симплексные отношения и запишем их в последнем столбце симплексной таблицы.
В качестве ведущей строки возьмем строку с минимальным симплексным отношением- это первая строка, так как
==200,
таким образом, ведущий элемент b11=2.
С этим ведущим элементом проведем одну итерацию метода Жордана- Гаусса, а именно:
все элементы ведущей строки разделим на ведущий элемент
(на 2);
в ведущем столбце все элементы кроме ведущего заменим нулями (запишем единичный столбец);
элементы b14=b15=0, следовательно, четвертый и пятый столбцы останутся без изменения;
все остальные элементы симплекс-таблицы преобразуем по правилу « определителя второго порядка »,
например, b22=(2.4 -1.3): 2=5/2.
Итак, получим новую симплекс-таблицу:
Базис
| х1
| х2
| х3
| х4
| х5
| аi0
| аi0 / аip
| x1
| 1
| 1/2
| 1/2
| 0
| 0
| 200
| 400
| x4
| 0
| 5/2
| -3/2
| 1
| 0
| 300
| 120
| x5
| 0
| 5/2
| -1/2
| 0
| 1
| 400
| 160
| f
| 0
| -10
| 30
| 0
| 0
| 1200
|
|
Таким образом, полагая свободные переменные х2=х3=0, получим новое опорное решение Х1=(200, 0, 0, 300, 400).
Значение целевой функции(см. в строке целевой функции) при опорном решении Х1 равно
f(X1)=1200.
Полученное решение Х1 снова не оптимально, так как в последней строке таблицы - в строке целевой функции
оценка 2=-10<0.
Но в столбце (втором) над этой оценкой есть положительные числа, следовательно, проведем еще одну итерацию с ведущим элементом b22=5/2
(выбрали ведущую строку – вторую - по минимальному симплексному отношению==120).
.
В результате получим симплексную таблицу
Базис
| х1
| х2
| х3
| х4
| х5
| bi0
| bi0 / bip
| x1
| 1
| 0
| 4/5
| -1/5
| 0
| 140
|
| x2
| 0
| 1
| -3/5
| 2/5
| 0
| 120
|
| x5
| 0
| 0
| 1
| -1
| 1
| 100
|
| f
| 0
| 0
| 4
| 0
| 0
| 13200
|
|
Новое опорное решение найдем при свободных переменных х3=х4=0, получим Х2=(140, 120, 0, 0, 100),
Значение целевой функции равно
f(X2)=13200.
Так как в последней строке целевой функции нет отрицательных оценок, то найденное опорное решение оптимально, причем
х1=140, х2=120, х3=х4=0, х5=100, fmax=13200.
Итак, задача решена, оптимальный план производства будет следующий:
нужно произвести 140 изделий типа А, 120 изделий типа В, которые обеспечат максимальную прибыль при их реализации в размере 13200 руб., причем будет экономия сырья третьего типа в количестве 100 кг.
0>0>0> |