
методы оптимальных решений / 404921108
.pdf
|
|
|
|
31 |
|
|
|
|
|
f ( X ) |
|
2x1 |
3x3 |
4x5 |
Max; |
||||
2x1 |
|
|
x3 |
x4 |
|
x5 |
|
2 |
|
x1 |
x2 |
4x3 |
|
2x5 |
|
1 |
|||
4x1 |
|
|
5x3 |
|
|
x5 |
x6 |
2 , |
|
x1 0, x2 |
0, x3 |
0, x4 |
0, x5 |
0, x6 0. |
|||||
Вычислим |
|
|
|
|
|
|
|
|
|
|
1 |
5 |
0 |
1 |
|
2 |
3 |
x5 |
|
Q B |
6 |
4 |
2 |
|
1 |
3 |
x |
||
|
|
|
|||||||
|
4 |
|
|
|
|
|
|
2 |
|
|
1 |
0 |
1 |
|
2 |
1 |
x |
||
|
|
|
|
|
|
|
|
|
3 . |
Все компоненты полученного вектора положительные, следовательно, план
X (0;3;1;0;3) является |
оптимальным |
планом новой задачи. Оптимальное |
значение целевой функции равно f Opt |
f ( X ) 2 0 3 1 4 3 9 . |
|
Если условие Q B |
0 не выполняется, то новую задачу придется решать |
ссамого начала.
6.13.Выпишем еще раз первую и последнюю симплекс-таблицы рассматриваемой задачи:
|
|
|
|
|
|
|
|
Таблица 1 |
||||
|
Б |
x1 |
x2 |
x3 |
x4 |
x5 |
x6 |
B |
|
|||
|
x4 |
2 |
0 |
1 |
1 |
1 |
|
0 |
|
2 |
|
|
|
x2 |
1 |
1 |
4 |
0 |
2 |
|
0 |
|
1 |
|
|
|
x6 |
4 |
0 |
5 |
0 |
1 |
|
1 |
|
0 |
|
|
|
f |
2 |
0 |
3 |
0 |
4 |
|
0 |
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Таблица 3 |
||||
|
Б |
x1 |
x2 |
x3 |
x4 |
x5 |
|
x6 |
|
B |
|
|
|
x5 |
7 / 2 |
0 |
0 |
5 / 4 |
1 |
1/ 4 |
5 / 2 |
|
|
||
|
x2 |
0 |
1 |
0 |
3/ 2 |
0 |
|
1/ 2 |
|
4 |
|
|
|
x3 |
3/ 2 |
0 |
1 |
1/ 4 |
0 |
|
1/ 4 |
|
1/ 2 |
|
|
|
f |
15 / 2 |
0 |
0 |
17 / 4 |
0 |
|
1/ 4 |
17 / 2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Четвертый, второй и шестой столбцы матрицы ограничений в первой таблице образуют единичную матрицу E . Оценки соответствующих этим столбцам переменных в последней таблице называются двойственными оценками оптимального базисного плана задачи ЛП в канонической базисной форме. Название связано с так называемой двойственной задачей ЛП, оптимальное решение которой совпадает с этими оценками. Двойственные оценки имеют специ-
альное обозначение yi* , причем номер двойственной оценки совпадает с номером ограничения в первой симплексной таблице.

