连锁店货物配送及选址优化-数学建模论文

1、.连锁店货物配送及选址优化摘 要梦想连锁是一家肉类食品加工与销售公司，降低运输成本、增大销售量对公司经营具有重要意义。在充分理解题意的基础上，通过对问题的深入分析，我们建立多个模型，对连锁店货物配送及选址进行了优化。在问题一中，要求设计运输成本最低的生产与配送方案，由于运输成本与路程有关，因而建立Floyd模型，运用MATLAB求解出两个生产基地分别到各个销售店的最短路径，再建立邻接矩阵，得到配送方案与生产方案，求解得到最低运输成本为10540.89元/天。方案见表1和表2。 在问题二中，为分析各城镇需求特征并预测全省猪肉需求量峰值时间，首先对各城镇猪肉需求量求平均值及方差，由于各城镇需求量之

2、间相差甚大，进一步求得变异系数，以变异系数为数据进行分析，有32个城镇猪肉需求量波动较大，其余城镇变化平稳。建立拟合模型，对全省总猪肉需求量进行拟合，运用MATLAB求得在2014年1月猪肉需求量达到峰值，进而对各城镇猪肉需求量进行拟合，求得各城镇在2014年1月猪肉需求量，排名前五的城镇为120、31、63、106、150，排名后五的城镇为102、84、74、30、143。在问题三中，要求设计使全省销售量达到最大的增设销售连锁店方案，将问题中的要求视为约束条件，建立非线性规划模型，以总销售量为目标函数，运用Lingo求解出全省新增设24家销售连锁店，全省销售量的最大值919424kg/日。增

3、设销售连锁店和原有销售连锁店销售量见表7和表8。在问题四中，要建立新的生产基地使运输成本最低，由于仍要求解最低运输成本，因而建立非线性规划模型，将具有销售店的城镇作为新增生产基地的选址地点，求解出每一个城镇若作为生产基地的运输成本，由于要求新增的每日产品生产达到250吨以上，因而对求解后的结果进行筛选，得到日产品生产在250吨以上的运输成本最低的新生产基地。新的生产基地选址在31号城镇，最低运输成本为1.5万元/天。本文建立多个模型，对梦想连锁公司的货物配送进行了优化，并且对模型的优缺点进行了合理的评价。关键词： Floyd；拟合；非线性规划；MATLAB.1. 问题重述梦想连锁是一家肉类食品

4、加工与销售公司，主营：鲜猪肉。公司在全省县级及以上城镇设立销售连锁店。通过建立数学模型求解以下问题：（1）目前公司现有2个生产基地、23家销售连锁店，生产基地设在120号和63号城镇，为23家连锁店提供鲜猪肉，连锁店的日销售量见附录1。若运输成本为0.45元/吨公里，设计生产与配送方案，使运输成本最低。（2）分析各城镇需求特征，并预测未来数年，何时全省鲜猪肉需求达到峰值，达到峰值时需求达到前5位和后5位的城镇是那些。（3）未来几年公司在全省的市场占有率可增至3成左右。公司产品的需求量与销售量并不完全一致，若在当地（同一城镇）购买，则这一部分需求量与销售量相同，若在不足10公里的其他城镇的销售连

5、锁店购买，则这一部分需求量只能实现一半，而在超过10公里的其他城镇的销售连锁店购买，销售量只能达到需求量的三成。公司决定在各城镇增设销售连锁店，基于现有条件、成本等的考虑，原有的23家销售连锁店销售能力可在现有销售量的基础上上浮20%，增设的销售连锁店销售能力控制在每日20吨至40吨内，并且要求增设的销售连锁店的销售量必须达到销售能力的下限。同一城镇可设立多个销售连锁店。设计增设销售连锁店方案，使全省销售量达到最大。（4）在增设销售连锁店的基础上，公司决定增加生产基地，地址设立在城镇所在地，每日产品生产必须达到250吨以上，在生产与销售各环节不能有产品积压。设计生产基地增设方案，使运输成本最低

