Алгоритм решения задач симплексным методом

8. После заполнения новой симплекс-таблицы алгоритм возвращается к 5-му пункту.

Конец работы алгоритма:

• либо когда в F-строке все коэффициенты (кроме С₀) будут больше нуля, тогда получаем оптимальное решение;
• либо когда существует столбец, у которого C_j<<0 и все коэффициенты A_ij<<0, в этом случае решения не существует.

 

 

2 Алгоритм решения задач симплексным методом

 

Формализованный алгоритм симплекс-метода состоит из двух основных этапов:

1) построение опорного плана;

2) построение оптимального плана.

Алгоритм симплекс-метода рассмотрим на примере реальной задачи.

Пусть предприятие может  изготавливать изделия четырех видов: U1, U2, U3, U4. Для изготовления каждого изделия требуется определенный вид оборудования. Известно, сколько времени требуется на изготовление каждого изделия на каждом оборудовании и какая прибыль может быть получена при реализации каждого изделия. Требуется так распределить изделия по видам оборудования, чтобы предприятие имело максимальную прибыль.

Таблица 2.1. Распределение изделий по видам оборудования.

	U 1	U2	U3	U4	bj
O1	3	5	2	7	15
O2	4	3	3	5	9
O3	5	6	4	8	30
сj	40	50	30	20

 

Где b_j - ресурсы оборудования O_i;

       a_ij - время изготовления j-ого изделия на i-ом оборудовании;

       c_j - прибыль от одного изделия U_j;

       x_j - количество j-ых изделий, которое необходимо выпустить на предприятии.

F(x)= 40x₁ + 50х₂ + 30х₃ + 20х₄→ max                                                               (2.1)

3x₁ + 5x₂ + 2x₃ + 7x₄ ≤ 15

4x₁ + 3x₂ + 3x₃ + 5x₄ ≤ 9                                                                                      (2.2)

5x₁ + 6x₂ + 4x₃ + 8x₄ ≤ 30

Система неравенств преобразуется к равенствам, если к каждому из неравенств ввести фиктивные изделия, остаточные переменные x5, x6, x7, которые используют оставшийся ресурс оборудования. При этом время изготовления фиктивных изделий на соответствующем оборудовании принимается равным:

a₁₅=1, a₂₅=0, a₃₅=0

a₁₆=0, a₂₆=1, a₃₆=0

a₁₇=0, a₂₇=0, a₃₇=1

 

Тогда система равенств и условие неотрицательности запишется в виде:

3*x₁ + 5*x₂ + 2*x₃ + 7*x₄ + 1*x₅ + 0*x₆ + 0*x₇ = 15

4*x₁ + 3*x₂ + 3*x₃ + 5*x₄ + 0*x₅ + 1*x₆ + 0*x₇ = 9                                           (2.3)                    

5*x₁ + 6*x₂ + 4*x₃ + 8*x₄ + 0*x₅ + 0*x₆ + 1*x₇ = 30

 

                                                                                                                             (2.4)

 

Целевая функция запишется следующим  образом:

                                       

                                                                                                                             (2.5)

При этом: с₅ = c₆ = c₇ = 0

Таким образом, задача распределения  продукции по видам оборудования сводится к нахождению таких значений переменных x_j (j от 1 до n), которые удовлетворяют системе равенств вида (2.3), дополнительным ограничениям (2.4) и максимизируют прибыль предприятия вида (2.5).

Решение общей распределительной  задачи выполняется в два этапа:

1-ый этап. Находим любое решение (как правило, неоптимальное), удовлетворяющее ограничениям (2.3) и (2.4) или убеждаемся, что решения не существует. Этот этап называется определением опорного плана (базиса).

2-ой этап. Производим последовательное улучшение данного плана до получения оптимального. В некоторых задачах опорный план не трудно определить, в противном случае используют специальные методы получения опорного плана (метод искусственного базиса).

В нашем случае базис выражается легко, его составляют дополнительно  введенные свободные переменные x₅, x₆, x₇. Выразим их через остальные и результат запишем в форме Таккера²:

x₅ = 15 - (3x₁ + 5x₂ + 2x₃ + 7x₄)

x₆ = 9 - (4x₁ + 3x₂ + 3x₃ + 5x₄)                         

x₇ = 30- (5x₁ + 6x₂ + 4x₃ + 8x₄)

В нашем случае x₅, x₆, x₇ называются базисными и соответственно x₁, x₂, x₃, x₄ небазисными переменными. Для получения опорного плана надо приравнять к нулю небазисные переменные, тогда получим:

Х = (0; 0; 0; 0; 15; 9; 30)

Опорный план получен, и первый этап решения задачи завершен.

На втором этапе (улучшение опорного плана) можно использовать несколько способов, однако на практике чаще всего используют симплекс-метод.

1) Выражение целевой  функции через небазисные переменные:

Y = 40x₁ + 50x₂ + 30x₃ + 20x₄

Записать целевую функцию в  форме Таккера:

Y = 0 - (-40x₁ - 50x₂ - 30x₃ - 20x₄)

2) Проверка базисного решения  на оптимальность. Если в полученном  выражении целевой функции, записанной  в форме Таккера, все коэффициенты  при небазисных переменных неотрицательны (≥ 0), то исходный базис является  оптимальным и находится соответствующее оптимальное базисное решение, минимизирующее целевую функцию. Для этого всем небазисным переменным присваивается значение равное 0, а все базисные переменные находят из системы ограничений, и одновременно определяется целевая функция. Если в полученном выражении целевой функции хотя бы один коэффициент при небазисных переменных отрицателен, то переходят к следующему этапу.

