Симплекс-метод с естественным базисом

Для применения симплекс-метода с естественным базисом КЗЛП должна содержать единичную подматрицу размером mxm – в этом случае очевиден начальный опорный план (неотрицательное базисное решение системы ограничений КЗЛП).

Для определенности предположим, что первые m векторов матрицы системы уравнений составляют единичную матрицу. Тогда первоначальный опорный план очевиден – (b₁, b₂,…, b_m,0,…,0).

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

Математический признак оптимальности состоит из следующих двух теорем:

1. Если для всех векторов А₁, А₂,…, А_n выполняется условие

где

то полученный опорный план является оптимальным.

2. Если для некоторого вектора, не входящего в базис, выполняется условие

,

то можно получить новый опорный план, для которого значение ЦФ будет больше исходного, при этом могут быть два случая:

- если все компоненты вектора, подлежащего вводу в базис, неположительны, то ЗЛП не имеет решения (конечного оптимума нет);

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

На основании признака оптимальности в базис вводится вектор Ак, давший минимальную отрицательную величину симплекс разности:

Чтобы выполнялось условие неотрицательности значений опорного плана, выводится из базиса вектор Аr, который дает минимальное положительное оценочное отношение

Строка Аr называется направляющей, столбец Ак и элемент ar к – направляющими.

Элементы направляющей строки в новой симплекс-таблице вычисляются по формулам:

а элементы i-й строки – по формулам:

Значения нового опорного плана рассчитываются по формулам:

для i = r;

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

Примечание. Для использования приведенной процедуры к минимизации линейной функции f(x₁,x₂,…, x_n) следует искать максимум – f(x₁,x₂,…, x_n), затем полученный максимум взять с противоположным знаком. Оптимальное решение то же.

Пример. Получить решение по модели:

Эта задача (модель) линейного программирования, приведем ее к каноническому виду путем введения дополнительных переменных x₃ и x₄:

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

Номер симплекс-таблицы	Базис		В план					Q
	А 3
А 4
	f(x)		-2	-3
I	А 2			1/3		1/3
А 4			2/3		-1/3
	f(x)		-1
II	А 2					1/2	-1/2
А 1					-1/2	3/2
	f(x)				1/2	3/2

В симплекс-таблице II получен оптимальный опорный план, поскольку все симплекс-разности (оценки) j. Оптимальные значения переменных равны: x₁*=75, x₂* =75 (основные переменные), x₃* =0, x₄* =0 (дополнительные переменные). Максимальное значение целевой функции равно 375.

Таким образом, в рассмотренной выше задаче об оптимальном использовании ограниченных ресурсов, оптимальная производственная программа состоит в выпуске 75 ед. изделий первого вида и 75 ед. изделий второго вида. С этой программой связана максимальная выручка от реализации готовой продукции – 375 у.е.