6、。（5）公司产品若采用载重1.5吨的小货车从生产基地运往销售连锁店，小货车在高速公路上限速100公里/小时，在普通公路上限速60公里/小时，销售连锁店需要的产品必须当日送达。假设：每日车辆使用时间不超过8小时，小货车装满或卸完1.5吨的货物均需要半小时，本市运输车辆行驶时间可忽略不计。在公司增设销售连锁店、增加生产基地后，为完成每日运输任务，请你为公司确定小货车的最小需求量，及各车辆的调运方案。2. 模型假设1) 假设每个销售店只有一辆配送车配送且载货量为无限大；2) 假设生产基地到连锁店的距离即为生产基地所在城镇到连锁店所在城镇的距离；3) 每个连锁店只有一个生产基地来配送；4) 设生产基地

7、生产量与配送销售量相等。3. 通用符号说明序号符号符号说明1各销售店与生产基地最短距离2城镇序号3各销售店日销售量4各城镇猪肉需求量平均值5各城镇猪肉需求量标准差6变异系数4. 问题一：模型的建立、求解4.1问题分析 本问题要求设计运输成本最低的生产与配送方案，首先拟根据附件中的全省交通网络数据画出全省交通网络拓扑图。由于问题本质上仍是最短路问题，拟建立Floyd模型。首先打算求解出所有城镇之间的最短路径，再建立0-1矩阵，选取距离销售连锁店最近的生产基地，即得到各个生产基地的配送方案，将距离加和，之后求解出运输成本即为最低运输成本。4.2模型准备根据附件全省交通网络数据，编写MATLAB程序

8、求解画出全省交通网络拓扑图，如图1。程序见附录1。图1 全省交通网络拓扑图 图1即为全省交通网络图。4.3模型建立建立Floyd模型对问题一进行求解。建模过程如下：Step1 求解最短路径全省交通网络共有154个位置点，现需要求从生产基地到销售连锁店的的最小运输成本。设邻接矩阵 （4-1）来存放各边长度，其中：； ，i，j之间不相连，在程序中以各边都不可能达到的充分无穷大的数代替。 是i，j之间的欧氏距离，i，j=1,2，n对于无向图是对称矩阵， 递推产生一个矩阵序列，其中，表示城镇到城镇路径上所经过的顶点序号不大于k的最短路径的长度。计算时用迭代公式： （4-2）k为迭代次数，i,j,k=1

9、,2,，n。最后，当k=n时，即是各城镇之间最短路径值。Step2 选取最近生产基地将具有销售店的城镇挑选出来，一共23个城镇。引入0-1变量令 （4-3）则运输成本为 （4-4）其中为各销售店日销售量，为各销售店与生产基地最短距离。4.4模型求解通过编写MATLAB程序得到各生产基地配送城镇分配方案。程序见附录2。表1 63号城镇生产基地分配方案城镇编号连锁店编号距离（km）日销售量（公斤）费用（元）106108.368481413.5505221116179.156103492.00860252220168.956375484.67531252417128.943251188.632773

10、2712135.19265563.265675314114.66239471235.593359633021733018028295064237.3118406.0526865719.0915570133.75408579828.1738759491.3284635表2 120号城镇生产基地分配方案城镇编号连锁店编号距离（km）日销售量（公斤）费用（元）19134.3114744891.1199881621103.6414783689.4495543413119.5445124.2606433611151.1911503782.61235654214110.589489472.18212994

11、15170.1712773978.1116345106263.7382231095.66229512010287330100325170123225.111808141.5772595141561.729258257.1316921451972.85396531299.924473根据表1和表2可得（1）63号城镇生产基地的货物配送至11个销售连锁店，配送成本为4008.86元/天，生产量为150084公斤/天。（2）120号城镇生产基地的货物配送至12个销售基地，配送成本为6532.03元/天，生产量为230208公斤/天。（3）两个生产基地总运输成本为10540.89元/天，即为最低运输成

12、本。具体路线如下表3 配送路线生产基地路线连锁销售店631011222427316363646579120116343642941061201201231411455. 问题二：模型的建立、求解5.1问题分析本问题要求分析各城镇需求特征并预测全省鲜猪肉何时达到峰值，以及达到峰值时排序前5位和后5位的城镇，在分析特征时，拟求解各城镇猪肉需求量的平均值、标准差，由于标准差不能完全反映需求特征，因而再求解变异系数，分析这三个数据体现各城镇需求量特征。拟建立拟合模型，对全省各月总需求量进行拟合，求解出何时达到峰值，再对每一个城镇需求量进行拟合，求解出在全省需求量达到峰值的月份各城镇猪肉需求量，之后对各