32
Из первой таблицы видно, что двойственные оценки образованы оценка-
ми переменных x4 , |
x2 и x6 |
(именно в этом порядке расположены столбцы еди- |
|||||||||
ничной |
матрицы в |
первой |
таблице). Из последней таблицы следует, что |
||||||||
y* |
4 |
17 / 4 |
, |
y* |
2 |
0 |
и |
y* |
6 |
1/ 4 |
. |
1 |
|
2 |
|
3 |
|
6.14. Двойственные оценки обладают следующими важными свойствами.
Двойственные оценки задачи ЛП в канонической базисной форме неотрицательные.
Утверждение является следствием условия оптимальности базисного плана.
Сумма произведений двойственных оценок и правых частей системы ограничений исходной задачи равна оптимальному значению целевой
функции: |
y*b y*b ... |
y* b |
f Opt |
|
|
|
|
|
|
|
|
|
|
|
|
||||
1 1 |
2 |
2 |
|
m |
m |
. |
|
|
|
|
|
|
|
|
|
|
|
|
|
Действительно, используя индексную строку последней симплексной |
|||||||||||||||||||
таблицы, получаем следующее выражение для целевой функции: |
|
|
|
||||||||||||||||
|
n |
6 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
f ( X ) |
|
x |
c |
|
x |
y* x |
x |
|
y* x |
|
|
x |
y* x |
f opt. |
|||||
|
|
k k |
0 |
|
1 1 |
2 2 |
3 3 |
1 4 |
|
5 5 |
3 6 |
|
|||||||
|
k 1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Первый базисный план имеет вид |
X 1 (0;b ;0;b ;0;b ) |
. Подставляя этот план в |
|||||||||||||||||
|
2 |
|
1 |
3 |
|||||||||||||||
целевую функцию и учитывая, что |
f ( X 1 ) |
0, получаем: |
|
|
|
|
|
|
|||||||||||
f ( X 1 ) |
y2*b2 |
|
y1*b1 |
y3*b3 |
f opt. |
0 |
|
y1*b1 |
y2*b2 |
y3*b3 |
|
f opt. . |
|||||||
В частности, с помощью этого свойства можно проверять правильность |
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
17 |
2 |
0 1 |
1 |
0 |
17 |
f Opt |
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
вычислений. Для задачи из пункта 6.2 имеем: 4 |
|
|
|
|
4 |
|
2 |
. |
|||||||||||
Из формулы |
f Opt |
y1*b1 |
y2*b2 ... |
ym* bm |
|
немедленно следует, что при |
изменении правой части ограничения с номером i на единицу оптималь-
ное значение целевой функции изменится на величину двойственной оцен-
ки |
y* |
|
|
|
|
|
|
i (конечно, при условии, что найденный базис остается оптимальным). |
|||||||
|
В рассматриваемом примере увеличение изменение b1 |
на единицу приве- |
|||||
дет к изменению оптимального значения на |
y* |
17 / 4 |
, а изменение |
b |
|||
1 |
|
2 не меня- |
|||||
ет оптимального значения, так как y2* 0 . |
|
|
|
|
|
||
|
Произведение двойственной оценки |
y* |
и значения соответствующей |
||||
|
i |
||||||
|
|
xOpt |
|
|
|
|
|
ей переменной |
k ( i ) в оптимальном плане равно нулю: |
|
|
y* |
xOpt |
0, i 1,..., m |
i |
k ( i ) |
. |
|
||
|
|

