| Продукция / Сырье | Нормы расхода сырья, кг/ед. | Объемы запасов сырья, кг | |
| П1 | П2 | ||
| С1 | 1 | 3 | 300 |
| С2 | 1 | 1 | 150 |
| Прибыль, у.е./ед. прод. | 2 | 3 | |
Составить план производства по критерию "максимум прибыли".
Решение. Введем необходимые обозначения. Обозначим объем производства продукции П1 через х1, продукции П2 через х2. Таким образом, формально (математически) план производства (производственная программа) - это вектор
. С учетом введенных обозначений математическая модель задачи по критерию "максимум прибыли" имеет вид

при ограничениях

Приведем эту ЗЛП к каноническому виду, введя дополнительные переменные х3 и х4:


или

КЗЛП имеет необходимое число единичных столбцов, т.е. обладает очевидным начальным опорным планом (0, 0, 300, 150). Решение осуществляется симплекс-методом с естественным базисом с оформлением расчетов в симплекс-таблицах (табл. 2.2).
В исходной симплекс-таблице с номером 0 строка оценок
определяется по приведенной выше формуле (2.22):

Исходный опорный план (0, 0, 300, 150) не является оптимальным, так как среди оценок ∆j имеются отрицательные. Переход к новому опорному плану осуществим, введя в базис вектор А2, имеющий минимальную отрицательную оценку. Определяем вектор, выходящий из базиса:

т.е. вектор А3 следует вывести из базиса. Разрешающим элементом является а12 = 3 (выделен рамкой). Переход к следующей симплекс-таблице осуществляем с помощью преобразований Жордана - Гаусса, при этом рекомендуется пользоваться двумя описанными выше вычислительными процедурами симплексного метода, включая правило прямоугольного треугольника.
Таблица 2.2
Решение ЗЛП в симплекс-таблицах
| Номер симплекс- таблицы | Базис | ci / cj | План В | 2 | 3 | 0 | 0 | Q |
| A1 | A2 | A3 | A4 | |||||
| 0 | А3 | 0 | 300 | 1 | 3 | 1 | 0 | 100 |
| А4 | 0 | 150 | 1 | 1 | 0 | 1 | 150 | |
| ∆j | - | 0 | -2 | -3 | 0 | 0 | - | |
| I | A2 | 3 | 100 | 1/3 | 1 | 1/3 | 0 | 300 |
| A4 | 0 | 50 | 2/3 | 0 | -1/3 | 1 | 75 | |
| ∆j | - | 300 | -1 | 0 | 1 | 0 | - | |
| II | A2 | 3 | 75 | 0 | 1 | 1/2 | -1/2 | |
| A1 | 2 | 75 | 1 | 0 | -1/2 | 3/2 | ||
| ∆j | - | 375 | 0 | 0 | 1/2 | 3/2 | - |
Второй опорный план (0, 100, 0, 50) не оптимальный; переход к следующему опорному плану осуществим, вводя в базис вектор A1 и выводя вектор А4.
В симплекс-таблице II получен оптимальный опорный план, поскольку все симплекс-разности (оценки)
Оптимальные значения переменных равны:
(основные переменные),
(дополнительные переменные).
Максимальное значение целевой функции равно 375.
Таким образом, в рассмотренной задаче об оптимальном использовании ограниченных ресурсов, оптимальная производственная программа состоит в выпуске 75 ед. продукции первого вида и 75 ед. продукции второго вида. С этой программой связана максимальная прибыль от реализации готовой продукции - 375 у.е.
Симплекс-метод с искусственным базисом (М-метод)
Применяется в тех случаях, когда затруднительно найти первоначальный опорный план исходной задачи ЛП, записанной в канонической форме.
М-метод заключается в применении правил симплекс-метода к так называемой М-задаче. Она получается из исходной добавлением к левой части системы уравнений в канонической форме исходной ЗЛП таких искусственных единичных векторов с соответствующими неотрицательными искусственными переменными, чтобы вновь полученная матрица содержала систему единичных линейно-независимых векторов. В линейную форму исходной задачи добавляется в случае ее максимизации слагаемое, представляющее собой произведение числа (-М) на сумму искусственных переменных, где М - достаточно большое положительное число.
В полученной задаче первоначальный опорный план очевиден. При применении к этой задаче симплекс-метода оценки ∆j теперь будет зависеть "от буквы М". Для сравнения оценок нужно помнить, что М - достаточно большое положительное число, поэтому из базиса будут выводиться в первую очередь искусственные переменные.
В процессе решения М-задачи следует вычеркивать в симплекс-таблице искусственные векторы по мере их выхода из базиса. Если все искусственные векторы вышли из базиса, то получаем исходную задачу. Если оптимальное решение М-задачи содержит искусственные векторы или М-задача неразрешима, то исходная задача также неразрешима.
Путем преобразований число вводимых переменных, составляющих искусственный базис, может быть уменьшено до одной.
Пример 2.8. Найти максимум целевой функции:
при условиях

Решение. Матрица условий содержит только один единичный вектор, добавим еще один искусственный вектор (искусственную неотрицательную переменную ух в первое ограничение):

