钢管订购和运输

1、钢管订购和运输要铺设一条A1TA2T2A15的输送天然气的主管道，如图一所示。经筛选后可以生产这种主管27或者建有施工公路）圆圈表示火车站，每段铁路、公路和管道旁的阿拉伯数字表示道钢管的钢厂有S,S,AS。图中粗线表示铁路，单细线表示公路，双细线表示要铺设的管道（假设沿管1道或者原来有公路里程（单位km）。为方便计，1km主管道钢管称为1单位钢管。一个钢厂如果承担制造这种钢管，至少需要生产500个单位。钢厂S在指定期限内能生产该钢管的最大数量为s个单位，钢管出厂销价1单位钢管为p万元，iii如下表：i1234567si80080010002000200020003000pi1601551551

2、601551501601单位钢管的铁路运价如下表：里程（km）W300301350351400401450451500运价（万元）2023262932里程（km）5016006017007018008019009011000运价（万元）37445055601000km以上每增加1至100km运价增加5万元。公路运输费用为1单位钢管每公里0.1万元（不足整公里部分按整公里计算）。钢管可由铁路、公路运往铺设地点（不只是运到点A,A,A,A，而是管道全线）。1215（1）请制定一个主管道钢管的订购和运输计划，使总费用最小（给出总费用）。（2）请就（1）的模型分析：哪个钢厂钢管的销价的变化对购运计划和

3、总费用影响最大，哪个钢厂钢管的产量的上限的变化对购运计划和总费用的影响最大，并给出相应的数字结果。摘要：本文建立一个钢管订购和运输模型，从钢厂到主管道结点的运费是影响总费用的重要因素。为使总费用最小，须使从钢厂到主管道结点的运费一一钢管运输费最小。对求网络中最短路径的Dijkstra算法进行改进，得到新的算法，可对含多种权重计算方式的网络进行搜索，得出最小费用路径（最短路径）。在此基础上，建立起描述总费用的函数，把钢管的订购和运输问题归结为在一定约束条件下求最小总费用的二次规划问题。对于问题（1）,运用Lingo软件包求出了较优的订购和运输计划（见表4,表5）,其最小费用为1278632万元。

4、对于问题（2）而言，可得出钢厂S钢厂销价变化对总费用影响最大，S,S钢厂钢管的销价的变656化对购运计划影响最大，S钢厂钢管的产量的上限的变化对总费用影响最大，购运计划影响较小。1关键词：穷举法非线性规划一、问题重述和分析要铺设一条ATATATA的输送天然气的主管道，如图一所示，经筛选后可以生产这种主管1215道的钢厂有S,S,A,S。图中粗线表示铁路，单细线表示公路，双细线表示要铺设的管道（假设沿管道127或者原来有公路，或者建有施工公路），圆圈表示火车站，每段铁路、公路和管道旁的阿拉伯数字表示里程（单位km）。为了方便，1km主管道称为1单位钢管。一个钢厂如果承担制造这种钢管，至少需要生产

5、500个单位。钢厂S在指定期限内能生产该钢管的i最大生产数量为s个单位，钢厂出厂销价为p万元，如下表：iii1234567si80080010002000200020003000160155155160155150160pi1单位钢管的铁路运价如下表：表2里程（km）<300301350351400401450451500运价（万元）2023262932里程（km）5016006017007018008019009011000运价（万元）37445055601000km以上每增加1至100km运价增加5万元。公路运输费用为1单位管道每公里0.1万元（不足整公里的按整公里计算）。管道可由铁

6、路、公路运往铺设地点（不只是运到点ATAtATA，而是管道全线）。1215要求：（1）制定一个主管道钢管的订购和运输计划，使总费用最小，并给出总费用。（2）就（1）的模型进行分析：哪个钢厂钢管的销价的变化对购运计划和总费用影响最大，哪个钢厂钢管的产量的上限的变化对购运计划和总费用的影响最大，并给出相应的数字结果。二、基本假设1在计算运费时，沿管道铺设路线上的公路与其它普通公路相同（1单位钢管每公里0.1万元）；2订购的钢管数量刚好等于需要铺设的钢管数量；3. 管道可由铁路、公路、管道全线运往铺设地点（不只是运到点A,A,A,A）；12154. 模型只考虑钢管销价费用和钢管从钢管厂运送到铺设点的