13、城镇需求量排序，得到排名前五和后五的城镇。5.2模型建立1.分析各城镇需求特征通过均值和方差以及变异系数的计算，可以看出每个镇对猪肉的需求量及其波动情况。 （5-1） （5-2） （5-3）2.拟合模型对同一月中每个城镇的猪肉需求量进行累加 （5-4）其中代表第个城镇；代表第月多项式拟合 设是一个定义在实轴上的的函数，是在点上的函数值。又设是某个空间函数。所谓曲线拟合就是要在函数空间中找到一个函数，使得 （5-5）其中是权重系数，特殊地，可取。就称为在点上对的拟合曲线。通常，我们要求函数空间是一个有限维的线性空间，即存在着一组线性无关的函数，使得 （5-6）于是，求拟合曲线就可以化简为求一组实

14、数a1,a2,，aN，使得下面目标函数达到最小值 （5-7）通过偏导数方法：，我们可以得到如下线性方程组 （5-8）其中于是，曲线拟合的系数可通过求解正规方程得到。只考虑多项式拟合，也即，i=2,3，N；并令。于是，我们有，（5-9）这时拟合曲线，也称为拟合多项式。预测出拟合多项式的峰值，就是全省鲜猪肉需求量的峰值。5.3模型求解（1）分析各城镇需求特征利用Excel求解各城镇需求量的平均值、标准差与变异系数，根据变异系数排序，选取变异系数大于10的结果，如表4。表4 各城镇需求量数值分析城镇平均值标准差变异系数城镇1501661.204833521.915921631.41791495城镇9

15、4107.889333333.2156979930.78682291城镇51147.789543.9006113629.70482434城镇40396.6691667112.730164228.41919002城镇681433.024167405.2057528.27626773城镇86866.582231.849296326.75445558城镇41456.5551667115.920003125.39014156城镇124922.238234.118759925.38593724城镇138577.2691667129.638570422.4572137城镇67187.037333341.3

16、140311522.08865493城镇47122.825.8765894221.07214122城镇1101445.038167283.434544519.6143293城镇2120.989522.8069130918.85032428城镇83659.1536667119.399144818.11400753城镇55176.738166731.8659530718.03003486城镇107967.2781667170.070706917.58239902城镇108406.907833370.734202117.38334736城镇61405.785666768.7538495516.943

17、39037城镇761478.902333241.315827316.31722541城镇109101.615166716.5450565116.28207388城镇11084.089333174.314077416.07930934城镇39358.74156.5170588315.75427922城镇1042141.91336.599305115.71491356城镇1012097.4935310.577001714.80705431城镇12999.274514.2830732514.38745423城镇1208497.0478331180.82756613.8969156城镇62185.47

18、2524.7900648513.36589783城镇10688.056833390.9460489413.21781058城镇69810.14599.2830762112.25497611城镇130290.785333332.8123266211.28403769城镇141834.427166788.5797937810.61564116城镇131147.789833315.5594801210.52811264分析表4，变异系数反映各城镇每个月份猪肉需求量的离散程序，这32个城镇变异系数大于10,因而这32个城镇猪肉需求量变化较大，波动较大；而其他城镇变异系数较小，各月份猪肉需求量变化平稳。

19、（2）猪肉需求量预测将每个月各个城镇需求数加和，对加和后的数据运用MATLAB进行拟合，程序见附录3。拟合后的公式为 得到的拟合趋势图如图2。图2 全省猪肉需求量变化趋势根据图2得，拟合趋势线为一元二次方程，在第60-90个月之间达到峰值，因而对之后30个月的猪肉需求量进行预测。 根据拟合后的公式预测出2012年12月之后30个月的月度总需求量，如表5。表5 预测未来30个月月度总需求量月份需求量/kg月份需求量/kg月份需求量/kg31119490.066541119826.34851119648.117332119546.847742119831.67852119601.996331195

20、98.483843119831.862853119550.729634119644.974744119826.902654119494.318135119686.320645119816.797255119432.761536119722.521346119801.546756119366.059837119753.576947119781212938119779.487348119755.610358119217.220939119800.252749119724.924459119135.083840119815.872940119689.093460119047

