数学建模送货线路设计问题答案仅供参考.docx
《数学建模送货线路设计问题答案仅供参考.docx》由会员分享，可在线阅读，更多相关《数学建模送货线路设计问题答案仅供参考.docx（58页珍藏版）》请在冰点文库上搜索。
装订线
第九届西安电子科技大学数学建模竞赛暨
全国大学生数学建模竞赛选拔赛题目
A（B）题
密封号
2010年5月4日
剪切线
密封号
2010年5月4日
通信工程学院第队
队员1
队员2
队员3
姓名
邓晓光
谭正中
刘春燕
班级
送货路线设计问题
1、问题重述
现今社会网络越来越普及，网购已成为一种常见的消费方式，随之物流行业也渐渐兴盛，每个送货员需要以最快的速度及时将货物送达，而且他们往往一人送多个地方，请设计方案使其耗时最少。
现有一快递公司，库房在图1中的O点，一送货员需将货物送至城市内多处，请设计送货方案，使所用时间最少。
该地形图的示意图见图1，各点连通信息见表3，假定送货员只能沿这些连通线路行走，而不能走其它任何路线。
各件货物的相关信息见表1，50个位置点的坐标见表2。
假定送货员最大载重50公斤，所带货物最大体积1立方米。
送货员的平均速度为24公里/小时。
假定每件货物交接花费3分钟，为简化起见，同一地点有多件货物也简单按照每件3分钟交接计算。
现在送货员要将100件货物送到50个地点。
请完成以下问题。
1.若将1~30号货物送到指定地点并返回。
设计最快完成路线与方式。
给出结果。
要求标出送货线路。
2.假定该送货员从早上8点上班开始送货，要将1~30号货物的送达时间不能超过指定时间，请设计最快完成路线与方式。
要求标出送货线路。
3.若不需要考虑所有货物送达时间限制（包括前30件货物），现在要将100件货物全部送到指定地点并返回。
设计最快完成路线与方式。
要求标出送货线路，给出送完所有快件的时间。
由于受重量和体积限制，送货员可中途返回取货。
可不考虑中午休息时间。
2、问题分析
送货路线问题可以理解为：
已知起点和终点的图的遍历问题的合理优化的路线设计。
图的遍历问题的指标：
路程和到达的时间，货物的质量和体积，以及最大可以负载的质量和体积。
在路线的安排问题中，考虑所走的路程的最短即为最合理的优化指标。
对于问题二要考虑到所到的点的时间的要求是否满足题意即采用多次分区域的假设模型从而找出最优的解
对于问题三则要考虑到体积和质量的双重影响，每次到达后找到达到最大的体积和质量的点然后返回，再依次分析各个步骤中可能存在的不合理因素达到模型的进一步合理优化得到最合理的解。
3、模型假设与符号说明
、模型的假设
（1）、到同一地点的货物要一次拿上，即不考虑再以后又经过时再带些货物
（2）、要求达到不超过的时间不包括此次在该点交易的时间。
（3）、所用的距离数据都精确到米而时间则精确到
（4）、同一地点有多件货物也简单按照每件3分钟交接计算。
、符号说明
其中i，j=1、2、3……50 并且 M=50kg V=1m3
4、模型的建立及求解
模型四
模型三
模型二
模型一
最短路径模型
图的遍历模型
多区域最短路
多阶段最短路
任意两点之间的最短路距离
由起始点遍历路径回到原点
多区域无返回起点的最短路
多阶段有返回起点的最短路
模型一
模型的建立
我们为了求出各个点的之间的最短的路径，使用Dijstra算法求解。
Dijkstra算法是图论中非常有名的一个算法。
图采用邻接矩阵的形式描述，w（i,j）表示结点i到结点j间的最短距离，如果没有直接连通，则为无穷大，计算机中可以用一个很大的数据代替（如matlab中的inf）。
但dijkstra算法只能求出从结点i到其它各结点的最短路径。
算法引入这样两个集合s和t，s是那些已经确定了到i结点的最短路径的结点，t为全集u和s的差集，即那些还未确定最短路径的结点。
而且s的初值是{i}，t的初值是u-{i}。
另外再引入一个标记数组d[n]，其中在某一步d[k]表示当前从i到k的较短路径，d[k]的初值为w（i，k）。
整个算法过程如下：
１、在t中选择一个d[k]最小的结点k，将k并入s，并从t中去掉，如果t为｛｝则转到３；
２、用k结点和t中其余结点进行一遍比较，如果d[i]>d[k]+m[k][i]，则用d[k]+m[k][i]取代原来的d[i]，重复１；
３、算法结束，此时d[k]中保存的就是从i到k结点的最短路径。
算法就以这样非常简单的形式完成了求解，时间复杂度是O（n^2），确定了从i到其余各结点的最短路径。
模型的求解
根据算法和相邻的点的距离可以用dijkstra求出任意两点的最短路径。
图1相邻的点的距离
使用循环的结构求出1-50各个点之间的最短距离。
程序1见附录
可以求出w和a
a为最短路径是的所过的的地点
如从O开始到其余50个点的a（0）=[074831511812141813131821122321024220291731190313025222623283138214036273437433841364140464240]
要从O点到16点则要先到23即0-23-16要从O点到23点则要先到17即0-17-23-16要从o点到17点则要先到21即0-21-17-23-16而O可以直接到21所以从0到16的最优路径是0-21-17-23-16最短的距离是w（0,16）=7493m
模型二
—对于问题一的求解
模型的建立
由前30件货物可以到达的地点可以知道i，j=13、14、16、17、18、21、23、24、26、27、31、32、34、36、38、39、40、42、43、45、49。
图2需要达到的点（红点标注的）
其中共经过21个点，运送30件货物
该30件货物m=<50kgv=m3,所以可以一次把货物携带进行运送。
由T与W关系可知要使所用的时间最小即所走的距离最短。
即
目标函数是：
T=W÷V+t0×30
约束条件是：
必须全部遍历回到0点
即求出从O出发遍历这图的21个点的并回到o的最短的距离
要距离最短则每一步也要最短，即从O开始找最短的点到达后继续找未遍历的最短的点则可求出最短的距离。
本题要求出回到O点则可以看到两个开始最短遍历的点在某点重合即可完成最短的遍历。
模型的求解
由图可以明显得出距离O最近的点是21点和26点。
由于32点到38点的距离小于32点到16点的距离为使从21点出来的线遍历右下的点完后再和26点出来的汇合则安排32点到35点断开。
有程序2（附录）可得：
0
1
31
12
13
2
32
13
14
3
34
14
16
4
36
15
17
5
38
16
18
6
39
17
21
7
40
18
23
8
42
19
24
9
43
20
26
10
45
21
27
11
49
22
遍历节点路线是：
0-2-4-40-45-49-42-43-38-36-39-27-31-26-0
最优的路线是：
0-2-3-0-45-42-49-42-43-38-36-27-39-27-31-26-0
总路程是：
W=53787m
最优时间是：
T=
模型三
—对于问题二的求解
模型的建立
由第一个模型建立的可以求出到达24时所用的时间是:
可知到24点的时间是：
t（24）=
由表可知必须在9点之前把货物送到24点即t（24）<1故模型一不适用于问题二的求解.
由下图3可知：
图3.考虑时间的点的位置
由于右边的点的地点需要的时间要比左边的早，所以先分两个阶段，即先走左边后走右边即先走圈内的元素由程序3（附录）可得：
从O出发经过13、18、24、26、27、31、34、39、40到达45
5
26
7
31
6
27
9
29
4
24
3
28
2
13
8
31
10
40
11
45
而到13点时必须在9点之前到达但>1，到45点时必须在9点半之前到达而>故分成两个阶段不成功，所以分四个阶段，求出各个阶段的最短距离和到达时的时间即可。
目标函数：
ti=Wi÷v+t0
约束条件是：
T到个点的时间最大值
模型的求解
①
④
③
②
图个阶段的圈图
对四个阶段分别求出到达的时间，由程序4（附录）可知
l分4个阶段
3
18
2
13
4
24
1.从0出发经过13、18到24。
满足t<1的条件
故路线为：
0-18-13-24
2
31
3
34
4
40
5
45
2.从24出发经过31、34、40到45。
满足t<
故路线为：
24-31-34-40-45
3
42
5
49
4
43
2
38
3.从45出发经过38、42、43到49。
满足t<
所以路线为：
45-42-49-43-38
4.从38出发经过14、16、17、21、23、26、27、32、36、39回到O。
10
36
8
27
11
39
7
26
5
21
4
17
6
23
9
32
3
16
2
14
满足t<4
故路线为：
38-36-27-39-27-3-32-16-14-21-0
所以总的遍历点顺序是：
0-4-40-45-42-49-43-38-36-27-39-26-2-14-0
总时间是T=
总距离是W=57912m
最优路线是：
0--49-42-43-38-36-27-39-27-3-32-23-16-14-21-0
到每个点的时间见附录
模型四
—对于问题三的求解
．模型的建立
本题中要遍历所有的50个点但由于M总=147kg，v总=m3而M<50kg,V<1m3故应该以M<50kg和V<1m3判断的标准到达的最远的点后返回。
目标函数：
W=150w（i,j）
约束条件：
M<50kg,V<1m3
模型的求解
由O开始逐渐依次找出最近的点后再找出离该点最近的点直到不满足约束条件。
见程序5（附录）
图5.改进后的遍历图
1第一阶段
2.第二阶段
3.第三阶段
4.第四阶段
模型的优化
由于总的m=148kgv=m3
所以最少要分四个阶段，但由于每次不可能刚好带满50kg而如果只要3次则最多只能带150kg只比原货物多2kg所以不可能是三次就把货物带完，最少要四次。
故只需要把上述的模型进行数据处理就好了。
过程如下：
1.由于到21点时M=49V=若走过14则M大于了50故直接从21点返回。
最优路线为：
0-26-3-38-35-32-23-17-21-0
走的距离W=27122m，花费的时间T=
2.若按程序给出的从13到8的路线是而当为13-11-12-8时更短故修改之；同时到达40后如果选择34则45的周围全被遍历过。
到45后M=,V=不满足要求，故从40到34后沿21-26返回。
最优路线为:
0--3-5-38-43-42-49-50-40-45-36-21-0
走得距离是：
W=83220，所用的时间是T=
3.当到达45点时若要去20点放货物的话则需要遍历许多已经遍历过的地点，故从45点沿36-21-0返回
最优路线为：
0-26-3-22-30-28-33-46-48-44-4-45-36-21-0
所走的距离为：
W=128970m，所用的时间是：
T=
4.只余下了5个点,所以由图可知
路线为：
0-26-3-22-20-2-5-2-4-3-8-12-13-18-o
总路程是：
W=171510m所用的时间是T=
由上面的四个阶段可以知道该问的最优路线是：
0-26-3-38-35-32-23---6-23-32-35-38-43-42-49-50-40-45-36-20-28-33-46-48-44-4-45-36-20-2-5-2-4-3-8-12-13-18-o
总路程是：
W=171510m
所用的时间是T=
5、模型的分析
①误差分析：
对于模型一是使用了精确地Dijkstra算法，故误差可以忽略不计
对于模型二假定了32到38点的断开存在一定的误差，但相对于断开其余的几点得到的数值要小，故该模型可以使用。
对于模型三，由于分区域的方法有很多，故不可避免的存在些许误差，但由于区域越多，路程越多，故选择分成4个区域最合适；分成的四个不同的时间的到达区域比较紧密故按照时间的不同划分了四个区域，从而大大的消除了误差，此模型可以使用。
对于模型四的误差比较大，由于未考虑货物的拆分可能会有一定的影响同时由于4个阶段的划分也是有一定的不确定性故误差存在。
对于该模型简化了考虑的条件，仅以M和V为判断标准，虽对准确性存在挑战，但该模型相对与其他的分类有明显的优越性。
故该模型适用于该问的求解。
②灵敏度分析
对于模型一、二、三，灵敏度很好，模型的准确性很高。
对于模型四由于质量和体积的制约，使其灵敏度不会很好，但准确性较高，因此模型可以使用。
6、模型评价、改进和推广
模型的评价
优点：
l充分利用了已知数据建立模型，使其具有很高的准确性和可行性
l使用了准确的算法和适当的假设，使模型的准确性和实用性到达统一
l运用功能强大的Matlab工具使数据处理误差达到最小
缺点
l由于数据较多，没法使用工具进行模型的验证，只能一步一步的精化模型
模型的改进
对于模型一和三主要是进行验证。
对于模型二断开的那个点可以去取别的点进行。
主要是模型四的改进，可以考虑到不同的地点送的货物进行拆分，从而渠道最优的解
模型的推广
可充分使用到图的遍历和最短路的一系列问题的求解中。
7、参考文献
FirstCourseinMathenmaticalModerling（ThirdEdition）
FrankMauriceWilliam
2.图论任韩。
3.数学建模案例选集姜启源谢金星
4.图论第3版德迪斯特尔 着
5.大学生数学建模竞赛辅导教材叶其效
6.基于matlab动态规划中最短路线的实现程序[J]电脑学习
施益昌郑贤斌李自立
7.物流配送问题的混沌优化算法研究中央民族大学学报（自然科学版）2009年11月第18卷第4期
8.Dijkstra算法在企业物流运输网络中的应用湖南农业大学学报（自然科学版）2005年8月第29卷4期
附录
附录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
50个位置点的坐标
位置点
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
相互到达信息
序号
位置点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
模型二中到达时的时间
点
到的时间
最大允许的时间
0
0
0
18
1
13
1
24
1
31
34
40
升级会员 