毕业设计数学建模论文供应链网络的建立与道路破坏问题文档格式.docx
《毕业设计数学建模论文供应链网络的建立与道路破坏问题文档格式.docx》由会员分享，可在线阅读，更多相关《毕业设计数学建模论文供应链网络的建立与道路破坏问题文档格式.docx（16页珍藏版）》请在冰豆网上搜索。
并非所有的道路都可以被破坏，可破坏的道路见表4。
当某条道路被破坏后，该条道路就不能再被使用，以前运输经过该道路的只有改道，但总是沿最短路运输。
如果破坏方选取的策略是使对方总费用增加25%，而每破坏一条道路都需要成本和代价，因此需要破坏最少的道路。
问破坏方选取哪几条线路进行破坏。
给出具体的破坏道路和总费用。
3.假定各道路能否被破坏具有随机性，当某条道路被破坏后，该条道路就不能再被使用，以前运输经过该道路的只有改道，但总是沿最短路运输。
由于破坏方选取一些边进行破坏时，这些边不一定被破坏，而是服从一定的概率分布。
设可破坏的边及各边破坏的概率见表4。
运输时产生的费用可按照各种情况下的平均费用来考虑。
如果破坏方选取的策略是使对方平均总费用增加最大。
给出具体的破坏道路和平均总费用。
9
三、问题的分析
由题意可知，问题一的就是为了建立一种模型，解决配送中心选址的问题，在供应点不确定的情况下，找出最优解，使得总的费用最少，此题经考虑供应点的固定建设费用和，一次运输费用，大大的降低了问题的复杂性；
问题二则仅需列出可能的道路破坏情况，计算出每一种情况下使对方增加的费用，然后进行组合，找出破坏最少道路使对方总费用增加25%的情况；
问题三则需用数学期望来计算其增加费用的平均值，然后结合概率论知识，计算出多条道路被破坏的情况下的平均增加费用。
。
四、建模过程
（1）问题一
1、模型假设：
1、在本题中，假设供应点的城市其供应量可以满足有需要的城市的需求。
2、本题在计算总费用过程中，只考虑固定建设费用和运输费用，不考虑其他费用。
3、运输费用与运输量成正比。
2、定义符号说明：
i,j——对应城市的编号；
——城市i被选中作为配送中心的固定建设费用；
——城市i、j之间的距离；
——城市j的需求量;
——城市i被选中作为配送中心的固定建设费用和为其供应的城市送货的运输费用之和。
3、模型建立：
由于本题要在49个城市中选取若干个作为配送中心，使得每一个城市都能满足其货物需求量，同时要使总体费用最小。
解决此问题需先确定出哪几个城市作为配送中心。
我们可以先假设在只能为相邻城市供应的情况下计算出
的大小，初步确定出几个满足条件的点作为配送中心，再根据各城市的基本建设费用和配送中心的方位区域、城市之间的距离精确确定出配送中心。
由公式
=
+
（1）
其中：
由于题目中给出的基本建设费为100000000的不能作为配送中心，所以i只可取1、3、4、5、7、10、11、12、13、14、15、16、17、18、19、20、22、23、24、25、26、27、28、30、40、43、44、45这二十八个值。
根据题目给出的数据由公式
（1）可求出表1中的数据：
Y1
1899905.0
Y3
1369900.0
Y4
643833.0
Y5
746290.0
Y7
792920.0
Y10
1742000.0
Y11
1984326.0
Y12
1225750.0
Y13
1698612.0
Y14
1970025.0
Y15
2359367.0
Y16
2150388.0
Y17
1601654.0
Y18
1616090.0
Y19
2112800.0
Y20
1241215.0
Y22
3556900.0
Y23
1912518.0
Y24
1834095.0
Y25
3197585.0
Y26
2056858.0
Y27
3604282.0
Y28
759004.0
Y30
719619.0
Y40
2341203.0
Y43
1343567.0
Y44
1223214.0
Y45
21204380
表1
对以上表格中的Y值排序，取其前五个值Y4、Y30、Y5、Y28、Y7结合49各城市的坐标图进行分析:
从图中和看出，4城市作为配送中心非常合适，由于4和5、30相邻，故不选5和30作为配送中心，城市28和7也非常符合作为配送中心。
又因为城市26的基本建设费用非常小，而且距离其周边城市非常远，以28号城市也可作为配送中心。
现在确定下了4、7、26、28这四个城市作为配送中心。
将上述配送中心相邻的城市都排除掉，只剩下1、2、9、10、11、12、13、14、15、17、18、19、20、21、22、23、24、25、32、33、34、35、36、37、38、43、44、45、48、49这几个城市需要考虑。
又这些城市中能作为配送中心的只有1、10、11、12、13、14、15、17、18、19、20、22、23、24、25、43、44、45这十八个点，且他们都分布在东南方区域。
不考虑已经确定作为配送中心的城市以及其周围相邻的城市，再次由公式1计算的表2中的数据：
1173051.0
22093095.0
3204730.0
9435870
1406984.0
1501143.0
2120438.0
991425.0
3138425.0
表2
再次对表2中的Y值由小到大进行排序，得到以下一组数：
44、12、43、20、1、17、23、18、13、10、24、14、11、19、45、15、25、22.。
对这些数据进行分析:
由于44号和43号城市基本固定建设费用过高，故不应选取44和43号城市作为配送中心，12号和20号城市符合要求，排除44号并且确定20号之后，使得Y45达到最小，故45号城市也应被确定为配送中心。
可在坐标图上可以看出，12号城市不足以辐射至整个东南沿海地区。
因此，需在东南方增加几个点为配送中心。
经坐标图中城市间距离分析可知，选取12号城市作为配送中心时，还需增加13号城市为配送中心可使费用更少。
在不选取12号城市为配送中心的情况下，可用11号城市代替12号和13号城市作为东南地区的配送中心。
故此处有两种情况，应分选4、7、12、13、20、23、26、28、45这九个个城市和选取4、7、11、20、23、26、28、45这八个城市作为配送中心的情行进行讨论：
一、在选取4、7、12、13、20、23、26、28、45九个城市的情况下，根据各个城市间的距离可确定各个供应点城市所供应的城市。
具体如下：
4——1、2、3、5、15、16、27、46、47、
7——6、8、39、40、41、42
12——9、10、14、38、43
13——11、32、36、37
20——19、21、25、24、33、35、34、48、49
23——22
26——无
28——29、30、31
45——17、18、44
利用matlab软件，使用公式1对选取的供应点和供应城市进行计算得：
=9618177（元）（matlab计算程序附于论文后程序1中）
二、选取4、7、11、20、23、26、28、45这八个城市作为配送中心的情况下，根据各个城市间的距离可确定各个供应点城市所供应的城市。
4——1、2、3、5、15、16、27、46、47
11——9、10、12、13、32、36、37、38、43
20——19、21、24、25、33、34、35、48、49
28——31、29、30
45——17、14、18、44
同理利用matlab软件，使用公式1对选取的供应点和供应城市进行计算得：
=9197077（元）（matlab计算程序附于程序2中）
对两种情况结果比较可知，选取八个城市作为供应点比选取九个点所用的费用更少，所以选取4、7、11、20、23、26、28、45这八个城市作为配送中心即为此题最优解。
最少费用为9197077元。
每一个供应点到需求点的连接图如下图1.（图1的画图程序附于程序3中）
（图1）
（2）、问题二
1、假设某条道路被破坏后，该条道路就不能再被使用；
2、假设运输路线被破坏后，新的运输线路总是沿着最短路径运输。
B：
城市的货物需求量；
W：
由于道路被破坏，更换线路后增加的费用；
ΔS：
更换线路后运送路程的增加量。
3、模型建立及说明：
总共有九条道路可以被破坏，但由于道路1道路2都与城市4相连，且相互之间对运输路线有影响，所以，道路1和道路2应该一起讨论。
同理，道路5和道路9也应该在一起讨论，所依，这就条道路总共应分11种情况惊醒讨论分析：
道路断开后比断开前增加的费用应使用如下公式计算
W=0.5×
B×
ΔS
（2）
（1）道路1破坏，道路2未被破坏：
4——5：
W=54635（4——5表示城市4到5的运输路线）
4——47：
变为（28—47）：
W=114359-77686=36673
W1=91308
（2）道路1未破坏，道路2被破坏；
4——1：
W=363440
4——2：
W=297070
4——3：
W=318450
4——15:
W=10220805
W2=1081168.5
（3）道路1,2同时被破坏
变为（7——1）：
W=695464-295680=399784
4——2：
W=311680
4——3：
W=318450
4——5：
变为（28——5）：
W=136030-59059=76935
4——15：
W=102208.5
变为（28——47）：
W=114359-77686=36673
W3=1245730.5
（4）道路3被破坏
7——40：
W4=99538
（5）道路4被破坏
11——10：
W=59592
11——12：
W=46508.5
11——38：
W=72675.5
11——43：
W=90534
W5=269310
（6）道路5被破坏，道路9未被破坏
20——19：
W=168990
20——33：
W=150715
20——34：
W=11825
20——35：
W=50095
W6=381625
（7）道路5未被破坏，道路9破坏
20——21：
W=56940
20——49：
W20075
W7=77015
（8）道路5和9同时被破坏：
20——19:
变为（45——19）：
W=506184-279030=227154
20——21:
变为（45_——21）：
W=145704-43680=102024
20——33:
变为（45——33）：
W=5085755-305986.5=202589
20——34:
变为（45——34）：
W=38995-23100=15895
20——35:
变为（45——35）：
W=164847.5-97510.5=67337
20——49:
变为（45——49）：
W=58795-22825=35970
W8=650969
（9）道路6被破坏，由于没有供应点与城市之间的鱼送路线经过该道路，所依，该道路破坏与否对运输费用无影响，即：
W9=0
（10）道路7被破坏：
45—17：
W=212670
45—14：
变为（11—14）：
W=158400-156420=1980
W10=212670+1980=214650
（11）道路8被破坏：
由于城市49只能通过道路8进行运输，如果道路8被破坏，则无法向49号城市运送货物，因此，这种情况不存在，即，道路8不可被破坏。
由以上11种情况可知，要使总费用增加25%，切破坏道路最少，则应该选取道路1、2、4、5、7、9或者2、3、4、5、7、9进行破坏。
1、道路1、2、4、5、7、9同时被破坏：
W总=W1+W2+W4+W5+W7+W9
=2380659.5
W总/
=2380659.5/9197077×
100%
=25.88%
2、道路2、3、4、5、7、9同时被破坏：
W总=W2+W3+W4+W5+W7+W9
=2315635.5
W总/
=2315635.5/9197077×
=25.18%
以上两种情况均满足题目要求。
（3）、问题三
1、模型假设
1、假设各条道路能否被破坏是随机的；
2、定义符号说明
Eζ（i，j，k…）——表示i，j，k…道路被破坏后费用增加的期望。
W（i，j，k…）——由于道路i，j，k…被破坏，更换线路后增加的费用；
P（i，j，k…）——i，j，k…道路被破坏的概率。
3、模型建立
由于该9条道路的破坏服从一定的该路分布，所以计算每条道路破坏时增加的费用应使用数学期望来计算其平均增加费用，由此题第二问可知：
1、仅1道路被破坏时增加的费用和此事件发生的概率为：
W1=91308，P
（1）=0。
6×
（1-0.7）=0.18
2、仅2道路被破坏时增加的费用和此事件发生的概率为：
W2=1081168.5，P
（2）=0.28
3、仅1和2道路被破坏时增加的费用和此事件发生的概率为：
W（1,2）=1245730.5，P（1,2）=0.42
4、仅3道路被破坏时增加的费用和此事件发生的概率为：
W3=99538，P（3）=0.45
5、仅4道路被破坏时增加的费用和此事件发生的概率为：
W4=269310，P（4）=0.5
6、仅5道路被破坏时增加的费用和此事件发生的概率为：
W5=381625，P（5）=0.22
7、仅6道路被破坏时增加的费用和此事件发生的概率为：
W6=0，P（5）=0.4
8、仅7道路被破坏时增加的费用和此事件发生的概率为：
W7=214650，P（7）=0.5
9、8道路被破坏：
此种情况无意义。
10、仅9道路被破坏时增加的费用和此事件发生的概率为：
W9=77015，P（9）=0.27
11、5和9道路同时被破坏时增加的费用和此事件发生的概率为：
W（5,9）=650969，P（5,9）=0.27
有上述条件可计算出以下道路破坏时的数学期望：
Eζ（1,2）=0.42*1245730.5+0.18*91308+0.28*1081168.5=942369.43
Eζ（3）=0.45*99538=44792.1
Eζ（4）=0.5*269310=134655
Eζ（5,9）=0.33*650969+0.22*381625+0.27*77015=319571.32
Eζ（6）=0
Eζ（7）=0.5*214650=107325
由公式
可知：
在一定范围内，破坏的道路越多，增加费用的期望越大，当破坏道路为1、2、3、4、5、7、9时，
Eζ（1,2,3,4,5,7,9）=1448712.85
当破坏道路为1,2,3,4,5,6,7,9时
Eζ（1,2,3,4,5,6,7,9）=1448712.85
当破坏道路为1,2,3,4,5,6,7,8,9时
Eζ（1,2,3,4,5,6,7,8,9）=0.4*Eζ（1,2,3,4,5,7,9）
=579485.14
即，破坏道路1,2,3,4,5,7,9或者1,2,3,4,5,6,7,9可以使对方平均总费用增加最大，并且增加费用为：
1448712.85元
程序1
272080+0.5*（480*1232+580*974+210*965+530*223+521*603+430*721+630*774+934*224+716*217）
ans=
1.7488e+06
>
457824+0.5*（330*715+230*989+717*583+429*317+347*204+1016*272）
1140106
370120+0.5*（440*1391+160*624+610*495+758*761+435*948）
1371644
580032+0.5*（745*677+490*246+266*174+592*568）
1.0838e+06
526680+0.5*（710*786+560*156+830*55+305*366+650*364+873*701+837*233+840*55+820*531）
1690637
684608+0.5*340*3257
1238298
196384+0.5*（230*194+500*151+234*1980）
488104
243808+0.5*（362*834+508*642+466*1150）
825780
1748800+1140106+1371644+1083800+1690637+1238298+31008+825780+488104
9618177
程序2
825780+0.5*632*495
982200
411616+0.5*（190*1391+279*624+439*487+745*636+1235*246+1011*174+153*568+877*761+604*948）
ans=
1877924
1748800+1140106+1877924+1690637+1238298+31008+
982200+488104
9197077
程序3
city_coordinate=[36392685
37122601
34882465
33262444
32382771
41962956
43123210
43863430
41771756
39181821
40611630
37801788
40291162
36761422
37152322
34292092
35071624
33941357
3439799
2935760
3140450
27691508
25451643
27781174
23701025
13041688
30072030
25622244
23812324
27882509
13323305
42631069
3538702
3470696
3526737
3928971
42011603
40162285
40892613
42962920
40953374
45122710
37512055
33341893
32291633
30542290
30892749
3044919
3053261
];
bestchrom=[47112023262845];
figure
（1）
title（'
×
î
Ó
Å
¹
æ
»
®
É
Ë
Í
Â
·
Ï
ß
'
）
cargox=city_coordinate（bestchrom,1）;
cargoy=city_coordinate（bestchrom,2）;
plot（cargox,cargoy,'
rs'
'
LineWidth'
2,...
'
MarkerEdgeColor'
r'
...
MarkerFaceColor'
b'
MarkerSize'
20）
holdon
plot（city_coordinate（:
1）,city_coordinate（:
2）,'
o'
k'
.
升级会员 