|
|
|
33 |
|
|
Действительно, если |
xOpt |
0 |
, то это переменная базисная и поэтому ее |
||
k ( i ) |
|
||||
|
|
|
|
||
оценка yi* 0 . Если же yi* |
|
|
|
Opt |
|
0 , то переменная xk ( i ) ‒ свободная, а так как рас- |
|||||
сматриваемый план базисный, то |
xOpt 0 |
. |
|||
k ( i ) |
|||||
|
|
|
|||
7. Решение задачи производственного планирования симплекс- |
|||||
методом |
|
|
|
|
|
7.1. Предприятие может выпускать 3 вида продукции П1 П3 , используя |
|||||
для этого 4 вида сырья М1 |
М 4 . Количество сырья каждого вида, необходимое |
для производства единицы продукции, запасы сырья и прибыль от реализации единицы продукции приведены в следующей таблице.
|
П1 |
П2 |
П3 |
Запасы |
М1 |
2 |
1 |
2 |
30 |
М 2 |
4 |
2 |
5 |
66 |
М3 |
2 |
3 |
2 |
50 |
М 4 |
2 |
4 |
1 |
48 |
Прибыль |
18 |
10 |
20 |
|
Требуется составить план производства, при выполнении которого прибыль от реализации произведенной продукции будет наибольшей.
Математическая модель этой задачи имеет следующий вид:
f ( X ) 18x1 |
10x2 |
20x3 |
Max; |
|
2x1 |
x2 |
2x3 |
30 |
|
4x1 |
2x2 |
5x3 |
66 |
|
2x1 |
3x2 |
2x3 |
50 |
|
2x1 |
4x2 |
x3 |
48, |
|
x1 0, x2 0, x3 0.
С помощью дополнительных (балансовых) переменных x4 , x5 , x6 , x7 преобразуем стандартную задачу к канонической базисной форме:
f ( X ) |
18x1 10x2 |
20x3 |
|
Max; |
||
2x1 |
x2 |
2x3 |
x4 |
|
|
30 |
4x1 |
2x2 |
5x3 |
|
x5 |
|
66 |
2x1 |
3x2 |
2x3 |
|
|
x6 |
50 |
2x1 |
4x2 |
x3 |
|
|
|
x7 48, |
x1 0, x2 |
0, x3 0, x4 0 x5 |
0, x6 0, x7 0. |

34
Здесь дополнительные переменные x4 , x5 , x6 , x7 играют роль базисных переменных, целевая функция имеет приведенный вид. Составляем первую симплексную таблицу.
|
|
|
|
|
|
|
|
Таблица 1 |
||
Б |
x1 |
x2 |
x3 |
x4 |
x5 |
x6 |
x7 |
B |
bi / aij |
|
x4 |
2 |
1 |
2 |
1 |
0 |
0 |
0 |
30 |
15 |
|
x5 |
4 |
2 |
5 |
0 |
1 |
0 |
0 |
66 |
66 / 5 |
|
x6 |
2 |
3 |
2 |
0 |
0 |
1 |
0 |
50 |
25 |
|
x7 |
2 |
4 |
1 |
0 |
0 |
0 |
1 |
48 |
|
|
f |
18 |
10 |
20 |
0 |
0 |
0 |
0 |
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
Базисные переменные |
x4 , x5 , x6 , x7 . |
Свободные |
переменные |
x1 , x2 , x3 . |
|||||||||||||||||||
Первый базисный план |
X 1 |
(0;0;0;30;66;50; 48) |
невырожденный и неопти- |
||||||||||||||||||||
мальный, так как есть отрицательные оценки 1 |
18, |
2 |
|
|
10 и |
3 |
20 . |
||||||||||||||||
Значение целевой функции f ( X 1 ) |
0. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
В качестве разрешающего выбираем третий столбец. В следующем плане |
|||||||||||||||||||||||
переменная x3 |
будет базисной. Минимальное отношение для элементов третьего |
||||||||||||||||||||||
столбца равно b2 / a23 |
66 / 5 . В следующем плане переменная x2 |
станет свобод- |
|||||||||||||||||||||
ной. Целевая функция увеличится на |
3 |
(b2 / a23 ) |
( |
20) |
66 / 5 |
|
264 . |
||||||||||||||||
Получаем вторую симплекс-таблицу: |
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Таблица 2 |
|||
|
Б |
|
x1 |
|
x2 |
|
|
x3 |
|
x4 |
|
x5 |
|
x6 |
|
|
x7 |
|
|
B |
bi / aij |
|
|
|
x4 |
|
2 / 5 |
1/ 5 |
|
0 |
|
1 |
|
2 / 5 |
|
0 |
|
|
0 |
|
18 / 5 |
|
9 |
|
|||
|
x3 |
|
4 / 5 |
2 / 5 |
|
1 |
|
0 |
|
1/ 5 |
|
0 |
|
|
0 |
|
66 / 5 |
33 / 2 |
|
||||
|
x6 |
|
2 / 5 |
|
11/ 5 |
|
0 |
|
0 |
|
2 / 5 |
|
1 |
|
|
0 |
|
118 / 5 |
|
59 |
|
||
|
x7 |
|
6 / 5 |
|
18 / 5 |
|
0 |
|
0 |
|
1/ 5 |
|
0 |
|
|
1 |
|
174 / 5 |
|
29 |
|
||
|
f |
|
2 |
|
2 |
|
|
0 |
|
0 |
|
4 |
|
0 |
|
|
0 |
|
|
264 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
Базисные переменные x3 , x4 , x6 , x7 . Свободные переменные x1, x2 , x5 . Вто- |
|||||||||||||||||||||||
рой базисный план X 2 |
(0;0;66 / 5;18 / 5;0;118 / 5;174 / 5) |
|
невырожденный и не- |
||||||||||||||||||||
оптимальный, так как есть отрицательные оценки |
1 |
|
|
2 и |
2 |
2 . Значение |
|||||||||||||||||
целевой функции f ( X 2 ) |
264 . |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
В качестве разрешающего столбца выбираем первый. В следующем плане |
|||||||||||||||||||||||
переменная x1 |
будет базисной. Минимальное отношение для элементов первого |