21、.8016由表5可得，第43个月即2014年1月全省鲜猪肉需求达到峰值，峰值为119831.8628kg/月。要求解在2014年1月各城镇猪肉需求量，采用同样的拟合方法对每一个城镇进行拟合，程序见附录4。得到每一个城镇猪肉需求量的拟合方程，见附录5。对应求解出2014年1月的猪肉需求量，将其进行排序，部分数据如表6。完整数据见附录6。表6 城镇猪肉需求量部分排名前五后五城镇猪肉需求量/kg城镇猪肉需求量/kg1208663.505929143118.8467228314641.77054130112.3482009633822.32188174107.0905441063560.9096438

22、498.253206191502289.52822910292.02227853根据表6可得排名第一和第二的城镇是公司生产基地所在地，符合常理。5.3模型检验5.3.1检验模型建立和方差SSE 该统计参数计算的是拟合数据和原始数据对应点的误差的平方和，计算公式如下： （5-10）SSR预测数据与原始数据均值之差的平方和，公式如下： （5-11）SST原始数据和均值之差的平方和，公式如下： （5-12）确定系数 （5-13）5.3.2模型检验求解得到，由于越接近1，说明模型对数据拟合的越好，因此本次建立的模型较好。6. 问题三：模型的建立、求解6.1问题分析本问题要求根据题目中的猪肉需求量变化、

23、销售连锁店的销售能力等条件设计增设连锁店方案使全省销售量达到最大，由于条件较多，拟建立非线性规划模型，将问题中的距离、成本、需求等条件作为约束条件，将全省销售量作为目标函数，进行非线性规划求解。6.2模型建立为公司设计增设销售连锁方案，使全省销售能力达到最大。根据题中已知条件，设表示新增的所有连锁店的销售能力；表示未来号城镇的需求；表示在不足10公里的其他城镇的销售连锁店购买，则这一部分需求量只能实现一半；表示超过10公里的其他城镇的销售连锁店购买，销售量只能达到需求量的三成；表示表示原有的23家店在新增设连锁店的销售能力。设增设的连锁店 （6-1）目标函数 （6-2）约束条件 （6-3）6.

24、3模型求解编写Lingo程序对模型进行求解，得到的增设销售店序号及各销售量如表6所示。程序见附录7，具体结果见附录7。表7 增设销售店及销售量城镇销售量/（kg/日）城镇销售量/（kg/日）6234537621909.5827491.810021208.51025228.810127736.11820673104204783140000110223893320304.6116203825021972.8120400005422079.7121400005625328.2123400006223053.512529010.76423658.115027089.3682512215422002.2

25、根据表7可得，一共需要增加24个销售连锁店，增加的销售连锁店的销售额为630570.8kg/日。得到的原销售店增长之后的销售量如表8所示。表8原有销售店及销售量城镇销售量/（kg/日）城镇销售量/（kg/日）城镇销售量/（kg/日）117692.83115474.17919577.51010177.234541.2945134.6117323.6361145910642360161761.1425330.212054147.422765063464891235370243901.264220814110270.82711118657035.51453832根据表8可得，原有销售店增长销售额之后

26、的总量为288853.2kg/日。根据上述表7和表8的销售连锁店方案，全省销售量达到最大，最大为919424kg/日。7. 问题四：模型的建立、求解7.1问题分析本问题要求建立新的生产基地且运输成本最低，拟建立非线性规划模型。将具有销售店的城镇作为新增生产基地的选址地点，求解出每一个城镇若作为生产基地的运输成本，由于要求新增的每日产品生产达到250吨以上，因而对求解后的结果进行筛选，得到日产品生产在250吨以上的运输成本最低的新生产基地。7.2模型建立总运输成本=单位运输成本*生产基地与连锁店的距离*运输重量，根据题意建立模型如下将总运输成本作为目标函数 （7-1） （7-2）其中，表示第i个

27、连锁店到距离它最近的生产基地的距离，表示第i个连锁店的销售量，分别表示为新建的，63镇的和120镇的生产基地的生产总量，表示各城镇的销售总量。7.4模型求解编写MATLAB程序对模型进行求解，得到的结果为在31号城镇建立新的生产基地，最低运输费用为1.5万元/天。因此，在31、63、120号城镇为生产基地，向45个城镇的连锁店配送货物，运输费用为1.5万元/天。8. 模型评价本文建立多个模型，对梦想连锁公司货物配送问题进行了求解，这些模型在求解过程中既有优点也有缺点。优点：（1）Floyd算法容易理解，可以算出任意两个节点之间的最短距离，求解最短距离是灵活可用，适合本文问题。（2）拟合模型较易

