Транспортная логистика
.pdfны в матричные. С помощью моделей производится оптимизация маршрутов транспортных потоков. С этой целью используется извест ная «транспортная задана линейного программирования»:
целевая функция
тп
;=1 7=1
при условиях
т |
п |
X ХЧ = B j ’ |
X ХУ = А‘ ’ ХУ - 0 ’ |
,=1 |
м |
где - стоимость перевозки груза из пункта / в пункт j на единицу перевозимого груза;
Ху - искомая величина грузопотока из пункта / в пунктj; Bj - количество получаемых грузов в пунктахj;
А, - количество отправления грузов из пунктов /;
1 |
, 2 |
пункты отправления груза; |
1 |
, 2„..J,...,п - пункты назначения груза. |
|
Кроме транспортных расходов в качестве критериев оптимиза ции могут быть использованы расстояния перемещения, объем транспортной работы и т. д.
8.2.1.Метод потенциалов для решения транспортной задачи
вматричной форме с ограничениями пропускной способности
Задача 1
1. Построить оптимальный план перевозок каменного угля с пя ти станций Aj (i = 1.. .5), обслуживающих шахты, до девяти крупных потребителей, имеющих подъездные пути В) (/' = 1,..., 9).
В задаче № 1 по своему варианту студенту необходимо привести все заполненные матрицы, полученные в процессе решения, с дан ными ресурсов станций отправления и потребностей станций на значения, а также заполнить матрицы с начальным и с оптималь ным планом перевозки.
2. Определить объем тонно-километровой работы для всех пла нов начального и оптимального планов перевозки груза.
280
Исходные данные
Данные о наличии ресурсов на станциях отправления А; приве дены в табл. 8.3; данные о размерах прибытия груза Bj на станциях назначения - в табл. 8.4. Расстояние перевозки от каждой ;-й стан ции отправления до каждой у'-й станции назначения указано в форме 1 в правом верхнем углу каждой клетки. В левом верхнем углу формы 1 указаны ограничения пропускной способности, если тако вые имеются. Матрица расстояний и ограничений пропускной спо собности принимается одинаковой для всех вариантов.
Таблица 8.3 Ресурсы станции отправления Аь тыс. т____________
Номер стан |
|
|
|
|
Варианты |
|
|
|
|
||
ции отправ |
|
2 |
|
|
5 |
6 |
7 |
8 |
9 |
0 |
|
ления («ро |
1 |
3 |
4 |
||||||||
ка матрицы) |
150 |
150 |
150 |
145 |
155 |
400 |
400 |
150 |
400 |
145 |
|
А\ |
|||||||||||
а 2 |
150 |
145 |
150 |
155 |
150 |
155 |
150 |
155 |
155 |
400 |
|
As |
145 |
155 |
155 |
150 |
400 |
145 |
150 |
150 |
150 |
150 |
|
л 4 |
155 |
400 |
145 |
150 |
150 |
150 |
155 |
145 |
150 |
150 |
|
А5 |
400 |
150 |
400 |
400 |
145 |
150 |
145 |
400 |
145 |
155 |
|
ИТОГО |
1000 |
1000 |
1000 |
1000 |
1000 |
1000 |
1000 |
1000 |
1000 |
1000 |
|
отправлено |
|||||||||||
|
|
|
|
|
|
|
|
|
|
||
Таблица 8.4 Объем прибытия груза В, на станции назначения, тыс. т
Номер стан |
|
|
|
|
Варианты |
|
|
|
|
||
ции отправ |
|
|
|
|
|
|
7 |
|
9 |
0 |
|
ления (стро |
1 |
2 |
3 |
4 |
5 |
6 |
8 |
||||
ка матрицы) |
100 |
90 |
100 |
100 |
100 |
85 |
85 |
115 |
95 |
80 |
|
в, |
|||||||||||
Вг |
100 |
90 |
100 |
100 |
100 |
85 |
100 |
120 |
80 |
115 |
|
в3 |
І 0 0 І |
110 |
95 |
100 |
100 |
80 |
120 |
80 |
100 |
80 |
|
В4 |
100 |
110 |
95 |
150 |
100 |
80 |
120 |
85 |
115 |
95 |
|
Bs |
100 |
100 |
105 |
90 |
100 |
100 |
80 |
100 |
80 |
100 |
|
в. |
150 |
150 |
155 |
90 |
І00 |
100 |
110 |
50 |
170 |
50 |
|
В-, |
100 |
100 |
100 |
110 |
90 |
120 |
80 |
165 |
85 |
155 |
|
В* |
100 |
100 |
100 |
110 |
90 |
120 |
80 |
120 |
50 |
120 |
|
Д, |
150 |
150 |
150 |
150 |
220 |
230 |
225 |
165 |
225 |
205 |
|
ИТОГО |
1000 |
1000 |
1000 |
1000 |
1000 |
1000 |
1000 |
1000 |
1000 |
1000 |
|
отправлено |
|||||||||||
|
|
|
|
|
|
|
|
|
|
||
281
Методические указания к выполнению задачи 1
1. Решение задачи начинают с составления исходной матрицы (форма 1). Для решения задачи в столбец Л/ заносят объемы ресур сов станций отправления, а в строку Bj - объемы прибытия грузов на станции назначения.
После записи ресурсов и потребностей груза по своему варианту на исходную матрицу для решения задачи строят начальный план любым известным способом (к примеру, методом минимального элемента).
После составления начального плана необходимо проверить ба ланс по строкам и столбцам матрицы. Число базисных клеток (ба зисной клеткой называется клетка, имеющая корреспонденцию, т. е. объем перевозки из i-й станции отправления на j -ю станцию назна чения), должно быть:
К - т + п ~ \ ,
где т - число строк матрицы; п - число столбцов.
Для рассматриваемой матрицы АТ = 9 + 5 - 1 = 13. Если число ба зисных клеток больше числа К, то начальный план составлен не правильно. Клетки, где величина перевозки равна пропускной спо собности, называются небазисными.
В форме 1 в базисных клетках корреспонденция перевозки отме чена жирным шрифтом.
Форма 1
Пример построения начального плана для задачи № 1 на матрице
X |
Л ,~ « 5 |
Вг ~ 100 |
В , - 120 |
Д ,” |
120 |
S i - 8 0 |
|
«SO |
B*>~ SO |
& »«225 |
|
||
|
80 |
40 |
90 |
|
103 |
150 |
30 50 |
?S |
80 |
|
90 |
|
|
«4,-iSO |
|
|
|
*000 |
|||||||||
|
|
|
|
1 0 * |
|
30 |
|
|
|
3 0 |
|||
|
|
|
|
|
|
|
|
|
|
||||
|
|
Ю |
30 |
43 |
|
40 |
25 |
65 |
|
*0 |
10 |
so |
ІШ5 |
|
|
т |
|
|
|
|
3S |
|
3C |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
4 » -1 5 5 |
10 |
20 |
35 |
75 |
|
160 |
90 |
80 |
70 |
40 |
|
60 |
1030 |
|
|
|
|
|
|
|
|
|
SB |
|
75 |
||
|
|
|
|
|
|
|
|
|
|
|
|
||
Л 4 - И 5 |
|
45 |
5 |
15 |
|
30 |
IW |
40 |
75 |
JO |
|
20 |
1070 |
|
|
t m |
|
|
|
|
|
|
|
|
43 |
||
|
|
|
|
|
|
|
|
|
|
|
|
||
Л . -»400 |
|
І5 |
15 25 |
10 |
20 |
*5 |
25 |
|
iso |
70 |
|
«5 |
, « , f |
|
|
|
130 |
20 |
4& |
e a |
5 8 |
|
|
U |
|||
|
|
|
|
|
|
) |
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
*5 |
101S |
1075 |
I05S |
1105 |
І030 |
lo s s |
1025 |
1070 |
1 |
0 ) |
|
||
|
|
|
|
.... |
|
|
|
|
|
|
|
|
|
282
Объем тонно-километровой работы равен 42725 т-км.
2. Оптимальный план перевозок на заданной матрице составим по методу потенциалов.
Любой допустимый план является оптимальным тогда, когда каждой строке и каждому столбцу могут быть присвоены некото рые числа U, и У„ называемые потенциалами и отвечающие сле дующим условиям:
|
Уг -и ,< с ц |
д л я |
лг„ = 0; |
(8.1) |
|
Vі - Ui = Су |
д л я |
0 < Хц < di{. |
(8.2) |
|
Vj -U j> с,, |
д л я |
Хц = dih |
(8.3) |
где Vj - |
потенциал j -го столбца; |
|
|
|
U, - |
потенциал /-й строки; |
|
|
|
c,/ - |
расстояние перевозки от і-го поставщика до /-го потребителя; |
|||
х,у - |
корреспонденция (размер перевозки) от г-го поставщика до |
|||
у'-го потребителя; |
|
|
|
|
d,j - |
величина пропускной способности. |
|
||
3. |
Присвоение потенциалов начинаем со строки, в которой среди |
|||
базисных клеток имеется максимальное расстояние. Этой строке присваиваем потенциал, равный 1000. Затем, используя условие оп тимальности (8.2), находим потенциалы остальных строк и столб цов следующим образом:
для у-го столбца
для /-Й строки
4. После присвоения всем строкам и столбцам потенциалов определяем, имеются ли нарушения неравенств (8.1) и (8.3), по форійуле
Ht/ = Vj - Ui - Cij.
Для свободных клеток нарушения являются положительными по своей величине; для клеток с поставкой, равной пропускной спо собности, - отрицательными.
283
5. Улучшение допустимого плана начинаем с клетки, имеющей
максимальное (по модулю) нарушение Z/™3*. Для этой клетки строим
замкнутый контур, в который входят только базисные клетки и вы бранная клетка с нарушением: из выбранной клетки с нарушением проводим ломаную линию, заканчивающуюся в той же клетке, двига ясь аналогично движению шахматной ладьи, направление движения при этом изменяем под прямым углом только в базисных клетках.
Следует заметить, что для каждой клетки с нарушением сущест вует только один контур улучшения плана. Нумерация клеток кон тура начинается с клетки с нарушением. Если клетка с нарушением свободная, то ей присваиваем № 1. Для клеток с поставками, рав ными пропускной способности, нумерация начинается с нуля. Да лее номера присваиваются по ходу контура. Число клеток в контуре всегда четное.
В найденном замкнутом контуре определяем корреспонденцию улучшения допустимого плана на данном этапе решения. Коррес понденция улучшения плана находится следующим образом:
* У = m in [х0четн, ( d ij - * „ ) _ , ] .
На величину хул изменяются все корреспонденции контура, на чиная с клегки с нарушением: уменьшаются корреспонденции, за писанные в четных клетках, и увеличиваются корреспонденции, записанные в нечетных клетках контура.
6. После пересмотра корреспонденции необходимо пересоста вить систему потенциалов всей матрицы и проверить соблюдение условия оптимальности (8.1) и (8.3). Если небазисные клетки удов летворяют этим условиям, то найдено оптимальное решение. Если имеются нарушения условий оптимальности, то расчет на матрице следует продолжить до тех пор, пока все клетки матрицы не будут удовлетворять условиям (8.1)-(8.3). Число нарушений и их величина всегда стремятся к нулю.
Используя вышеизложенный алгоритм и произведя необходимое количество итераций, получим оптимальный план перевозок (в ле вом верхнем углу каждой клетки находится величина нарушения, если таковая имеется), (табл. 8.5—8.12).
284
Оттенком I |
I отмечены четные клетки замкнутого контура; |
1 \ |
I - нечетные; |
[ |
I - клетки с максимальным (по модулю) на |
рушением;
■небазисные клетки.
Таблица 8.5 Построение оптимального плана перевозок. Шаг 1
|
f t - |
*5 |
f t |
-100 |
f t * 120 |
f t - |
120 |
By-9 0 |
f t - 110 |
|
X |
|
|
|
|
|
|
|
|
|
|
J50 |
|
90 |
35 |
40 |
90 |
|
|
|
150 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
m |
l |
|
|
І Ш |
Лгш150 |
|
10 |
|
30 |
45 |
|
|
|
|
65 |
ts |
|
|
|
|
I T |
|
|
|
|
|
|
|
|
|
|
|
ш |
і і |
|
||
|
|
|
|
|
|
|
|
|
||
Aj * 155 |
25 |
20 |
)0 |
35 |
75 |
|
l«0 |
|
90 |
SO |
|
|
|
|
|
|
|
|
|
|
|
Л -М 5 |
|
45 |
|
5 |
35 |
5 |
30 |
|
no |
40 |
|
|
|
m |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
400 |
55 |
15 |
45 |
25 |
10 |
|
|
|
|
SO |
|
|
|
|
m |
|
|
w |
m |
* • |
|
|
|
|
|
|
few rS S ? |
|||||
|
|
|
|
|
|
& ш .< д |
|
|||
К |
1075 |
1075 |
1015 |
1105 |
|
10Э0 |
1015 |
|||
* i 8
75
И
TO
75
20
9«
1025
ft-to |
f t -225 |
Ц |
|
SO |
■ y ? : M |
1000 |
|
|
|||
30 |
30 |
(065 |
|
|
|
||
40 |
«0 |
юзо |
|
M |
75 |
||
|
|||
30 |
20 |
1070 |
|
|
45 |
||
|
|
||
70 |
|
1005 |
|
|
|
||
1070 |
1090 |
|
Таблица 8.6 Построение оптимального плана перевозок. Шаг 2
|
S.-85 |
B-i« 100 |
5, = 120 |
84»120 |
Bs-m |
tf6~!I0 |
|
«,-eo |
tf,-S0 |
Й ,-225 |
|
|
|||||
,4,“ t50 |
<0 |
35 |
40 |
90 |
|
|
|
150 |
|
|
|
5 |
75 |
80 |
, / / / / / / |
1000 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
30 |
|
f |
|
P |
25 |
H |
65 |
|
|
|
30 |
30 |
|
|
|
|
|
45 |
|
10 |
|
|
|
|
|
|||||||
Ai=150 |
|
|
|
|
|
|
|
|
1065 |
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
25 20 i 10 |
35 |
75 |
|
|
160 |
90 |
30 |
80 |
Я |
70 |
40 |
«о |
|
|
||
155 |
|
|
|
|
1030 |
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
80 |
IS |
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
A -145 |
45 |
|
5 |
35 |
|
5 |
30 |
110 |
30 |
40 |
|
|
75 |
30 |
20 |
1070 |
|
|
106 |
|
|
|
|
|
1 |
|
|
|
|
|
. 45 . |
|
|||
|
|
|
|
|
|
|
|
|
|
i . |
. |
. |
|||||
* =400 |
j j j j j j l |
|
25 |
10 |
1 |
|
Щ |
25 |
|
80 |
|
|
20 ! |
70 |
|
1005 |
|
|
I |
120 |
|
Ш |
W |
(0 |
|
! |
|
я |
|
ш |
|
||||
|
|
|
|
|
|
|
|
|
|
||||||||
Kr |
1075 |
j 1075 |
1015 |
j |
1105 |
10» |
1 10*5 |
1 |
|
г |
|
|
|
|
|||
1025 1 |
1070 |
1090 |
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
285
Таблица 8.7 Построение оптимального плана перевозок. Шаг 3
Таблица 8.8
286
Таблица 8.9 Построение оптимального плана перевозок. Шаг 5
Таблица 8.10 Построение оптимального плана перевозок. Шаг 6
287
Таблица 8.11 Построение оптимального плана перевозок. Шаг 7
Таблица 8.12
Итоговый оптимальный план перевозок
288
8.2.2.Метод потенциалов для решения транспортной задачи
всетевой форме
Задача 2
Построить оптимальный план перевозки груза на сети (рис. 8.4) от трех станций отправления до девяти станций назначения.
Рис. 8.4. Сеть транспортных узлов
Для данной задачи студент приводит две сети: с начальным пла ном перевозок и с оптимальным планом перевозок.
Исходные данные
Размеры груза со станции отправления приведены в табл. 8.13, а прибытия на станцию назначения - в табл. 8.14. Для всех вариантов задан один и тот же полигон (см. рис. 8.4), на котором производится решение задачи. В знаменателе для звеньев 1-7, 1-11, 2-9, 3-10 зада ны ограничения пропускной способности.
Методические указания к выполнению задачи 2
1. Решение задачи в сетевой форме начинают с составления на чального плана, который не допускает встречных перевозок на уча стках заданного полигона. Начальный (или любой допустимый)
289