Получим следующую М-задачу: найти максимум целевой функции
- Му1 при условиях

М-задачу решаем симплекс-методом. Начальный опорный план (0, 0, 6, 8), решение проводим в симплекс-таблицах (табл. 2.3).
В начальной таблице наименьшее ∆j соответствует вектору А1 - он вводится в базис, а искусственный вектор P1 из базиса выводится, так как ему отвечает наименьшее Q. Столбец, соответствующий Р1, из дальнейших симплексных таблиц вычеркивается.
Таблица 2.3
Решение М-задачи в симплекс-таблицах
| Номер симплекс- таблицы | Базис | ci / cj | План В | 3 | 2 | 1 | -M | Q |
| A1 | A2 | A3 | P1 | |||||
| 0 | P1 | -M | 8 | 2 | 1 | 0 | 1 | 4 |
| A3 | 1 | 6 | 1 | 1 | 1 | 0 | 6 | |
| ∆j | - | -8М+6 | -2М-2 | -М-1 | 0 | 0 | - | |
| 1 | A1 | 3 | 4 | 1 | 0,5 | 0 | - | |
| A3 | 1 | 2 | 0 | 0,5 | 1 | - | ||
| ∆j | - | 14 | 0 | 0 | 0 | - |
Полученный новый опорный план является опорным планом исходной задачи. Для него все ∆j > 0, поэтому он является и оптимальным. Таким образом, получен оптимальный план исходной задачи (4, 0, 2) и максимальное значение целевой функции max
.
Пример 2.9. Решить ЗЛП:

Решение. Приведем ЗЛП к каноническому виду, перейдя к задаче "на максимум":

Для нахождения опорного плана переходим к М-задаче:

Дальнейшее решение проводим в симплекс-таблицах (табл. 2.4).
Таблица 2.4
Решение ЗЛП М-методом
| Номер симплекс-таблицы | Базис | ci / cj | В | -10 | 5 | 0 | 0 | 0 | -M | -M | Q |
| A1 | A2 | A3 | A4 | A5 | P1 | P2 | |||||
| 0 | P1 | -M | 3 | 2 | -1 | -1 | 0 | 0 | 1 | 0 | 3/2 |
| P2 | -M | 2 | 1 | 1 | 0 | -1 | 0 | 1 | 1 | 2 | |
| A5 | 0 | 1 | -1 | -2 | 0 | 0 | 0 | 0 | 0 | - | |
| - | ∆j | - | -5 М | -3 М+ + 10 | -5 | М | М | 0 | - | ||
| I | A1 | -10 | 3/2 | 1 | -1/2 | -1/2 | 0 | 0 | 0 | - | |
| P2 | -M | 1/2 | 0 | 3/2 | 1/2 | -1 | 0 | 1 | 1/3 | ||
| A5 | 0 | 5/2 | 0 | -5/2 | -1/2 | 0 | 1 | 0 | - | ||
| - | ∆j | - | -M/2-15 | 0 | -3М/2 | -М/2+5 | М | 0 | 0 | - | |
| II | A1 | -10 | 5/3 | 1 | 0 | -1/3 | -1/3 | 0 | - | ||
| A2 | 5 | 1/3 | 0 | 1 | 1/3 | -2/3 | 0 | - | |||
| A5 | 0 | 10/3 | 0 | 0 | 0 | -5/3 | 1 | - | |||
| - | ∆j | - | -15 | 0 | 0 | 5 | 0 | 0 | - |
В симплекс-таблице II получен опорный план исходной ЗЛП; поскольку все симплекс-разности
, то этот план является и оптимальным, т.е.
(исходные переменные),
(дополнительные переменные), при этом 
При рассмотрении графического метода выделялись три особых случая решения ЗЛП. В симплекс-методе эти случаи определяются следующим образом.
1. Если найден оптимальный план и оценки всех свободных переменных строго больше нуля, то оптимальный план является единственным', если оценки некоторых свободных переменных в оптимальном плане равны нулю, то этот план будет неединственным, так как ввод этих переменных в базис не нарушает критерия оптимальности и не меняет оптимальное значение целевой функции. В соответствии с этим оптимальный план в табл. 2.2 является единственным, а в табл. 2.3 и 2.4 - несдинственным (первый особый случай).
2. Если в процессе решения ЗЛП М-методом искусственные переменные не выводятся из базиса, это является свидетельством того, что область определения исходной ЗЛП является пустым множеством: в этом случае ЗЛП не имеет решения ввиду противоречивости системы ограничений (второй особый случай).
3. Если в направляющем столбце все элементы ajk неположительны (см. 2.23), то это свидетельствует о незамкнутости области определения ЗЛП; в этом случае ЗЛП не имеет решения ввиду неограниченности целевой функции (третий особый случай).
Для автоматизации решения задач линейного программирования могут быть использованы стандартные офисные средства Microsoft Excel - надстройка Поиск решения, использующая симплексный метод (линейная оптимизация с помощью надстройки Поиск решения подробно рассмотрена, например, в литературе).
Однако для корректного и эффективного использования программных средств необходимо знать основы линейного программирования, изложенные выше в данной главе.
Дата: 2019-03-05, просмотров: 302.