2979
.pdfm |
|
|
aij yi cj, |
j 1, n. |
(13) |
i 1
Лемма 1. Если X – некоторый опорный план исходной задачи (9)-(11), а Y – произвольный план двойственной задачи (12)-(13), то значение целевой функции исходной задачи на плане X всегда не превосходит значения целевой функции на плане Y, т.е.
F(X ) F * (Y ).
Лемма 2. Если F(X *) F * (Y*) для некоторых планов
X* и Y* задач (9)-(11), то X* - оптимальный план исходной задачи, а Y* - оптимальный план двойственной задачи.
Теоремы двойственности.
Теорема 1. Если одна из пары двойственных задач (9)- (11) или (12)-(13) имеет оптимальный план, то и другая имеет оптимальный план и значения целевых функций задач при их оптимальных планах равны между собой, т.е. Fmax F *min .
Если же целевая функция одной из пары двойственных задач не ограничена (для исходной задачи (9)-(11) – сверху, для двойственной (12)-(13) – снизу), то другая задача вообще не имеет планов.
Теорема 2. |
План |
X * x* , x* ,..., x* |
задачи (9)-(11) и |
||
план Y * y* , y* ,..., y* |
|
1 2 |
n |
|
|
задачи (12)-(13) являются оптималь- |
|||||
1 2 |
m |
|
|
|
|
ными планами этих задач тогда и только тогда, когда j 1, n выполняется равенство
|
m |
|
|
|
aij |
yi* c j x*j |
0. |
i 1 |
|
|
Геометрическая интерпретация двойственных за-
дач.
Если число переменных в прямой и двойственной задачах, образующих данную пару, равно двум, то, используя геометрическую интерпретацию ЗЛП, можно легко найти
141
решение данной пары задач. При этом имеет место один из следующих трёх взаимно исключающих друг друга случаев:
1)обе задачи имеют планы;
2)планы имеет только одна задача;
3)для каждой задачи двойственной пары множество планов пусто.
Пример.
Для задачи: F 2x1 7x2 max
2x1 3x2 14 |
|
|
x1 x2 8 |
|
x1 , x2 0
составить двойственную задачу и найти решение обеих задач.
Решение.
Двойственная задача имеет вид:
F * 14y 8y |
2 |
min |
|||
|
|
1 |
|
|
|
2 y1 |
y2 |
2 |
|
||
|
3y1 |
y2 |
7 |
|
|
|
|
y1 , y2 0
Как в исходной, так и в двойственной задаче, число неизвестных равно двум. Следовательно, их решение можно найти, используя геометрическую интерпретацию ЗЛП.
Как видно из представленных ниже рисунков, максимальное значение целевая функция исходной задачи принимает в точке В(2;6). Следовательно, X*=(2;6) является оптимальным планом, при котором Fmax 2 2 7 6 46 .
Минимальное значение целевая функция двойственной задачи принимает в точке D(1;4). Значит, Y*=(1;4) является оптимальным планом двойственной задачи, при котором
Fmin* 14 1 8 4 46.
142
y2 I
F1*
4 |
F * |
|
0 |
4
0 1 II
y2 |
I |
|
F1*
4 |
F * |
|
0 |
0 1 II
I : 2 y1 y2 |
2 |
n |
||
II : |
3y1 y2 7 |
|||
|
F1* : 14 y1 8y2 46
F0* : 14 y1 8y2 56
|
|
|
y1 |
|
I : 2 y1 y2 |
2 |
n |
||
II : |
3y1 y2 7 |
|||
|
F1* : 14 y1 8y2 46
F0* : 14 y1 8y2 56
4
y1
Таким образом, значения целевых функций исходной и двойственной задач при их оптимальных планах равны между собой. Из этих рисунков также видно, что при всяком плане исходной задачи значение целевой функции не превос-
143
ходит 46, а значение целевой функции двойственной задачи при любом её плане не меньше 46.
Таким образом, при любом плане исходной задачи значение целевой функции не превосходит значения целевой функции двойственной задачи при её произвольном плане.
Соответствие между первоначальными переменными одной из двойственных задач и дополнительными переменными другой задачи.
С помощью теорем двойственности можно, решив симплекс-методом исходную задачу, найти оптимум целевой функции и оптимальное решение двойственной задачи.
Если среди векторов P1 , P2 , ..., Pn , составленных из ко-
эффициентов при неизвестных в системе уравнений (10), имеется m единичных, то компоненты оптимального плана двойственной задачи совпадают с соответствующими элементами (m+1)-ой строки столбцов единичных векторов в последней симплекс-таблице, если данный коэффициент c j 0 ,
и равны сумме соответствующего элемента этой строки и c j , если c j 0 .
Переменные исходной задачи
|
Первоначальные |
|
Дополнительные |
||||
x1 |
x2 |
...x j |
...xn |
xn 1 |
xn 2 |
...xn j |
...xn m |
|
|
|
|
|
|
|
|
ym 1 |
ym 2 |
...ym j |
...ym n |
y1 |
y2 |
...y j |
...ym |
Дополнительные Первоначальные
Переменные двойственной задачи
Если пара двойственных задач симметричная, т.е. система ограничений исходной задачи содержит неравенства
144
вида ≤, то компоненты оптимального плана двойственной задачи совпадают с соответствующими числами (m+1)-ой строки последней симплекс-таблицы решения исходной задачи.
Указанные числа стоят в столбцах векторов, соответствующих дополнительным переменным.
Пример. Для ЗЛП:
F x1 x2 max
2x1 |
|
x2 |
2 |
|
|
x1 |
2x2 |
2 |
|
|
||||
|
x1 |
|
x2 |
5 |
|
x1 , x2 0
составить двойственную задачу и найти решения обеих задач симплекс-методом.
Решение.
Запишем исходную ЗЛП в каноническом виде.
|
F x1 x2 max |
|||
2x1 x2 x3 2 |
||||
|
x1 x2 x4 2 |
|||
|
||||
|
x1 x2 x5 5 |
|||
|
||||
|
|
|
|
|
x j |
0, |
j 1,5. |
Двойственная ЗЛП:
F* 2y1 2y2 5y3 min
2 y1 y2 y3 1 |
|||||
|
y1 |
2 y2 |
y3 1 |
||
|
|||||
|
y j |
0, |
j 1,3. |
||
|
|
|
|
|
|
Решим исходную ЗЛП симплекс-методом.
145
|
i |
Базис |
|
Сб |
|
|
|
P0 |
|
|
1 |
|
-1 |
|
0 |
|
0 |
|
|
0 |
|
|
||||||
|
|
|
|
|
|
|
P1 |
|
P2 |
|
P3 |
|
|
P4 |
|
|
P5 |
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
1 |
|
P3 |
|
0 |
|
|
|
|
2 |
|
|
-2 |
|
1 |
|
1 |
|
0 |
|
|
0 |
|
|
||||
|
2 |
|
P4 |
|
0 |
|
|
|
|
2 |
|
|
1 |
|
-2 |
|
0 |
|
1 |
|
|
0 |
* |
|
||||
|
3 |
|
P5 |
|
0 |
|
|
|
|
5 |
|
|
1 |
|
1 |
|
0 |
|
0 |
|
|
1 |
|
|
||||
|
4 |
|
|
|
|
|
|
|
|
|
0 |
|
|
-1 |
|
1 |
|
0 |
|
0 |
|
|
0 |
|
|
|||
|
1 |
|
P3 |
|
0 |
|
|
|
|
6 |
|
|
0 |
|
-3 |
|
1 |
|
2 |
|
|
0 |
|
|
||||
|
2 |
|
P1 |
|
1 |
|
|
|
|
2 |
|
|
1 |
|
-2 |
|
0 |
|
1 |
|
|
0 |
** |
|
||||
|
3 |
|
P5 |
|
0 |
|
|
|
|
3 |
|
|
0 |
|
3 |
|
0 |
|
-1 |
|
|
1 |
|
|
||||
|
4 |
|
|
|
|
|
|
|
|
|
2 |
|
|
0 |
|
-1 |
|
0 |
|
1 |
|
|
0 |
|
|
|||
|
1 |
|
P3 |
|
0 |
|
|
|
|
9 |
|
|
0 |
|
0 |
|
1 |
|
1 |
|
|
1 |
|
|
||||
|
2 |
|
P1 |
|
1 |
|
|
|
|
4 |
|
|
1 |
|
0 |
|
0 |
|
-2/3 |
|
|
2/3 |
|
|
||||
|
3 |
|
P2 |
|
-1 |
|
|
|
1 |
|
|
0 |
|
1 |
|
0 |
|
-1/3 |
|
|
1/3 |
|
|
|||||
|
4 |
|
|
|
|
|
|
|
|
|
3 |
|
|
0 |
|
0 |
|
0 |
|
2/3 |
|
|
1/3 |
|
|
|||
|
|
|
|
|
2 |
|
|
5 |
|
2 |
|
|
|
|
|
|
|
3 |
3 |
|
|
|
||||||
где * 0 |
min |
|
|
, |
|
|
|
|
|
2 |
и ** 0 |
min |
|
|
|
|
1. |
|
|
|||||||||
1 |
|
1 |
|
3 |
|
|
||||||||||||||||||||||
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
3 |
|
|
|
|||||||||
X * 4; 1 ; Fmax |
|
3. |
|
Для |
того, |
чтобы найти |
оптимальный |
план двойственной задачи, установим соответствие между основными и дополнительными (базисными) переменными прямой и двойственной задач.
основные |
|
|
|
|
|
дополнительные |
|
|
x1 , x2 |
x3 x4 x5 |
Прямая задача |
|
|
|
|
|
|
|
|
y4 |
y5 |
y1 |
y2 |
y3 |
Двойственная задача |
|
|
|
|||
дополнительные |
основные |
|
|
||
Компоненты оптимального плана двойственной задачи |
|||||
расположены |
в |
последней |
строке |
последней симплекс- |
таблицы в столбце векторов P3 , P4 , P5 , образующих единичный базис исходной ЗЛП.
Y*= (0, 2/3, 1/3); F *min 2 0 2 23 5 13 3 Fmax .
146
5. Транспортная задача (ТЗ)
Имеются m пунктов отправления (поставщиков) грузов А1, А2 , ..., Аm , на которых сосредоточены запасы какого-
либо однородного груза в |
объёмах соответственно |
||
|
|
|
|
a1 , a2 , ..., am . Величины ai , i 1, m, |
определяют максимально |
возможные размеры поставок груза с пунктов отправления.
m
Суммарный запас груза у поставщиков составляет ai .
i 1
Кроме того, имеется n пунктов назначения B1 , B2 , ..., Bn , которые подали заявки на грузы соответственно
n
b1 , b2 , ..., bn . Суммарная величина заявок составляет bj .
j 1
Стоимость перевозки одной единицы груза от поставщика Ai к потребителю B j обозначим через сij (транспорт-
ный тариф), через xij – количество единиц груза, перевозимого от поставщика Ai к потребителю B j .
Необходимо составить оптимальный план, т.е. найти такие значения объёмов перевозок грузов xij , i 1, m ,
j 1, n , от поставщиков Ai к потребителям B j , чтобы вывез-
ти все грузы от поставщиков, удовлетворить заявки каждого потребителя, обеспечить минимальные транспортные расходы на перевозку груза.
Все исходные данные ТЗ можно записать в виде таблицы, которая называется транспортной.
147
Пункты от- |
|
Пункты назначения B j |
|
Запасы |
||||
правления |
|
|
||||||
|
|
|
|
… |
|
ai |
||
Ai |
B1 |
B2 |
… B j |
Bn |
||||
|
|
|
|
|
|
|
||
A1 |
с11 |
с12 |
... |
с1 j |
… |
с1n |
a1 |
|
x11 |
x12 |
x |
x |
|||||
|
|
|
|
|||||
|
|
|
|
1 j |
|
1n |
|
|
A2 |
c21 |
c22 |
… |
c2 j |
… |
c2n |
a2 |
|
x21 |
x22 |
x2 j |
x2n |
|||||
|
|
|
|
|||||
… |
… |
… |
… |
… |
… |
… |
… |
|
|
|
|
|
|
||||
Ai |
ci1 |
ci 2 |
… |
cij |
… |
cin |
ai |
|
xi1 |
xi 2 |
xij |
xin |
|||||
|
|
|
|
|||||
… |
… |
… |
… |
… |
… |
… |
… |
|
|
|
|
|
|
||||
Am |
cm1 |
cm2 |
… |
cmj |
… |
cmn |
am |
|
xm1 |
xm2 |
xmj |
xmn |
|||||
|
|
|
|
|||||
|
|
|
|
|
|
|
m |
|
Заявки |
|
|
|
|
|
|
ai |
|
b1 |
b2 |
… |
b j |
… bn |
i 1 |
|||
|
||||||||
b j |
n |
|||||||
|
|
|
|
|
|
bj |
||
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
j 1 |
|
Математически постановка ТЗ заключается в опреде- |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
лении плана перевозок - |
|
матрицы |
X |
xij |
, i 1, m, |
j 1, n, |
|||||||
которая удовлетворяет условиям: |
|
|
|
|
|
|
|
|
|||||
n |
|
|
|
|
|
|
|
|
|
|
|
|
|
xij |
ai , |
i 1, m ; |
|
|
|
|
(1) |
||||||
j 1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
m |
|
|
|
|
|
|
|
|
|
|
|
|
|
xij |
bj , |
j 1, n ; |
|
|
|
|
(2) |
||||||
i 1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
xij 0, |
i 1, m, |
j 1, n ; |
(3) |
||||||||||
|
|
148 |
|
|
|
|
|
|
|
|
|
|
|
и обеспечивает min целевой функции
m |
n |
|
F cij xij min . |
(4) |
|
i 1 |
j 1 |
|
Для решения ТЗ необходимо и достаточно, чтобы суммарные запасы груза в пунктах отправления были равны сумме заявок пунктов назначения:
m
ai
i 1
n |
|
bj . |
(5) |
j 1
|
Модель ТЗ, удовлетворяющая условию (5), называется |
||||||||
закрытой, в противном случае – открытой. |
|
|
|
||||||
|
В |
случае |
превышения |
запаса |
|
над |
заявками |
||
m |
n |
|
|
|
|
|
|
|
|
ai |
bj |
вводится фиктивный (n+1)-ый пункт назначения |
|||||||
i 1 |
j 1 |
|
|
|
|
|
|
|
|
|
|
|
m |
n |
|
|
|
|
|
с потребностью bn 1 ai |
bj |
и соответствующие тари- |
|||||||
|
|
|
i 1 |
j 1 |
|
|
|
|
|
|
|
|
|
||||||
фы считаются равными нулю: ci, n 1 0, i 1, m . |
|
||||||||
|
|
|
m |
n |
|
|
|
|
|
|
При условии ai |
bj вводится фиктивный (m+1)- |
|||||||
|
|
|
i 1 |
j 1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
n |
m |
ый пункт отправления с запасом груза am 1 bj |
ai и |
||||||||
|
|
|
|
|
|
|
|
j 1 |
i 1 |
соответствующие |
тарифы |
принимаются |
равными нулю: |
cm 1, j 0, j 1, n .
Теорема. Число линейно независимых уравнений в системе (1) – (2) (другими словами, ранг системы (1) – (2)) равно m+n-1. Если в опорном плане перевозок число отличных от нуля компонент равно m+n-1, то он называется невы-
рожденным, если меньше – то вырожденным.
149
Методы определения опорного плана ТЗ.
Метод северо-западного направления (северозападного угла).
На каждом шаге рассматривают первый из оставшихся пунктов отправления и первый из оставшихся пунктов назначения. Заполнение клеток условий начинается с левой верхней клетки для неизвестного x11 (северо-западный угол) и
заканчивается для неизвестного xmn , т.е. идёт как бы по диагонали таблицы с северо-запада на юго-восток.
Пример. На 3 базы A1, A2 , A3 поступил однородный
груз в количествах, соответственно равных 140, 180, 160 единиц. Этот груз требуется перевезти в 5 пунктов назначения B1, B2 , B3 , B4 , B5 соответственно в количествах 60, 70, 120,
130 и 100 единиц. Тарифы перевозок единицы груза с каждого из пунктов отправления в соответствующие пункты назначения указаны в таблице.
B j |
B1 |
B2 |
B3 |
B4 |
B5 |
|
ai |
Ai |
|
||||||
|
|
|
|
|
|
|
|
A1 |
2 |
3 |
4 |
2 |
|
4 |
140 |
|
|
|
|
|
|
||
A2 |
8 |
4 |
1 |
4 |
|
1 |
180 |
|
|
|
|
|
|
||
A3 |
9 |
7 |
3 |
7 |
|
2 |
160 |
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
bj |
60 |
70 |
120 |
130 |
100 |
|
480 |
|
|
|
|
|
|
|
|
Найти план перевозок ТЗ методом северо-западного направления.
Решение.
Здесь число пунктов отправления m=3, а число пунктов назначения n=5. Следовательно, опорный план задачи
150