数学建模送货路线设计问题.docx
《数学建模送货路线设计问题.docx》由会员分享,可在线阅读,更多相关《数学建模送货路线设计问题.docx(43页珍藏版)》请在冰豆网上搜索。
数学建模送货路线设计问题
送货路线设计问题
摘要:
本文主要讨论的是送货路线的设计问题。
总体的解题思路是将问题中的地点、路线分别抽象成数学中的点、线,然后利用图论的相关知识理论来考虑这些问题。
最后,设计方法程序,并利用Matlab运行,解决问题。
问题一要求根据1-30号货物设计一条最快的送货路线,由于货物的总质量mzong和总体积vzong(mzong=;vzong=)均未超出最大限度50和1,所以,该问题可转化成求最短路问题。
解决方法:
首先,写出每个点的带权邻接矩阵;然后,运用Floyd求任意两点间的最短距离;最后,用H圈构造运算法,并通过矩阵翻转的二边逐次修正法,得到最短距离和最快完成路线图,如下:
o→18→13→24→31→27→39→34→40→45→49→42→43→36→38→32→23→16→14→17→21→26→o
lucheng=+004米t=lucheng/1000*v+t*21/60=小时
问题二设计一条路线,要求在时间允许的条件下,使总路程最小。
解决思路是利用问题一中的方法,结合每个货物的时间限制,最终得到路线图,如下:
o→18→13→24→31→27→39→34→40→45→49→42→43→38→36→32→23→16→14→17→21→26→o
lucheng2=+004t2=lucheng2/1000*v+t*21/60=小时
问题三将1-100号货物全部送到指定地点,mzong=148,vzong=,显然不能一次性送到。
解题思想是根据仓库到各个点的最小距离将地点分为三部分,分别派送。
分完组后在利用第一问的思想给予优化求出最佳的H圈.得到的送货路线分别为:
第一组路线:
o→26→31→27→39→27→36→45→40→47→40→50→49→42→43→38→35→32→23→17→21→o;
第二组路线:
o→26→31→34→40→37→41→44→48→46→33→28→30→22→20→22→29→25→19→24→31→26→o;
第三组路线:
o→21→17→23→16→14→9→10→7→1→6→1→8→3→4→2→5→15→12→11→13→1811→o。
送货时间为:
t3=lucheng/1000*v+t*100/60=小时
关键词:
图论带权邻接矩阵Floyd算法最优Hamilton圈二边逐次修正
一、问题重述
现今社会网络越来越普及,网购已成为一种常见的消费方式,随之物流行业也渐渐兴盛,每个送货员需要以最快的速度及时将货物送达,而且他们往往一人送多个地方,请设计方案使其耗时最少。
现有一快递公司,库房在图1中的O点,一送货员需将货物送至城市内多处,请设计送货方案,使所用时间最少。
该地形图的示意图见图1,各点连通信息见表3,假定送货员只能沿这些连通线路行走,而不能走其它任何路线。
各件货物的相关信息见表1,50个位置点的坐标见表2。
假定送货员最大载重50公斤,所带货物最大体积1立方米。
送货员的平均速度为24公里/小时。
假定每件货物交接花费3分钟,为简化起见,同一地点有多件货物也简单按照每件3分钟交接计算。
现在送货员要将100件货物送到50个地点。
请完成以下问题。
1.若将1~30号货物送到指定地点并返回。
设计最快完成路线与方式。
给出结果。
要求标出送货线路。
2.假定该送货员从早上8点上班开始送货,要将1~30号货物的送达时间不能超过指定时间,请设计最快完成路线与方式。
要求标出送货线路。
3.若不需要考虑所有货物送达时间限制(包括前30件货物),现在要将100件货物全部送到指定地点并返回。
设计最快完成路线与方式。
要求标出送货线路,给出送完所有快件的时间。
由于受重量和体积限制,送货员可中途返回取货。
可不考虑中午休息时间。
以上各问尽可能给出模型与算法。
图1快递公司送货地点示意图
O点为快递公司地点,O点坐标(11000,8250),单位:
米
二、模型假设
1.将仓库视为第51个点,参与计算。
2.送货员在路上无特殊情况,不会因抛锚等现象而耽误时间;
3.同一地点要送多件货物,那么这些物品在同一次中运送;
4.要求到达的时间不包括此次在该点交接的时间;
5.送货员只沿着已知的路线行走;
6.道路是双向的,无单向路线;
7.送货员取货的时间不计。
三、符号说明
1问中涉及到的符号
a各货物号信息(货物号、运送地点、重量、体积和最晚时间)矩阵
b50个位置点的坐标矩阵
c互通点信息矩阵
d任意两相通两点间距离
e对应两相通两点间距离
e1对e进行去重后得到的矩阵
f带权邻接矩阵
D任意两点间最小距离矩阵
u初始H圈
mzong货物的总质量
vzong货物的总体积
luxian最短路线
lucheng最小路程
t1最短时间
t货物交接时所需时间(3分钟)
v送货员的行驶速度(24千米每小时)
2问中涉及到的符号
luxian2最短路线
lucheng2最小路程
t2最短时间
3问中涉及到的符号
luxian3最短路线
lucheng3最小路程
t3最短时间
D3分组矩阵
四、问题的分析与模型的建立
将快递网图中,每个投递点看作图中的一个节点,各投点之间的公路看作图中对应节点间的边,各条路的长度(或行驶时间)看作对应边上的权,所给快递网就转化为加权网络图,问题就转化为在给定的加权网络图中寻找从给定点0出发,行遍所有顶点至少一次再回到O点,使得总权(路程或时间)最小,此即最佳推销员回路问题。
1)问题一是需将30个货物送达21个固定点并返回,O点和另外21个点构成了一个典型的最短路问题。
即先利用Floyd计算两点间的最短距离,再随机构造哈密顿圈,利用优化算法对此H圈优化,使H圈的权最小。
2)问题二本小问是在一问的基础上加入时间的限制,解题思想是以第一问的过程为基础,从随机产生的H圈中选出符合时间要求的多条路线,再从中学出事的路程权重最小的路线。
并检验其是否符合时间的要求。
3)问题三主要是对路线的分组,分组后检验,调整使得每组货物质量小于50kg,体积小于1m3,然后利用问题一,解出每组的最佳H圈。
五、模型的分析与求解
由附录1给定的数据知,前30号货物由于货物的总质量mzong和总体积vzong分别为和均未超出最大限度50和1,显然送货员能够一次带上所有货物到达各送货点,且货物要送达总共为21个,如下:
13,14,16,17,18,21,23,24,26,27,31,32,34,36,38,39,40,42,43,45,49
本模型运用图论中Floyd算法与最佳H圈中的相关结论,建立了关于该类问题的优化模型,将出发点O和21个送货点结合起来构造完备加权图。
用矩阵翻转来实现二边逐次修正,求最佳哈密尔顿圈(H圈)。
由完备加权图,确定初始H圈,列出该初始H圈加点序边框的距离矩阵,然后用二边逐次修正法对矩阵进行“翻转”,就可得到近似最优解的距离矩阵,从而确定近似最佳H圈。
由于用矩阵翻转方法来实现二边逐次修正法的结果与初始圈有关,故为了的到得到较优的计算结果,在用MATLAB编程时,随机搜索出200个初始H圈。
在所有H圈中,找出权最小的一个,即要找的最佳H圈的近似解。
最佳H圈的近似解min{H0,H1,H2,…H99}
送货路线:
o→18→13→24→31→27→39→34→40→45→49→42→43→36→38→32→23→16→14→17→21→26→o
送货时间:
lucheng=+004米t=lucheng/24000+3*21/60=小时
本小问是在一问的基础上加入时间的限制,解题思想是以第一问的过程为基础,从随机产生的H圈中选出符合时间要求的多条路线,即选择符合每个点时间要求的最佳H圈。
为了更有针对性,可将一问的最佳路线作为初始的H圈进行计算。
得到结果,如下:
o→18→13→24→31→27→39→34→40→45→49→42→43→38→36→32→23→16→14→17→21→26→o
lucheng2=+004t2=lucheng2/24000+3*21/60=小时
现根据距离分组,在调整,然后求解。
51号到各个地点的最小距离如下:
12345678910
1006816296104671400416563113628100850977758092
11121314151617181920
696567525295509411558749336212182696813417
21222324252627282930
17971191853954709893413923997142231082013205
31323334353637383940
29296707155495254762446778975621457776885
41424344454647484950
115779751883313943786014312921615806117229928
0→26→31→27→39→27→36→45→40→47→40→50→49→42→43→38→35→32→23→
17→21→0;
0→26→31→34→40→37→41→44→48→46→33→28→30→22→20→22→29→25→19→
24→31→26→0;
0→21→17→23→16→14→9→10→7→1→6→1→8→3→4→2→5→15→12→11→13→18
11→0。
计算三个区域各自送货员走的总路程:
123
计算时间:
++/24000+3/60*100=小时
六、模型的不足及改进的方向
不足:
由于数据量大,且最佳H圈与原始圈的选取有关,只能去近似最佳圈,因此对于第二问随机性很强,只能多设置一下循环次数,以求精确。
第三问的手动画图、分组比较麻烦,要尝试多次才能找出符合要求的点。
参考文献
【1】赵静、但琦,数学建模与数学实验(第3版)高等教育出版社
【2】姜启源、谢金星、叶俊,数学模型,北京:
高等教育出版社,2003
相关程序数据
图1快递公司送货地点示意图
O点为快递公司地点,O点坐标(11000,8250),单位:
米
表1各货物号信息表
货物号
送达地点
重量(公斤)
体积(立方米)
不超过时间
1
13
9:
00
2
18
9:
00
3
31
9:
30
4
26
12:
00
5
21
12:
00
6
14
12:
00
7
17
12:
00
8
23
12:
00
9
32
12:
00
10
38
10:
15
11
45
9:
30
12
43
10:
15
13
39
12:
00
14
45
9:
30
15
42
10:
15
16
43
10:
15
17
32
12:
00
18
36
12:
00
19
27
12:
00
20
24
9:
00
21
31
9:
30
22
27
12:
00
23
26
12:
00
24
34
9:
30
25
40
9:
30
26
45
9:
30
27
49
10:
15
28
32
12:
00
29
23
12:
00
30
16
12:
00
31
1
32
2
33
3
34
4
35
5
36
6
37
7
38
8
39
9
40
10
41
11
42
12
43
13
44
14
45
15
46
16
47
17
48
18
49
19
50
20
51
21
52
22
53
23
54
24
55
25
56
26
57
27
58
28
59
29
60
30
61
31
62
32
63
33
64
34
65
35
66
36
67
37
68
38
69
39
70
40
71
41
72
42
73
43
74
44
75
45
76
46
77
47
78
48
79
49
80
50
81
25
82
46
83
32
84
23
85
20
86
25
87
19
88
41
89
46
90
37
91
32
92
33
93
36
94
38
95
17
96
11
97
15
98
12
99
10
100
7
表250个位置点的坐标
位置点
X坐标(米)
Y坐标(米)
1
9185
500
2
1445
560
3
7270
570
4
3735
670
5
2620
995
6
10080
1435
7
10025
2280
8
7160
2525
9
13845
2680
10
11935
3050
11
7850
3545
12
6585
4185
13
7630
5200
14
13405
5325
15
2125
5975
16
15365
7045
17
14165
7385
18
8825
8075
19
5855
8165
20
780
8355
21
12770
8560
22
2200
8835
23
14765
9055
24
7790
9330
25
4435
9525
26
10860
9635
27
10385
10500
28
565
9765
29
2580
9865
30
1565
9955
31
9395
10100
32
14835
10365
33
1250
10900
34
7280
11065
35
15305
11375
36
12390
11415
37
6410
11510
38
13915
11610
39
9510
12050
40
8345
12300
41
4930
13650
42
13265
14145
43
14180
14215
44
3030
15060
45
10915
14235
46
2330
14500
47
7735
14550
48
885
14880
49
11575
15160
50
8010
15325
表3相互到达信息
序号
位置点1
位置点2
1
1
3
2
1
8
3
2
20
4
2
4
5
3
8
6
3
4
7
4
2
8
5
15
9
5
2
10
6
1
11
7
18
12
7
1
13
8
12
14
9
14
15
9
10
16
10
18
17
10
7
18
11
12
19
12
13
20
12
25
21
12
15
22
13
18
23
13
19
24
13
11
25
14
18
26
14
16
27
14
17
28
14
21
29
15
22
30
15
25
31
16
23
32
17
23
33
18
31
34
19
24
35
20
22
36
21
26
37
21
36
38
21
17
39
22
30
40
23
17
41
24
31
42
25
41
43
25
19
44
25
29
45
27
31
46
28
33
47
29
22
48
30
28
49
30
41
50
31
26
51
31
34
52
32
35
53
32
23
54
33
46
55
33
28
56
34
40
57
35
38
58
36
45
59
36
27
60
37
40
61
38
36
62
39
27
63
40
34
64
40
45
65
41
44
66
41
37
67
41
46
68
42
43
69
42
49
70
43
38
71
44
48
72
44
50
73
45
50
74
45
42
75
46
48
76
47
40
77
48
44
78
49
50
79
49
42
80
50
40
81
O
18
82
O
21
83
O
26
程序
问题一的程序
1.%作图,标号,标距离
clc;
a=[%货物信息数据
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
310
320
330
340
350
360
370
380
390
400
410
420
430
440
450
460
470
480
490
500
510
520
530
540
550
560
570
580
590
600
610
620
630
640
650
660
670
680
690
700
710
720
730
740
750
760
770
780
790
800
810
820
830
840
850
860
870
880
890
900
910
920
930
940
950
960
970
980
990
1000
];
b=[%货物坐标数据
19185500
21445560
37270570
43735670
52620995
6100801435
7100252280
871602525
9138452680
10119353050
1178503545
1265854185
1376305200
14134055325
1521255975
16153657045
17141657385
1888258075
1958558165
207808355
21127708560
2222008835
23147659055
2477909330
2544359525
26108609635
271038510500
285659765
2925809865
3015659955
31939510100
321483510365
33125010900
34728011065
351530511375
361239011415
37641011510
381391511610
39951012050
40834512300
41493013650
421326514145
431418014215
44303015060
451091514235
46233014500
47773514550
4888514880
491157515160
50801015325
51110008250
];
c=[%连通数据
113
218
3220
424
538
634
742
8515
952
1061
11718
1271
13812
14914
15910
161018
17107
181112
191213
201225
211215
221318
231319
241311
251418
261416
271417
281421
291522
301525
311623
321723
331831
341924
352022
362126
372136
382117
392230
402317
412431
422541
432519
442529
452731
462833
472922
483028
493041
503126
513134
523235
533223
543346
553328
563440
573538
583645
593627
603740
61383