7、钢管运费，而不考虑其它费用，如不计换车、转站的时间和费用，不计装卸费用等；5. 不计运输时由于运输工具出现故障等意外事故引起工期延误造成损失；6. 销售价和运输价不受市场价格变化的影响。三、符号说明S：i第i钢管厂s：i表示S的最大生产能力iA：j表示需要铺设管道路径上的车站x：ij从所有S运往A的钢管数ijc：ij表示单位钢管从S地运往A地的最小费用ijp：i从S订购钢管的单位价格iQ：订购的所有钢管全部运到A(j=1,2,A,15)点的总运费jT：当钢管从钢厂S.运到点A后，钢管向A的左右两边运输(铺设)管道的运输费用ijjZ：用于订购和运输的总费用y：运到A地向左铺设的数目jjz：运到A

8、地向右铺设的数目jjd：单位钢管1公里的公路运输费用A:表示A和A之间需要铺设的管道长度j,j+1jj+1四、模型的建立与求解问题一1、模型的建立钢管的订购和运输方案是直接影响工程费用的主要原因，因此，选取费用最小的路线运送货物，合理的订购计划是决定该工程费用的重要因素，首先利用图论的方法，来确定从钢管生产厂家到施工结点的费用最小路线，然后建立工程费用的优化模型，从中优化出最佳购运方案。对本问题而言，实际上是一个要求制定订购和运输计划，使总费用最小的优化问题。本模型的总费用包括钢管的销价和运输总的费用。首先，向某厂订购钢管，然后将在每个厂订购的钢管运往需要铺设的全路段。欲解决本问题可以按以下方

9、案进行思考：首先，需要确定将货物从i地运往j地的最优路线(费用最小)；然后，求出向每个钢管厂的订购计划，并确定出运输计划；最后计算将运往j地的钢管铺到各个管道上的运输费用，我们不妨假设运往以j为终点的钢管只铺到与j点相邻的两段管道上。因此，本问题可以按以下步骤求解。第一步：确定从i地到j地的最优路径，从而确定出单位钢管从i地运往j地的最小运费。s(i=1,2,A7)表示钢管厂S(i=1,2,A7)的最大生产能力，A(j=h2,A,15)表示需要铺设钢管路iij径上的车站。假设从S运往A的钢管用于铺设A点左右侧的钢管数为x单位，单位产品从S到A地ijji,jij的运费为f万元，用c表示单位钢管从

10、S地运往A地的最小费用，贝y：i,ji,jijc=minF(1)ijij第二步：建立从S厂运送x单位钢管到A点的运费的模型：ii,jj用q表示订购的所有钢管全部运到A(j=1,2,A,15)点的总运费，贝y：jQ=玆7xc；(2)ijijj=1i=1第三步：将运到A处的钢管铺到相邻两段路上的运输费用j对于运到A的钢管，它向左运输的总量y，它向左运输的总费用为：jjyxd+(y一l)xd+(y一2)xdLlxdjjj=0.1x(1+2+L+y)=0.05y(y+1)(万元);jjj同理它向右运输的总费用为：=0.05z(z+1)jj6)用T表示当钢管从钢厂S运到点A后，钢管向A的左右两边运输(铺

11、设)管道的运输费用，得：ijjT=0.05艺G+yLjj=1)y+(1+zjj3)y和z之间存在的关系：jj£x=y+z;(j=1,2,L,15)】ijjj(4)1=1z+y=A；(j=1,2,L,14)jj+1j,j+1A表示A和A之间需要铺设的管道长度)j,j+1jj+1第四步：建立订购费用的模型设W表示订购管道的总费用，则可建立如下模型：W=工工pxii,ji=1j=15)又因为一个钢厂如果承担制造钢管任务，至少需要生产500个单位，钢厂S在指定期限内最大生产量i为s.个单位，故500<£x<s或瓦x=0，用Z表示订购和运输的总费用，由(2)、(3)、(4

12、)、iijiijj=2j=2(5)，本问题可建立如下的非线性规划模型：目标函数m1nZ=W+Q+T=£7£15(p+c)x+0.05iij1=1j=1ijj=1+y)y+(1+z)zjjjj约束条件工x=y+z；(j=1,2,L,15)1jjj1=1z+y=A;(j=1,2,L,14)<jj+1j,j+1500x<s或£x=0;(i=1,2,L,7)ijijj=2j=2x>0i=1,L,7,j=2,L,15ij（A表示A和A之间需要铺设的管道长度）j,j+1jj+12、模型的求解（1）首先求解cij由于钢管从钢厂S运到运输点A要通过铁路和公路运输

