数学建模论文-送货线路设计问题.pdf

优秀送货路线设计问题1、问题重述现今社会网络越来越普及，网购已成为一种常见的消费方式，随之物流行业 也渐渐兴盛，每个送货员需要以最快的速度及时将货物送达，而且他们往往一人 送多个地方，请设计方案使其耗时最少。现有一快递公司，库房在图1中的0点，一送货员需将货物送至城市内多处，请设计送货方案，使所用时间最少。该地形图的示意图见图1,各点连通信息见 表3,假定送货员只能沿这些连通线路行走，而不能走其它任何路线。各件货物 的相关信息见表1,50个位置点的坐标见表2。假定送货员最大载重50公斤，所带货物最大体积1立方米。送货员的平均 速度为24公里/小时。假定每件货物交接花费3分钟，为简化起见，同一地点有 多件货物也简单按照每件3分钟交接计算。现在送货员要将100件货物送到50个地点。请完成以下问题。1.若将130号货物送到指定地点并返回。设计最快完成路线与方式。给出 结果。要求标出送货线路。2.假定该送货员从早上8点上班开始送货，要将130号货物的送达时间不 能超过指定时间，请设计最快完成路线与方式。要求标出送货线路。3.若不需要考虑所有货物送达时间限制（包括前30件货物），现在要将100 件货物全部送到指定地点并返回。设计最快完成路线与方式。要求标出送货线路，给出送完所有快件的时间。由于受重量和体积限制，送货员可中途返回取货。可 不考虑中午休息时间。2、问题分析送货路线问题可以理解为：已知起点和终点的图的遍历问题的合理优化的路 线设计。图的遍历问题的指标：路程和到达的时间，货物的质量和体积，以及最大可 以负载的质量和体积。在路线的安排问题中，考虑所走的路程的最短即为最合理 的优化指标。对于问题二要考虑到所到的点的时间的要求是否满足题意即采用多次分区 域的假设模型从而找出最优的解对于问题三则要考虑到体积和质量的双重影响，每次到达后找到达到最大的 体积和质量的点然后返回，再依次分析各个步骤中可能存在的不合理因素达到模7型的进一步合理优化得到最合理的解。3、模型假设与符号说明3.1、模型的假设(1)、到同一地点的货物要一次拿上，即不考虑再以后又经过时再带些货物(2)、要求达到不超过的时间不包括此次在该点交易的时间。(3)、所用的距离数据都精确到米而时间则精确到O.OOOlh(4)、同一地点有多件货物也简单按照每件3分钟交接计算。3.2、符号说明其中 i,j=l、2、3.50 并且 M=50kg V=lm3符号符号说明b j送货点的标号w从。开始回到。的总的路程T从。开始回到。的总的时间W(i,j)从于始到j的最短路程V送货员的平均速度ti从。开始到i的所用最短的时间到达1可以芾的最多的货物的质量4到达i可以带的最大的货物的体积M一次可以携带的最多货物的总高质量V一次可以携带的最多货物的总高质量4、模型的建立及求解7馍型一1.1模型的建立我们为了求出各个点的之间的最短的路径，使用D ijstra算法求解。D ijkstra算法是图论中非常有名的一个算法。图采用邻接矩阵的形式描述，w(ij)表示结点i到结点j间的最短距离，如 果没有直接连通，则为无穷大，计算机中可以用一个很大的数据代替(如matlab 中的inf)o但dijkstra算法只能求出从结点i到其它各结点的最短路径。算法引入这样两 个集合s和t,s是那些已经确定了到i结点的最短路径的结点，t为全集u和s 的差集，即那些还未确定最短路径的结点。而且s的初值是i,t的初值是u-i。另外再引入一个标记数组dn,其中在某一步dk表示当前从i到k的较短路径,dk的初值为w(i,k)o整个算法过程如下：1、在t中选择一个dk最小的结点k,将k并入s,并从t中去掉，如果t 为 则转到3；2、用k结点和t中其余结点进行一遍比较，如果didk+mki,则用 dk+mki取代原来的di,重复1；3、算法结束，此时dk中保存的就是从i到k结点的最短路径。算法就以这样非常简单的形式完成了求解，时间复杂度是0(M2),确定了从7i到其余各结点的最短路径。1.2模型的求解根据算法和相邻的点的距离可以用dijkstra求出任意两点的最短路径。图1相邻的点的距离使用循环的结构求出1-50各个点之间的最短距离。程序1见附录2.1 可以求出w和aa为最短路径是的所过的的地点如从 0 开始到其余 50 个点的 a(O)=0 7 4 8 3 15 118 12 14 18 13 13 18 21 12 23 21 0 24 22 0 29 17 31 19 0 31 30 25 22 26 23 28 31 38 21 40 36 27 34 37 43 38 41 36 41 40 46 42 40要从0点到16点则要先到23即0-23-16要从0点到23点则要先到17即 0-17-23-16要从o点到17点则要先到21即0-21-17-23-16而0可以直接到21所 以从0到16的最优路径是0-21-17-23-16最短的距离是w(0,16)=7 493m演型二一对于问题一的求解2.1 模型的建立由前30件货物可以到达的地点可以知道i,j=13、14、16、17、18、21、23、24、26、27、31、32、34、36、38、39、40、42、43、45、49o73500图2需要达到的点（红点标注的）其中共经过21个点，运送30件货物该30件货物Em=47.3kg 50kg v=0.837 1m3,所以可以一次把货物携带进行运送。由T与W关系可知要使所用的时间最小即所走的距离最短。即目标函数是：T=W4-V+tox3O约束条件是：必须全部遍历回到。点即求出从O出发遍历这图的21个点的并回到。的最短的距离 要距离最短则每一步也要最短，即从O开始找最短的点到达后继续找未遍历的 最短的点则可求出最短的距离。本题要求出回到。点则可以看到两个开始最短遍历的点在某点重合即可完 成最短的遍历。2.2 模型的求解由图可以明显得出距离O最近的点是21点和26点 由于32点到38点的距 离小于32点到16点的距离为使从21点出来的线遍历右下的点完后再和26点出 来的汇合则安排32点到35点断开。有程序2(附录2.2)可得:013112132321314334141643615175381618639172174018238421924943202610452127114922顺序为：Columns 1 through 131 7 5 8 13 4 3 6 2 9 14 18 21Columns 14 through 2322 19 20 16 15 17 11 12 10 1总的路程为：5.37 87 e+004总的时间是：3.7411遍历节点路线是：0-21-17-23-32-16-14-18-13-24-34-40-45-49-42-43-38-36-39-27-31-26-0最优的路线是：0-21-17-23-32-23-16-14-21-18-13-19-24-31-34-40-45-42-49-42-43-38-36-27-39-27-31-26-0总路程是：W=537 87 m最优时间是：T=3.7 411h模型三对于问题二的求解3.1 模型的建立由第一个模型建立的可以求出到达24时所用的时间是:到达24时所走的距离4.051le+004到24所用的日寸间2.0880可知到24点的时间是：t(24)=2.08807由表2.1可知必须在9点之前把货物送到24点即t(24)1,到45点时必须在9点半之前到达而1.74121.5 故分成两个阶段不成功，所以分四个阶段，求出各个阶段的最短距离和到达时的时间即可。目标函数:ti=Wi+v+2 t0约束条件是：T到个点的时间最大值3.2模型的求解对四个阶段分别求出到达的时间，由程序4（附录2.4）可知分4个阶段1.从0出发经过13、18到24。满足t l的条件顺序为：1 3总的路程为：1.1010。+004总的时间是：0.60873180.09092130.27 064240.55872故路线为：0-18-13-242.从24出发经过31、34、40到45。满足t 1.51 2 3 4 5总的路程为：1.9963。+0042310.7 3293340.92974401.047 75451.2317总的时同是：1.331|故路线为：24-31-34-40-45431.7 32243.从45出发经过38、42、43到49。2381.8913满足t 2.25顺序为：1 3 5 4 2总的路程为：2.9793e+004 总的时间是：1.9413所以路线为：45-42-49-43-384.从 38 出发经过 14、16、17、21、23、26、27、32、36、39 回至U 0。满足t 410362.00548272.147 211392.32147262.55405212.7 4544172.87 146232.99539323.15003163.44202143.6007顺序为：1 10 8 11 75469321总的路程为：5.7912e+004总的时间是：I3.9130故路线为：38-36-27-39-27-31-26-21-17-23-32-16-14-21-0所以总的遍历点顺序是：0-18-13-24-31-34-40-45-42-49-43-38-36-27-39-26-21-17-23-32-16-14-0 总时间是T=3.9130h总距离是 W=57 912m最优路线是：0-18-13-19-24-31-34-40-45-42-49-42-43-38-36-27-39-27-31-26-21-17-2 3-32-23-16-14-21-0到每个点的时间见附录1.4模型四一对于问题三的求解4.1.模型的建立本题中要遍历所有的50个点但由于M曰二147 kg,vm二2.8m?而心7M 50kg,VGm，故应该以M 50kg和VVlm，判断的标准到达的最远的点后返回。目标函数：WN7w(i,j)约束条件：M50kg,Vlm34.2模型的求解由0开始逐渐依次找出最近的点后再找出离该点最近的点直到不满足约 束条件。见程序5(附录2.5)图5.改进后的遍历图1第一阶段顺序为:0 26 31 27总的路程为：2.7122e+004总的时间是：1.73012.第二阶段咂序为：Columns 1 through 13 0 18 13 12Columns 14,tHroigR 21 16 43 42 49息的路程为：8.4850e-F004总的时同是：4.535439 36 38 35 32 23 17 2111 8 3 1 6 7 10 950 40 34 03.第三阶段顺序为：Columns 1 through 130 24 19 25 29 22 30 28 33 46 48 44 41Columns 14 through 19 37 47 45 20 15 0总的路程为：1.5936。+005 总的时间是：7.53994.第四阶段顺序为：0 4 2 5 0息的路程为；1.9347e+005 总的时间是：8.21124.3模型的优化由于总的m=148kg Sv=2.8m3所以最少要分四个阶段，但由于每次不可能刚好带满50kg而如果只要3次则 最多只能带150kg只比原货物多2kg所以不可能是三次就把货物带完，最少要四 次。故只需要把上述的模型进行数据处理就好了。过程如下：1.由于到21点时M=49 V=0.87 57若走过14则M大于了 50故直接从21点返回。第一次待货物 顺序为：0 26 31 27 39 36 38 35 32 23 17 21 0总的路程为：2.7122e+004 总的时间是：1.7301最优路线为：0-26-31-27-39-27-36-38-35-32-23-17-21-0走的距离W=27 122m,花费的时间T=l.73012.若按程序给出的从13到8的路线是13-12-11-12-8而当为13-11-12-8时更短 故修改之；同时到达40后如果选择34则45的周围全被遍历过。到45后M=46.83,V=l.0247不满足要求，故从40到34后沿21-26返回。第2次带货物 顺序为：Columns 1 through 130 18 13 12 11 8 3 1 6 7 10 9 14Columns 14 through 21 16 43 42 49 50 40 34 0总的路程为：8.4850e+004 总的时间是：4.5354最优路线为：0-18-13-11-12-8-3-1-6-1-7-10-9-14-16-23-32-35-38-43-42-49-50-40-45-36-21-0走得距离是：W=83220,所用的时间是T=4.46753.当到达45点时若要去20点放货物的话则需要遍历许多已经遍历过的地点,故从45点沿36-21-0返回第3次带货物 顺序为：Columns 1 through 13 0 24 19 25 29 22 30 28 33 46 48 44 41Columns 14 through 19 37 47 45 20 15 0总的路程为:1.5936e+005总的时间是：7.539目最优路线为：0-26-31-24-19-25-29-22-30-28-33-46-48-44-41-37-40-47-40-45-36-21-0所走的距离为：W=12897 0m,所用的时间是：T=6.12384.只余下了 5个点,所以由图可知第4次带货物顺序为：0 4 2 5 0总的路程为：1.9347e+005总的时间是：8.2112路线为：0-26-31-24-19-25-15-22-20-2-5-2-4-3-8-12-13-18-O总路程是：W=17 1510m所用的时间是T=7.3964由上面的四个阶段可以知道该问的最优路线是：0-26-31-27-39-27-36-38-35-32-23-17-21-0-18-13-11-12-8-3-1-6-1-7-10-9-14-16-23-32-35-38-43-42-49-50-40-45-36-21-0-26-31-24-19-25-29-22-30-2 8-33-46-48-44-41-37-40-47-40-45-36-21-0-26-31-24-19-25-15-22-20-2-5-2-4-3-8-12-13-18-0总路程是：W=17 1510m所用的时间是T=7.39645、模型的分析误差分析：对于模型一是使用了精确地D ijkstra算法，故误差可以忽略不计对于模型二假定了 32到38点的断开存在一定的误差，但相对于断开其余 的几点得到的数值要小，故该模型可以使用。对于模型三，由于分区域的方法有很多，故不可避免的存在些许误差，但由于区域越多，路程越多，故选择分成4个区域最合适；分成的四个不同 7的时间的到达区域比较紧密故按照时间的不同划分了四个区域，从而大大的 消除了误差，此模型可以使用。对于模型四的误差比较大，由于未考虑货物的拆分可能会有一定的影响同 时由于4个阶段的划分也是有一定的不确定性故误差存在。对于该模型简化 了考虑的条件，仅以M和V为判断标准，虽对准确性存在挑战，但该模型相 对与其他的分类有明显的优越性。故该模型适用于该问的求解。灵敏度分析对于模型一、二、三，灵敏度很好，模型的准确性很高。对于模型四由于质量和体积的制约，使其灵敏度不会很好，但准确性较 高，因此模型可以使用。6、模型评价、改进和推广6.1模型的评价优点：充分利用了已知数据建立模型，使其具有很高的准确性和可行性 使用了准确的算法和适当的假设，使模型的准确性和实用性到达统 运用功能强大的Matlab工具使数据处理误差达到最小缺点 由于数据较多，没法使用工具进行模型的验证，只能一步一步的精 化模型6.2模型的改进对于模型一和三主要是进行验证。对于模型二断开的那个点可以去取别的点进行。主要是模型四的改进，可以考虑到不同的地点送的货物进行拆分，从而渠 道最优的解6.3模型的推广可充分使用到图的遍历和最短路的一系列问题的求解中。7、参考文献1.A F irst Course in Mathenmatical Moderling(Third Edition)F rank R.Giordiano Maurice D.weir William P.F ox2.图论 任韩。3.数学建模案例选集姜启源谢金星4.图论第3版 德迪斯特尔著5.大学生数学建模竞赛辅导教材叶其效6.基于matlab动态规划中最短路线的实现程序J电脑学习 施益昌郑贤斌李自立7.物流配送问题的混沌优化算法研究 中央民族大学学报（自然科学版）2009年11月第18卷第4期8.D ijkstra算法在企业物流运输网络中的应用湖南农业大学学报（自然科学版）2005年8月第29卷4期附录附录1.、表格1.1各货物号信息表货物号送达地点重量（公斤）体积（立方米）不超过时间1132.500.03169：002180.500.03549：003311.180.02409：304261.560.035012：005212.150.030512：006141.7 20.010012：007171.380.010912：008231.400.042612：009320.7 00.048112：0010381.330.021910：1511451.100.02879：3012430.950.022810：1513392.560.059512：0014452.280.03019：3015422.850.019010：1516431.7 00.07 8210：1517320.250.041212：0018361.7 90.018412：0019272.450.044512：0020242.930.04209：0021310.800.01089：3022272.250.001812：0023261.570.021012：0024342.800.01039：3025401.140.01559：3026450.680.03829：3027491.350.014410：1528320.520.002012：0029232.910.048712：0030161.200.042912：003111.260.02503221.150.05017Z3331.630.04833441.230.00063551.410.03873660.540.00673770.7 00.01293880.7 60.03463992.140.008740101.070.012441111.370.051042122.390.042843130.990.004844141.660.049145150.450.020946162.040.009847171.950.032448182.120.055449193.870.026250202.010.032451211.380.041952220.390.000153231.660.050254241.240.053455252.410.001256261.260.005957270.420.022458281.7 20.058059291.340.037 260300.060.040261310.600.027 462322.190.050363331.890.049464341.810.032565351.000.005566361.240.017 767372.510.036168382.040.011069391.070.04407 0400.490.03297 1410.510.00947 2421.380.045573431.310.01217 4441.260.00057 5450.980.041371.2 50个位置点的坐标76461.350.024177472.120.023078480.540.054279491.010.056680501.120.028481250.7 90.001182462.120.049283322.7 70.003484232.290.005485200.210.049086251.290.008887191.120.024988410.900.003889462.380.043490371.420.002091321.010.030092332.510.013393361.170.002094381.820.030895170.330.034596110.300.017 297154.430.053698120.240.005699101.380.017 510071.980.0493位置点X坐标（米）Y坐标（米）1918550021445560372705704373567052620995610080143571002522808716025259138452680101193530501178503545126585418513763052001413405532515212559751615365704571.3相互到达信息17141657385188825807519585581652078083552112770856022220088352314765905524779093302544359525261086096352710385105002856597652925809865301565995531939510100321483510365331250109003472801106535153051137536123901141537641011510381391511610399510120504083451230041493013650421326514145431418014215443030150604510915142354623301450047773514550488851488049115751516050801015325序号位置点1位置点211321832204245386347742851595210611171812711381214914159101610181710718111219121320122521121522131823131924131125141826141627141728142129152230152531162332172333183134192435202236212637213638211739223040231741243142254143251944252945273146283347292248302849304175031265131345232355332235433465533285634405735385836455936276037406138366239276340346440456541446641376741466842436942497043387144487244507345507445427546487647407748447849507949428050408101882O21830261.4模型二中到达时的时问点到的时间最大允许的时间000180.09091130.27 061240.55871310.7 3291.5340.92971.5401.047 71.57451.23171.5421.42972.25491.56182.25431.7 3222.25381.89132.25362.00544272.147 24392.32144262.5544212.7 4544172.87 144232.99534323.154163.4424143.6007403.9031附录2、MATLAB程序代码2.1、D ijstra 求解clc clear all a=11000 8250;9185 500;1445 560;7270 570;3735 670;2620 995;100801435;10025 2280;7160 2525;13845 2680;11935 3050;7850 3545;6585 4185;7630 5200;13405 5325;2125 5975;15365 7045;14165 7385;8825 8075;5855 8165;780 8355;12770 8560;2200 8835;14765 9055;7790 9330;4435 9525;108609635;10385 10500;565 9765;2580 9865;1565 9955;9395 10100;1483510365;1250 10900;7280 11065;15305 11375;12390 11415;6410 11510;13915 11610;9510 12050;8345 12300;4930 13650;13265 14145;14180 14215;3030 15060;10915 14235;2330 14500;7735 14550;885 14880;11575 15160;8010 15325;%a是各个点的坐标 for i=l:51for j=l:51t=a(i,:)-a(j,:);c(i,j)=sqrt(t(1)A2+t(2)A2);%两点之间的直线距离 endenda=1 18;10 21;15 17;24 35;32 45;413;1 8;2 20;2 4;3 8;3 4;4 2;5 15;5 2;6 1;7 18;7 1;8 12;9 14;97;11 12;12 13;12 25;1222;15 31;25 23;33 44;4125;16 41;25 46;33 37;4123;1719;25 28;34 46;4223;18 29;2740;3543;4215;1331;19 31;28 38;36 49;4318;13 24;20 33;29 45;36 38;4419;1322;21 22;3027;3748;4411;14 26;21 28;30 40;38 50;4518;14 36;21 41;31 36;39 50;4516;14 17;22 26;31 27;40 42;4610;1017;1430;2334;3234;4048;4740;4844;4950;4942;5040;018;0 21;0 26;%通路表b=zeros(51);for i=l:83b(a(i,l)+l,a(i,2)+l)=l;b(a(i,2)+lza(i,l)+l)=l;enda=b.*c;7for i=l:51for j=l:51if a(i,j)=0 a(i,j)=inf;endif i=ja(i,j)=0;endendendw=a;for p=l:51n=size(w,1);wl=w(p,:);for i=l:nl(i)=wl(i);z(i)=1；ends=;s(D=l;u=s(1);k=l；while k1(u)+w(u,i)1(i)=1(u)+w(u,i);z(i)=u;endendendend11=1；for i=l:nfor j-1:k if i=s(j)ll(i)=ll(i);else11(i)=inf;endendendlv=inf;for i=l:nif 11(i)lv7lv=ll(i);v=i;endends(k+1)=v;k=k+l;u=s(k);endif p=la=l;t=z;elsea=a;l;t=t;z;endendfor i=l:51a(i,i)=inf;%把相同的点赋值为无穷大endsave w.txt a-ascii;%保存最小星巨离save t.txt t-ascii;%保存最小路径经过的点2.2、问题一得求解clear allclcformat shortw=.数据太多省略.;pl=7；p2=10;sum=0;w(:,1)=inf;w(:,pl)=inf;w(:,p2)=inf;w(13,16)=inf;w(16,13)=inf;xl=1,pl;x2=P2,l;for i=l:15si,tl=min(w(pl,:);s2,t2=min(w(p2,:);sum=sum+sl+s2;w(:,tl)=inf;w(:,t2)=inf;Pl=tl；p2=t2;if tl=9|t2=97disp(到达2 4时所走的距离，)disp(sum)T=sum/1000/24+3*i/60;disp。到24所用的时间，)disp(T)endif tl=t2xl=xl,t1;x=xl,x2;break;endxl=xl,tl;x2=t2,x2;x=xl,x2;enddisp顺序为：)disp(x)disp总的路程为：，)disp(sum)T=sum/1000/24+3*30/60;disp总的时间是：Idisp(T)2.3、问题二的2阶段求解clear allclcformat shortw=.数据太多省略.;P=l；x=l；sum=0;v=w;w(:,p)=inf;for i=l:10s,t=min(w(p,:);sum-sum+s;T=sum/1000/24+3*i/60;disp(t,T)w(:,t)=inf;P=t；x=x;t;enddisp顺序为：，)disp(x)disp。总的路程为：，)disp(sum)7T=sum/1000/24+3*30/60;disp。总的时间是：，)disp(T)2.4问 题二4阶段的解法clcclear allw=inf 5295.49 2182.03 4709.245295.49 inf 3113.46 5714.342182.03 3113.46 inf 3883.844709.24 5714.34 3883.84 inf；disp第一个区域，)P=l；x=l；sum=0;v=w;T=0;w(:,p)=inf;for i=l:3s,t=min(w(p,:);sum=sum+s;T=s/1000/24+T;disp(t,T)T=T+3/60;w(:,t)=inf;P=t；x=x;t;end disp顺序为，)disp(x*)dispx总路程是：，)disp(sum)disp总时间是)disp(T)disp第二个区域，)w=inf 1780.15 4104.9 5735.68 8234.281780.15 inf 2324.75 3955.53 6454.13 4104.9 2324.75 inf 1630.78 4847.795735.68 3955.53 1630.78 inf 3217.018234.28 6454.13 4847.79 3217.01 inf；P=l；X=1;7W(:,p)=inf;T=0.6087;for i=l:4 s,t=min(w(p,:);sum=sum+s;T=T+3/60;T=s/1000/24+T;disp(t,T)if(i=D；T=T+3/60;endif(i=4);T=T+3/60*2;end w(:,t)=inf;P=t；x=x;t;enddisp顺序为Idisp(x*)dispx总路程是：，)disp(sum)disp总时间是)disp(T)disp第三个区域，)w=inf 4719.88 2351.72 3269.39 4323.1 4719.88 inf 3536.11 2618.44 5507.49 2351.72 3536.11 inf 917.67 1971.38 3269.39 2618.44 917.67 inf 2889.05 4323.1 5507.49 1971.38 2889.05 inf；x=1 3 5 4 2;T=1.3317;for i=l:4 m=i;s=w(x(i),x(i+1);sum=sum+s;if(i=4)m=m+l;endT=s/1000/24+T;disp(x(i+l),T)T=T+3/60;enddisp(顺序为)7disp(x)disp(x总路程是：，)disp(sum)disp总时间是)disp(T)disp第四个区域，)W=.数据太多省略.;p=l;x=l；v=w;w(:,p)=inf;w(:,12)=inf;T=1.9413;for i=l:10s,t=min(w(p,:);sum=sum+s;T=s/1000/24+T;disp(t,T)T=T+3/60;w(:,t)=inf;P=t；x=x;t;if i=2T=T+3/60;endif i=4T=T+3/60;endif i=7T=T+3/60;endif i=8T=T+3/60*2;enddisp(p,sum)endsum=sum+v(t,12);disp。顺序是：,)disp(x,1)disp总距离是：,)disp(sum)T=sum/1000/24+3*30/60;disp总时间是：,)disp(T)2.5、问题3的初步设定7clcclear allw=.；i=l；while i50)|(V1)break;endn=i;x=x;t-1;%disp(t-l,M,V)a(:,t)=inf;endsum=sum+v(p,1);7disp顺序为：，)disp(x 0)disp总路程是:)disp(sum)T=sum/1000/24+3*i/60;disp所用时间是：，)disp(T)disp第二阶段，)P=l；x=0;M=0;v=o；a(:,p)=inf;for i=l:50s,t=min(a(p,:);M=M+w(t-l,2);V=V+w(t-l,3);sum=sum+s;if(M50)|(V1)break;endn=n+l;P=t；x=x;t-1;%disp(t-1,M,V)a(:,t)=inf;enddisp顺序为：,)disp(x,0)disp(，总路程是：，)disp(sum)T=sum/1000/24+3*i/60;disp所用时间是：I disp(T)disp(，第三阶段，)P=l；x=0;M=0;V=0;a(:,p)=inf;for i=l:50-ns,t=min(a(p,:);M=M+w(t-1,2);V=V+w(t-1,3);sum=sum+s;if(M50)|(V1)7break;endP=t;x=x;t-1;%disp(t-1,M,V)a(:,t)=inf;n=n+l;endsum=sum+v(p,1);disp L顺序为：，)disp(x,0)disp(，总路程是：，)disp(sum)T=sum/1000/24+3*i/60;disp所用时间是：Idisp(T)disp(，第四阶段，)P=l；x=0;M=0;v=o；a(:,p)=inf;for i=l:50-ns,t=min(a(p,:);M=M+w(t-l,2);V=V+w(t-l,3);sum=sum+s;if(M50)|(V1)break;endP=t；x=x;t-1;a(:,t)=inf;n=n+l;endsum=sum+v(p,1);disp顺序为：)disp(x,0)disp(总、品巨离是：)；disp(sum)T=sum/1000/24+3*i/60;disp总时间是:);disp(T)2.6、问题三的优化clcclear allload w.txt7a=w;load x.txtw=x;P=l；x=0;M=0;V=0;sum=0;v=a;a(:,p)=inf;for i=l:50s,t=min(a(p,:);M=M+w(t-l,2);V=V+w(t-l,3);sum=sum+s;P=t；if(M50)|(V1)break;endn=i；x=x;t-1;disp(t-l,M,V)a(:,t)=inf;endsum=sum+v(p,1);disp顺序为：)disp(x,0)disp(，总路程是：，)disp(sum)T=sum/1000/24+3*i/60;disp所用时间是：，)disp(T)disp(，第二阶段，)P=l；x=0;M-0；v=o；a(:,p)=inf;for i=l:50s,