Жолобов Ввведение в Математическое 2008
.pdf
|
15 |
20 |
25 |
15 |
10 |
10 |
20 |
|
50 |
|
|
|
|
10 |
|
|
|
|
|
|
|
|
|
|
5 |
5 |
1 |
4 |
2 |
1 |
15 |
|
20 |
|
30 |
|
50 |
|
|
E32=0 |
|
|
|
|
|
|
|
3 |
5 |
5 |
1 |
4 |
5 |
2 |
|
|
|
|
|
|
|
40 |
|
|
|
|
|
|
|
|
|
|
20 |
7 |
5 |
1 |
3 |
5 |
4 |
|
|
100 |
|
20 |
|
|
|
|
|
|
|
|
|
|
|
|
20 |
1 |
2 |
3 |
4 |
6 |
6 |
|
|
40 |
30 |
|
|
|
100 |
|
|
|
|
|
|
|
|
|
15 |
20 |
25 |
15 |
10 |
10 |
20 |
|
50 |
|
|
|
|
10 |
|
|
|
|
|
|
|
|
|
|
5 |
5 |
1 |
4 |
2 |
1 |
15 |
|
20 |
|
30 |
|
50 |
|
|
E34=2 |
|
|
|
|
|
|
|
3 |
5 |
5 |
1 |
4 |
5 |
2 |
|
|
|
|
|
|
|
40 |
|
|
|
|
|
|
|
|
|
|
20 |
7 |
5 |
1 |
3 |
5 |
4 |
|
|
100 |
|
20 |
|
|
|
|
|
|
|
|
|
|
|
|
20 |
1 |
2 |
3 |
4 |
6 |
6 |
|
|
40 |
30 |
|
|
|
100 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
15 |
20 |
25 |
15 |
10 |
10 |
20 |
|
50 |
|
|
|
|
10 |
|
|
|
|
|
|
|
|
|
|
5 |
5 |
1 |
4 |
2 |
1 |
15 |
|
20 |
|
30 |
|
50 |
|
|
E46=-2 |
|
|
|
|
|
|
|
3 |
5 |
5 |
1 |
4 |
5 |
2 |
|
|
|
|
|
|
|
40 |
|
|
|
|
|
|
|
|
|
|
20 |
7 |
5 |
1 |
3 |
5 |
4 |
|
|
100 |
|
20 |
|
|
|
|
|
|
|
|
|
|
|
|
20 |
1 |
2 |
3 |
4 |
6 |
6 |
|
|
40 |
30 |
|
|
|
100 |
|
|
|
|
|
|
|
|
2.2.13. Дополнительные ограничения в постановке транспортной задачи
При постановке и нахождении решения конкретных транспортных задач часто бывает необходимо учитывать дополнительные ограничения, которые еще не рассматривались.
261
Рассмотрим подробно некоторые из таких ограничений.
1. Запрещена перевозка из пункта Ai в пункт Bj . В этом случае вводится очень высокий тариф (М) перевозки из пункта Ai в пункт Bj единицы продукции:
cij=M .
2. Из пункта Ai в пункт Bj требуется обязательно перевести точно dij единиц продукции. В клетку на пересечении i-й строки и j-го столбца вносится величина dij. Корректируется запас ai и потреб-
ность bj : ai ai dij , bj bj dij . В дальнейшем эта клетка счита-
ется свободной, а для того чтобы соответствующая переменная не попала в состав базисных переменных оптимального решения, этой клетке приписывается очень большой тариф: cij=M .
3. Из пункта Ai в пункт Bj требуется перевести не менее ij единиц продукции. Считаем, что из пункта Ai в пункт Bj уже перевезено ij единиц продукции. Уменьшаем запас ai и потребность bj на ве-
личину ij : |
ai ai dij , |
bj |
bj ij . Далее задача решается |
обычным методом, после чего корректируется полученное решение. Смысл этой корректировки заключается в следующем. Если в
оптимальном решении переменная xij принимает значение xij , то в окончательном решении ей приписывается значение xij xij ij .
4. Из пункта Ai в пункт Bj требуется перевести не более ij единиц продукции. В исходную таблицу вводится дополнительный
столбец Bj (рис. 2.11):
262
|
j |
|
|
|
|
j |
j |
|
||
|
|
|
|
|
|
|
|
|
|
|
|
c1j |
|
|
|
|
c1j |
c1j |
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
cij |
|
|
|
cij |
M |
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
cmj |
|
|
|
|
cmj |
cmj |
|
||
|
|
|
|
|
|
|
|
|
||
|
bj |
|
таблица |
|
|
ij |
bj- ij |
|
||
|
Старая |
|
|
|
Новая |
таблица |
|
Рис.2.11. Моделирование дополнительного ограничения
В окончательном решении принимается xij xij xij , i 1, m .
Естественно, что открывается возможность комбинировать
эти приемы и решать задачи с достаточно «серьезными» ограничениями20.
Пример 2.12
Построить транспортную модель с дополнительными ограничениями.
|
1 |
2 |
3 |
4 |
5 |
|
|
|
|
|
|
|
|
1 |
5 |
3 |
2 |
4 |
8 |
160 |
|
|
|
|
|
110 |
|
|
|
|
|
|
|
|
2 |
7 |
6 |
5 |
3 |
1 |
90 |
|
|
|
|
|
||
|
|
|
|
|
|
|
3 |
8 |
9 |
4 |
5 |
2 |
140 |
|
|
|
|
|
80 |
|
|
|
|
|
|
|
|
|
90 |
60 |
80 |
70 |
90 |
|
|
|
10 |
|
|
30 |
|
Дополнительные ограничения:
Из A1 в B2 должно быть перевезено не менее 50 ед. (x12 50); Из A3 в B5 должно быть перевезено не менее 60 ед. (x35 60); Из A2 в B4 должно быть перевезено не более 40 ед. (x24 40).
Учтем теперь ограничение x24 40. Для этого нужно ввести но-
20 В частности, достаточно просто учитывается ограничение dн xij dв.
263
вый столбец |
4 . |
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
2 |
3 |
|
4 |
|
5 |
|
6= 4 |
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
5 |
3 |
2 |
|
4 |
|
8 |
|
4 |
|
110 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
7 |
6 |
5 |
|
3 |
|
1 |
|
M |
|
90 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
3 |
|
8 |
9 |
4 |
|
5 |
|
2 |
|
5 |
|
80 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
90 |
10 |
80 |
|
40 |
|
30 |
|
70-40 |
|
|
|
|
|
|
|
|
30 |
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Теперь обычным способом решаем эту задачу (находим исходное опорное решение, а далее, например, используем метод потенциалов).
Однако помним, что полученные значения x12 и x35 нужно будет увеличить соответственно на 50 и 60 единиц.
2.2.14. Транспортная модель с промежуточными пунктами
Транспортная модель с промежуточными пунктами соответствует реальной ситуации, когда между исходными и конечными пунктами перевозок имеются промежуточные пункты для временного хранения грузов (транзитные пункты). Это – более общая, чем обычная транспортная модель, где перевозки осуществляются непосредственно между пунктами отправления и назначения.
Проведем сквозную нумерацию всех пунктов (исходных, ко-
нечных, транзитных): 1, 2,…, R.
Транспортную модель представим орграфом, вершинам которого соответствуют пункты 1, 2,…, R.
В том случае, если существует возможность перевозки про-
дукции непосредственно из |
k в l, эти вершины свяжем взве- |
||
шенной дугой (рис. 2.12), |
направленной из k |
в сторону |
l |
( k,l {1, 2,..., R}, k l ). В качестве веса этой дуги |
примем сk,l |
– |
|
стоимость перевозки единицы продукции из k в l. |
|
||
k |
l |
|
|
сk,l
Рис.2.12. Дуга между промежуточными пунктами
264
Модель предполагает возможность сосредоточения запасов продукции в любом пункте. Кроме того, каждый пункт может иметь собственную потребность в определенном количестве продукции. Этот факт будем отображать на графе следующим образом
(рис. 2.13):
ak |
k |
сk,l |
l |
bl |
|
|
Рис.2.13. Запасы и потребности в модели
Здесь ak – количество единиц продукции, сосредоточенной в пункте k ; bl – количество единиц продукции, составляющее по-
требность пункта l .
В рассматриваемой модели перевозки транзитом могут осуществляться через любые пункты (в соответствии с направлением взвешенных дуг на орграфе), даже через некоторые пункты назначения.
Поэтому все множество пунктов 1, 2,…, R можно разбить на три класса.
1.Пункты, которым соответствуют как входящие, так и исходящие взвешенные дуги, назовем транзитными пунктами (ТП).
2.Пункты, которым соответствуют только исходящие взвешенные дуги, назовем истинными пунктами отправления (ИПО).
3.Пункты, которым соответствуют только входящие взвешенные дуги, назовем истинными пунктами назначения (ИПН).
Представленную таким орграфом транспортную модель с промежуточными пунктами преобразуем в обычную (закрытую) транспортную модель с помощью введения, т.н. буфера.
Объем буфера (B) должен быть таким, чтобы вместить объем всего предложения (или спроса):
B = Общий объем предложения (спроса).
Объемы спроса (предложения) перечисленных выше трех подмножеств пунктов определяются следующим образом:
объем предложения ИПО = объем исходного предложения;
объем спроса ИПН = объем исходного спроса;
объем предложения ТП = объем исходного предложения + B;
265
объем спроса ТП = объем исходного спроса + B.
При построении закрытой транспортной модели в качестве поставщиков принимаются все ИПО и все ТП. В качестве потребителей – все ТП и все ИПН.
Для того чтобы запретить перевозки между пунктами, не связанными в исходном орграфе взвешенными дугами, соответствующим маршрутам (запрещенным), приписывается достаточно высокая стоимость (M). Маршрутам, соответствующим петлям (из одного ТП в тот же самый ТП) приписывается нулевая стоимость.
Ниже рассматривается пример решения транспортной задачи с промежуточными пунктами.
Пример 2.13
|
|
|
|
T3 |
800 |
|
|
|
|
8 |
|
1000 |
P1 |
3 |
T1 |
5 |
|
|
|
4 |
|
6 |
|
|
|
|
|
|
|
|
|
|
7 |
T4 |
900 |
|
|
|
|
||
|
|
2 |
|
4 |
|
1200 |
P2 |
5 |
T2 |
3 |
|
|
|
|
|
9 |
|
|
|
|
|
D1 |
500 |
Здесь P1 и P2 – ИПО; T1, T2, T3, T4 – ТП; D1 – ИПН.
B = 1000+1200=800+900+500=2200.
Строим транспортную таблицу и решаем задачу:
266
|
|
T1 |
|
T2 |
|
T3 |
T4 |
D1 |
|
|
|
|
P1 |
3 |
|
4 |
|
M |
M |
M |
|
1000 |
|
|
|
|
1000 |
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
P2 |
2 |
|
5 |
|
M |
M |
M |
|
1200 |
|
|
1200 |
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
T1 |
0 |
|
7 |
|
8 |
6 |
M |
|
2200 |
|
|
1000 |
|
|
|
800 |
400 |
|
|
|
||
|
|
|
|
|
|
|
|
|
|||
|
T2 |
M |
|
0 |
|
M |
4 |
9 |
|
2200 |
|
|
|
|
1200 |
|
|
1000 |
|
|
|
||
|
|
|
|
|
|
|
|
|
|
||
|
T3 |
M |
|
M |
|
0 |
5 |
M |
|
2200 |
|
|
|
|
|
|
2200 |
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
T4 |
M |
|
M |
|
M |
0 |
3 |
|
2200 |
|
|
|
|
|
|
|
1700 |
500 |
|
|
||
|
|
|
|
|
|
|
|
|
|
||
|
|
2200 |
|
2200 |
|
3000 |
3100 |
500 |
|
11000 |
|
Решение задачи на графе: |
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
T3 |
800 |
||
|
|
|
|
|
|
|
800 |
|
|
|
|
1000 |
P1 |
|
|
T1 |
400 |
|
|
|
|
||
|
|
1000 |
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
T4 |
900 |
||
|
|
1200 |
|
|
|
1000 |
|
|
|
|
|
1200 |
P2 |
|
|
T2 |
|
500 |
|
||||
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
D1 |
500 |
Контрольные вопросы и задачи к разделу 2.2
1.Сравните метод Фогеля и метод плавающих зон с точки зрения вычислительных затрат на решение задачи и степени приближенности получаемого приближенного решения к оптимальному.
2.Дана транспортная задача. Эта задача решается методом Фогеля.
Имеем значение ЦФ Zисх. Введем в эту задачу дополнительное ограничение (например, запретим какой-либо маршрут). Получим значение ЦФ Zдоп. Можно ли утверждать следующее и почему:
а)Zисх Zдоп б)Zисх Zдоп?
267
3.Возможен ли в методе потенциалов случай зацикливания, аналогичный тому, который возможен в симплекс-методе в случае вырожденности задачи?
4.Докажите утверждение, что косвенная стоимость – это коэффициент, с которым свободная переменная входит в целевую функцию, если все базисные переменные выразить через свободные.
5.Каким образом можно привести открытую транспортную задачу с промежуточными пунктами к закрытой?
6.Решите задачи методами северо-западного угла, наименьшей стоимости, Фогеля и методом плавающих зон:
|
1 |
2 |
3 |
4 |
|
|
1 |
11 |
10 |
6 |
17 |
55 |
|
17 |
11 |
16 |
7 |
|||
|
40 |
|||||
|
9 |
14 |
15 |
9 |
||
3 |
80 |
|||||
50 |
40 |
45 |
40 |
|||
|
|
|||||
|
1 |
2 |
3 |
4 |
|
|
1 |
16 |
6 |
15 |
17 |
45 |
|
10 |
7 |
6 |
6 |
|||
2 |
50 |
|||||
10 |
7 |
9 |
10 |
|||
3 |
115 |
|||||
55 |
55 |
45 |
55 |
|||
|
|
|||||
|
1 |
2 |
3 |
4 |
|
|
1 |
11 |
7 |
12 |
9 |
60 |
|
12 |
15 |
16 |
17 |
|||
2 |
45 |
|||||
15 |
12 |
11 |
14 |
|||
3 |
115 |
|||||
60 |
60 |
55 |
45 |
|||
|
|
7. В процессе решения задачи с помощью метода потенциалов была получена следующая транспортная таблица:
|
1 |
2 |
|
3 |
4 |
|
1 |
16 |
15 |
|
12 |
|
11 |
|
20 |
|
25 |
10 |
||
|
|
|
||||
2 |
8 |
15 |
|
16 |
|
16 |
40 |
20 |
|
|
|
|
|
|
|
|
|
|
||
3 |
17 |
6 |
|
12 |
|
7 |
|
|
|
|
|
60 |
|
|
|
|
|
|
|
|
|
|
|
268 |
|
|
Продолжите решение задачи.
8. Решите задачу распределительным методом, применив метод наименьшей стоимости для поиска исходного опороного решения:
|
|
1 |
2 |
3 |
4 |
|
|
|
1 |
16 |
10 |
14 |
9 |
60 |
|
|
|
|
|
|
|
||
|
2 |
7 |
16 |
6 |
15 |
45 |
|
|
|
|
|
|
|
||
|
3 |
15 |
8 |
12 |
17 |
115 |
|
|
|
|
|
|
|
||
|
|
60 |
60 |
55 |
45 |
|
|
9. Постройте закрытую транспортную модель на основании сле- |
|||||||
дующей таблицы: |
|
|
|
|
|
|
|
с дополнительными ограничениями |
|
|
|||||
|
|
1 |
2 |
3 |
4 |
|
|
|
1 |
11 |
11 |
16 |
7 |
55 |
|
|
|
|
|
|
|
||
|
2 |
9 |
10 |
10 |
15 |
60 |
|
|
|
|
|
|
|
||
|
3 |
16 |
15 |
15 |
10 |
105 |
|
|
|
|
|
|
|
||
|
|
60 |
55 |
50 |
55 |
|
|
а) запрещена перевозка из пункта A1 в пункт B3; |
|||||||
б) из пункта A2 в пункт B4 необходимо перевезти ровно 30 |
|||||||
единиц продукции; |
|
|
|
|
|
|
|
в) из пункта A1 в пункт B2 необходимо перевезти не менее 45 |
|||||||
единиц продукции; |
|
|
|
|
|
|
|
г) из пункта A3 в пункт B1 можно перевезти не более 5 еди- |
ниц продукции.
10. Решите транспортную задачу с промежуточными пунктами:
|
300 |
|
|
|
4 |
5 |
|
|
|
|
|
|
|
3 |
|
200 |
4 |
7 |
150 |
|
|
|
|
|
5 |
4 |
|
|
2 |
3 |
|
|
|
||
|
|
|
|
|
300 |
|
350 |
269
2.3. Метод отсечения
2.3.1.Основные понятия метода отсечения
Всвоей основе метод отсечения использует тот факт, что, если все вершины многогранного множества D имеют целочисленные координаты, то и оптимальное опорное решение соответствующей задачи ЛП будет удовлетворять требованию целочисленности.
Рассмотрим ЛЦП-задачу:
|
|
n |
|
|
|
|
|
|
|
|
|
|
|
|
cj xj max |
|
|
|
(2.23) |
||||||
|
|
j 1 |
|
|
|
|
|
|
ЛП |
|
|
|
|
|
n |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|||
|
aij xj ai , |
i 1, m |
|
|
ЛЦП |
(2.24) |
||||||
D j 1 |
|
|
|
|
|
|
|
|
|
|
||
Dц |
|
|
|
|
|
|
|
|
|
|
|
|
xj 0, |
j 1,n |
|
|
|
(2.25) |
|||||||
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
xj целые, |
j 1,n |
|
|
|
(2.26) |
||||||
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
Допустимое множество этой задачи, определенное ограниче-
ниями (2.24)-(2.26), обозначим Dц.
Очевидно, Dц – множество целых точек выпуклого многогранного множества D, определенного ограничениями (2.24)-(2.25).
С целью упрощения последующих выкладок будем считать, что D – замкнутое ограниченное множество. То есть Dц – конечное множество.
Обозначим Dц* выпуклую линейную оболочку множества Dц.
Множество |
D* |
называется выпуклой линейной оболочкой множе- |
|||
|
ц |
|
|
|
|
ства D {x1 |
, x2 ,..., xN }, если содержит все возможные выпуклые |
||||
ц |
|
|
|
|
|
линейные комбинации вида: |
|
|
|
||
|
|
N |
|
|
N |
|
|
i xi , |
i 0 ( i |
|
), i 1. |
|
|
1, N |
|||
|
|
i 1 |
|
|
i 1 |
|
|
|
270 |
|
|