13、，而铁路运输费用是分段函数，与全程运输ij总距离有关。又由于钢厂S直接与铁路相连，所以可先求出钢厂S到铁路与公路相交点b的最短路径。iij依据钢管的铁路运价表，算出钢厂S到铁路与公路相交点b的最小铁路运输费用，并把费用作为边权赋ij给从钢厂S到b的边。再将与b相连的公路、运输点A及其与之相连的要铺设管道的线路（也是公路）ijji添加到图上，根据单位钢管在公路上的运价规定，得出每一段公路的运费，并把此费用作为边权赋给相应的边。这样就转换为以单位钢管的运输费用为权的赋权图，再利用E.W.Dijkstra的最短路算法计算出一个minZ=工25i=1j=1(p+c)xiijij+0.05国(1+y)y

14、+G+z)zLjjjj-j=1单位钢管从钢厂运到工地的最少费用系数阵c，MATLAB程序见附录一。ij表3表中对应的值为对应两点最优路径单位钢管运输费用S1S2S3S4S5S6S7Al170.7215.7230.7260.7255.7265.7275.7A2160.3205.3220.3250.3245.3255.3265.3A3140.2190.2200.2235.2225.2235.2245.2A498.6171.6181.6216.6206.6216.6226.6A538111121156146156166A620.595.5105.5140.5130.5140.5150.5A73.18

15、696131121131141A821.271.286.2116.2111.2121.2131.2A964.2114.248.284.279.284.299.2A10921428262576276A11961468651335166A121061569661514556A13121.2171.2111.276.271.226.238.2A1412817811883731126A151421921329787282（2）根据以上结果,继续求解非线性规划模型：工x=y+z；(j=1,2,L,15)ijjji=1z+y=A;(j=1,2,L,14)s.t十j吕+1500正x<s或艺x=0;(i

16、=1,2,L,7)ijiijj=2j=2x>0i=1,L,7,j=2,L,15ij由于不能直接处理约束条件：500瓦x<s或ijj=22x=0，我们可先将此条件改为瓦iijj=2x<s，得ijij=2minZ=工艺(p+c)xiijiji=1j=1到如下模型：+0.05国G+y)y+G+z)zLjjjj-j=1工x=y+z；(j=1,2,L,15)ijjji=1z+y=A;(j=1,2,L,14)st<jj+1j,j+1为x<s;(i=1,2,L,7)ijij=2x>0i=1,L,7,j=2,L,15ij用LINGO求解，程序见附录二。分析结果后发现购运方案

17、中钢厂S的生产量不足500单位，下面我7们采用不让钢厂S生产和要求钢厂S的产量不小于500个单位两种方法计算：771）不让钢厂S生产，程序见附录三。7计算结果：Z=1278632（万元）（此时每个钢厂的产量都满足条件）.12）要求钢厂S的产量不小于500个单位，程序见附录四。7计算结果：Z=1285281（万元）（此时每个钢厂的产量都满足条件）.2比较这两种情况，得最优解为，minZ=min（Z,Z）=Z=1278632（万元）。121所以根据上述的模型，得运输总费用最小为：1278632（万元）。具体的购运计划和铺设方案如表4，表5：表4问题一的订购和调运方案订购量A2A3A4A5A6A7A

18、8A9A10AiiA12A13A14A15S1800004029520026500000000S80017900321003000000000S31000003360000664000000S4000000000000000S10150508920000004150000S6155600000000351086333621165S7000000000000000表5问题一的铺设方案YzA10.0000000.000000A2104.000075.00000A3226.0000282.0000A4468.00000.000000A5606.00009.500000A6184.500015.500

19、00A7189.500076.00000A8125.0000175.0000A9505.0000159.0000A10321.000030.00000A11270.0000145.0000A1275.0000011.00000A13199.0000134.0000A14286.0000335.0000A15165.00000.000000问题二：针对问题一的求解模型，讨论钢厂钢管的销售价格变化对购运计划和总费用影响及钢厂钢管产量的上限变化对购运计划和总费用的影响定义方案中运往各点A的运输量的变化量的绝对值之和称为运输方案变化量i1、讨论钢厂钢管的销售价格变化对购运计划和总费用的影响当钢厂钢管销

20、售价格变化时，会对购运计划和总费用造成影响。为了更好地观察每一个钢厂钢管销售价格所造成的影响，采用比较法，即每次只让一个钢厂钢管的销售价格发生相同的变化，其余钢厂钢管的销售价格不发生变化。我们将各个钢厂单位钢管的销价分别增加1万元和减少1万元，借助LINGO软件得出相应的总费用、运输方案、订购方案变化情况如表6、表7所示表6各个钢厂单位钢管的销价分别增加1万元钢厂总费用总费用变化量运输方案变化量订购方案变化量S1127943280000S2127943280000S31279632100000S41278632000S5127963910074030S612798341202712712S71

21、278632000表7各个钢厂单位钢管的销价分别减少1万元钢厂总费用总费用变化量运输方案变化量订购方案变化量S1127783280000S2127783280000S31277632100000S41278632000S512772631369712712S6127706815644030S71278632000由上述表格观察分析可得：S钢厂销价变化对总费用影响最大，S,S钢厂钢管的销价的变化对购656运计划影响最大。2、讨论钢厂钢管产量的上限的变化对购运计划和总费用的影响同样采用比较法，即每次只让一个钢厂钢管产量的上限的发生相同的变化，其余钢厂钢管产量的上限不发生变化。将各个钢厂的产量的上限

22、分别增加100个单位和减少100个单位，分别计算，得到购运计划和总费用变化情况如表8、表9所示表8各个钢厂钢管的产量的上限分别增加100个单位钢厂总费用总费用变化量运输方案变化量订购方案变化量S1126833210300218200S212751323500404200S3127613225001786200S41278632000S51278632000S6127863208440S71278632000表9各个钢厂钢管的产量的上限分别减少100个单位钢厂总费用总费用变化量运输方案变化量订购方案变化量S1128893210300260200S2128213235001244200S31281

23、1322500200200S41278632000S51278632000S61278632000S71278632000由上述表格观察分析可得：S1钢厂钢管的产量的上限的变化对总费用影响最大，购运计划影响较小。五、结果分析由于总费用由订购费用和运输费两部分组成，运输费又由一般线路上的运输费和铺设管道上的运输费组成。利用求网络中最短路径的Dijkstra算法，进行改进得到新的算法，可对含多种权重计算方式的网络进行搜索，得出最小费用路径（最短路径），算出两点之间的最优路径，进而根据非线性规划，借助于Lingo软件求解即可求出相应的结果。六、模型的评价及改进1优点：1) 本问题中运用了求网络中最短

24、路径的Dijkstra算法，进行改进得到新的算法，可对含多种权重计算方式的网络进行搜索，算出两点之间的最优路径，计算结果准确。2) 本问题构造出的模型算法较简单，也可以运用相应的其他编程软件来得到比较满意的结果。2缺点：1) 由于本问题有现成的比较先进的解法，但由于缺乏基本的数学软件资料和相应的算法资料，不能将其准确求解。2) 于在求解最短路时，我们用人工计算容易将问题复杂化，同时，容易出错。3) 作为图论问题的技术而言，求解过程较难，且不易求出最优解3模型改进本问题中需要用现成的先进的网络流求解，但本问题没有求解，实际上本问题中加权穷举法，可采用相应的网络法来代替。在求解问题三时可用最小树求

25、解，同时，我们可将本问题运用于时间的变化等范围的推广。参考文献：1 甘应爱，田丰等等.运行学.清华大学出版社,北京,1994。2 袁亚湘.孙文瑜著.最优化理论与方法.科学出版社，北京，1997.3 徐俊明著.图论及其应用.中国科学技术大学出版社，合肥，1997.附录一S1-S7与A1-A15间的最小费用矩阵A=;fori=1:39forj=1:39A(i,j)=inf;endendforj=1:39A(j,j)=0;endA(1,2)=104;A(2,1)=104;A(2,3)=301;A(3,2)=301;A(3,4)=750;A(4,3)=750;A(4,5)=606;A(5,4)=606

26、;A(5,6)=194;A(6,5)=194;A(6,7)=205;A(7,6)=205;A(7,8)=201;A(8,7)=201;A(8,9)=680;A(9,8)=680;A(9,10)=480;A(10,9)=480;A(10,11)=300;A(11,10)=300;A(11,12)=220;A(12,11)=220;A(12,13)=210;A(13,12)=210;A(13,14)=420;A(14,13)=420;A(14,15)=500;A(15,14)=500;A(15,22)=20;A(22,15)=20;A(15,39)=20;A(39,15)=20;A(14,38)=

27、30;A(38,14)=30;A(14,21)=110;A(21,14)=110;A(13,37)=62;A(37,13)=62;A(12,35)=10;A(35,12)=10;A(11,34)=10;A(34,11)=10;A(10,32)=70;A(32,10)=70;A(9,31)=42;A(31,9)=42;A(8,30)=12;A(30,8)=12;A(7,16)=31;A(16,7)=31;A(7,29)=10;A(29,7)=10;A(6,28)=5;A(28,6)=5;A(5,27)=10;A(27,5)=10;A(4,26)=600;A(26,4)=600;A(3,24)=2

