数学建模实验二.doc

钢管订购和运输 一、实验目的 学习数学规划建模方法，并能够利用lindo或lingo 完成计算，并完成一篇完整论文，对问题结果做简单分析。 二、实验内容 原始问题(2000-B-1) 有7个生产厂，可以生产输送天然气主管道的钢管Sii=1,2,3,4,5,6,7。要沿着Aii=1…15的主管道铺设, 如题图一所示。图中粗线表示铁路，单细线表示公路，双细线表示要铺设的管道(假设沿管道或者原来有公路，或者建有施工公路)，圆圈表示火车站，每段铁路、公路和管道旁的阿拉伯数字表示里程(单位km)。为方便计，1km主管道钢管称为1单位钢管。 一个钢厂如果承担制造这种钢管，至少需要生产500个单位。钢厂Si在指定期限内能生产该钢管的最 大数量为Si个单位，钢管出厂销价1单位钢管为Pi万元，如下表： 表 SEQ 图表 \* ARABIC 1钢铁厂的最大产量和单位钢管的价格 表 SEQ 图表 \* ARABIC 2 1单位钢管的铁路运价如下表 里程(km) ≤300 301～350 351～400 401～450 451～500 运价(万) 20 23 26 29 32 里程(km) 501～600 601～700 701～800 801～900 901～1000 运价(万) 37 44 50 55 60 1000km以上每增加1至100km运价增加5万元。公路运输费用为1单位钢管每公里0.1万元（不足整公里部分按整公里计算）。钢管可由铁路、公路运往铺设地点（不只是运到点，而是Aii=1…15管道全线）。 请制定一个主管道钢管的订购和运输计划，使总费用最小（给出总费用)？ Floyd算法介绍 通过一个图的权值矩阵求出它的每两点间的最短路径矩阵。从图的带权邻接矩阵A=[a(i,j)]n*n开始，递归地进行n次更新，最终得到D(n)，矩阵的i行j列元素是i号顶点到j号顶点的最短路径长度，称D(n)为图的距离矩阵，同时还可引入一个后继结点矩阵path来记录两点间的最短路径。算法过程为： 从任意一条单边路径开始。所有两点之间的距离是边的权，如果两点之间没有边相连，则权为无穷大。 对于每一对顶点 u 和 v，看看是否存在一个顶点 w 使得从 u 到 w 再到 v 比已知的路径更短。如果是更新它。 把图用邻接距阵G表示出来，如果从Vi到Vj有路可达，则G[i,j]=d，d表示该路的长度；否则G[i,j]=无穷大。 定义一个距阵D用来记录所插入点的信息，D[i,j]表示从Vi到Vj需要经过的点，初始化D[i,j]=j。 把各个顶点插入图中，比较插点后的距离与原来的距离，G[i,j] = min( G[i,j], G[i,k]+G[k,j] )，如果G[i,j]的值变小，则D[i,j]=k。 在G中包含有两点之间最短道路的信息，而在D中则包含了最短通路径的信息。 比如，要寻找从V5到V1的路径。根据D，假如D(5,1)=3则说明从V5到V1经过V3，路径为{V5,V3,V1}，如果D(5,3)=3，说明V5与V3直接相连，如果D(3,1)=1，说明V3与V1直接相连。 问题分析 首先，所有钢管必须运到天然气主管道铺设路线上的节点Aii=1…15，然后才能向左或右铺设。必须求出每个钢管厂Sii=1,2,3,4,5,6,7到每个节点Aii=1…15的每单位钢管的最小运输费用。 对最小运费的求解，我们采用Floyd算法。先求出铁路网上钢管厂到铁路上任意两点Vi，Vj的最短路线的长度Lij,用matlab求得Lij对应的铁路单位运费Dij；同理用Floyd 算法求出公路网上的任意两点Vj，Vk 的最短公路路线的长度Ljk，结果乘以0.1得到公路运费D1jk。Cik=min(Dij+D1jk)，j表示所有运输中转点，于是就得到从某钢厂到某铺设点运输单位钢管的最少运输费用。 每个铺设点分别向R.L两个方向展开，通过Lingo编程求出最小铺设费用。运输费用加上购买费用再加上铺设费用就是我们所要求的总费用。 4、模型的假设与符号说明 1) 基本假设 eq \o\ac(○,1)钢管在运输中由铁路运转为公路运时不计中转费用； eq \o\ac(○,2)所需钢管均由Sii=1,2,3,4,5,6,7 钢厂提供； eq \o\ac(○,3)假设运送的钢管路途中没有损耗。 2) 符号说明 Si:钢厂Si的最大生产能力; Pi:钢厂Si 的出厂钢管单位价格(单位: 万元) ; d:公路上一单位钢管的每公里运费(d = 0. 1 万元) ; Dij铁路网上两点间的单位钢管最少运输费用； D1jk:公路网上两点间的单位钢管最少运输费用； e:铁路上一单位钢管的运费(分段函数见表1) ; Cij:1单位钢管从钢厂Si运到Aj的最小费用(单位: 万元) ; bj:从Aj 到Aj+1之间