钢管的订购和运输模型程序
19569-数学建模-钢管订购和运输 (1)
钢管订购和运输张伟 丁林阁 邓小涛 指导教师：数模组 海军航空工程学院摘要 本模型研究了管道铺设过程中钢管的订购和运输问题，它通过图论和非线性规划的知识建立。
模型使总费用达到最小，很好地解决了向哪个钢厂定货，定货多少，如何运输的问题，并且可以推广到更一般的网络。
同时针对模型中涉及的变量多、求解复杂这一问题，我们对模型进行了适当的简化,大大减少了变量的个数,从而减少了计算量。
一、问题重述要铺设一条1521A A A →→→ 的天然气主管道，可以生产这种主管道钢管的钢厂有721,,S S S 七家。
钢厂的位置，管道的铺设路线，以及从钢厂到铺设地的运输网络（运输网络包括沿管道的公路）均已知。
每个钢厂的钢管价格及其生产能力不全一样，且一个钢厂若要生产这种钢管，至少需要生产500个单位（1千米钢管记为1个单位）。
铁路的运价和公路的运价不一样。
要求在这种情况下，（1）制定一个钢管的订购和运输计划，使总费用最小，并给出总费用。
(2)分析：哪个钢厂钢管的销价的变化对购运计划和总费用影响最大，哪个钢厂钢管的产量的上限的变化对购运计划和总费用的影响最大，给出相应的数字结果。
(3)如果铺设的管道不是一条线，而是一个树形图，铁路、公路和管道构成网络，就这种更一般的情形给出一种解决办法，并对图二按（1）的要求给出模型和结果。
二、问题的假设在问题所给条件成立的前提下，我们进一步作如下假设： 1. 假设公路运输费用不是整公里的按整公里计算是合理的。
2. 假设沿管道的公路（施工公路）运输费用也为每公里0.1万元（不足整公里部分按整公里计算）。
3. 假设不考虑铁路、公路及施工公路的运输能力限制。
4. 假设运输费用为单程运输的费用，即从出发点到目的地的单程费用，不考虑空车返回的费用。
5. 假设运输费用已包含装卸费用。
关于假设的一点说明：根据上述假设我们认为在铺设管道的过程中每隔一公里，卸下一单位钢管供工人铺设是合理的。
三、符号约定i S ：生产主管道钢管的钢厂 ；j A ：管道节点 ；1,+j j l ：从j A 到1+j A 铺设钢管的路段长度（单位：公里，14,...,1=j ）； i s : 钢厂i S 在指定期限内生产钢管的最大数量（单位：单位钢管）； i P : 钢厂i S 单位钢管的出厂价格（单位：万元）； ij x ：从 钢厂i S 运到j A 的钢管数量（单位：单位钢管）； ij c ：表示1单位钢管从 钢厂i S 到j A 的最小费用（单位：万元）； j X ：运到j A 的钢管总数（单位：单位钢管）； j L ： 从j A 往左铺设的钢管总数（单位：单位钢管），j L 为j X 的一部分； j R ： 从j A 往右铺设的钢管总数（单位：单位钢管），这里j j j L X R -=； 其中 15,...,1;7,...,1==j i 四、问题分析本问题分两部分：一部分是图论中的最短路径的问题：确定1单位钢管从 钢厂i S 到j A 的最小费用；另一部分是非线性规划问题：求总的最小费用。
钢管订购与运输问题一的数学模型与求解
钢管订购与运输问题一的数学模型与求解
钢管订购与运输问题是一种组合优化问题，它涉及到钢管的订购和运输，旨在找到最佳的订购和运输方案，以最小的成本获得最大的收益。
这个问题通常可以用数学模型来表示。
设 n 个工地需要订购 m 根钢管，钢管订购和运输费用分别为
c1(订购费用)、c2(运输费用),订购钢管的最早时间 t0 为早订购时间，最迟时间为 t1 为晚订购时间，运输时间不计费用。
则钢管订购与运输问题的数学模型可以表示为:
minimize Σi=1~n c1(t1-t0) + Σj=i+1~n c2(t2-t1)
subject to:
t1≤t0
t2≥t1
t1+t2≤t0+30
x1=1, x2=1, ..., xnm=1
其中，x1、x2、...、xnm 是订购钢管的数量，1 表示订购，0 表示不订购。
通过这个数学模型，我们可以制定出钢管订购与运输问题的求解方法，以找到最佳的订购和运输方案。
在实际问题中，我们通常需要对求解结果进行评估和优化，以便找到更加优秀的方案。
因此，钢管订购与运输问题的数学模型和求解方法只是问题的第一步，实际应用中还需要进行进一步的分析和优化。
钢管的订购及运输优化方案
钢管的订购及运输优化方案钢管是一种常见的工业材料，主要用于建筑、桥梁、机器制造和能源开采等领域。
订购和运输钢管需要考虑多方面因素，如规格、数量、质量、运输距离、运输方式等。
本文将介绍一些钢管订购及运输的优化方案。
一、钢管订购方案1. 确定钢管规格和数量在订购钢管前，首先需要了解工程或项目的具体需求，确定钢管的规格和数量。
不同的工程或项目需要的钢管规格和数量可能会有所不同，选择合适的规格和数量可以避免浪费和损失。
2. 寻找可靠的供应商选择可靠的供应商可以确保钢管的质量和供应稳定性。
可以通过市场调研、参加行业展会或咨询同行业的项目经理、工程师等人员来寻找可靠的供应商。
3. 确定采购合同和交付方式在确定供应商后，需要签订采购合同并确定交付方式。
采购合同要明确规定钢管的规格、数量、价格和交付日期等具体条款，避免误解和纠纷。
交付方式可以选择集装箱运输、散装运输或其他方式，根据具体情况灵活选择。
4. 质量控制为确保钢管的质量，采购方可以要求供应商提供产品质量证明、实际样品或第三方检测报告。
在收到钢管后，可以进行抽检或全检，检查钢管的尺寸、表面状态、壁厚和材质等指标，避免存在不合格品质的钢管进入工程或项目。
二、钢管运输方案1. 选择合适的运输方式钢管的运输可以选择公路运输、铁路运输、水路运输或航空运输等方式。
具体选择哪种方式需要综合考虑运输距离、运输量、运输时间、运输成本及货物安全等各方面因素。
2. 管理运输过程在钢管运输过程中，需要对货车、火车、船舶或飞机等交通工具进行监控，确保运输过程中货物的安全。
可以使用GPS或其他定位技术实时掌握货物的位置和状态，及时处理运输中遇到的问题和风险。
3. 管理卸货和储存在将钢管卸货到工厂、工地或仓库后，需要将其储存到指定位置并标记钢管的规格、数量等信息。
可以采用RFID等智能化技术对钢管进行管理，便于日后的存储和使用。
4. 管理短途运输在项目工期中，可能需要短途运输钢管到具体施工位置。
钢管的订购和运输模型程序
模型一的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;endendendA(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);endendT=A;m=1;while m<=39for i=1:39for j=1:39if T(i,j)>T(i,m)+T(m,j);T(i,j)=T(i,m)+T(m,j);endendendm=m+1;endT2=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&T(i,j)<=300 T2(i,j)=20;elseif T(i,j)>300&T(i,j)<=350 T2(i,j)=23;elseif T(i,j)>350&T(i,j)<=400 T2(i,j)=26;elseif T(i,j)>400&T(i,j)<=450 T2(i,j)=29;elseif T(i,j)>450&T(i,j)<=500 T2(i,j)=32;elseif T(i,j)>500&T(i,j)<=600 T2(i,j)=37;elseif T(i,j)>600&T(i,j)<=700 T2(i,j)=44;elseif T(i,j)>700&T(i,j)<=800 T2(i,j)=50;elseif T(i,j)>800&T(i,j)<=900 T2(i,j)=55;elseif T(i,j)>900&T(i,j)<=1000 T2(i,j)=60;else T2(i,j)=60+ceil((T(i,j)-1000)/100)*5;endendendB=ones(39,39);for i=1:39;for j=1:39;if i==j B(i,j)=0;else B(i,j)=inf;endendendB(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);endendR=B;g=1;while g<=39;for i=1:39;for j=1:39;if R(i,j)>R(i,g)+R(g,j);R(i,j)=R(i,g)+R(g,j);endendendg=g+1;endR2=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);endendendf=1;while f<=39for i=1:39for j=1:39if w(i,j)>w(i,f)+w(f,j)w(i,j)=w(i,f)+w(f,j);endendendf=f+1;endQ=ones(39,39);for i=1:7for j=25:39cost(i,j)=w(i,j)endend模型一的matlab结果cost =Columns 25 through 39170.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.0000215.7000 205.3000 190.2000 171.6000 111.0000 95.5000 86.0000 71.2000 114.2000 142.0000 146.0000 156.0000 171.2000 178.0000 192.0000230.7000 220.3000 200.2000 181.6000 121.0000 105.5000 96.0000 86.2000 48.2000 82.0000 86.0000 96.0000 111.2000 118.0000 132.0000260.7000 250.3000 235.2000 216.6000 156.0000 140.5000 131.0000 116.2000 84.2000 62.0000 51.0000 61.0000 76.2000 83.0000 97.0000255.7000 245.3000 225.2000 206.6000 146.0000 130.5000 121.0000 111.2000 79.2000 57.0000 33.0000 51.0000 71.2000 73.0000 87.0000265.7000 255.3000 235.2000 216.6000 156.0000 140.5000 131.0000 121.2000 84.2000 62.0000 51.0000 45.0000 26.2000 11.0000 28.0000275.7000 265.3000 245.2000 226.6000 166.0000 150.5000 141.0000 131.2000 99.2000 76.0000 66.0000 56.0000 38.2000 26.0000 2.0000模型一的lingo程序model:!钢管购买与运输铺优化设问;sets:!cj表示厂家;md表示目的地；cj/1..7/:p,s,t;md/1..15/:y,z,D;link(cj,md):cost,x;endsets!这里是数据;data:s=800 800 1000 2000 2000 2000 3000;D=104,301,750,606,194,205,201,680,480,300,220,210,420,500,0;p=160,155,155,160,155,150,160;cost=170.7 160.3 140.2 98.6 38.0 20.5 3.1 21.2 64.2 92.0 96.0 106.0 121.2 128.0 142.0215.7 205.3 190.2 171.6 111.0 95.5 86.0 71.2 114.2 142.0 146.0 156.0 171.2 178.0 192.0230.7 220.3 200.2 181.6 121.0 105.5 96.0 86.2 48.2 82.0 86.0 96.0 111.2 118.0 132.0260.7 250.3 235.2 216.6 156.0 140.5 131.0 116.2 84.2 62.0 51.0 61.0 76.2 83.0 97.0255.7 245.3 225.2 206.6 146.0 130.5 121.0 111.2 79.2 57.0 33.0 51.0 71.2 73.0 87.0265.7 255.3 235.2 216.6 156.0 140.5 131.0 121.2 84.2 62.0 51.0 45.0 26.2 11.0 28.0275.7 265.3 245.2 226.6 166.0 150.5 141.0 131.2 99.2 77.0 66.0 56.0 38.2 26.0 2.0;enddata!目标函数;min=w;w=@sum(link(i,j):(p(i)+cost(i,j))*x(i,j))+0.1/2*@sum(md(j):y(j)^2+y(j)+z(j)^2+z(j));!约束条件;@for(cj(i):@sum(md(j):x(i,j))>=500*t(i);s(i)*t(i)>=@sum(md(j):x(i,j));@bin(t(i)));@for(md(j):@sum(cj(i):x(i,j))=y(j)+z(j));@for(md(j)|j#ne#15:D(j)=y(j)+z(j+1));z(15)=0;y(1)=0;@gin(@sum(link(i,j):x(i,j)));end铺设方案Variable Value Reduced CostY( 1) 0.000000 0.000000Y( 2) 75.00000 0.000000Y( 3) 282.0000 0.000000Y( 4) 0.000000 0.000000Y( 5) 9.500000 0.000000Y( 6) 15.50000 0.000000Y( 7) 76.00000 0.000000Y( 8) 175.0000 0.000000Y( 9) 159.0000 0.000000Y( 10) 30.00000 0.000000Y( 11) 145.0000 0.000000Y( 12) 11.00000 0.000000Y( 13) 134.0000 0.000000Y( 14) 500.0000 0.000000Y( 15) 0.000000 141.0500Z( 1) 0.000000 405.7500Z( 2) 104.0000 0.000000Z( 3) 226.0000 0.000000Z( 4) 468.0000 0.000000Z( 5) 606.0000 0.000000Z( 6) 184.5000 0.000000Z( 7) 189.5000 0.000000Z( 8) 125.0000 0.000000Z( 9) 505.0000 0.000000Z( 10) 321.0000 0.000000Z( 11) 270.0000 0.000000Z( 12) 75.00000 0.000000Z( 13) 199.0000 0.000000Z( 14) 286.0000 0.000000Z( 15) 0.000000 0.000000 厂家分配方案X( 1, 1) 0.000000 28.00000X( 1, 2) 0.000000 28.00000X( 1, 3) 0.000000 23.00000X( 1, 4) 250.4441 0.000000X( 1, 5) 84.05588 0.000000X( 1, 7) 265.5000 0.000000 X( 1, 8) 0.000000 23.00000 X( 1, 9) 0.000000 99.00000 X( 1, 10) 0.000000 143.0000 X( 1, 11) 0.000000 171.0000 X( 1, 12) 0.000000 174.0000 X( 1, 13) 0.000000 208.0000 X( 1, 14) 0.000000 230.0000 X( 1, 15) 0.000000 264.0000 X( 2, 1) 0.000000 0.000000 X( 2, 2) 179.0000 0.000000 X( 2, 3) 321.0000 0.000000 X( 2, 4) 0.000000 0.000000 X( 2, 5) 0.000000 0.000000 X( 2, 6) 0.000000 2.000000 X( 2, 7) 0.000000 9.900000 X( 2, 8) 300.0000 0.000000 X( 2, 9) 0.000000 76.00000 X( 2, 10) 0.000000 120.0000 X( 2, 11) 0.000000 148.0000 X( 2, 12) 0.000000 151.0000X( 2, 14) 0.000000 207.0000 X( 2, 15) 0.000000 241.0000 X( 3, 1) 0.000000 5.000000 X( 3, 2) 0.000000 5.000000 X( 3, 3) 147.2817 0.000000 X( 3, 4) 119.7061 0.000000 X( 3, 5) 69.01229 0.000000 X( 3, 6) 0.000000 2.000000 X( 3, 7) 0.000000 9.900000 X( 3, 8) 0.000000 5.000000 X( 3, 9) 664.0000 0.000000 X( 3, 10) 0.000000 50.00000 X( 3, 11) 0.000000 78.00000 X( 3, 12) 0.000000 81.00000 X( 3, 13) 0.000000 115.0000 X( 3, 14) 0.000000 137.0000 X( 3, 15) 0.000000 171.0000 X( 4, 1) 0.000000 5.000000 X( 4, 2) 0.000000 5.000000 X( 4, 3) 0.000000 5.000000 X( 4, 4) 0.000000 5.000000X( 4, 6) 0.000000 7.000000 X( 4, 7) 0.000000 14.90000 X( 4, 8) 0.000000 5.000000 X( 4, 9) 0.000000 6.000000 X( 4, 10) 0.000000 0.000000 X( 4, 11) 0.000000 13.00000 X( 4, 12) 0.000000 16.00000 X( 4, 13) 0.000000 50.00000 X( 4, 14) 0.000000 72.00000 X( 4, 15) 0.000000 106.0000 X( 5, 1) 0.000000 5.000000 X( 5, 2) 0.000000 5.000000 X( 5, 3) 39.71834 0.000000 X( 5, 4) 97.84982 0.000000 X( 5, 5) 462.4318 0.000000 X( 5, 6) 0.000000 2.000000 X( 5, 7) 0.000000 9.900000 X( 5, 8) 0.000000 5.000000 X( 5, 9) 0.000000 6.000000 X( 5, 10) 262.7391 0.000000 X( 5, 11) 415.0000 0.000000X( 5, 13) 0.000000 50.00000 X( 5, 14) 0.000000 67.00000 X( 5, 15) 0.000000 101.0000 X( 6, 1) 0.000000 10.00000 X( 6, 2) 0.000000 10.00000 X( 6, 3) 0.000000 5.000000 X( 6, 4) 0.000000 5.000000 X( 6, 5) 0.000000 5.000000 X( 6, 6) 0.000000 7.000000 X( 6, 7) 0.000000 14.90000 X( 6, 8) 0.000000 10.00000 X( 6, 9) 0.000000 6.000000 X( 6, 10) 88.26093 0.000000 X( 6, 11) 0.000000 13.00000 X( 6, 12) 86.00000 0.000000 X( 6, 13) 333.0000 0.000000 X( 6, 14) 786.0000 0.000000 X( 6, 15) 0.000000 37.00000 X( 7, 1) 0.000000 9.000000 X( 7, 2) 0.000000 9.000000 X( 7, 3) 0.000000 4.000000X( 7, 5) 0.000000 4.000000X( 7, 6) 0.000000 6.000000X( 7, 7) 0.000000 13.90000X( 7, 8) 0.000000 9.000000X( 7, 9) 0.000000 10.00000X( 7, 10) 0.000000 4.000000X( 7, 11) 0.000000 17.00000X( 7, 12) 0.000000 0.000000X( 7, 13) 0.000000 1.000000X( 7, 14) 0.000000 4.000000X( 7, 15) 0.000000 0.000000 最少总费用Local optimal solution found.Objective value: 1281354.Extended solver steps: 0Total solver iterations: 36模型三的matlab程序A=ones(39,39);for i=1:39;for j=1:39;if i==j A(i,j)=0;else A(i,j)=inf;endendendA(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);endendT=A;m=1;while m<=39for i=1:39for j=1:39if T(i,j)>T(i,m)+T(m,j);T(i,j)=T(i,m)+T(m,j);endendendm=m+1;endT2=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&T(i,j)<=300 T2(i,j)=20;elseif T(i,j)>300&T(i,j)<=350 T2(i,j)=23;elseif T(i,j)>350&T(i,j)<=400 T2(i,j)=26;elseif T(i,j)>400&T(i,j)<=450 T2(i,j)=29;elseif T(i,j)>450&T(i,j)<=500 T2(i,j)=32;elseif T(i,j)>500&T(i,j)<=600 T2(i,j)=37;elseif T(i,j)>600&T(i,j)<=700 T2(i,j)=44;elseif T(i,j)>700&T(i,j)<=800 T2(i,j)=50;elseif T(i,j)>800&T(i,j)<=900 T2(i,j)=55;elseif T(i,j)>900&T(i,j)<=1000 T2(i,j)=60;else T2(i,j)=60+ceil((T(i,j)-1000)/100)*5;endendB=ones(39,39);for i=1:39;for j=1:39;if i==j B(i,j)=0;else B(i,j)=inf;endendendB(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(18,34)=70;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);endendg=1;while g<=39;for i=1:39;for j=1:39;if R(i,j)>R(i,g)+R(g,j);R(i,j)=R(i,g)+R(g,j);endendendg=g+1;endR2=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);endendendf=1;while f<=39for i=1:39for j=1:39if w(i,j)>w(i,f)+w(f,j)w(i,j)=w(i,f)+w(f,j);endendendf=f+1;endQ=ones(39,39);for i=1:7for j=19:39cost(i,j)=w(i,j)endend模型三的matlab结果Columns 19 through 3995.0000 100.0000 105.0000 115.0000 125.0000 140.0000 170.7000160.3000 140.2000 98.6000 38.0000 20.5000 3.1000 21.2000 89.2000 92.0000 122.0000 106.0000 121.2000 128.0000 142.0000145.0000 150.0000 155.0000 165.0000 175.0000 190.0000 215.7000 205.3000 190.2000 171.6000 111.0000 95.5000 86.0000 71.2000 139.2000 142.0000 172.0000 156.0000 171.2000 178.0000 192.000080.0000 90.0000 95.0000 105.0000 115.0000 130.0000 230.7000 220.3000 200.2000 181.6000 121.0000 105.5000 96.0000 86.2000 130.0000 82.0000 112.0000 96.0000 111.2000 118.0000 132.000044.0000 55.0000 60.0000 70.0000 80.0000 95.0000 260.7000 250.3000 235.2000 216.6000 156.0000 140.5000 131.0000 116.2000 110.0000 62.0000 83.0000 61.0000 76.2000 83.0000 97.000037.0000 50.0000 50.0000 65.0000 70.0000 85.0000 255.7000 245.3000 225.2000 206.6000 146.0000 130.5000 121.0000 111.2000 105.0000 57.0000 73.0000 51.0000 71.2000 73.0000 87.000050.0000 37.0000 44.0000 20.0000 14.0000 26.0000 265.7000 255.3000 235.2000 216.6000 156.0000 140.5000 131.0000 121.2000 110.0000 62.0000 67.0000 45.0000 26.2000 11.0000 28.000060.0000 50.0000 55.0000 32.0000 23.0000 4.0000 275.7000265.3000 245.2000 226.6000 166.0000 150.5000 141.0000 131.2000 124.0000 76.0000 78.0000 56.0000 38.2000 26.0000 2.0000模型三的lingo程序model:!钢管购买与运输铺优化设问;sets:cj/1..7/:p,s,t;md/1..21/:y,z,v,D;link(cj,md):cost,x;endsets!这里是数据;data:s=800 800 1000 2000 2000 2000 3000;D=104 301 750 606 194 205 201 680 480 300 220 210 420 500 0 0 0 0 0 0 0;p=168 155 155 160 155 150 160;cost= 95 100 105 115 125 140 170.7 160.3 140.2 98.6 38 20.5 3.1 21.2 89.2 92 122 106 121.2 128 142145 150 155 165 175 190 215.7 205.3 190.2 171.6111 95.5 86 71.2 139.2 142 172 156 171.2 178 19280 90 95 105 115 130 230.7 220.3 200.2 181.6 121 105.5 96 86.2 130 82 112 96 111.2 118 13244 55 60 70 80 95 260.7 250.3 235.2 216.6 156 140.5 131 116.2 110 62 83 61 76.2 83 9737 50 50 65 70 85 255.7 245.3 225.2 206.6 146 130.5 121 111.2 105 57 73 51 71.2 73 8750 37 44 20 14 26 265.7 255.3 235.2 216.6 156 140.5 131 121.2 110 62 67 45 26.2 11 2860 50 55 32 23 4 275.7 265.3 245.2 226.6 166 150.5 141 131.2 124 76 78 56 38.2 26 2;!目标函数;min=@sum(link(i,j):(p(i)+cost(i,j))*x(i,j))+0.05*(@sum(md(j)|j#ge#2#and#j#le#21:y(j)^2+y(j))+ @sum(md(j)|j#le#14:z(j)^2+z(j))+@sum(md(j)|j#eq#9#or#j#eq#11#or#j#eq#17:v(j)^2+v(j))+@ sum(md(j)|j#eq#17#or#j#eq#19#or#j#eq#20:z(j)^2+z(j)));@for(cj(i):@sum(md(j):x(i,j))>=500*t(i));@for(cj(i):@sum(md(j):x(i,j))<=s(i)*t(i));@for(cj(i):@bin(t(i)));@for(md(j)|j#ne#9#and#j#ne#11#and#j#ne#17:v(j)=0); @for(md(j):@sum(cj(i):x(i,j))=y(j)+z(j)+v(j));@for(md(j)|j#lt#15:D(j)=y(j)+z(j+1));v(9)+y(16)=42;v(11)+v(17)=10;y(17)+z(18)=130;z(17)+y(19)=190;z(19)+y(20)=260;z(20)+y(21)=100;@gin(@sum(link(i,j):x(i,j)));endWelcome To Download !!!欢迎您的下载，资料仅供参考！。
2000年数学建模B题钢管订购和运输
关于下面3个问题（可以是其中某个小问题），试分别建立模型。
包括给出问题分析和建模思路、模型假设、变量说明、模型建立。
不需要求解。
1 B 题 钢管订购和运输要铺设一条1521A A A →→→ 的输送天然气的主管道, 如图一所示(见反面)。
经筛选后可以生产这种主管道钢管的钢厂有721,,S S S 。
图中粗线表示铁路，单细线表示公路，双细线表示要铺设的管道(假设沿管道或者原来有公路，或者建有施工公路)，圆圈表示火车站，每段铁路、公路和管道旁的阿拉伯数字表示里程(单位km)。
为方便计，1km 主管道钢管称为1单位钢管。
一个钢厂如果承担制造这种钢管，至少需要生产500个单位。
钢厂i S 在指定期限内能生产该钢管的最大数量为i s 个单位，钢管出厂销价1单位钢管为i p 万元，如下表：1单位钢管的铁路运价如下表：1000km 以上每增加1至100km 运价增加5公路运输费用为1单位钢管每公里0.1万元（不足整公里部分按整公里计算）。
钢管可由铁路、公路运往铺设地点（不只是运到点1521,,,A A A ，而是管道全线）。
（1）请制定一个主管道钢管的订购和运输计划，使总费用最小（给出总费用)。
（2）请就（1）的模型分析：哪个钢厂钢管的销价的变化对购运计划和总费用影响最大，哪个钢厂钢管的产量的上限的变化对购运计划和总费用的影响最大，并给出相应的数字结果。
（3）如果要铺设的管道不是一条线，而是一个树形图，铁路、公路和管道构成网络，请就这种更一般的情形给出一种解决办法，并对图二按（1）的要求给出模型和结果。
77问题分析问题一，首先，所有钢管必须运到天然气主管道铺设路线上的节点1521A A A →→→ ，然后才能向左或右铺设。
必须求出每个钢管厂721,,S S S 到每个节点1521A A A →→→ 的每单位钢管的最小运输费用。
问题二，通过问题一里面Lingo 编程运行得出的结果，分析哪个钢厂钢管的销价的变化对购运计划和总费用影响最大，哪个钢厂钢管的产量的上限的变化对购运计划和总费用的影响最大。
钢管订购与运输的优化模型
钢管订购与运输的优化模型钢管订购与运输是现代经济中的一个重要问题。
钢管是建筑、制造、输送等多个领域必不可少的材料，一般情况下我们需要从厂家或供应商那里订购所需钢管，并通过运输将其送到指定地点。
在订购和运输的过程中，我们需要考虑许多因素，如运输距离、交通方式、需求量、时间限制、价格等等。
针对这些问题，我们需要使用优化模型来提高订购和运输的效率和经济性。
一、钢管需求模型在实际工作中，我们需要尽可能准确地了解钢管的需求情况。
这样才能更好地制定订购和运输计划。
钢管需求模型是一个重要的决策工具，可以帮助我们进行有效的决策。
其主要内容如下：（一）需求量需求量是指市场中对钢管的总需求量。
建立需求量模型需要考虑市场状况、产品质量、价格、季节等因素。
我们可以通过市场调研、历史销售数据等途径进行预测。
需求结构是指不同规格的钢管在市场中的占比情况。
了解需求结构可以帮助我们更好地制定订购方案，避免过度订购或订购不足情况的发生。
需求时间是指市场上对钢管的需求时间分布。
了解需求时间可以帮助我们更好地制定订购和运输计划，减少废弃和过度库存，提高物流效率。
在了解钢管需求的情况之后，我们需要制定订购计划。
为了提高采购效率，我们需要采用优化模型。
该模型的主要内容包括：订购量是指在一定时间内企业需要订购的钢管数量。
订购量的大小直接影响企业的成本和库存水平。
因此，我们需要根据实际需求，结合采购成本、库存水平等因素进行考虑，制定出合理的订购量。
（二）订购频率订购频率是指企业在一段时间内订购钢管的次数。
频繁而杂乱的订购计划会耗费大量的人力、物力和财力，同时也增加了库存和物流的费用。
因此，我们需要根据实际情况，制定出合适的订购频率。
（三）订购价格订购价格是采购者与供应商之间协商的价格。
采购者需要确保订购价格与采购成本相符，同时也要考虑到供应商的利润。
因此，合理的订购价格既要考虑到采购方的利益，也要考虑到供应商的利益。
订购钢管需要通过运输将其送到指定地点。
钢管订购与运输的优化模型
论文题目：钢管订购与运输的优化模型27队队长：杨璐学号：******** 专业：信计队员：高春妮学号：******** 专业：数应队员：贺瑞瑞学号：******** 专业：计科2012 年 07 月 19日钢管订购与运输的优化模型摘 要 本文讨论了在铺设天然气管道的过程中如何合理订购与运输钢管以使总费用最小的优化问题。
问题一是在一定约束条件下以钢管订购和运输的总费用为目标函数的非线性优化问题。
总费用由订购钢管的总费用、从钢厂到站点运输钢管的总费用及从站点开始铺设钢管的总费用三部分组成。
订购钢管的总费用和从钢厂到站点运输钢管的总费用分别通过在各厂购买量与各厂出厂销价和在各厂购买量与从各钢厂到各站点运输单位钢管的最小费用的线性运算得到。
从站点开始铺设钢管的总费用通过等差数列求和得到。
在求从钢厂到站点运输钢管的总费用时，关键是采用弗洛伊德算法，用MATLAB 软件编程求出单位钢管从各钢厂运往各站点最小运输费用（见表3）。
约束条件可由题目相应已知条件给出，故可建立钢管订购与运输的优化模型一。
利用LINGO 软件编程求解出此模型，得到钢管订购和运输的最小总费用为1280837万元，并经整理分析给出钢管订购与运输方案（见表4）和钢管铺设方案（见表5）。
问题二是对问题一中模型的灵敏度分析，通过控制变量的方法即每次只让一家钢厂的销价或生产上线发生变化并且每次的变化是相同的，分别得出：6S 钢厂钢管的销价变化对总费用影响最大，5S 和6S 钢厂钢管的销价变化对购运计划影响最大；1S 钢厂钢管的产量上限变化对总费用影响最大，3S 钢厂钢管的产量上限变化对购运计划影响最大。
问题三是对问题一的推广，站点向左右两边铺设变为向三个方向铺设，在问题一中模型基础上增加一些支路变量和约束条件，同时在目标函数中增加相应的铺设费用，由此建立了钢管订购与运输的优化模型二。
利用LINGO 软件编程求解得到钢管订购和运输的最小总费用为1414800万元，并给出钢管订购与运输方案（见表11）和钢管铺设方案（见表12）。
钢管的订购及运输优化方案
钢管的订购及运输优化方案一、钢管的订购1.1 需求分析在进行钢管订购之前，首先需要进行需求分析。
钢管的形状、尺寸、材质、用途、数量等方面的需求都需要进行仔细的分析，以确保最终的采购结果符合要求。
1.2 供应商选择在确定钢管的需求后，需要寻找合适的供应商进行订购。
选择供应商要综合考虑各方面因素，包括价格、质量、交期、可靠性等方面。
可以通过询价、对比、调查等手段来寻找合适的供应商。
1.3 合同签订在选择好供应商之后，需要签订合同。
合同中要明确钢管数量、型号、质量标准、交货期限、运输方式、付款方式等内容，规范供应商的行为，确保合同执行的顺利进行。
二、钢管的运输2.1 运输方式选择根据钢管的数量、尺寸、重量、运输距离、要求到达时间等因素，选择合适的运输方式。
常见的运输方式包括铁路、公路、水路、空运等。
2.2 包装方式选择钢管在运输过程中需要进行包装，以保证其不受损坏。
包装方式应根据钢管的特点进行选择，常用的包装方式包括裸装、编织袋、塑料薄膜、木箱等。
2.3 运输路线优化在确定运输方式和包装方式之后，应针对具体的运输路线进行优化。
优化的原则包括缩短运输时间、降低运输成本、提高运输效率等方面。
2.4 运输管理在钢管运输过程中，需要进行运输管理。
管理内容包括钢管的装车、卸货、运输途中的安全监管等方面。
同时，应建立健全的运输记录管理体系，确保运输全程可追溯。
三、钢管订购及运输优化方案为了更好的优化钢管的订购及运输过程，可采取以下措施：3.1 制定钢管需求分析标准建立钢管需求分析的标准化体系，规范钢管订购的流程和细节。
该标准应涵盖钢管的形状、材质、性能、用途、数量等方面，确保符合实际需求。
3.2 建立供应商评价体系通过建立供应商评价体系，以价格、品质、信誉、交期等为考核指标，对供应商进行评价和排名，选用优质和稳定的供应商，确保采购的钢管质量和交货期的稳定。
3.3 采用智能供应链管理系统建立智能供应链管理系统，通过物流信息技术支持物流实时监控、自动化分配、预警预测、异常处理等功能，实现钢管订购及运输全流程的可视化和管理。
钢材采购及运输节点流程
钢材采购及运输节点流程1. 采购流程1.1 采购需求确认- 由采购部门根据相关项目需求，确认钢材的采购需求。
- 确定所需钢材的规格、种类、数量等详细信息。
1.2 供应商选择- 采购部门根据公司设定的供应商管理规则，选择合适的供应商。
- 对供应商的信誉度、交货期、价格等进行综合考量。
1.3 询价与比较- 采购部门向选定的供应商发送询价单，要求供应商提供钢材的报价、交货期等信息。
- 根据供应商的报价和其他相关因素进行比较，选择最合适的供应商。
1.4 合同签订- 采购部门与选定的供应商进行合同谈判，明确双方的权责义务和交易条件。
- 确定具体的采购合同，包括钢材规格、价格、数量、交货期等内容。
1.5 采购订单发布- 采购部门根据合同内容生成采购订单，并将之发送给供应商。
- 采购订单包括采购合同的详细信息以及采购部门的收货地址、联系人等信息。
1.6 供应商确认与安排生产- 供应商收到采购订单后，进行订单确认，并进行生产安排。
- 进行生产时，供应商需严格按照合同约定的规格和质量标准进行生产。
1.7 发货准备- 供应商在完成钢材生产后，进行发货准备工作。
- 检查钢材的质量并进行包装，确保货物在运输过程中不受损。
1.8 收货与验收- 采购部门根据合同约定的交货期，安排接收钢材货物。
- 对收到的钢材进行验收，检查质量和数量是否符合合同要求。
1.9 支付结算- 采购部门根据验收结果和合同约定的支付方式，进行支付结算。
- 将货款支付给供应商，并完成相关的财务记录。
2. 运输流程2.1 运输计划制定- 采购部门根据采购订单、交货期等信息制定运输计划。
- 确定运输的起点、终点、路线、运输方式等，并预估运输时间。
2.2 货物装载- 采购部门与供应商协商确定装载时间和地点。
- 货物按照规定的包装要求进行装载，并制定装载清单。
2.3 运输过程监控- 采购部门与物流公司保持沟通，监控运输过程。
- 确保货物按照计划按时到达目的地。
钢管订购与运输论文
钢管订购和运输优化模型关键词：最优化，Floyd 算法；钢管订购与运输；lingo ；摘要本文研究钢管订购与运输问题，该问题给定钢厂生产能力、不同路径长度下运输价格等限制，需要制定使总费用最小的订购运输计划，并分析在生产能力、单位订购价格变化时，对订购及运输计划的影响。
对问题一，我们先将整个网络图简化，利用Floyd 算法通过MATLAB 编程计算出铁路网上，任意两点的最短距离，从而可以知道从每个钢厂S 运往铁路网上各点个单位运输价格，结合公路运输路线，计算出钢管从钢厂i 运往节点j 运输单价ij c 。
设运往节点j ，刚才分别向左右两个方向运输的钢管量分别为j y 、j z ，那么铺设管道过程中的运输费也很容易表示。
我们就可以建立最小费用的优化模型。
再利用lingo 进行求解就可以得出最终结果。
最小总费用为：1278632万元，最终的订购计划为： S1 S2 S3 S4 S5 S6 S7 订购量800800100013211251问题二通过对模型参数的灵敏度分析，我们能得出模型的实际合理性，其中钢厂的产量影响通过问题一的LINGO 求解分析得出答案，而销价对结果的影响要通过自行改变钢厂的销价来实现。
分析结果显示：S1的产量变化对总费用的影响是最大的；S6钢厂销价的变化对总费用的影响是最大的，S4钢厂的销价变化对总费用的影响是最小的。
问题三，相对于问题一管道的铺设路径稍微复杂了一点，我们在问题一的基础上增设管道向第三方运输的运量j m ，再修改问题一的模型，得到本问题模型，再利用lingo 求解，最小总费用为：1404949万元，最终订购计划为： S1 S2 S3 S4 S5S6S7订购量80080010001475.086 1827.914 0最后，我们对模型进行了评价和改进，使模型更加完整和优化。
AbstractIn this paper, we study on the problem of pipe Ordering and Transporting,which Limit steel production capacity and transportation unit price under a different path length, and develop a plan of Pipe ordering and shipping making the minimumcost ,then analyze the production capacity of the unit subscription price changes, see which impacts of the ordering and transportation plans.The first problem, we first simplify the entire network diagram the use MATLAB programs and Floyd algorithm to calculate the shortest distance rail network, between any two points, which can be shipped to know from each steel S-line rail transportation unit price points, combined with road transport routes, calculate the transport unit from steel mills shipped node. Let destined for a node, the amount of steel were just about transport in both directions, respectively, then the process of laying the pipeline transportation costs can easily expressed. We can create the minimum cost optimization model. Reuse lingo for solving the final result can be drawn. The minimum total cost: 12,786,320,000 yuan, the final subscription plan is: S1 S2 S3 S4 S5 S6 S7 subscribe 800 800 1000 0 1321 1251 0 Second problem by sensitivity analysis of the model parameters, we can drawpractical rationality of the model, which yield a steel through problem solving LINGO analysis of the answer, but the selling price impact on the outcome of the change through their own mills the selling price to achieve. Analysis showed: the impact on the total production cost changes S1 is the largest; affect S6 steel selling price changes on the total cost is greatest impact on the selling price changes S4 mills total cost is minimal.Question three, relative to the laying of a pipeline path problem is slightly morecomplicated, we added the issue on the basis of a third-party pipeline to transport the volume, and then modify the problem a model to get the problem model, re-use lingo solved The minimum total cost: 14,049,490,000 yuan, the final subscription plan is:S1S2 S3 S4 S5S6S7subscribe 80080010001475.086 1827.914 0Finally, the model was evaluated and improved to make the model more complete and optimized.一、问题重述1.1问题提出要铺设一条1521A A A →→→ 的输送天然气的主管道,如题中图一所示，图中包括可以生产这种主管道钢管的钢厂有721,,S S S ，钢厂生产的钢管需要通过铁路运输网、公路运输网运到管道铺设处。
钢管订购和运输计划
钢管订购和运输计划一、引言本文档旨在详细描述钢管订购和运输计划的各个方面，包括订购过程、运输方式、时间安排等内容。
钢管作为建筑、工程和制造业的重要材料之一，对于项目的顺利进行具有重要意义。
因此，钢管的订购和运输需得到合理安排和重视。
二、钢管订购2.1 计算需求量在进行钢管订购之前，首先需要计算所需的钢管数量。
这一计算通常由项目负责人、工程师或建筑师来完成。
计算需求量时，需要考虑以下因素：•项目规模和要求•钢管的类型和规格•使用钢管的位置和用途2.2 选择供应商选择合适的供应商是钢管订购过程中的关键步骤。
在选择供应商时，应考虑以下几个方面：•供应商的信誉和声誉•产品质量和性能•价格和交货时间2.3 发出订单一旦选择了合适的供应商，就需要发出订单。
订单应包括以下信息：•钢管的规格和数量•交货日期和地点•付款方式和条款•其他特殊要求三、钢管运输3.1 运输方式钢管的运输方式多种多样，常见的有以下几种：•公路运输：适合短程或小批量运输，成本较低。
•铁路运输：适合远距离和大批量运输，安全可靠。
•水运：适合长距离和大宗运输，成本相对较低。
•空运：适合迫切需要和紧急情况下的运输，费用较高。
3.2 运输安排在确定运输方式后，需要进行具体的运输安排。
主要包括以下几个方面：•运输时间表：明确每次运输的时间，确保与工程进度相匹配。
•运输车辆或船舶：根据货物的规模和距离选择合适的运输工具。
•路线规划：选择最优的运输路线，考虑效率和安全性。
3.3 运输风险和控制在钢管运输过程中，存在着一定的风险，如交通事故、货物丢失或损坏等。
为了减少这些风险，可以采取以下措施：•选择可靠的运输公司或车队，避免使用低质量的运输工具。
•对货物进行包装和固定，确保在运输过程中不会受到损坏。
•跟踪和监控货物的运输情况，及时处理可能出现的问题。
四、总结本文档详细介绍了钢管订购和运输计划的各个方面。
钢管作为重要的建筑材料，其订购和运输对于项目的进展具有重要意义。
(交通运输)钢管订购和运输计划精编
（交通运输）钢管订购和运输计划（交通运输）钢管订购和运输计划钢管的订购和运输计划摘要在钢管的订购和运输计划中，在第壹问中用最短路算法，求解出每个钢厂到站点的最小费用（包括运输费和出厂销售价），考虑到在铺设时管道要沿铺设路线离散地卸货，即运货到A j后，仍要在铺设路线上运输，因为不足整公里部分要按照整公里计算，所以我们认为沿管道路线每铺设1公里就要卸下1单位钢管，因此从某点A j向左铺设或向右铺设y时，此段运费应为：点A j向右铺设z j，从A j+1向左铺设y j+1，为了保证合拢，则z j+y j+1=a j，在这些条件之下，利用软件，求解出总费用最小。
分析模型的销售价灵敏度的时候，将各个钢厂单位钢管的销售价分别增加和减少若干万元，再用求解第壹问题的模型，见总费用的变化大小，变化大的就是影响结果比较大的；用同样的方法能够分析生产上限的灵敏度。
第三问得时候，我们利用求解第壹问的方式来求解问题。
关键字：最短路算法，，分别改变同样的条件来对比壹，问题重述（略）二，符号说明：a ij站点A j至A j+1的里程（铺设管道需要的钢管量）s i s i钢厂的最大生产量x ij从钢厂s i到A j的钢管数量c ij从钢厂s i运往A j的单位钢材费用最短路，即亮点运输单位钢材所需的最少费用，包括运输费和出厂销价y j A j点往左铺设的钢管数量zjA j点往右铺设钢管的数量f总费用三，问题分析：(1)对问题壹的分析：从钢厂s i向点A j运输钢管时，为了降低费用，应该走费用最小的路径，从壹个工厂s i到壹个点A j的路线且不唯壹，需要从中找出费用最短的路，相应的最小费用为c ij，包括运输费和销售费。
从图我们能够见到，七个钢材厂要到A1这点必须要经过A2，所以在考虑最低费用路径的时候，能够把A1和A2见做壹个点来考虑，。
根据图，我们由最短路问题的算法。
例：从s1到最短的铁路为：2902km，根据1单位钢管的铁路运价表，可知铁路花费为：60+5*20=160万元，公路运费为3*0.1=0.3万元，且且s1钢厂出厂1单位刚窜为160万元，所以，总费用=铁路运费+公路运费+销售价即=320.3（万元）;用同样的方法，我们能够得到A j的最小费用（单位：万元）：A2 A3 A4 A5 A6 A7 A8 A9A10A11A12A13A14A15S 1 320.3300.2258.6198180.5163.1181.2224.2252256266281.2288302S 2 360.3345.2326.6266250.5241226.2269.2297301311326.2333347S37353327262524202324252627283 5.3 5.2 6.6 6 0.5 1 1.2 3.2 7 1 1 6.2 3 7S 4 410.3395.2376.6316300.5291276.2244.2222211221236.2243257S 5 400.3380.2361.6301285.5276266.2234.2212188206226.2228242S 6 405.3385.2366.6306290.5281271.2234.2212201195176.2161178S 7 425.3405.2386.6326310.5301291.2259.2236226216198.2186162在铺设时管道要沿铺设路线离散地卸货，即运货到A j后，仍要在铺设路线上运输，因为不足整公里部分要按照整公里计算，所以我们认为沿管道路线每铺设1公里就要卸下1单位钢管，因此从某点A j向左铺设或向右铺设y时，此段运费应为：设从点A j向右铺设z j，从A j+1向左铺设y j+1，为了保证合拢，则z j+y j+1=a j，j=1,2…15.问题的实质是确定从钢厂向运输钢管的数量，以及从A j向左，右铺设的里程(km)数，使总费用最小。
钢管的订购和运输(模型)
图2的网络 , 请就这种更一般的情形给出一种解决方法,并对图2按(1)的要求
给出模型和结果.
图1
图2
最短路线问题算例
2-6-10
6-9-10
600
300
300
2
6
9-10
100
1-4-10 650
1
4-6-12,7
yij 0 i 1,2,7; j 1,2,,14,15
s4 410.3 395.2 376.6 316 300.5 291 276.2 244.2 222 211 221 236.2 243 257
s5 400.3 380.2 361.6 301 285.5 276 266.2 234.2 212 188 206 226.2 228 242
s6 400.3 385.2 366.6 306 290.5 278.1 266.2 234.2 211 197 187 166.2 161 178
cij : 单位钢管从Si 运到Aj所需的最小购运费用
xi : Si是否承担制造任务
yij : Si 运抵Aj的钢管数量(不含路过的部分)
zj
:
运到Aj沿Aj
Aj
铺
1
设
的
量
： 单位钢管每公里路程的运费
7 15
14
min
i 1
cij yij
j 1
2
j 1
z j (z j
1) (w j, j1 z j )(w j, j1 z j
公路运输费用为1单位钢管每公里0.1万元(不足整公里部分按1公里计).
钢管可由铁路公路运往铺设地点不只是运到各Aj而是运到管道全线
钢管订购和运输的规划模型
钢管订购和运输的规划模型陈丹妮摘要：本文就天然气管道钢管的订购和运输问题，建立了使订购和运输总费用最小的优化模型.我们把计算分为订货和铁路，公路费用的计算及管道上运输费用的计算两个部分.对第一部分的计算，我们采用了增减约束条件的方法，避免了求解一组多分支规划的繁重的计算.对第二部分的计算，我们综合各种可能情况作出比较，从而使计算简化，并求出了最优的钢管订购和运输计划.对于第二问，我们把每个钢厂的销价及生产上限在一定范围内浮动，观察比较得出钢厂3S 钢管的销价的变化对购运计划和总费用影响最大，钢厂1S 钢管的产量的上限的变化对购运计划和总费用的影响最大.在第三问中运用第一问的方法建立模型，同样求出了铁路，公路和管道构成网络时总费用最小的钢管订购和运输计划.一 题的重述要铺设一条1521...A A A →→→的输送天然气的主管道.经筛选后可以生产这种主管道钢管的钢厂有721,...,,S S S .连接钢厂i S （i=1,…,15）和)15,...,1(=j A j 的有铁路和公路.沿管道或者原来有公路，或者建有施工公路.一个钢厂如果承担制造这种钢管，至少需要生产500个单位.已知钢厂i S 在指定期限内能生产该钢管的最大数量，钢管出厂销价及1单位钢管的铁路运价和公路运输费.钢管不只是运到点,,...,,1521A A A 而是管道全线.问如何制定一个主管道钢管的订购和运输计划，使总费用最小；哪个钢厂的销价变化对购运计划和总费用影响最大，哪个钢厂钢管的产量的上限的变化对可以计划和总费用的影响最大；如果要铺设的管道不是一条线，而是一个树形图，铁路，公路和管道构成网络，如何建立相应的模型和如何求解.7二 本假设与符号约定1) 1km 主管道钢管称为1单位钢管； 2) 假设在钢厂i S 的订购货量为i x 个单位；3) 对于图一，铁路和公路相交的车站从左到右分别记为1721,...,,t t t ；4) 对于图二，铁路和公路相交的车站或者铁路和管道相交的车站从左到右分别记为1821,...,t t t ；5) 假设钢厂i S 流经j t 站的钢管量为j i x ,个单位； 6) 假设j A 处的到货量为j a ；7) 假设1单位钢管从钢厂i S 运到j A 的运价为j i k ,；8) 钢厂i S 在指定期限内能生产该钢管的最大数量是i S 个单位； 9) 钢管出厂销价1单位钢管为i P 万元；10) 假设铁路运费是整段计算的（从货物上车到下车一次性收费），二不是分段计算； 11) 沿管道公路的运费计算与其他公路一致，且不考虑流量限制的问题.三 问题的分析从图上可以看出，各钢厂订购的钢管必先经铁路或公路运往主管道与公路的各节点iA 上再沿主管道进行运输和铺设.因此，我们可以把运输的总费用分为在非管道（铁路或公路）7上的运输费用和主管道上的运输费用两部分来计算.对于非管道上的运输.由于钢厂承担制造钢管后至少生产500个单位，所以对于每一个钢厂来说，订购量要么为0，要么就大于或等于500个单位，这就构成一组62个的多分支线性规划问题，计算将非常复杂.但我们可以采用如下办法简化计算：对所有钢厂的产量先不设下限进行求解，若解出来的订购量都符合不小于500个单位的情况则为可行解，若解出来的订购量中有不为0的，但小于500个单位，则在约束条件中加进这个订购量的下限进行求解，直至得出符合条件的最优解.对于管道上的钢管运输铺设的费用则比较复杂，钢管从一个i A 点出发，可以单纯沿管道公路进行运输，也可以一边运输一边铺设，要使运输费最优是类似一次规划的非线性规划问题，由于变量多，计算量大，因此要进行一定的简化.我们现证明一重要结论：当管道上各节点的钢管量等于与节点相连接的两边管道总长度的一半时，管道上钢管的运输费最小.设运价为y ，运量为x,y 是x 的函数，并且有1.0-==-dxdyk （其中路程单位为km ）.的钢管从其中一端点出发，y-x 的关系如图所示： ys x < s x >运费g 即是图中阴影部分的面积.当x<s 时，])([2122x s x k dy y gso-+==⎰，当x=s 时，221)(kS dx s x k dy y g s s =-==⎰⎰， 当x>s 时，)2(210s x ks dy y g s-==⎰， 容易看出，当x>s 时，对g 来求导有：)24(2)]1)((22[2's x kx s x k g -=--+=，推出2s x =为稳定点.在[0,s]区间上，,41)2(,21)(,21)0(222ks s g ks s g ks g ===所以当=x 2s时，费用是最小的，由此方法我们计算出管线上的最小运输费t=61593.275万元.四 模型的建立和求解1,通过上面的分析，我们首先先令各钢厂订购的钢管运往各节点的铁（公）路运费和订购费最优，然后我们把各钢厂订购的钢管分成17份分别运往与公路相连的火车站.由于铁路轨道成树状分布，所以这样的最优路线是确定的.通过对图一的分析，我们发现，15141413111098,,,A A A A A A A A ---- 这四段管道路有这样的情况：1单位钢管从这些管道路之一运过的运费，比从连结该管道路两端点的最短的公（铁）路线运过的运费要高.也就是说，与其将钢管经过这些管道路运输，不如发生“倒运”.因此，这些管道路左右两边的钢管存货应该要满足两边管道铺设的需要，而不应该经这四段管道路进行货物调送.根据前面的假设，我们列出如下以铁（公）路运费和订购费为目标函数的线性规划：∑∑∑===+=17171,,71min j i i i j i ji i x p x kfs.t)1)....(7,...,2,1(171,==∑=i x xi j ji )2....(0,,≥j i i x x∑==712,)3....(i ji a x∑==7132,)4....(i i a x∑==7143,)5....(i i a x∑==7154,)6....(i i a x∑==7165,)7....(i i a x∑==+7177,6,)8....()(i i i a x x∑==7188,)9....(i i a x∑==7199,)10....(i i a x∑==711010,)11....(i i a x∑==711111,)12....(i i a x∑==711212,)13....(i i a x∑==711313,)14....(i i a x∑==+711415,14,)15....()(i i i a x x∑==+711517,16,)16....()(i i i a x x)17....(236182≥∑=i ia)18....(2130159≥∑=i ia)19....(3521102≥∑=i ia)20....(13501511≥∑=i ia)21....(4251132≥∑=i i a )22....(5001514≥∑=i ia)23....(4671142≥∑=i ia)24....(517171≥∑=i ix)25....(i i s x ≤由于只对i x 作非负限制时，计算出7x 低于下界500，所以需另加约束条件)26....(5007≥x重新求解得：f=1015556,500 ,740 ,1331 ,0 ,1000 ,800 ,8007654321=======x x x x x x x .这样，我们得到各节点的钢管量，然后一边运输一边铺设这些钢管，求出所需运费为 q=366409.05万元，所以这样的运输方案得到的总费用为m=1381965.05万元.对于这个方案我们还要进行调整.由上面的讨论我们知道，当管道上各节点的钢管量刚好等于与节点相连接的两边管道总长度的一半时，在管道上的运输费用最小.我们把（151413121110987654321,,,,,,,,,,,,,,a a a a a a a a a a a a a a a ）＝( 0, 254.5, 525.5, 678, 400, 199.5, 203, 440.5, 580, 390, 260, 215, 315, 460,250 ) 作为约束条件加进上述规划中，解得m ’=13066563万元，可见这样的运输方案更优.我们可以再考虑把各钢厂的钢管运到各节点后，再通过运输调整到运输最小时分钢量分布.调整的运输费用为∑∑===151,,152j j i ji i a hz .其中j i h ,是1单位钢管从i A 到j A 沿管线运输的价格，j i a ,是从i A 到j A 沿管线的钢管运输量. 因此我们又得到如下规划t z f F ++=min 其中，∑∑∑=-=+=17171,,71j i i i j i ji i x p x kfs.t （加上一规划约束条件中的（1）至（16）及（25））∑==1511,52i i a∑==1512,5.202i i a∑==1513,5.525i i a∑==1514,678i i a∑==1515,400i i a∑==1516,5.199i i a∑==1517,203i i a∑==1518,5.440i i a∑==1519,580i i a∑==15110,390i i a∑==15111,260i i a∑==15112,215i i a∑==15113,315i i a∑==15114,460i i a∑==15115,250i i a用Maple 软件解得：F=1203697.575,500 ,890 ,1181 ,0 ,1000 ,800 ,8007654321=======x x x x x x x经过比较，我们认为这个订购和运输的方案是最优的.由此可得详细的订运方案如下： （1）7个钢厂的订购量分别为500,890,1181,0,1000,800,8007654321=======x x x x x x x .（2）钢厂1S 的钢管分3批运输，第一批197个单位运往5A ，第二批400个单位运到6A ，第三批203个单位直接沿公路运到7A .钢厂2S 的钢管分两批运输，第一批359.5个单位运往4A ，另一批经8t 运到8A .钢厂5S 的钢管也分两批运输，第一批420个单位经, ,78t t 456 , ,t t t运到4A ，另一批580个单位经9t 运往9A .钢厂5S 的钢管分六批运输，第一批45个单位运往3A ,第二批166个单位运往4A ，第三批61个单位运往5A ，第四批199.5个单位经6t 运到7A ，第五批390单位运到10A ，第六批260个单位运往11A .钢厂6S 的钢管分三批运输，第一批115个单位往12A ，第二批315个单位往13A ，第三批420个单位直接沿公路运往14A .钢厂7S 的钢管全部直接运到15A .2，根据我们建立的模型，保持其它条件不变，令各个钢厂的钢管销价i p 上浮或下降5％，可得到总费用的变化幅度和运购计划的变化情况，如下表：从上表比较可得，钢厂i S 的钢管销价的变化对总费用及购运计划影响最大.用同样的方法，保持其他条件不变，令各个钢厂钢管产量的上限上浮或下降10%，得同样由上表可得出，钢厂1S 的钢管产量的上限的变化对总费用及购运计划的影响最大. 3、如图二，要铺设的管道不是一条线，而是一个树形图，铁路、公路和管道构成网络.管道运输最优时各个i A 的存钢量为()()50,180,225,65,430,622,250,460,315,215,0,390,0,5.440,203,5.199,400,678,5.525,5.254,0,,,,,,,,,,,,,,,,,,,,212019181716151413121110987654321==a a a a a a a a a a a a a a a a a a a a a L对于212019,,A A A 这三个兼为火（汽）车中转站的点，我们把它们一分为二看待.以19A 为例，一方面看成1219,A t 的货物由此经过，一方面看成19A ，其钢管量为与之相连接的两段管道总长度的一半，并且钢管直接从它运到主管道.根据第一问的做法，先把i A 的钢管量预置成L 的数量值，这样沿主管道的运输费用就能降到最低，在此基础上对各钢厂的定货量及其分流方式进行调配.然后，使用第一问的方法列出线性规划如下：()∑∑∑∑====≥==+=211,,7117171,,0,7,...,2,1..min j j i i i j i i j i ii j i j i x x i x x t s x p x k f3712,7121,a x a xi i i i ==∑∑== ∑∑====7154,7143,i i i i a x a x∑∑===+=7177,6,7165,)(i i i i i a x x a x∑∑====71169,7188,i i i i a x a x∑∑====711711,711010,i i i i a x a x∑∑====711213,711812,i i i i a x a x∑∑===+=711518,17,711314,)(i i i i i a x x a x∑∑====+711919,711416,15,)(i i i i i a x a x x∑∑====712121,712020,i i i i a x a xi i i is x x≤=∑= 715903初次求解结果，50007<<x ，因此我们加入约束条件5007≥x ，再次求解：500,1260,1543,0,1000,800,800,14503577654321========x x x x x x x f最小费用 475.1518014=+=e f m （万元），其中e 是当节点上的钢管量取自中的数值时，仿照问题一中的计算方法所得出的管道上的运输费用，475.67647=e 万元.五 模型的优缺点分析及其推广我们建立的模型具有较强的可行性和可操作性，并且具有相当的实际意义.虽然我们未能对多个分支规划组逐个进行求解从而得出最优解，但我们对模型进行了适当的近似简化处理，减少了计算量和计算难度，最后得出可行解.我们建立模型的方法和思想对其它类似题材也适用，在建筑运输方面适用性较强，并可以推广到社会生活中相关的多个领域中去.对于类似的问题，对模型的决策性因素加以具体对照分析即可.参考文献：[1]程里文，吴江，张玉林，运筹学模型与方法教程，清华大学出版社，北京，2000[2]L.库珀，U.N勃哈特，L.J勒布朗（美），运筹学模型概论，上海科学技术出版社，上海，1987[3]刘宝碇，赵瑞清，随机规划与模糊规划，清华大学出版社，北京，1998[4]李世奇，杜慧琴，Maple计算机代数系统应用及程序设计，重庆大学出版社，1999.（文章编辑：黄绮玲\颜学友）接83页。
钢管订购和运输的数学模型
１钢管订购和运输的数学模型摘要： 本文先对钢管订购和运输问题做了深入的分析，通过对问题的简化和等价转换,将问题归结为一非线性规划模型，利用软件LINGO 和LINDO 对问题1和3都作出最优解（分别为：127.966亿元与140.5170 亿元），在解1时给出简化模型和算法（解为：130.057亿元）。
在解问题3时，充分考虑了网络的特性,简化了算法。
由于本题是铁路,公路混合网,本文提出等价转换方法将之变为纯公路网,运用固有最段路径算法简化了计算过程.1 问题的提出计划铺设一条输送天然气的主管道，已知有五个生产主管道钢管的钢厂和铁路，公路混合的交通运输网。
试根据钢厂的位置距离销价生产能力以及运输费用等情况制定钢管订购和运输的最佳方案，使总费用最少。
已知运输网（略）及以下数据：（注：为方便计，1km 主管道钢管称为1单位钢管）钢厂i S 在指定期限内能生产该钢管的最大数量为i s 个单位，钢管出厂销价1单位钢管为i p 万元，如下表：1单位钢管的铁路运价如下表：公路运输费用为1单位钢管每公里0.1万元（不足整公里部分按整公里计算）。
2 问题假设H1：运输通路畅通无阻，即任一厂的钢管可到达网络上任一点。
H2：运输费用只与里程和所经过的线路有关，不考虑在铁路和公路之间转换时所增加的额外费用。
H3：总费用只包括运输费和所用钢管的总价格。
H4：钢管必须运到管道全线，设堆放点之间最小距离为1km。
H5：公路运输时不考虑空车来回开的费用3符号说明A i 节点iB i 铁路公路交点C ji厂j到节点i的单位费用(包括销价)S j厂j2*r i运到i点的钢量r i - w i ，r i+w i , r1i，r2i节点i向两侧铺运的距离，r3i为往第三方铺运的距离x ij厂j向i点供的钢量D 问题1中的管道总长d j j厂的销价v i 节点A i与节点A i+1之间的距离4问题分析4．1 问题1的分析1.将运输费用分成两部分：①在管道通路上的运费f1,简称铺运费，这种运输方式称铺运。
钢管订购和运输计划
钢管订购和运输计划钢管是一种广泛应用的建筑材料，其在建筑、制造和其他行业中都有大量的应用需求。
然而，钢管的订购和运输计划并不是一项简单的任务，必须考虑到许多因素，如尺寸、数量、质量、交货时间、成本等因素。
在本文中，将讨论钢管订购和运输计划的关键方面，以便您能够更好地制定计划，满足您的需求。
订购钢管的步骤首先，订购钢管之前应该考虑到需要采购的数量、尺寸、材质、规格等因素，这些都是影响钢管选择和订购的重要因素。
在订购时，应该选择信誉良好、质量可靠的钢管供应商，以确保产品质量。
同时，还需要考虑到交货时间、运输方式、支付条款等因素。
其次，订购钢管之前要确定好计划，包括订购的数量、交货时间、付款方式等，以及产品的相关信息，如材质、规格等。
这些都可以通过与供应商之间进行沟通和协商来确定，以确保您获得最满意的订购计划。
最后，在订购钢管之前，还需要对供应商的信誉度和产品质量进行评估，包括查看其资质、生产能力、技术水平、售后服务等。
只有找到可靠、有信誉的供应商，才能保证订购到质量优良的钢管，满足您的需求。
钢管运输计划运输是钢管订购过程中不可或缺的一环，因为它可以直接影响到订购的成本和交货时间。
以下是一些关于钢管运输计划的关键方面：运送距离：决定了运输成本。
运输方式：运输方式包括海运、铁路运输、公路运输等，选择最佳的运输方式将有助于节约成本和缩短交货时间。
包装方式：钢管通常需要加强包装，以保护产品在运输过程中不被破坏。
运输包装应该符合国际标准，并保证产品不容易损坏。
运输时间：运输时间是关键因素之一，因为它可以影响到交货时间，因此需要尽可能找到最快的运输方式。
交货方式：交货方式包括门到门、门到港口、港口到港口等，选择最适合的交货方式可以有助于节约成本和提高效率。
结论钢管订购和运输计划需要综合考虑许多因素，包括产品尺寸、数量、质量、交货时间、成本、运输方式等。
只有充分了解这些关键因素，才能成功制定好的钢管订购和运输计划。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
模型一的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;endendendA(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);endendT=A;m=1;while m<=39for i=1:39for j=1:39if T(i,j)>T(i,m)+T(m,j);T(i,j)=T(i,m)+T(m,j);endendendm=m+1;endT2=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&T(i,j)<=300 T2(i,j)=20;elseif T(i,j)>300&T(i,j)<=350 T2(i,j)=23;elseif T(i,j)>350&T(i,j)<=400 T2(i,j)=26;elseif T(i,j)>400&T(i,j)<=450 T2(i,j)=29;elseif T(i,j)>450&T(i,j)<=500 T2(i,j)=32;elseif T(i,j)>500&T(i,j)<=600 T2(i,j)=37;elseif T(i,j)>600&T(i,j)<=700 T2(i,j)=44;elseif T(i,j)>700&T(i,j)<=800 T2(i,j)=50;elseif T(i,j)>800&T(i,j)<=900 T2(i,j)=55;elseif T(i,j)>900&T(i,j)<=1000 T2(i,j)=60;else T2(i,j)=60+ceil((T(i,j)-1000)/100)*5;endendendB=ones(39,39);for i=1:39;for j=1:39;if i==j B(i,j)=0;else B(i,j)=inf;endendendB(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);endendR=B;g=1;while g<=39;for i=1:39;for j=1:39;if R(i,j)>R(i,g)+R(g,j);R(i,j)=R(i,g)+R(g,j);endendendg=g+1;endR2=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);endendendf=1;while f<=39for i=1:39for j=1:39if w(i,j)>w(i,f)+w(f,j)w(i,j)=w(i,f)+w(f,j);endendendf=f+1;endQ=ones(39,39);for i=1:7for j=25:39cost(i,j)=w(i,j)endend模型一的matlab结果cost =Columns 25 through 39170.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.0000215.7000 205.3000 190.2000 171.6000 111.0000 95.5000 86.0000 71.2000 114.2000 142.0000 146.0000 156.0000 171.2000 178.0000 192.0000230.7000 220.3000 200.2000 181.6000 121.0000 105.5000 96.0000 86.2000 48.2000 82.0000 86.0000 96.0000 111.2000 118.0000 132.0000260.7000 250.3000 235.2000 216.6000 156.0000 140.5000 131.0000 116.2000 84.2000 62.0000 51.0000 61.0000 76.2000 83.0000 97.0000255.7000 245.3000 225.2000 206.6000 146.0000 130.5000 121.0000 111.2000 79.2000 57.0000 33.0000 51.0000 71.2000 73.0000 87.0000265.7000 255.3000 235.2000 216.6000 156.0000 140.5000 131.0000 121.2000 84.2000 62.0000 51.0000 45.0000 26.2000 11.0000 28.0000275.7000 265.3000 245.2000 226.6000 166.0000 150.5000 141.0000 131.2000 99.2000 76.0000 66.0000 56.0000 38.2000 26.0000 2.0000模型一的lingo程序model:!钢管购买与运输铺优化设问;sets:!cj表示厂家;md表示目的地；cj/1..7/:p,s,t;md/1..15/:y,z,D;link(cj,md):cost,x;endsets!这里是数据;data:s=800 800 1000 2000 2000 2000 3000;D=104,301,750,606,194,205,201,680,480,300,220,210,420,500,0;p=160,155,155,160,155,150,160;cost=170.7 160.3 140.2 98.6 38.0 20.5 3.1 21.2 64.2 92.0 96.0 106.0 121.2 128.0 142.0215.7 205.3 190.2 171.6 111.0 95.5 86.0 71.2 114.2 142.0 146.0 156.0 171.2 178.0 192.0230.7 220.3 200.2 181.6 121.0 105.5 96.0 86.2 48.2 82.0 86.0 96.0 111.2 118.0 132.0260.7 250.3 235.2 216.6 156.0 140.5 131.0 116.2 84.2 62.0 51.0 61.0 76.2 83.0 97.0255.7 245.3 225.2 206.6 146.0 130.5 121.0 111.2 79.2 57.0 33.0 51.0 71.2 73.0 87.0265.7 255.3 235.2 216.6 156.0 140.5 131.0 121.2 84.2 62.0 51.0 45.0 26.2 11.0 28.0275.7 265.3 245.2 226.6 166.0 150.5 141.0 131.2 99.2 77.0 66.0 56.0 38.2 26.0 2.0;enddata!目标函数;min=w;w=@sum(link(i,j):(p(i)+cost(i,j))*x(i,j))+0.1/2*@sum(md(j):y(j)^2+y(j)+z(j)^2+z(j));!约束条件;@for(cj(i):@sum(md(j):x(i,j))>=500*t(i);s(i)*t(i)>=@sum(md(j):x(i,j));@bin(t(i)));@for(md(j):@sum(cj(i):x(i,j))=y(j)+z(j));@for(md(j)|j#ne#15:D(j)=y(j)+z(j+1));z(15)=0;y(1)=0;@gin(@sum(link(i,j):x(i,j)));end铺设方案 Variable Value Reduced CostY( 1) 0.000000 0.000000 Y( 2) 75.00000 0.000000 Y( 3) 282.0000 0.000000 Y( 4) 0.000000 0.000000 Y( 5) 9.500000 0.000000 Y( 6) 15.50000 0.000000 Y( 7) 76.00000 0.000000 Y( 8) 175.0000 0.000000 Y( 9) 159.0000 0.000000 Y( 10) 30.00000 0.000000 Y( 11) 145.0000 0.000000 Y( 12) 11.00000 0.000000 Y( 13) 134.0000 0.000000 Y( 14) 500.0000 0.000000 Y( 15) 0.000000 141.0500 Z( 1) 0.000000 405.7500 Z( 2) 104.0000 0.000000 Z( 3) 226.0000 0.000000 Z( 4) 468.0000 0.000000 Z( 5) 606.0000 0.000000 Z( 6) 184.5000 0.000000 Z( 7) 189.5000 0.000000 Z( 8) 125.0000 0.000000 Z( 9) 505.0000 0.000000 Z( 10) 321.0000 0.000000 Z( 11) 270.0000 0.000000 Z( 12) 75.00000 0.000000 Z( 13) 199.0000 0.000000 Z( 14) 286.0000 0.000000 Z( 15) 0.000000 0.000000 厂家分配方案X( 1, 1) 0.000000 28.00000 X( 1, 2) 0.000000 28.00000 X( 1, 3) 0.000000 23.00000 X( 1, 4) 250.4441 0.000000 X( 1, 5) 84.05588 0.000000 X( 1, 6) 200.0000 0.000000 X( 1, 7) 265.5000 0.000000 X( 1, 8) 0.000000 23.00000 X( 1, 9) 0.000000 99.00000 X( 1, 10) 0.000000 143.0000X( 1, 12) 0.000000 174.0000 X( 1, 13) 0.000000 208.0000 X( 1, 14) 0.000000 230.0000 X( 1, 15) 0.000000 264.0000 X( 2, 1) 0.000000 0.000000 X( 2, 2) 179.0000 0.000000 X( 2, 3) 321.0000 0.000000 X( 2, 4) 0.000000 0.000000 X( 2, 5) 0.000000 0.000000 X( 2, 6) 0.000000 2.000000 X( 2, 7) 0.000000 9.900000 X( 2, 8) 300.0000 0.000000 X( 2, 9) 0.000000 76.00000 X( 2, 10) 0.000000 120.0000 X( 2, 11) 0.000000 148.0000 X( 2, 12) 0.000000 151.0000 X( 2, 13) 0.000000 185.0000 X( 2, 14) 0.000000 207.0000 X( 2, 15) 0.000000 241.0000 X( 3, 1) 0.000000 5.000000 X( 3, 2) 0.000000 5.000000 X( 3, 3) 147.2817 0.000000 X( 3, 4) 119.7061 0.000000 X( 3, 5) 69.01229 0.000000 X( 3, 6) 0.000000 2.000000 X( 3, 7) 0.000000 9.900000 X( 3, 8) 0.000000 5.000000 X( 3, 9) 664.0000 0.000000 X( 3, 10) 0.000000 50.00000 X( 3, 11) 0.000000 78.00000 X( 3, 12) 0.000000 81.00000 X( 3, 13) 0.000000 115.0000 X( 3, 14) 0.000000 137.0000 X( 3, 15) 0.000000 171.0000 X( 4, 1) 0.000000 5.000000 X( 4, 2) 0.000000 5.000000 X( 4, 3) 0.000000 5.000000 X( 4, 4) 0.000000 5.000000 X( 4, 5) 0.000000 5.000000 X( 4, 6) 0.000000 7.000000 X( 4, 7) 0.000000 14.90000 X( 4, 8) 0.000000 5.000000 X( 4, 9) 0.000000 6.000000X( 4, 11) 0.000000 13.00000 X( 4, 12) 0.000000 16.00000 X( 4, 13) 0.000000 50.00000 X( 4, 14) 0.000000 72.00000 X( 4, 15) 0.000000 106.0000 X( 5, 1) 0.000000 5.000000 X( 5, 2) 0.000000 5.000000 X( 5, 3) 39.71834 0.000000 X( 5, 4) 97.84982 0.000000 X( 5, 5) 462.4318 0.000000 X( 5, 6) 0.000000 2.000000 X( 5, 7) 0.000000 9.900000 X( 5, 8) 0.000000 5.000000 X( 5, 9) 0.000000 6.000000 X( 5, 10) 262.7391 0.000000 X( 5, 11) 415.0000 0.000000 X( 5, 12) 0.000000 11.00000 X( 5, 13) 0.000000 50.00000 X( 5, 14) 0.000000 67.00000 X( 5, 15) 0.000000 101.0000 X( 6, 1) 0.000000 10.00000 X( 6, 2) 0.000000 10.00000 X( 6, 3) 0.000000 5.000000 X( 6, 4) 0.000000 5.000000 X( 6, 5) 0.000000 5.000000 X( 6, 6) 0.000000 7.000000 X( 6, 7) 0.000000 14.90000 X( 6, 8) 0.000000 10.00000 X( 6, 9) 0.000000 6.000000 X( 6, 10) 88.26093 0.000000 X( 6, 11) 0.000000 13.00000 X( 6, 12) 86.00000 0.000000 X( 6, 13) 333.0000 0.000000 X( 6, 14) 786.0000 0.000000 X( 6, 15) 0.000000 37.00000 X( 7, 1) 0.000000 9.000000 X( 7, 2) 0.000000 9.000000 X( 7, 3) 0.000000 4.000000 X( 7, 4) 0.000000 4.000000 X( 7, 5) 0.000000 4.000000 X( 7, 6) 0.000000 6.000000 X( 7, 7) 0.000000 13.90000 X( 7, 8) 0.000000 9.000000X( 7, 10) 0.000000 4.000000X( 7, 11) 0.000000 17.00000X( 7, 12) 0.000000 0.000000X( 7, 13) 0.000000 1.000000X( 7, 14) 0.000000 4.000000X( 7, 15) 0.000000 0.000000最少总费用Local optimal solution found.Objective value: 1281354.Extended solver steps: 0Total solver iterations: 36模型三的matlab程序A=ones(39,39);for i=1:39;for j=1:39;if i==j A(i,j)=0;else A(i,j)=inf;endendendA(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);endendT=A;m=1;while m<=39for i=1:39for j=1:39if T(i,j)>T(i,m)+T(m,j);T(i,j)=T(i,m)+T(m,j);endendendm=m+1;endT2=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&T(i,j)<=300 T2(i,j)=20;elseif T(i,j)>300&T(i,j)<=350 T2(i,j)=23;elseif T(i,j)>350&T(i,j)<=400 T2(i,j)=26;elseif T(i,j)>400&T(i,j)<=450 T2(i,j)=29;elseif T(i,j)>450&T(i,j)<=500 T2(i,j)=32;elseif T(i,j)>500&T(i,j)<=600 T2(i,j)=37;elseif T(i,j)>600&T(i,j)<=700 T2(i,j)=44;elseif T(i,j)>700&T(i,j)<=800 T2(i,j)=50;elseif T(i,j)>800&T(i,j)<=900 T2(i,j)=55;elseif T(i,j)>900&T(i,j)<=1000 T2(i,j)=60;else T2(i,j)=60+ceil((T(i,j)-1000)/100)*5;endendendB=ones(39,39);for i=1:39;for j=1:39;if i==j B(i,j)=0;else B(i,j)=inf;endendendB(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(18,34)=70;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);endendR=B;g=1;while g<=39;for i=1:39;for j=1:39;if R(i,j)>R(i,g)+R(g,j);R(i,j)=R(i,g)+R(g,j);endendendg=g+1;endR2=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);endendendf=1;while f<=39for i=1:39for j=1:39if w(i,j)>w(i,f)+w(f,j)w(i,j)=w(i,f)+w(f,j);endendendf=f+1;endQ=ones(39,39);for i=1:7for j=19:39cost(i,j)=w(i,j)endend模型三的matlab结果Columns 19 through 3995.0000 100.0000 105.0000 115.0000 125.0000 140.0000 170.7000 160.3000 140.2000 98.6000 38.0000 20.5000 3.1000 21.2000 89.2000 92.0000 122.0000 106.0000 121.2000 128.0000 142.0000145.0000 150.0000 155.0000 165.0000 175.0000 190.0000 215.7000 205.3000 190.2000 171.6000 111.0000 95.5000 86.0000 71.2000 139.2000 142.0000172.0000 156.0000 171.2000 178.0000 192.000080.0000 90.0000 95.0000 105.0000 115.0000 130.0000 230.7000 220.3000 200.2000 181.6000 121.0000 105.5000 96.0000 86.2000 130.0000 82.0000 112.0000 96.0000 111.2000 118.0000 132.000044.0000 55.0000 60.0000 70.0000 80.0000 95.0000 260.7000 250.3000 235.2000 216.6000 156.0000 140.5000 131.0000 116.2000 110.0000 62.0000 83.0000 61.0000 76.2000 83.0000 97.000037.0000 50.0000 50.0000 65.0000 70.0000 85.0000 255.7000 245.3000 225.2000 206.6000 146.0000 130.5000 121.0000 111.2000 105.0000 57.0000 73.0000 51.0000 71.2000 73.0000 87.000050.0000 37.0000 44.0000 20.0000 14.0000 26.0000 265.7000 255.3000 235.2000 216.6000 156.0000 140.5000 131.0000 121.2000 110.0000 62.0000 67.0000 45.0000 26.2000 11.0000 28.000060.0000 50.0000 55.0000 32.0000 23.0000 4.0000 275.7000 265.3000 245.2000 226.6000 166.0000 150.5000 141.0000 131.2000 124.0000 76.0000 78.0000 56.0000 38.2000 26.0000 2.0000模型三的lingo程序model:!钢管购买与运输铺优化设问;sets:cj/1..7/:p,s,t;md/1..21/:y,z,v,D;link(cj,md):cost,x;endsets!这里是数据;data:s=800 800 1000 2000 2000 2000 3000;D=104 301 750 606 194 205 201 680 480 300 220 210 420 500 0 0 0 0 0 0 0;p=168 155 155 160 155 150 160;cost= 95 100 105 115 125 140 170.7 160.3 140.2 98.6 38 20.5 3.1 21.2 89.2 92 122 106 121.2 128 142145 150 155 165 175 190 215.7 205.3 190.2 171.6 111 95.5 86 71.2 139.2 142 172 156 171.2 178 19280 90 95 105 115 130 230.7 220.3 200.2 181.6 121 105.5 96 86.2 130 82 112 96 111.2 118 13244 55 60 70 80 95 260.7 250.3 235.2216.6 156 140.5 131 116.2 110 62 83 61 76.2 83 9737 50 50 65 70 85 255.7 245.3 225.2 206.6 146 130.5 121 111.2 105 57 73 51 71.2 73 8750 37 44 20 14 26 265.7 255.3 235.2 216.6 156 140.5 131 121.2 110 62 67 45 26.2 11 2860 50 55 32 23 4 275.7 265.3 245.2 226.6 166 150.5 141 131.2 124 76 78 56 38.2 26 2;!目标函数;min=@sum(link(i,j):(p(i)+cost(i,j))*x(i,j))+0.05*(@sum(md(j)|j#ge#2#a nd#j#le#21:y(j)^2+y(j))+@sum(md(j)|j#le#14:z(j)^2+z(j))+@sum(md(j)|j# eq#9#or#j#eq#11#or#j#eq#17:v(j)^2+v(j))+@sum(md(j)|j#eq#17#or#j#eq#19 #or#j#eq#20:z(j)^2+z(j)));@for(cj(i):@sum(md(j):x(i,j))>=500*t(i));@for(cj(i):@sum(md(j):x(i,j))<=s(i)*t(i));@for(cj(i):@bin(t(i)));@for(md(j)|j#ne#9#and#j#ne#11#and#j#ne#17:v(j)=0);@for(md(j):@sum(cj(i):x(i,j))=y(j)+z(j)+v(j));@for(md(j)|j#lt#15:D(j)=y(j)+z(j+1));v(9)+y(16)=42;v(11)+v(17)=10;y(17)+z(18)=130;z(17)+y(19)=190;z(19)+y(20)=260;z(20)+y(21)=100;@gin(@sum(link(i,j):x(i,j)));end。