Лабораторная работа №2. Симплексный метод определения оптимального плана

Симплексный метод определения оптимального плана

Данный метод (предложен Канторовичем в 1939г., Данцигом в 1949г.) наиболее эффективен при решении сложных задач и представляет итерационный (шаговый) процесс, начинающийся с 0-ого (опорного) решения (вершины n-мерного многогранника). Далее в поисках оптимального варианта плана предполагается движение по угловым точкам (вершинам многогранника) до тех пор, пока значения целевой функции не достигнет максимальной (минимальной) величины. Рассмотрим алгоритм симплексного метода на примере задачи планирования товарооборота при ограниченных ресурсах сырья.

Предприятие реализует n товарных групп, располагая m ограниченными материально-денежными ресурсами b_i ≥0 (1 ≤ i ≤ m). Известны расходы ресурсов каждого i- вида на производство и реализацию единицы товара каждой группы, представленные в виде матрицы (a _ij) и прибыль, получаемая предприятием от реализации единицы товара j-группы, входящая в целевую функцию Z(X). Метод линейного программирования не отличается от системы (1)-(2):

Z(X) = с₁ Х₁ + с₂ Х₂ + с₃ Х₃ + … +с_n Х_n →max(min) (1)

a₁₁ X₁ + a₁₂ X₂ +…a_1n X _n ≤ b_1,

а₂₁ X₁ + a₂₂ X₂ +…a_2n X _n ≤ b₂ (2) ……………………………………..

a_m1 X₁ + a_m2 X₂ +…a _mn X _n ≤ b _m,

X₁≥0 X₂≥0 X₃≥0 …X _n ≥0

Этапы решения поставленной задачи симплексным методом включают:

1) Составление нулевого опорного плана. Вводим новые неотрицательные (базисные) переменные, благодаря которым система неравенств (2) становится системой уравнений:

a₁₁ X₁ + a₁₂ X₂ +…a_1n X _n + X_n+1 = b₁

a₂₁ X₁ + a₂₂ X₂ +…a_2n X _n + X_n+2 = b₂ (3) ……………………………………..

a_m1 X₁ + a_m2 X₂ +…a _mn X _n + X_n+m = b _m,

Если принимать вводимые переменные за векторы-столбцы, то они представляют собой единичные (базисные!) векторы. Отметим, что базисные переменные имеют простой физический смысл – это остаток конкретного ресурса на складе при заданном плане выпуска продукции, поэтому данный базис называют естественным. Решаем систему (3) относительно базисных переменных:

X_n+1 = b_1, -a₁₁ X₁ - a₁₂ X₂ -…a_1n X _n

X_n+2 = b₂ - a₂₁ X₁ - a₂₂ X₂ -…a_2n X _n (4) ………………………………………..

X_n+m = b _m, - a_m1 X₁ + a_m2 X₂ +…a _mn X _n

Целевую функцию перепишем в виде

Z(X)=0-(-с₁Х₁-с₂Х₂-с₃Х₃-…-с_nХ_n) (5) Полагая, что искомые основные переменные Х₁= X₂=X₃=…=X_n =0, получаем нулевой опорный план Х =(0,0, …0,b₁,b_2,b_3…b_m), при котором Z(X) = 0 (все ресурсы на складе, ничего не производится!). Заносим план в симплексную таблицу.

2) Из отрицательных коэффициентов индексной строки выбираем наибольший по абсолютной величине, что определяет ведущий столбец и показывает – какая переменная на следующей итерации (шаге) перейдет из основных (свободных) в базисные (фактически выбирается товарная группа, чья реализация приносит максимальный доход). Затем запасы сырья b_i делим на соответствующие коэффициенты затрат, результаты заносим в таблицу и определяем минимальное значение Q_min (выбирается ресурс, чей запас наиболее сильно ограничивает выпуск выбранной товарной группы). Это значение выделяет ведущую строку и переменную Х_i, которая при следующем шаге (итерации) выйдет из базиса и станет свободной.

3) Переход к новому плану осуществляется в результате пересчета симплексной таблицы методом Жордана-Гаусса. Сначала заменим в базисе Х_j на Х_i ведущего столбца. Разделим все элементы ведущей строки на разрешающий элемент (РЭ), в результате чего на месте РЭ в ведущей строке будет 1. Так как Х_i стал базисным, то остальные его коэффициенты должны быть равны 0. Новые элементы этого плана находятся по правилу прямоугольника

НЭ=СЭ – (А*В)/РЭ (6)

Оценка полученного плана производится по коэффициентам индексной строки: если они все положительны, то план оптимален, если нет, то план можно улучшить, производя следующую итерацию (шаг).

Пример №2.

На приобретение оборудования для производственного участка выделено 20тыс. руб. Оборудование может быть размещено на площади, не превышающей 72 кв.м. Можно заказать оборудование двух типов: типа А, требующие производственную площадь 6кв.м и дающие 6 тыс.ед. продукции в смену (цена 5000 руб.) и типа В, требующие площадь 12 кв.м и дающие 3тыс.ед., (цена 2000 руб.). Каков оптимальный план приобретения оборудования, обеспечивающий максимальную производительность участка?

Обозначим количество приобретаемого оборудования типа А и В через Х₁ и Х₂ соответственно.

Производительность участка (целевая функция): Z(X) =6Х₁+3Х₂.

Основные ограничения связаны