28、求解，较好的拟合猪肉需求量并预测未来最大值。（3）非线性规划求解简单。缺点：拟合模型没有对曲线进行三次及以上方程拟合，可能二次方程不是最优解。参考文献1 郝自军,何尚录，最短路问题的Floyd算法的若干讨论，重庆工学院学报(自然科学)，2008年22（5）：156-159。2 周先曙，最短路径问题及其解法研究，电脑知识与技术，2010年6（6）：1403-1405。3司守奎，数学建模算法与应用，国防工业出版社，2015年。附录附录1 网络拓扑图MATLAB程序n=154;a=zeros(n); x=137.5167376210263284278.529529936241039527723520

29、0167225400414424411420403376380377374370368374365356358354357359347261270313324333334323312317316315.5316318291.5284281299302305281274273.5267270275276270276288293.5296303298.5293296276250247262269276282282273276138.5155163173183238213189210201150135143142140121109138.5167168184194200212220246246252

30、257255.5249247254251.5259261255240241236235232235.5245225225219219228.5213211208.5206206200200200200206194184184193.5193184184192.5192192170168167;y=462399400390445409425382444443408.5520496465466.5462443447422400396401404406404424424423427.5431448450459495513528553537.5514511.5511511497497498451448

31、444434412415425421434451457458.5448444446440422419415405403392.5387386378376372352350384399397402398.5386.5389361378396390364370382412413433434400395407414430432441442442435440442442443443444455458460.5466464469472477478470494495514514496487486.5474457.54514514624724814874965075155145074974844664664

32、63.5475475484484496.5496.5507514516.5507495.5;scatter(x,y,g+);title(网络拓扑图);hold on;a(1,100)=21.14;a(1,16)=30.67;a(2,85)=9.91;a(2,90)=28.11;a(3,24)=6.2;a(5,60)=4.37;a(5,62)=27.87;a(6,63)=12.94;a(6,51)=9.73;a(6,66)=7.61;a(7,54)=23.47;a(8,68)=5.55;a(8,67)=11.45;a(8,71)=6.79;a(9,48)=16.58;a(9,55)=7.67;a

33、(10,32)=9.3;a(10,31)=6.3;a(11,19)=15.41;a(11,22)=13.35;a(13,45)=35.52;a(13,119)=23.69;a(14,126)=13.6;a(14,109)=15.04;a(15,140)=18.62;a(15,105)=25.31;a(15,141)=6.41;a(15,142)=6.47;a(16,143)=18.06;a(17,108)=22.58;a(18,19)=29.33;a(18,31)=37.59;a(19,20)=26.51;a(20,21)=14.78;a(21,3)=35.53;a(22,20)=4.27;a

34、(22,21)=10.68;a(23,11)=8.89;a(23,21)=11.99;a(24,25)=4.83;a(24,27)=18.78;a(25,26)=20.64;a(26,27)=3.03;a(27,28)=4.44;a(29,28)=5.03;a(29,10)=17.27;a(30,29)=7.33;a(30,23)=40.53;a(30,31)=20.47;a(31,32)=9.67;a(31,33)=13.24;a(33,12)=76.92;a(33,34)=36.59;a(34,35)=18.78;a(35,36)=15.49;a(36,12)=38.84;a(36,37)

35、=27.98;a(38,37)=90.93;a(38,39)=25.43;a(39,40)=43.56;a(39,13)=20.83;a(40,41)=11.33;a(40,45)=14.35;a(41,42)=9.87;a(42,35)=25.12;a(42,43)=15.01;a(43,34)=21.62;a(44,41)=14.64;a(44,43)=11.72;a(45,44)=11.17;a(45,46)=51.68;a(46,47)=3.22;a(47,48)=4.14;a(48,49)=10.81;a(50,49)=23.17;a(50,10)=57.96;a(51,50)=27

36、.73;a(51,52)=12.84;a(53,52)=5.02;a(53,7)=5.03;a(54,49)=17.73;a(54,9)=10.45;a(55,47)=15.19;a(56,46)=13.5;a(56,55)=6.92;a(56,57)=25.9;a(57,55)=23.85;a(57,13)=38.19;a(57,58)=12.78;a(58,59)=4.07;a(59,9)=25.52;a(60,58)=7.59;a(60,61)=7.15;a(61,59)=5.7;a(61,7)=18.16;a(62,7)=4.66;a(62,63)=3.36;a(63,53)=5.45

37、;a(64,63)=7.31;a(64,65)=11.78;a(65,6)=9.07;a(66,67)=12.05;a(67,68)=6.16;a(68,69)=7.3;a(69,3)=76.68;a(70,69)=9.38;a(70,8)=5.45;a(71,72)=5.45;a(72,70)=6.96;a(72,73)=29.86;a(73,74)=26.56;a(74,75)=34.86;a(75,76)=21.38;a(76,77)=7.95;a(77,78)=9.21;a(78,65)=3.17;a(78,79)=7.16;a(79,66)=7.62;a(80,79)=12.75;a

38、(80,67)=14.25;a(80,71)=15.47;a(81,76)=15.25;a(81,77)=9.3;a(81,80)=9.41;a(81,82)=30.09;a(82,72)=23.74;a(82,86)=113.17;a(74,86)=81.41;a(83,86)=39.54;a(83,84)=24.8;a(84,2)=12.84;a(85,86)=28.31;a(85,87)=30.43;a(86,87)=12.68;a(87,88)=58.27;a(88,75)=9.85;a(88,89)=40.2;a(89,4)=23.38;a(90,91)=31.41;a(91,92)

39、=9.6;a(93,84)=6.62;a(93,95)=10.2;a(94,93)=16.53;a(94,83)=18.09;a(94,95)=14.94;a(96,95)=7.47;a(96,97)=17.32;a(97,98)=19.92;a(97,102)=29.66;a(98,96)=28;a(98,99)=15.04;a(101,97)=12.44;a(101,16)=20.64;a(101,100)=30.36;a(101,102)=7.75;a(102,90)=33.26;a(102,103)=17.53;a(103,104)=10.44;a(104,105)=6.46;a(10

40、5,92)=8.67;a(105,106)=12.93;a(106,91)=10.96;a(106,107)=8.08;a(106,130)=21.63;a(107,17)=5.23;a(108,5)=17.39;a(108,76)=48.23;a(109,108)=11.91;a(109,110)=6.83;a(110,111)=5.68;a(111,5)=17.73;a(111,112)=6.21;a(112,113)=7.15;a(112,118)=7.28;a(113,110)=6.81;a(114,113)=5.9;a(114,115)=8.03;a(115,112)=6.61;a(

41、115,116)=5.61;a(116,117)=8.18;a(117,118)=8.86;a(119,117)=16.69;a(120,119)=15.82;a(120,116)=22.06;a(121,120)=20.07;a(121,39)=30.16;a(121,38)=32.14;a(121,122)=5.09;a(122,123)=18.49;a(123,120)=5.11;a(124,123)=10.36;a(124,125)=3.77;a(125,120)=10.51;a(125,126)=15.96;a(126,116)=7.82;a(126,114)=5.81;a(127,

42、14)=13.22;a(127,128)=6.79;a(128,17)=8.21;a(128,109)=22.99;a(129,127)=9.05;a(129,128)=6.04;a(129,107)=8.68;a(130,127)=8;a(130,131)=14.78;a(131,14)=10.17;a(131,125)=16.78;a(132,131)=18.62;a(133,124)=22.71;a(133,132)=6.53;a(134,133)=10.1;a(134,123)=28.59;a(135,134)=12.24;a(135,136)=8.4;a(136,122)=31.92

43、;a(137,136)=6.66;a(137,138)=7.31;a(138,135)=6.04;a(138,139)=10.87;a(139,134)=9.1;a(140,139)=13.46;a(141,132)=18.2;a(142,104)=24.54;a(143,142)=10.98;a(143,103)=24.92;a(143,144)=11.95;a(144,145)=9.64;a(144,147)=9.02;a(145,142)=9.39;a(145,146)=9.18;a(146,140)=7.51;a(147,146)=9.33;a(147,148)=13.55;a(148

44、,149)=9.12;a(149,146)=13.22;a(149,139)=7.65;a(150,138)=8.77;a(150,149)=10.79;a(150,151)=7.65;a(151,137)=8.18;a(152,99)=100.69;a(152,38)=94.21;a(152,151)=23.56;a(152,153)=9.88;a(153,150)=24.11;a(153,154)=12.38;a(154,16)=34.66;a(154,148)=18.15;for i=1:nfor j=1:nif a(i,j)=0plot(x(i);x(j),y(i);y(j),b-);

