ФГБОУ ВПО
ДВГУПС
Кафедра «Управление эксплуатационной
работой».
Математическое моделирование систем и
процессов
Гр.232
Выполнил: Илларионова Е.В
Проверил: Широков А.П.
2013
Транспортная задача закрытого типа, представленная в матричной форме, с ограничениями пропускной способности.
Условие: Для заданного варианта транспортной задачи в матричной форме с ограничениями пропускной способности необходимо найти оптимальный план, при котором будет выполняться условие наименьшего суммарного пробега порожних вагонов. Для этого необходимо:
1) Составить математическую модель задачи;
2) Составить начальный план, проверить по условию вырождения, рассчитать суммарные вагоно-километры порожнего пробега;
) Решить задачу методом потенциалов, рассчитать суммарные вагоно-километры порожнего пробега оптимального плана;
) Сравнить начальный и оптимальный варианты.
Исходные данные транспортной задачи с ограничениями
|
Станция отправления |
Ресурсы, тыс.т |
Станция назначения(потребности) |
|||||||||
|
|
|
В1 |
В2 |
В3 |
В4 |
В5 |
В6 |
В7 |
В8 |
В9 |
|
|
|
|
Объем потребления, тыс.т |
|||||||||
|
|
|
18 |
20 |
12 |
15 |
10 |
15 |
20 |
25 |
15 |
|
|
А1 |
40 |
25 10 |
30 7 |
45 |
60 |
35 |
55 |
60 |
90 |
75 |
|
|
А2 |
25 |
45 |
58 |
46 |
75 |
48 |
54 |
65 |
58 |
95 |
|
|
А3 |
25 |
35 8 |
35 |
40 |
65 |
40 |
80 |
35 |
65 |
55 |
|
|
А4 |
15 |
48 |
50 |
30 6 |
55 |
45 |
45 |
55 |
50 |
40 |
|
|
А5 |
45 |
40 |
52 |
38 |
40 |
70 |
37 8 |
30 10 |
75 |
45 |
|
математическое моделирование транспортная задача
Транспортная задача является закрытой, так как суммарные запасы ресурсов
равны суммарным потребностям:
+25+25+15+45=18+20+12+15+10+15+20+25+15
=150
Математическая модель задачи:
. Целевая функция:
2. Система ограничения
3. Условие не отрицательности
Для решения задачи необходимо составить опорный план. Исходный опорный план составляется с помощью метода наименьшего критерия в столбце.
Начальный опорный план, построенный методом наименьшего критерия в столбце.
|
Станция отправления |
Ресурсы, тыс.т |
Станция назначения(потребности) |
||||||||||
|
|
|
В1 |
В2 |
В3 |
В4 |
В5 |
В6 |
В7 |
В8 |
В9 |
||
|
|
|
Объем потребления, тыс.т |
||||||||||
|
|
|
18 |
20 |
12 |
15 |
10 |
15 |
20 |
25 |
15 |
||
|
А1 |
40 |
25 10 10 |
30 7 7 |
45 |
60 |
35 10 |
55 |
60 4 |
90 |
75 9 |
||
|
А2 |
25 |
45 |
58 |
46 |
75 |
48 |
54 |
65 |
58 25 |
95 |
||
|
А3 |
25 |
35 8 8 |
35 13 |
40 |
65 |
40 |
80 |
35 4 |
65 |
55 |
||
|
А4 |
15 |
48 |
50 |
30 6 6 |
55 |
45 |
45 7 |
55 2 |
50 |
40 |
||
|
А5 |
45 |
40 |
52 |
38 6 |
40 15 |
70 |
37 8 8 |
30 10 10 |
75 |
45 6 |
||
Значение целевой функции составит:
F=25*10+30*7+35*10+60*4+75*9+58*25+35*8+35*13+35*4+30*6+45*7+55*2+38*6+40*15+37*8+30*10+45*6=6349 (ваг-км)
Скорректированный начальный план, построенный методом наименьшего критерия в столбце
|
Станция отправления |
Ресурсы тыс.т |
Станция назначения (потребности) |
Ui |
|||||||||
|
|
|
В1 |
В2 |
В3 |
В4 |
В5 |
В6 |
В7 |
В8 |
В9 |
|
|
|
|
|
Недостаток порожних вагонов |
|
|
|
|||||||
|
|
|
18 |
20 |
12 |
15 |
10 |
15 |
20 |
25 |
15 |
|
|
|
А1 |
40 |
25 10 10 |
30 7 7 |
45Н23 |
60Н10 |
35 10 |
55 |
60 4 |
90 |
75 9 |
50 |
|
|
А2 |
25 |
45 Н3 |
58 |
46 0 |
75 |
48 |
54 |
65 |
58 25 |
95 |
72 |
|
|
А3 |
25 |
35 8 8 |
35 13 |
40 Н3 |
65 |
40 |
80 |
35 4 |
65 |
55 |
75 |
|
|
А4 |
15 |
48Н17 |
50 Н5 |
30 6 6 |
55Н10 |
45 |
45 7 |
55 2 |
50Н25 |
40Н30 |
55 |
|
|
А5 |
45 |
40 0 |
52 |
38 6 |
40 15 |
70 |
37Н17 8 8 |
30 10 10 |
75 |
45 6 |
80 |
|
|
Vj |
120 |
110 |
118 |
120 |
85 |
100 |
110 |
130 |
125 |
|
||
1.F=25*10+30*7+35*10+60*4+75*9+46*0+58*25+35*8+35*13+35*4+30*6+45*7+55*2+40*0+38*6+40*15+37*8+30*10+45*6=6349
(ваг-км)
.Проверка исходного опорного плана на условие «вырождения».
Кзб≤m+n-1
Где Кзб - число занятых базисных клеток; m- число строк
матрицы (пунктов отправления); n-число столбцов (пунктов назначения).
≤ 5+9-1
≤13
Условие вырождения выполняется, но план является «вырожденным». Для устранения «вырождения» назначаются фиктивные перевозки х51=0,х23=0.
. Потенциалы столбцов определяются следующим образом:
Vi=Ui+cij
Где cij- критерий расстояния в заданной клетке.
Потенциалы строк определяются через занятые клетки, связанных со столбцами, получившими потенциал по формуле:
i=Vj
- cij
Проверка на оптимальность. План считается оптимальным, если соблюдается следующие условия:
i - Ui ≤ cij, при xij=0 (клетка свободна)
Vi
- Ui = cij, при
(клетка
базисная).i - Ui ≥cij, при xij=
(насыщенная клетка)
Формальное правило улучшения плана:
а) начиная с клетки с нарушением, двигаясь по горизонталям и вертикалям ходом «шахматной ладьи», строят замкнутый контур с вершинами в базисных клетках;
б) начиная с клетки с нарушением, нумеруют вершины контура (направление обхода контура значения не имеет);
в) в четный вершинах находится минимальная перевозка.
г) для балансировки матрицы в нечетных вершинах контура найденное значение прибавляется к значениям перевозок в этих клетках (с учетом возможных ограничений в этих клетках), в четных вершинах - вычитается из значений перевозок.
Формальное правило улучшения плана при нарушении условия оптимальности для насыщенных клеток:
а) начиная с насыщенной клетки с нарушением, двигаясь по горизонталям и вертикалям (ходом шахматной ладьи) строят замкнутый контур с вершинами в базисных клетках.
б) начиная с клетки с нарушением, нумеруют вершины контура (направление обхода контура значения не имеет);
в) в нечетных вершинах находится минимальная перевозка;
г) для балансировки матрицы в четных вершинах контура найденное значение прибавляется к значениям перевозок в этих клетках (с учетом возможных ограничений в этих клетках), в нечетных вершинах-вычитается из значений перевозок.
|
Станция отправления |
Ресурсы тыс.т |
Станция назначения (потребности) |
Ui |
|||||||||
|
|
|
В1 |
В2 |
В3 |
В4 |
В5 |
В6 |
В7 |
В8 |
В9 |
|
|
|
|
|
Недостаток порожних вагонов |
|
|
|
|||||||
|
|
|
18 |
20 |
12 |
15 |
10 |
15 |
20 |
25 |
15 |
|
|
|
А1 |
40 |
25 10 10 |
30 7 7 |
45Н23 |
60Н10 |
35 10 |
55Н25 |
60 6 |
90 |
75 7 |
20 |
|
|
А2 |
25 |
45 Н3 |
58 |
46 0 |
75 |
48 |
54Н4 |
65 |
58 25 |
95 |
||
|
А3 |
25 |
35 8 8 |
35 13 |
40 |
65 |
40 |
80 |
35 4 |
65 |
55 |
45 |
|
|
А4 |
15 |
48 |
50 |
30 6 6 |
55 |
45 |
45 7 |
55 |
50 |
40 2 |
55 |
|
|
А5 |
45 |
40 0 |
52 |
38 6 |
40 15 |
70 |
37 8 8 |
30 10 10 |
75 |
45 6 |
50 |
|
|
Vj |
90 |
80 |
88 |
90 |
55 |
100 |
80 |
100 |
95 |
|
||
Скорректированный план перевозок (первая итерация)
F=25*10+30*7+35*10+60*6+75*7+46*0+58*25+35*8+35*13+35*4+30*6+45*7+40*2+40*0+38*6+40*15+37*8+30*10+45*6=6289
(ваг-км)
Проверка на условие «вырождения»:
13≤5+9-1
=13
Скорректированный план перевозок (вторая итерация)
|
Станция отправления |
Ресурсы тыс.т |
Станция назначения (потребности) |
Ui |
|||||||||
|
|
|
В1 |
В2 |
В3 |
В4 |
В5 |
В6 |
В7 |
В8 |
В9 |
|
|
|
|
|
Недостаток порожних вагонов |
|
|
|
|||||||
|
|
|
18 |
20 |
12 |
15 |
10 |
15 |
20 |
25 |
15 |
|
|
|
А1 |
40 |
25 10 10 |
30 7 7 |
45Н20 |
60Н7 |
35 10 |
55 7 |
60 6 |
90 |
75 |
50 |
|
|
А2 |
25 |
45 Н3 |
58 |
46 0 |
75 |
48 |
54 |
65 |
58 25 |
95 |
69 |
|
|
А3 |
25 |
35 8 8 |
35 13 |
40 0 |
65 |
40 |
80 |
35 4 |
65 |
55 |
75 |
|
|
А4 |
15 |
48 |
50 |
30 6 6 |
55 |
45 |
45 |
55 |
50 |
40 9 |
82 |
|
|
А5 |
45 |
40 0 |
52 |
38 6 |
40 15 |
70 |
37 8 8 |
30 10 10 |
75 |
45 6 |
77 |
|
|
Vj |
117 |
110 |
115 |
117 |
85 |
105 |
110 |
127 |
122 |
|
||