|
|
|
|
|
|
35 |
|
|
|
|
|
|
|
|
|
столбца равно b1 / a11 |
9 . В следующем плане переменная x4 |
станет свободной. |
|||||||||||||
Целевая функция увеличится на |
1 |
(b1 / a11 ) |
( |
2) 9 |
18 . |
|
|
||||||||
Переходим к третьей симплекс-таблице: |
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
Таблица 3 |
||
|
Б |
x1 |
|
x2 |
x3 |
x4 |
|
x5 |
x6 |
|
x7 |
|
B |
bi / aij |
|
|
x1 |
1 |
|
1/ 2 |
0 |
5 / 2 |
|
1 |
0 |
|
0 |
|
9 |
18 |
|
|
x3 |
0 |
|
0 |
1 |
2 |
|
1 |
0 |
|
0 |
|
6 |
|
|
|
x6 |
0 |
|
2 |
0 |
1 |
|
0 |
1 |
|
0 |
|
20 |
10 |
|
|
x7 |
0 |
|
3 |
0 |
3 |
|
1 |
0 |
|
1 |
|
24 |
8 |
|
|
f |
0 |
|
1 |
0 |
5 |
|
2 |
0 |
|
0 |
|
282 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Базисные переменные x1, x3 , x6 , x7 . Свободные переменные x2 , x4 , x5 . Тре- |
||
тий базисный план X 3 (9;0;6;0;0; 20; 24) |
невырожденный и неоптимальный, |
|
так как |
есть отрицательная оценка 2 |
1. Значение целевой функции |
f ( X 3 ) |
282 . |
|
В качестве разрешающего выбираем второй столбец. В следующем плане |
переменная x2 будет базисной. Минимальное отношение для элементов второго
столбца равно b2 / a42 |
8 . В следующем плане переменная x7 станет свободной. |
|||||||||||||||||||||||||||
Значение целевой функции увеличится на |
|
2 |
(b4 / a42 ) |
( 1) |
8 |
|
8 . |
|
|
|||||||||||||||||||
|
Переходим к четвертой таблице: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Таблица 4 |
|
|
|||
|
|
|
Б |
|
|
x1 |
|
x2 |
|
|
x3 |
|
x4 |
|
|
|
x5 |
|
x6 |
x7 |
|
B |
|
|
bi / aij |
|
|
|
|
|
|
x1 |
|
|
1 |
|
0 |
|
|
|
0 |
|
3 |
|
|
|
7 / 6 |
|
0 |
1/ 6 |
|
5 |
|
|
|
|
|
|
|
|
x3 |
|
|
0 |
|
0 |
|
|
|
1 |
|
2 |
|
|
1 |
|
0 |
0 |
|
6 |
|
|
|
|
|
|
|
|
|
x6 |
|
|
0 |
|
0 |
|
|
|
0 |
|
1 |
|
|
|
2 / 3 |
|
1 |
2 / 3 |
|
4 |
|
|
|
|
|
|
|
|
x2 |
|
|
0 |
|
1 |
|
|
|
0 |
|
1 |
|
|
1/ 3 |
|
0 |
1/ 3 |
|
8 |
|
|
|
|
|
|
|
|
|
f |
|
|
0 |
|
0 |
|
|
|
0 |
|
4 |
|
|
|
7 / 3 |
|
0 |
1/ 3 |
|
290 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
|
Базисные переменные x1, x2 , x3 , x6 . Свободные переменные x4 , x5 , x7 . Чет- |
|||||||||||||||||||||||||||
вертый базисный план X 4 |
(5;8;6;0;0; 4;0) невырожденный, оптимальный и |
|||||||||||||||||||||||||||
единственный, |
так как оценки всех свободных переменных положительные. |
|||||||||||||||||||||||||||
Значение целевой функции |
f ( X |
4 ) 290 |
. |
Двойственные оценки |
y* |
4 |
4 |
, |
||||||||||||||||||||
|
|
|
|
|
1 |
|
|
|||||||||||||||||||||
y* |
5 |
7 / 3 |
, |
y* |
0 |
, |
y* |
7 |
1/ 3 |
. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
2 |
|
3 6 |
|
4 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Отбрасывая дополнительные переменные, получаем решение исходной задачи:
X Opt (5;8;6), f Opt 290 .