45、 text(x(i)+x(j)/2,(y(i)+y(j)/2,num2str(a(i,j),Fontsize,6); %显示边的权值hold on; %带箭头的连线endtext(x(i),y(i),num2str(i),Fontsize,6,color,r); %显示点的序号hold on;endhold on;end附录2 Floyd算法MATLAB程序clear;clc;n=154; a=zeros(n);a(1,100)=21.14;a(1,16)=30.67;a(2,85)=9.91;a(2,90)=28.11;a(3,24)=6.2;a(5,60)=4.37;a(5,62)=27.

46、87;a(6,63)=12.94;a(6,51)=9.73;a(6,66)=7.61;a(7,54)=23.47;a(8,68)=5.55;a(8,67)=11.45;a(8,71)=6.79;a(9,48)=16.58;a(9,55)=7.67;a(10,32)=9.3;a(10,31)=6.3;a(11,19)=15.41;a(11,22)=13.35;a(13,45)=35.52;a(13,119)=23.69;a(14,126)=13.6;a(14,109)=15.04;a(15,140)=18.62;a(15,105)=25.31;a(15,141)=6.41;a(15,142)=6

47、.47;a(16,143)=18.06;a(17,108)=22.58;a(18,19)=29.33;a(18,31)=37.59;a(19,20)=26.51;a(20,21)=14.78;a(21,3)=35.53;a(22,20)=4.27;a(22,21)=10.68;a(23,11)=8.89;a(23,21)=11.99;a(24,25)=4.83;a(24,27)=18.78;a(25,26)=20.64;a(26,27)=3.03;a(27,28)=4.44;a(29,28)=5.03;a(29,10)=17.27;a(30,29)=7.33;a(30,23)=40.53;a(

48、30,31)=20.47;a(31,32)=9.67;a(31,33)=13.24;a(33,12)=76.92;a(33,34)=36.59;a(34,35)=18.78;a(35,36)=15.49;a(36,12)=38.84;a(36,37)=27.98;a(38,37)=90.93;a(38,39)=25.43;a(39,40)=43.56;a(39,13)=20.83;a(40,41)=11.33;a(40,45)=14.35;a(41,42)=9.87;a(42,35)=25.12;a(42,43)=15.01;a(43,34)=21.62;a(44,41)=14.64;a(44

49、,43)=11.72;a(45,44)=11.17;a(45,46)=51.68;a(46,47)=3.22;a(47,48)=4.14;a(48,49)=10.81;a(50,49)=23.17;a(50,10)=57.96;a(51,50)=27.73;a(51,52)=12.84;a(53,52)=5.02;a(53,7)=5.03;a(54,49)=17.73;a(54,9)=10.45;a(55,47)=15.19;a(56,46)=13.5;a(56,55)=6.92;a(56,57)=25.9;a(57,55)=23.85;a(57,13)=38.19;a(57,58)=12.7

50、8;a(58,59)=4.07;a(59,9)=25.52;a(60,58)=7.59;a(60,61)=7.15;a(61,59)=5.7;a(61,7)=18.16;a(62,7)=4.66;a(62,63)=3.36;a(63,53)=5.45;a(64,63)=7.31;a(64,65)=11.78;a(65,6)=9.07;a(66,67)=12.05;a(67,68)=6.16;a(68,69)=7.3;a(69,3)=76.68;a(70,69)=9.38;a(70,8)=5.45;a(71,72)=5.45;a(72,70)=6.96;a(72,73)=29.86;a(73,74)=26.56;a(74,75)=34.86;a(75,76)=21.38;a(76,77)=7.95;a(77,78)=9.21;a(78,65)=3.17;a(78,79)=7.16;a(79,66)=7.62;a(80,79)=12.75;a(80,67)=14.25;a(80,71)=15.47;a(81,76)=15.25;a(81,77)=9.3;a(8