3) Проверка задачи на наличие  решения. Если при какой-либо  небазисной переменной, имеющей  отрицательный коэффициент целевой функции, окажется, что столбец коэффициентов при этой же переменной в системе ограничений состоит из одних неположительных чисел, то максимум целевой функции равен бесконечности, и говорят, что целевая функция не ограничена.

4) Выбор из небазисных переменных той, которая способна при введении её в базис увеличить значение целевой функции. Наиболее часто на практике в качестве такой переменной выбирают ту, которой соответствует наибольший отрицательный коэффициент целевой функции. Следовательно, в базис вводится x₂.

5) Определение, какая из базисных  переменных должна быть выведена  из базиса. Для всех положительных  коэффициентов столбца при вводимой  переменной в системе ограничений  находят отношение свободного  члена соответствующему коэффициенту столбца. Минимальное из полученных соотношений указывает строку выводимой переменной. 15/5; 9/3; 30/6 - в соответствии с этими отношениями из базиса выводим х₅.

6) Выражение вводимой в базис  переменной через выводимую и  другие небазисные переменные.

x₂ = 3 - (3/5x₁ + 2/5x₃ + 7/5x₄ + 1/5x₅)

7) Выражение остальных базисных  переменных х₆, х₇ и целевой функции через новые небазисные переменные х₁, х₃, х₄, х₅:

x₆ = 0 - (11/5x₁ + 9/5x₃ + 4/5x₄ - 3/5x₅)

x₇ = 12 - (1/5x₁ + 8/5x₃ - 2/5x₄ - 6/5x₅)

Y = 150 - (-10x₁ - 10x₃ + 50x₄ + 10x₅)

8) Повторение операций  пунктов 2) - 7) до тех пор, пока не будет найдено оптимальное решение.

Для того, чтобы решить задачу симплексным методом, необходимо выполнить следующее:

1) Привести задачу к каноническому виду;

2) Найти начальное опорное решение с "единичным базисом" (если опорное решение отсутствует, то задача не имеет решения ввиду несовместимости системы ограничений);

3) Вычислить оценки разложений векторов по базису опорного решения и заполнить таблицу симплексного метода;

4) Если выполняется признак единственности оптимального решения, то решение задачи заканчивается;

5) Если выполняется условие существования множества оптимальных решений, то путем простого перебора находят все оптимальные решения.

 

Первый этап. Дана целевая функция, для которой необходимо найти наибольшее или наименьшее значение, если значения всех неизвестных неотрицательные:

Строится исходная оптимизационная модель. Далее исходная матрица условий преобразуется в приведенную каноническую форму, которая среди всех других канонических форм выделяется тем, что:

а) правые части условий (свободные члены bi) являются величинами неотрицательными;

б) сами условия являются равенствами;

в) матрица условий  содержит полную единичную подматрицу.

Если свободные члены  отрицательные, то обе части неравенства  умножаются на – 1, а знак неравенства  меняется на противоположный. Для преобразования неравенств в равенства вводятся дополнительные переменные, которые, обычно, обозначают объем недоиспользованных ресурсов. В этом и состоит их экономический смысл.

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

В оптимальном решении  задачи все искусственные переменные (ИП) должны быть равными нулю. Для этого вводят искусственные переменные в целевую функцию задачи с большими отрицательными коэффициентами (-М) при решении задачи на max, и с большими положительными коэффициентами (+М), когда задача решается на min. В этом случае даже незначительное ненулевое значение искусственной переменной будет резко уменьшать (увеличивать) значение целевой функции. Обычно М в 1000 раз должно быть больше, чем значения коэффициентов при основных переменных.

Второй этап. Строится исходная симплекс-таблица и отыскивается некоторое начальное базисное решение. Множество переменных, образующих единичную подматрицу, принимается за начальное базисное решение. Значения этих переменных равны свободным членам. Все остальные небазисные переменные равны нулю.

Третий этап. Проверка базисного решения на оптимальность осуществляется при помощи специальных оценок коэффициентов целевой функции. Если все оценки коэффициентов целевой функции отрицательны или равны нулю, то имеющееся базисное решение – оптимальное. Если хотя бы одна оценка коэффициента целевой функции больше нуля, то имеющееся базисное решение не является оптимальным и должно быть улучшено.

Четвертый этап. Переход к новому базисному решению. Очевидно, что в оптимальный план должна быть введена такая переменная, которая в наибольшей степени увеличивает целевую функцию. При решении задач на максимум прибыли в оптимальный план вводится продукция, производство которой наиболее выгодно. Это определяется по максимальному положительному значению оценки коэффициента целевой функции.

Столбец симплексной  таблицы с этим номером на данной итерации называется генеральным столбцом.

Далее, если хотя бы один элемент генерального столбца а_ij строго положителен, то отыскивается генеральная строка (в противном случае задача не имеет оптимального решения).

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

Элемент симплексной  таблицы, находящийся на пересечении  генеральных столбца и строки, называется генеральным элементом.

Затем все элементы генеральной  строки (включая свободный член), делятся на генеральный элемент. В результате этой операции генеральный  элемент становится равным единице. Далее необходимо, чтобы все другие элементы генерального столбца стали бы равны нулю, т.е. генеральный столбец должен стать единичным. Все строки (кроме генеральной) преобразуются следующим образом. Полученные элементы новой строки умножаются на соответствующий элемент генерального столбца и полученное произведение вычитается из элементов старой строки.