36
Максимальная прибыль от реализации возникнет при производстве пяти единиц первой, восьми – второй и шести единиц третьей продукции и составит 290 условных единиц.
7.2. В оптимальном плане преобразованной задачи дополнительные пе-
ременные x4 x5 x7 0 . Следовательно, на оптимальном плане исходной задачи первое, второе и четвертое ограничения выполняются как равенства. Дополнительная переменная x6 = 4 – положительная. На оптимальном плане третье ограничение исходной задачи выполняется как строгое неравенство.
При выполнении оптимального плана ресурсы (запасы) первого, второго и четвертого сырья расходуются полностью. Такие ресурсы называются дефицитными. Третий ресурс оказывается недефицитным, так как после выполнения оптимального плана останется его запас в количестве x6 = 4 единиц.
7.3. Нумерация двойственных оценок |
y* |
4 |
4 |
, |
y* |
5 |
7 / 3 |
, |
1 |
|
2 |
|
y* |
6 |
0 |
, |
y* |
7 |
1/ 3 |
однозначно соответствует нумерации ресурсов. В оп- |
3 |
|
4 |
|
тимальном плане преобразованной задачи дополнительная переменная x6 4 – базисная и ее оценка (третья двойственная оценка) равна нулю. При этом третий ресурс является недефицитным.
В оптимальном плане недефицитным ресурсам отвечают нулевые двойственные оценки. Положительным двойственным оценкам отвечают дефицитные ресурсы.
Отметим, что дефицитному ресурсу может отвечать нулевая двойствен-
ная оценка. В рассматриваемом примере двойственные оценки y1* , y2* , y4* , отвечающие дефицитным ресурсам, оказались положительными. Это обстоятельство позволяет рассматривать двойственные оценки как некоторую ценность ресурсов, используемых в производстве. При этом ценность недефицитных ресурсов нулевая – они и так имеются в избыточном количестве.
8.Метод искусственного базиса
8.1.С помощью симплекс-метода задача ЛП в канонической форме может быть преобразована к канонической базисной форме. Потребность в таком преобразовании возникает, в частности, при переходе от задачи в стандартной форме к канонической форме, если среди правых частей системы ограничений имеются отрицательные значения.
Для примера рассмотрим задачу (I) в стандартной форме:
f ( X ) |
2x1 |
3x2 |
x3 x4 |
Max; |
2x1 |
x2 |
x3 |
2x4 |
7 |
x1 |
x2 |
x3 |
x4 |
3 |
2x1 |
2x2 |
x3 |
3x4 |
2, |
x1 |
0, x2 |
0, x3 |
0, x4 |
0. |