с денежными средствами: 5Х₁+2Х₂ ≤ 20,

с площадью производственного участка: 6Х₁+12Х₂ ≤ 72.

Вводим новые базисные переменные Х₃ (остаток денежных средств после закупки оборудования) и Х₄ (остаток площадей после размещения оборудования) и перепишем ограничения в виде системы уравнений:

5X₁+2Х₂+X₃=20 (X₃=20 – 5X₁ - 2X₂)

6Х₁+12Х₂+X₄ = 72 (X₄=72 – 6X₁ – 12X₂)

При этом функция цели: Z(X) =6Х₁+3Х₂+0Х₃+0Х₄.

Составляем опорный (0-ой) план: Х= (0, 0, 20, 72), т.е. пока ничего не приобретено (деньги не потрачены, площади пустуют). Составляем симплексную таблицу

Очевидно, что ведущий столбец соответствует Х₁, так как имеет самый большой индекс 6. Находим минимальное значение Q_min=4 (самое жесткое ограничение ресурса), определяя ведущую строку, показывающую, что из базисных переменных выводится Х₃, а вместо нее вводится Х₁. Пересчитываем элементы ведущей строки, разделив их на 5, а по формуле (6) определяем элементы второй и индексной строк. Целевая функция для 1-ого плана равна Z(X) =6*4+3*0=24.Однако, один из коэффициентов индексной строки для столбца Х₂ остается отрицательным -0,6, следовательно данный план не оптимален, и требуется еще одна итерация (шаг) для его улучшения. Выбираем ведущим 2-ой столбец и по минимальному значению Q_min=5 определяем ведущую строку с базисной переменной Х₄. Выполнив те же преобразования, получаем 2-ой план, который будет оптимальным, так как все индексные коэффициенты положительны. Проанализируем полученные результаты. При оптимальном решении целевая функция имеет максимальное значение 27тыс.руб., при этом оба ресурса выведены из базиса, следовательно израсходованы полностью! Убедимся в этом: 5*2+2*5=20 тыс.руб., 6*2+12*5=72 кв.м. Искомое решение Х= (2; 5;0;0;).Так бывает далеко не всегда.

Пример №3.

Коммерческое предприятие реализует три вида продукции А,В, и С. Нормативы затрат материально-денежных ресурсов на реализацию единицы продукции (или 1 тыс. товарооборота) представлены в таблице.

В последней строке задан доход от продажи товара каждой группы. Необходимо определить объем и структуру товарооборота, при котором доход предприятия будет максимальным. Математическая модель задачи: определение вектора Х=(Х₁; Х₂; Х₃), где Х – количество проданных изделий А, В,С, который удовлетворяет неравенствам

0,1X₁ + 0,2X₂ +0,4 X₃ ≤ 1100_,

0,05X₁ + 0,02X₂ +0,02 X₃ ≤ 120,

3X₁+1X₂+2X₃≤8000. X₁≥0 X₂≥0 X₃≥0

и обеспечивает максимальное значение целевой функции Z(X)=3Х₁+5Х₂+4Х₃.

Для построения опорного 0-ого плана вводим дополнительно базисные переменные Х₄, Х₅,Х₆ (определяем объемы неиспользованных ресурсов), преобразовывая систему неравенств в систему уравнений,

0,1X₁ + 0,2X₂ +0,4 X₃+Х₄= 1100_,

0,05X₁ + 0,02X₂ +0,02 X₃+Х₅= 120,

3X₁+1X₂+2X₃+Х₆=8000.

решаем полученную систему относительно базисных переменных Х₄, Х₅, Х₆:

Х₄=1100-0,1X₁ - 0,2X₂ -0,4 X_3,

Х₅= 120-0,05X₁ - 0,02X₂ -0,02 X_3,

Х₆=8000-3X₁ -X₂-2X_3.

Функцию цели запишем в виде Z(X)=0 – (- 3Х₁ – 5Х₂ – 4Х₃) и полагая, что основные переменные Х₁=Х₂=Х₃=0, получим опорный план Х=(0; 0; 0; 1100; 120; 8000), при котором целевая функция Z(X)=0 (товары не продаются, ресурсы не используются, доход нулевой). Опорный план запишем в симплексную таблицу:

Естественно, 0-ой план не оптимален, так как коэффициенты в индексной строке отрицательные -3, -5, -4.Ведущий столбец содержит Х₂, а ведущая стока содержит Х₄, РЭ=0,2 находится на их пересечении. Формируем 1-ый план аналогично предыдущему примеру, в индексной строке только коэффициент при Х₁ отрицателен -0,5, поэтому ведущий столбец содержит Х₁. Строка, содержащая Х₅, становится ведущей, и после итерации получаем 2-ой план, все коэффициенты которого в индексной строке положительны. Это означает, что план Х= (250; 5375; 0; 0; 0; 1875) оптимален, и целевая функция при данных значениях Х имеет максимальное значение Z(X) =27625 тыс.руб.

Проанализируем полученное решение: для получения максимальной прибыли необходимо продать 250 ед. товаров группы А, 5375 ед. товаров группы В, товары группы С вообще не реализуются! В оптимальном плане среди базисных переменных находится Х₆. Это указывает на то, что 1875м² складских помещений не используются, то есть данный ресурс в исходном объеме избыточен, и данная его часть (1875м²) может быть применена в других целях.

Поиск по сайту: