matlab钢管的订购和运输,钢管的订购和运输模型程序

钢管的订购和运输模型程序

模型一的matlab程序 !A为邻接矩阵；T，T2分别为铁路的矩阵和铁路的费用矩阵；R，R2分别为公路的矩阵和公路的费用矩阵；w,Q为总费用矩阵；cost为S到A的各项费用矩阵 A=ones(39,39); for i=1:39; for j=1:39; if i==j A(i,j)=0; else A(i,j)=inf; end end end A(8,10)=450;A(9,10)=80;A(10,11)=1150;A(11,16)=1100;A(12,13)=306; A(13,14)=195;A(1,14)=20;A(1,16)=202;A(2,16)=1200;A(16,17)=720; A(3,17)=690;A(17,18)=520;A(18,19)=170;A(4,19)=690;A(19,20)=160; A(5,15)=462;A(15,19)=88;A(20,21)=70;A(20,22)=320;A(22,23)=160; A(6,23)=70;A(23,24)=290;A(7,24)=30; for j=1:39; for i=1:j-1; A(j,i)=A(i,j); end end T=A; m=1; while mT(i,m)+T(m,j); T(i,j)=T(i,m)+T(m,j); end end end m=m+1; end T2=ones(39,39); for i=1:39; for j=1:39; if T(i,j)==0 T2(i,j)=T(i,j); elseif T(i,j)>0 else T2(i,j)=60+ceil((T(i,j)-1000)/100)*5; end end end B=ones(39,39); for i=1:39; for j=1:39; if i==j B(i,j)=0; else B(i,j)=inf; end end end B(8,26)=3;B(9,27)=2;B(11,28)=600;B(12,29)=10;B(13,30)=5; B(14,31)=10;B(1,31)=31;B(16,32)=12;B(17,33)=42;B(18,34)=70; B(15,35)=10;B(21,36)=10;B(22,37)=62;B(6,38)=110;B(23,38)=30; B(24,39)=20;B(7,39)=20; B(25,26)=104;B(26,27)=301;B(27,28)=750;B(28,29)=606;B(29,30)=194; B(30,31)=205;B(31,32)=201;B(32,33)=680;B(33,34)=480;B(34,35)=300; B(35,36)=220;B(36,37)=210;B(37,38)=420;B(38,39)=500; for j=1:39; for i=1:j-1; B(j,i)=B(i,j); end end R=B; g=1; while gR(i,g)+R(g,j); R(i,j)=R(i,g)+R(g,j); end end end g=g+1; end R2=0.1.*R; w=ones(39,39); for i=1:39; for j=1:39; if T2(i,j)>R2(i,j) w(i,j)=R2(i,j); else w(i,j)=T2(i,j); end end end f=1; while fw(i,f)+w(f,j) w(i,j)=w(i,f)+w(f,j); end end end f=f+1; end Q=ones(39,39); for i=1:7 for j=25:39 cost(i,j)=w(i,j) end end 模型一的matlab结果 cost = Columns 25 through 39 170.7000 160.3000 140.2000 98.6000 38.0000 20.5000 3.1000 21.2000 64.2000 92.0000 96.0000 106.0000 121.2000 128.0000 142.0000 215.7000 205.3000 190.2000 171.6000 111.0000 95.5000 86.0000 71.2000 114.2000 142.0000 146.0000 156.0000