37
С помощью неотрицательных балансовых переменных x5, x6, x7 неравен-
ства преобразуются в уравнения. Получаем новую задачу (I I) в канонической форме следующего вида:
f ( X ) 2x1 |
3x2 |
x3 |
x4 |
|
Max; |
|
||||
2x1 |
x2 |
x3 |
|
2x4 |
x5 |
|
|
|
7 |
|
x1 |
x2 |
x3 |
|
x4 |
|
|
x6 |
|
3 |
|
2x1 |
2x2 |
x3 |
|
3x4 |
|
|
|
x7 |
2, |
|
x1 0, x2 |
0, x3 |
0, x4 |
0, x5 |
0, x6 |
0, x7 |
0. |
||||
|
|
2 |
1 |
1 |
2 |
1 |
0 |
0 |
|
|
|
|
7 |
|
|||||||
|
|
1 |
1 |
1 |
1 |
0 |
1 |
0 |
3 |
|
Расширенная матрица |
2 |
2 |
1 |
3 |
0 |
0 |
1 |
2 |
системы ограниче- |
|
|
|
|
|
|
|
|
|
ний задачи (II) содержит единичную матрицу E = (A5 A6 A7), но правые части уравнений не удовлетворяют условию неотрицательности.
После умножения первого и третьего уравнений на минус единицу правые части системы ограничений становятся положительными и задача (II) при-
нимает вид: |
|
|
|
|
|
|
|
f ( X ) |
2x1 |
3x2 |
x3 x4 |
|
Max; |
||
2x1 |
x2 |
x3 |
2x4 |
x5 |
|
|
7 |
x1 |
x2 |
x3 |
x4 |
|
x6 |
|
3 |
2x1 |
2x2 |
x3 |
3x4 |
|
|
x7 |
2, |
x1 0, x2 |
0, x3 |
0, x4 |
0, x5 |
0, x6 |
0, x7 0. |
Правые части этой задачи неотрицательные, но при этом в матрице системы ограничений пропадает базис из единичных векторов. Расширенная матрица этой системы ограничений имеет вид:
2 |
1 |
1 |
2 |
1 |
0 |
0 |
|
7 |
|||||||
1 |
1 |
1 |
1 |
0 |
1 |
0 |
3 |
2 |
2 |
1 |
3 |
0 |
0 |
1 |
2 |
В этой матрице имеется только один единичный столбец, а именно A6.
Введем в первое и третье уравнения неотрицательные искусственные переменные x8 x9 и так, чтобы появились дополнительные единичные столбцы:
2x1 |
x2 |
x3 |
2x4 |
x5 |
|
x8 |
|
7 |
x1 |
x2 |
x3 |
x4 |
|
x6 |
|
|
3 |
2x1 |
2x2 |
x3 |
3x4 |
|
|
x7 |
x9 |
2, |
x1 |
0, x2 |
0, x3 |
0, x4 |
0, x5 0, x6 |
0, x7 |
0, x8 |
0, x9 0. |
Расширенная матрица получившейся системы ограничений имеет следующий вид:

