2000年数学建模题钢管订购和运输

1钢管订购和运输摘要本文根据问题的条件和要求，建立两个模型，两个模型均为单目标非线性规划模型，并通过求解这两个模型，完整地解决了问题。由于铁路运输费用函数具有不可加性，不能直接应用现有的最短路算法来求解铁路和公路交通网中任意两点间最小费用路问题。本文采用了一种分步递推算法，巧妙解决了这一问题。在单目标非线性规划模型中，将管道铺设分为两个过程。先将钢管从钢管厂运到管道与道路交叉口，再从交叉口铺设到管道线上。这样，总的运输费用就化为两个过程的运输费用之和。本模型是以总费用为目标函数的非线性规划模型，利用Lingo软件，求出问题一的最优解为1278632万元。对于问题二通过对模型1的灵敏度分析，确定了5S钢厂的销价的变化对购运计划和总费用的影响最大，确定S1钢厂的生产上限的变化对物运计划和总费用的影响最大。问题三模型的建立原理和问题一的相同，利用Lingo软件，求得最优解为1407149万元.关键词：Floyd算法单目标非线性规划灵敏度分析2问题重述有7个生产厂，可以生产输送天然气主管道的钢管721,,SSS。要沿着1521AAA的主管道铺设,如题图一所示。图中粗线表示铁路，单细线表示公路，双细线表示要铺设的管道(假设沿管道或者原来有公路，或者建有施工公路)，圆圈表示火车站，每段铁路、公路和管道旁的阿拉伯数字表示里程(单位km)。为方便计，1km主管道钢管称为1单位钢管。一个钢厂如果承担制造这种钢管，至少需要生产500个单位。钢厂iS在指定期限内能生产该钢管的最大数量为is个单位，钢管出厂销价1单位钢管为ip万元，如下表：123456780080010002000200020003000160155155160155150160iisip1单位钢管的铁路运价如下表：里程(km)≤300301～350351～400401～450451～500运价(万元)2023262932里程(km)501～600601～700701～800801～900901～1000运价(万元)37445055601000km以上每增加1至100km运价增加5公路运输费用为1单位钢管每公里0.1万元（不足整公里部分按整公里计算）。钢管可由铁路、公路运往铺设地点（不只是运到点，而是管道全线）。（1）请制定一个主管道钢管的订购和运输计划，使总费用最小（给出总费用)。（2）请就（1）的模型分析：哪个钢厂钢管的销价的变化对购运计划和总费用影响最大，哪个钢厂钢管的产量的上限的变化对购运计划和总费用的影响最大，并给出相应的数字结果。（3）如果要铺设的管道不是一条线，而是一个树形图，铁路、公路和管道构成网络，请就这种更一般的情形给出一种解决办法，并对题图二按（1）的要求给出模型和结果。1521,,,AAA3问题分析问题一，首先，所有钢管必须运到天然气主管道铺设路线上的节点1521AAA，然后才能向左或右铺设。必须求出每个钢管厂721,,SSS到每个节点1521AAA的每单位钢管的最小运输费用。对最小运费的求解，我们采用Floyd算法,先求出铁路网上钢管厂到铁路上任意两点iV，jV的最短路线的长度ijL,用matlab求得ijL对应的铁路单位运费ijD；同理用Floyd算法求出公路网上的任意两点jV，kV的最短公路路线的长度jkL，结果乘以0.1得到公路运费jkD1。)1min(jkijikDDC，j表示所有运输中转点，于是就得到从某钢厂到某铺设点运输单位钢管的最少运输费用。（具体算法及程序见附录）每个铺设点分别向y,z两个方向展开，通过Lingo编程求出最小铺设费用。运输费用加上购买费用再加上铺设费用就是我们所要求的总费用。问题二，通过问题一里面Lingo编程运行得出的结果，分析哪个钢厂钢管的销价的变化对购运计划和总费用影响最大，哪个钢厂钢管的产量的上限的变化对购运计划和总费用的影响最大。问题三，利用同问题一一样的方法，从而可求出某钢厂到某某铺设点运输单位钢管的最少运输费用。（具体算法及程序见附录）模型的假设与符号说明1)基本假设:○1要铺设的管道侧有公路，可运送所需钢管。○2钢管在运输中由铁路运转为公路运时不计中转（换车）费用；○3所需钢管均由)7,...,1(iSi钢厂提供；④假设运送的钢管路途中没有损耗。2)符号说明:iS:钢厂iS的最大生产能力;4ip:钢厂iS的出厂钢管单位价格(单位:万元);d:公路上一单位钢管的每公里运费(d=0.1万元);e:铁路上一单位钢管的运费(分段函数见表1);ijc:1单位钢管从钢厂iS运到jA的最小费用(单位:万元);jb:从jA到1jA之间的距离(单位:千米);ijx:钢厂iS运到jA的钢管数;yj:运到jA地的钢管向左铺设的数目;jZ：运到jA地的钢管向右铺设的数目;it:=不提供钢管；，钢厂提供钢管；，钢厂iiS0S1W:所求钢管订购、运输的总费用(单位:万元);模型的建立与求解问题一的模型：针对题图一，我们采用Floyd算法，用matlab编程求出单位钢管从iS运输到jA的最小运输费用，具体数据如下表1：表1单位钢管从iS运输到jA的最小运输费用（单位：万元）S1S2S3S4S5S6S7A1170.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.6A538.0111.0121.0156.0146.0156.0166.0A620.595.5105.5140.5130.5140.5150.5A73.186.096.0131.0121.0131.0141.0A821.271.286.2116.2111.2121.2131.2A964.2114.248.284.279.284.299.2A1092.0142.082.062.057.062.077.0A1196.0146.086.051.033.051.066.0A12106.0156.096.061.051.045.056.0A13121.2171.2111.276.271.226.238.2A14128.0178.0118.083.073.011.026.0A15142.0192.0132.097.087.028.02.0对表1的数据进行分析，我们得到一个非线性规划模型：5目标函数是总费用W,它包含三项:钢管出厂总价Q,运输费P,及铺设费T.即W=Q+P+T其中ijijixpQ71151，ijijijxcP71151，铺设费T可以如下来确定：jA开始从左右两个方向铺设，jy与zj单位长钢管的费用为(1)12...2jjjyydddyd与(1)2jjzzd故1511122jjjjjyyzzTd目标函数为：715715151111111min22jjjjiijijijijijjyyzzWpxcxd约束条件为:①生产能力的限制:iijijitsxt151500，)7,...,1(i)10(或it②运到jA的钢管用完:71ijjjixyz，)15,...,1(j③jA与j-1A之间的钢管:j1jjzyb，)14,...,1(j④变量非负性限制:0,0,0jijjxyz,)15,...,1,7,...,1(ji⑤运到jA的钢管整数限制:Nxij模型一715715151111111min22jjjjiijijijijijjyyzzWpxcxds.t.iijijitsxt151500，)7,...,1(i)10(或it671ijjjixyz，)15,...,1(jj1jjzyb，)14,...,1(j1y=0,15z=00,0,0jijjxyz,)15,...,1,7,...,1(jiit=0或1(i=1,..,7)d=0.05;根据模型二编写Lingo程序，程序运行后，得到最优最小费用为1282142W万元。问题二的模型对模型1的灵敏度分析（1）确定哪个钢厂的销价的变化对购运计划和总费用的影响最大我们假设该钢厂的销价变化在10%ip万元以内，这是较为合理的，将目标函数的w表示为ip的函数：127127127w=fp,,...,),......ppfffwpppppp（因此在销价的变化量相同时，ifp越大，则ip的变化对w的变化影响越大。由模型Obj1计算得到的数据可以知道62000fp单位是最大的，所以S6的销价变化对购运计算和总费用的影响最大，我们可以通过简单的分析来证明：由于S6提供的数量最大，销价只要很小变化的，就会引起总费用的很大变化，同时，当价格越来越高，由于P6和1516ixj互为消长的关系，当S6越来越小，它在总需求中占的份额减少，影响减弱，S6下降的速度也将放慢。除了销价的升高，我们还必须考虑销价的降低，此时应尽量满足提供量最少的点S5，当价格越来越低，由互为消长的关系，S5点的提供量将增加，它在总需求中占的份额增加，影响增强，对于S5上升的速度将放慢。2）确定哪个钢厂的生产上限的变化对物运计划和总费用的影响最大由于S1是A1到A8的最优首选，因此若S1与其他Si同时扩大相同的S容量，则S1会更优，所以推断S1应为影响最大者。由最小费用矩阵C可以知道，Ai(i=1,…,8)所需的钢管量最好都能由S1提供，则此时S1达到最大需求量，在模型Obj1的条件下，S1为25367单位，而S1的上限为800单位，考虑到实际钢厂的投入与产出，在很短的时期内生产要达到原来的3倍，不符合实际意义，所以考虑S1在10%范围内变化。同理对于A0点，最优为S3全部提供，即S3应提供634单位，对于A12,A13,A14,A15,由S6全部提供为最优，即S6应提供1205单位，A11,A10由S5全部提供为最优，即S5应提供796单位。利用计算机模拟，得到5个供货钢厂分别扩建1%，2%，4%，6%，8%，10%时的成本的增长率，见表。可以看出，相同的S下A1产生的增长率最大，符合上述分析。一旦工厂扩建范围超过最大需求量，则不再会使目标函数优化，则此时增长率为0。即是上图中S5,S6的情况。而对于S1，一旦S12536,则其增长率也为0。（S1的数字结果见表2）△表2某个Si在变化S的情况下目标函数减小量及减小的比率△sS1S2S3S5S6△z(万元)zz(%)△z（万元)zz(%)△z（万元)zz(%)△z(万元)zz(%)△z(万元)zz(%)1%8720.0683280.0253100.02400002%17440.1366560.0516200.04800004%34880.27213120.10212400.09600006%52320.40819680.15318600.14500008%69760.54426240.20424800.193000010%87200.68532800.25631000.2420000问题三的模型题图二为树形图，采用Floyd算法，用matlab编程求出单位钢管从iS运输到jA的最小运输费用，具体数据如下表2：表2单位钢管从iS运输到jA的最小运输费用（单位：万元）S1S2S3S4S5S6S7A1170.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.6A538.0111.0121.0156.0146.0156.0166.