28、;A(24,3)=2;A(2,23)=3;A(23,2)=3;A(23,25)=450;A(25,23)=450;A(24,25)=80;A(25,24)=80;A(25,26)=1150;A(26,25)=1150;A(26,30)=1100;A(30,26)=1100;A(17,30)=1200;A(30,17)=1200;A(30,16)=202;A(16,30)=202;A(16,29)=20;A(29,16)=20;A(29,28)=195;A(28,29)=195;A(28,27)=306;A(27,28)=306;A(30,31)=720;A(31,30)=720;A(31,18

29、)=690;A(18,31)=690;A(31,32)=520;A(32,31)=520;A(32,33)=170;A(33,32)=170;A(33,19)=690;A(19,33)=690;A(33,36)=160;A(36,33)=160;A(33,34)=88;A(34,33)=88;A(34,20)=462;A(20,34)=462;A(33,36)=160;A(36,33)=160;A(36,35)=70;A(35,36)=70;A(36,37)=320;A(37,36)=320;A(37,38)=160;A(38,37)=160;A(38,21)=70;A(21,38)=70;A

30、(38,39)=290;A(39,38)=290;A(39,22)=30;A(22,39)=30;A;H=;HH=;forsw=16:39L=;fori=1:40h(i)=0;endfori=16:39L(i)=A(sw,i);endT=L(16);s=16;fori=17:39ifT>L(i);T=L(i);s=i;endendh(s)=1;forw=16:38fori=16:39whileh(i)=1i=i+1;endifi=40break;endL1(i)=A(s,i);L3(i)=T+L1(i);ifL3(i)<L(i)L(i)=L3(i);endendi=16;while