38
2 |
1 |
1 |
2 |
1 |
0 |
0 |
1 |
0 |
|
7 |
|||||||||
1 |
1 |
1 |
1 |
0 |
1 |
0 |
0 |
0 |
3 |
2 |
2 |
1 |
3 |
0 |
0 |
1 |
0 |
1 |
2 |
Она содержит единичную матрицу E = (A8 A6 A9) и неотрицательный вектор правых частей. Базис {A8 A6 A9}, построенный с помощью искусственных переменных x8 и x9, называется искусственным базисом. Полученная система ограничений соответствует канонической базисной форме.
8.2. В отличие от балансовых переменных, которые вводятся в целевую функцию с нулевыми коэффициентами, использование искусственных переменных (искусственного базиса) предполагает построение вспомогательной це-
левой функции f ( X ) и решение вспомогательной задачи f (X ) Max для построенной системы ограничений. В результате решения вспомогательной задачи либо возникает первый опорный план канонической задачи (II), либо выясняется, что исходная задача не имеет ни одного опорного плана.
Существует два способа построения вспомогательной целевой функции. В методе искусственной целевой функции вспомогательная целевая функция
имеет вид f (X , X ) xk1 xk2 ... xkp , где ( X , X ) ‒ план вспомогательной задачи, в котором X (x1, x2 ,..., xn ) ‒ переменные исходной канонической за-
дачи, а X (xk1 ,..., xk p ) ‒ искусственные переменные. В данном примере
f (X , X ) f (x1, x2 ,..., x8 , x9 ) |
x8 |
x9 . |
|
|
В методе штрафных функций вспомогательная целевая функция имеет |
||
вид |
f (X , X ) f (X ) M xk1 |
M |
xk2 ... M xkp , где M обозначает доста- |
точно большое положительное число (значительно превосходящее модули всех чисел задачи). Вспомогательная задача с такой целевой функцией обычно называется M-задачей.
И в первом, и во втором случае целевая функция выражается через свободные переменные, в результате вспомогательная задача принимает каноническую базисную форму. Отметим, что все искусственные переменные на первом шаге входят в число базисных переменных. Вспомогательные задачи решаются симплекс-методом. Возможны три случая.
Если в процессе решения вспомогательной задачи возникает базисный план, в котором все искусственные переменные становятся свободными, то исходная задача принимает каноническую базисную форму, и решение вспомогательной задачи можно не продолжать.

|
|
39 |
|
|
|
|
Для метода искусственной функции приведенное условие означает, что |
||||
найден оптимальный |
план вспомогательной задачи |
вида ( X * ,O) , |
где |
||
O |
(0,...,0) ‒ вектор искусственных переменных, |
а вспомогательная целевая |
|||
функция достигает своего оптимального значения |
f Opt |
f (X * ,O) 0 . |
При |
||
этом план X * является базисным для исходной задачи, а значение исходной це- |
|||||
левой функции равно |
f ( X * ) . |
|
|
|
|
|
Для M - задачи возникает базисный план такого же вида ( X * ,O) , где план |
||||
X * |
является базисным для исходной задачи, а значения исходной и вспомога- |
||||
тельной целевых функций совпадают: f (X * ) f (X * ,O) . |
|
|
Если вспомогательная задача имеет оптимальный план, в котором хотя бы одна искусственная переменная больше нуля, то исходная задача не имеет решения из-за несовместности системы ограничений.
Если M-задача не имеет решения из-за неограниченности вспомогатель-
ной целевой функции на множестве допустимых планов, то и исходная задача не имеет решения или из-за неограниченности исходной целевой функции на множестве допустимых планов, или из-за несовместности системы ограничений.
8.3. Пример построения канонической базисной формы методом искусственной функции. После введения балансовых и искусственных переменных задача (I) из 8.1 принимает вид:
f ( X ) |
2x1 |
3x2 |
x3 |
x4 |
Max; |
|
|
2x1 |
x2 |
x3 |
2x4 |
x5 |
|
x8 |
7 |
x1 |
x2 |
x3 |
x4 |
|
x6 |
|
3 |
2x1 |
2x2 |
x3 |
3x4 |
|
|
x7 |
x9 2, |
x1 |
0, x2 |
0, x3 |
0, x4 |
0, x5 0, x6 |
0, x7 |
0, x8 |
0, x9 |
0. |
|||
Базисные |
переменные |
x8, |
x6 |
и x9. Вспомогательную |
целевую |
функцию |
|||||
f (X , X ) |
x8 x9 |
выражаем через свободные переменные. Из первого и третьего |
|||||||||
уравнений системы ограничений находим x8 7 |
2x1 |
x2 |
x3 |
2x4 |
x5 , |
||||||
x9 2 2x1 2x2 |
x3 |
3x4 |
x7 , и f (X , X ) 4x1 |
x2 |
2x3 |
x4 |
x5 |
x7 9. |
Получаем вспомогательную задачу в канонической базисной форме:

