Метод искусственного базиса
Симплексный метод решения задач базируется на введении дополнительных переменных, позволяющих образовать единичную матрицу, в которую не допускаются отрицательные и другие числа, кроме нуля и единицы. Наличие единичной матрицы является необходимым условием при решении задач симплексным методом.
Если
же ограничения задачи заданы в виде
неравенств вида
или уравнений
то невозможно сразу получить начальное базисное решение, если матрица, составленная из коэффициентов при неизвестных системы ограничений, не позволяет образовать единичную матрицу. Причем уравнения отражают жесткие условия ограничений по ресурсам, не допускающих никаких отклонений. Для соблюдения равенств вводятся искусственные переменные yi, равные нулю. Векторы искусственных переменных образуют необходимую для решения единичную матрицу. Такой базис называется искусственным, а метод решения называется методом искусственного базиса. Причем искусственные переменные не имеют отношения к содержанию поставленной задачи, однако они позволяют построить стартовую точку, а процесс оптимизации вынуждает эти переменные принимать нулевые значения и обеспечить допустимость оптимального решения.
Преобразование разнородных ограничений, представляющих собой смесь уравнений и неравенств разного вида, заключается в образовании базиса решения путем одновременного введения свободных и искусственных переменных, что придает симплексному методу большую гибкость.
Задача.
Определить
минимальное и максимальное значение
целевой функции
при следующих смешанных условиях -
ограничениях:
Запишем систему ограничений в виде равенств, для чего введем дополнительные переменные x3 и x4, в результате получим следующую систему:
Затем введем искусственные переменные y1 и y2 в первое и второе уравнения:
Для постановки задачи на минимум целевую функцию запишем в виде:
Для решения этой же задачи на максимум целевую функцию следует записать иначе:
За
использование искусственных переменных,
вводимых в целевую функцию, накладывается
так называемый “штраф” величиной М,
очень большое положительное число,
которое обычно не задается (М
).
С целью формулировки задачи для решения ее в табличной форме воспользуемся выражениями из системы уравнений для искусственных переменных:
которые подставим в целевую функцию:
Таким образом, стартовая точка решения определяется x1 = x2 = x3 =0, x4 = 0, y1= 3, y2 = 8, следовательно, уравнение целевой функции для симплексной таблицы будет иметь такой вид:
Последовательность решения задачи симплексным методом представлена в таблице 10.
Поскольку
задача решается на минимум, то ведущий
столбец выбирают по максимальному
положительному числу в индексной строке
в плане I
(9М-3),
а все остальные преобразования проводятся
по стандартной схеме до тех пор, пока
не получатся в индексной строке только
неположительные числа (
),
а в оптимальном решении должны
отсутствовать положительные значения
искусственных переменных y1
и y2.
Оптимальному решению соответствует
точка с координатами x1
=11/3
и x2
=1/3,
где
а
y1=0
и y2=0.
Таблица 10
План |
Базисные переменные |
Значения базисных переменных |
x1 |
x2 |
x3 |
x4 |
y1 |
y2 |
|
I |
y1 |
3 |
2 |
1 |
0 |
0 |
1 |
0 |
3 |
y2 |
8 |
4 |
8 |
-1 |
0 |
0 |
1 |
1 |
|
x4 |
5 |
1 |
2 |
0 |
1 |
0 |
0 |
2,5 |
|
|
11M |
(6M-1) |
(9M-3) |
-M |
0 |
0 |
0 |
|
|
II |
y1 |
2 |
1,5 |
0 |
1/8 |
0 |
1 |
-1/8 |
11/3 |
x1 |
1 |
1/2 |
1 |
-1/8 |
0 |
0 |
1/8 |
2 |
|
x4 |
3 |
0 |
0 |
1/4 |
1 |
0 |
-1/4 |
- |
|
|
3+2М |
(0,5+1,5M) |
0 |
M-3/8 |
0 |
0 |
3-9M/8 |
|
|
III |
x1 |
11/3 |
1 |
0 |
1/12 |
0 |
2/3 |
-1/12 |
|
x2 |
1/3 |
0 |
1 |
-1/6 |
0 |
-1/3 |
1/6 |
|
|
x4 |
3 |
0 |
0 |
1/4 |
1 |
0 |
-1/4 |
|
|
|
21/3 |
0 |
0 |
-5/12 |
0 |
(-1/3-M) |
3-9M/8 |
|
Для решения этой же задачи на максимум целевую функцию следует записать иначе:
Затем, подставив выражение для искусственных переменных, получим:
Преобразуем это выражение в удобную форму для записи в симплексную таблицу:
Полученные выражения вместе с условиями ограничения записываем для преобразования в симплексную таблицу 11.
Преобразования проводим до тех пор, пока все значения в индексной строке не будут положительными.
В плане IV получили x1=1/3, x2=21/3, а Fmax(-X4)=71/3.
Таблица 11
План |
Базис |
Значения |
x1
|
x2
|
x3
|
x4
|
y1
|
y2
|
|
I |
y1 |
3 |
2 |
1 |
0 |
0 |
1 |
0 |
3 |
y2 |
8 |
4 |
8 |
-1 |
0 |
0 |
1 |
1 |
|
x4 |
5 |
1 |
2 |
0 |
1 |
0 |
0 |
2,5 |
|
|
-11M |
-6M-1 |
-9M-3 |
M |
0 |
0 |
0 |
|
|
II |
y1 |
2 |
1,5 |
1 |
1/8 |
0 |
1 |
-1/8 |
1 |
x2 |
1 |
1/2 |
0 |
-1/8 |
0 |
0 |
1/8 |
2 |
|
x4 |
3 |
0 |
0 |
¼ |
1 |
0 |
-1/4 |
- |
|
|
3-2M |
0,5-1,5M |
0 |
-M-3/8 |
0 |
0 |
3+9M/8 |
|
|
III |
x1 |
11/3 |
1 |
0 |
1/12 |
0 |
2/3 |
-1/12 |
16 |
x2 |
1/3 |
0 |
1 |
-1/6 |
0 |
-1/3 |
1/6 |
- |
|
x4 |
3 |
0 |
0 |
1/4 |
1 |
0 |
-1/4 |
12 |
|
|
21/3 |
0 |
0 |
-5/12 |
0 |
-1/3+M |
5/12+M |
|
|
IV |
x1 |
1/3 |
1 |
0 |
0 |
-1/3 |
2/3 |
0 |
|
x2 |
21/3 |
0 |
1 |
0 |
2/3 |
-1/3 |
0 |
|
|
x3 |
12 |
0 |
0 |
0 |
4 |
0 |
-1 |
|
|
|
71/3 |
0 |
0 |
1 |
12/3 |
M-1/3 |
M |
|