31、h(i)=1i=i+1;endifi=40T=L(39);else%设最小数_T=L(i);s=i;endfori=16:39%比较得最小数_whileh(i)=1i=i+1;endifi=40break;endifT>L(i);T=L(i);s=i;endendh(s)=1;endH(sw,:)=L;endfori=16:39%将距离H转化为费用HH_forj=16:39HH(i,j)=ch741(H(i,j);%ch741函数是将距离转化为费用.endendfori=l:15%将公路的费用也输入费用矩阵HH中forj=1:39HH(i,j)=A(i,j)*0.1;endendfori

32、=16:39forj=1:15HH(i,j)=A(i,j)*0.1;endendHH;%计算任意两点间的最小费用，并将数据输入矩阵日日日中_HHH=;forsw=16:22;%计算S1S7与其他点的最小费用_L=;fori=1:40h(i)=0;endfori=1:39L(i)=HH(sw,i);endT=L(1);s=1;%T为最小数的暂时空间fori=2:39ifT>L(i);T=L(i);%找最小数_s=i;%S为已用过的标志endendh(s)=1;forw=1:38%W为循环次数_fori=1:39%换数whileh(i)=1i=i+1;endifi=40break;endL1

33、(i)=HH(s,i);L3(i)=T+L1(i);ifL3(i)<L(i)L(i)=L3(i);endendi=1;whileh(i)=1i=i+1;endifi=40T=L(39);else%设最小数_T=L(i);s=i;endfori=1:39%比较得最小数whileh(i)=1i=i+1;endifi=40break;end%T;_%LL(i);ifT>L(i);T=L(i);s=i;endendh(s)=1;endHHH(sw-15,:)=L;endforj=1:7fori=1:15c(j,i)=HHH(j,i);endenddisp('sl-S7与A1-A15

34、间的最小费用矩阵')c将距离转化为费用函数是:functionf=ch741(ll)rr=ll;ifrr=0real1=0;elseifrr<300.5&rr>0real1=20;elseifrr>300.5&rr<350.5real1=23;elseifrr>350.5&rr<400.5real1=26;elseifrr>400.5&rr<450.5real1=29;elseifrr>450.5&rr<500.5real1=32;elseifrr>500.5&rr<

35、600.5real1=37;elseifrr>600.5&rr<700.5real1=44;elseifrr>700.5&rr<800.5real1=50;elsereal1=fix(rr-800.5)/100)*5+55;endf=real1;附录二model:sets:gch/1.7/:p,s;gd/1.15/:A,y,z;links(gch,gd):x,c;endsetsdata:p=160155155160155150160;s=80080010002000200020003000;c=170.7160.3140.298.63820.53.121

36、.264.29296106121.2128142215.7 205.3190.2171.611195.58671.2114.2142146156171.2178192230.7 220.3200.2181.6121105.59686.248.2828696111.2118132260.7 250.3235.2216.6156140.5131116.284.262516176.28397255.7 245.3225.2206.6146130.5121111.279.257335171.27387265.7 255.3235.2216.6156140.5131121.284.262514526.2

37、1128275.7 265.3245.2226.6166150.5141131.299.276665638.2262;enddatamin=W+Q+T;W=sum(links(i,j):p(i)*x(i,j);Q=sum(links(i,j):c(i,j)*x(i,j);T=sum(gd(j):(1+y(j)*y(j)+(1+z(j)*z(j)*0.05;z(1)+y(2)=104;z(2)+y(3)=301;z(3)+y(4)=750;z(4)+y(5)=606;z(5)+y(6)=194;z(6)+y(7)=205;z(7)+y(8)=201;z(8)+y(9)=680;z(9)+y(10

38、)=480;z(10)+y(11)=300;z(11)+y(12)=220;z(12)+y(13)=210;z(13)+y(14)=420;z(14)+y(15)=500;sum(gch(i):x(i,1);sum(gch(i):x(i,2);sum(gch(i):x(i,3);sum(gch(i):x(i,4);sum(gch(i):x(i,5);sum(gch(i):x(i,6);sum(gch(i):x(i,7);sum(gch(i):x(i,8);sum(gch(i):x(i,9);=sum(gch(i)=sum(gch(i)=sum(gch(i)=sum(gch(i)=sum(gch

39、(i)=sum(gch(i)y(1)+z(1)y(2)+z(2)y(3)+z(3)y(4)+z(4)y(5)+z(5)y(6)+z(6)y(7)+z(7)y(8)+z(8)y(9)+z(9)y(10)y(11)y(12)y(13)y(14)+z(14)y(15)+z(15)+z+z+z+z(10)(11)(12)(13):x(i,10);:x(i,11);:x(i,12);:x(i,13);:x(i,14);:x(i,15);for(gch(i):sum(gd(j):x(i,j)<=s(i);end附录三model:sets:gch/1.7/:p,s;gd/1.15/:A,y,z;lin

40、ks(gch,gd):x,c;endsetsdata:p=160155155160155150160;s=80080010002000200020003000;c=170.7160.3140.298215.7205.3190.2171.6111230.7220.3200.2181.6121260.7250.3235.2216.6156255.7245.3225.2206.6146265.7255.3235.2216.6156275.7265.3245.2226.6166enddatamin=W+Q+TfW.63820.53.121.264.29296106121.212895.58671.21

41、14.2142146156171.21781929686.248.2828696111.2118132131116.284.262516176.28397121111.279.257335171.27387131121.284.262514526.21128141131.299.276665638.2262;105.5140.5130.5140.5150.5=sum(links(i,j):p(i)*x(i,j);=sum(links(i,j):c(i,j)*x(i,j);=sum(gd(j):(1+y(j)*y(j)+(1+z(j)*z(j)*0.05;(1)+y(2)=104;(2)+y(3

42、)=301;(3)+y(4)=750;(4)+y(5)=606;142z(5)+y(6)=194;z(6)+y(7)=205;z(7)+y(8)=201;z(8)+y(9)=680;z(9)+y(10)=480;z(10)+y(11)=300;z(11)+y(12)=220;z(12)+y(13)=210;z(13)+y(14)=420;z(14)+y(15)=500;y(1)+z(1)=sum(gch(i):x(i,1);y(2)+z(2)=sum(gch(i):x(i,2);y(3)+z(3)=sum(gch(i):x(i,3);y(4)+z(4)=sum(gch(i):x(i,4);y(

43、5)+z(5)=sum(gch(i):x(i,5);y(6)+z(6)=sum(gch(i):x(i,6);y(7)+z(7)=sum(gch(i):x(i,7);y(8)+z(8)=sum(gch(i):x(i,8);y(9)+z(9)=sum(gch(i):x(i,9);y(10)+z(10)=sum(gch(i):x(i,10);y(11)+z(11)=sum(gch(i):x(i,11);y(12)+z(12)=sum(gch(i):x(i,12);y(13)+z(13)=sum(gch(i):x(i,13);y(14)+z(14)=sum(gch(i):x(i,14);y(15)+z(1