|
|
|
|
40 |
|
|
|
|
|
f ( X , X ) |
4x1 |
x2 2x3 |
x4 |
x5 x7 |
9 |
Max; |
|
||
2x1 |
x2 |
x3 |
2x4 |
x5 |
|
|
x8 |
|
7 |
x1 |
x2 |
x3 |
x4 |
|
x6 |
|
|
|
3 |
2x1 |
2x2 |
x3 |
3x4 |
|
|
x7 |
|
x9 |
2, |
x1 |
0, x2 |
0, x3 |
0, x4 |
0, x5 |
0, x6 |
|
0, x7 |
0, x8 |
0, x9 0. |
8.4. Симплексная таблица для вспомогательной задачи дополняется еще одной индексной строкой – для исходной целевой функции. Этот прием позволяет автоматически получать значения исходной целевой функции в процессе решения вспомогательной задачи.
Таблица 1
Б |
x1 |
x2 |
x3 |
x4 |
x5 |
x6 |
x7 |
x8 |
x9 |
B |
bi / aij |
x8 |
2 |
1 |
1 |
2 |
1 |
0 |
0 |
1 |
0 |
7 |
7 / 2 |
|
|
|
|
|
|
|
|
|
|
|
|
x6 |
1 |
1 |
1 |
1 |
0 |
1 |
0 |
0 |
0 |
3 |
3 |
x9 |
2 |
2 |
1 |
3 |
0 |
0 |
1 |
0 |
1 |
2 |
1 |
|
|
|
|
|
|
|
|
|
|
|
|
f |
2 |
3 |
1 |
1 |
0 |
0 |
0 |
0 |
0 |
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
f |
4 |
1 |
2 |
1 |
1 |
0 |
1 |
0 |
0 |
9 |
|
|
|
|
|
|
|
|
|
|
|
|
|
Оценки вспомогательной задачи находятся в последней строке. Первый базисный план X 1 (0,0,0,0,0,3,0,7, 2) невырожденный и неоптимальный, так
как оценка 1 4 0 . Разрешающий столбец – первый, разрешающая строка – третья. Искусственная базисная переменная х9 в следующей таблице становится свободной, поэтому элементы девятого столбца можно не вычислять. Индексная (предпоследняя) строка исходной задачи преобразуется при переходе к следующему базисному плану по стандартному правилу. При этом происходит пересчет целевой функции исходной задачи. Получаем следующую симплексную таблицу.
Таблица 2
Б |
x1 |
x2 |
x3 |
x4 |
x5 |
x6 |
x7 |
x8 |
x9 |
B |
bi / aij |
x8 |
0 |
3 |
0 |
5 |
1 |
0 |
1 |
1 |
|
5 |
1 |
|
|
|
|
|
|
|
|
|
|
|
|
x6 |
0 |
0 |
3 / 2 |
1/ 2 |
0 |
1 |
1/ 2 |
0 |
|
2 |
4 |
|
|
|
|
|
|
|
|
|
|
|
|
x1 |
1 |
1 |
1/ 2 |
3 / 2 |
0 |
0 |
1/ 2 |
0 |
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
f |
0 |
1 |
2 |
2 |
0 |
0 |
1 |
0 |
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
f |
0 |
3 |
0 |
5 |
1 |
0 |
1 |
0 |
|
5 |
|
|
|
|
|
|
|
|
|
|
|
|
|
План X 2 (1, 0, 0, 0, 0, 2, 0,5, 0) |
не является оптимальным, так как есть |
|
отрицательные оценки 2 |
3 и 4 |
5 . Анализ минимальных отношений |
элементов для второго и четвертого столбца показывает, что в любом случае