数学建模选修论文资料Word文件下载.docx

1、关键字： 改进的最短路径 光的传播 Matlab 数学模型一、问题的重述某油田计划在铁路线一侧建造两家炼油厂，同时在铁路线上增建一个车站，用来运送成品油。由于这种模式具有一定的普遍性，油田设计院希望建立管线建设费用最省的一般数学模型与方法。利用模型分析管线布置和管线费用的情况，具体问题如下：1. 针对两炼油厂到铁路线距离和两炼油厂间距离的各种不同情形，提出你的设计方案。在方案设计时，若有共用管线，应考虑共用管线费用与非共用管线费用相同或不同的情形。2. 设计院目前需对一更为复杂的情形进行具体的设计。两炼油厂的具体位置由附图所示，其中A厂位于郊区（图中的I区域），B厂位于城区（图中的II区域），

2、两个区域的分界线用图中的虚线表示。图中各字母表示的距离（单位：千米）分别为a = 5，b = 8，c = 15，l = 20。 若所有管线的铺设费用均为每千米7.2万元。 铺设在城区的管线还需增加拆迁和工程补偿等附加费用，为对此项附加费用进行估计，聘请三家工程咨询公司（其中公司一具有甲级资质，公司二和公司三具有乙级资质）进行了估算。估算结果如下表所示：工程咨询公司公司一公司二公司三 附加费用（万元/千米）212420请为设计院给出管线布置方案及相应的费用。3. 在该实际问题中，为进一步节省费用，可以根据炼油厂的生产能力，选用相适应的油管。这时的管线铺设费用将分别降为输送A厂成品油的每千米5.6

3、万元，输送B厂成品油的每千米6.0万元，共用管线费用为每千米7.2万元，拆迁等附加费用同上。请给出管线最佳布置方案及相应的费用。二、模型假设1、管道均以直线段铺设，不考虑地形影响。2、不考虑管道的接头处费用。3、不考虑施工之中的意外情况，所有工作均可顺利进行。4、共用管线的价格如果和非公用管线不一致，则共用管线价格大于任意一条非公用管线价格，小于两条非公用管线价格之和。三、符号说明h：共用管道的高度（问题一中b）h1:共用管道高度h2:管线与分界线的交点到B厂与铁路平行线的距离w：方案的经费a：A厂到铁路的距离b：B厂到铁路的距离c:A厂到城郊分界线的距离l:A、B两厂之间的铁路长度x：A厂离

4、共用管道的距离（问题一中的c）y：共用管道的高度（问题一中的c）m：共用管道的费用（问题一）n：非共用管道费用（问题一）y1：为o点的纵坐标y2：为o1点的纵坐标x1：为o点的横坐标x2：为o1点的横坐标L: 为管线总长度（问题一中的b）四、问题分析要考虑有和没有共用管线，还要考虑共用管线与非共用管线费用相同和不同两种情况。同时还要考虑两个工厂是否在铁路的同一侧，如果两个工厂在铁路的同一侧那么一定没有共用管线。 不在铁路的同一侧那么就要考虑有和没有共用管线这个问题。计算共用管线的长度时，用光学原理，把一个工厂当作光源发射一束光经过一个平面的反射通过另一个工厂，这样能够保证路线最短。这个平面与铁

5、路的距离即为共用管线的长度。同时与这个平面的交点就是两厂的管线的交点。当共用管线与非共用管线费用不相同时可以通过建立方程组来解答。当共用管线与非共用管线费用不相同时要建立方程组来计算其最小费用从而来确定方案的可行性，共用管线与非共用管线长度作为变量来控制总费用，那么我们就可以列出一个方程组，从而在变量的约束条件下可以确定最小费用。把这个问题分两部分来考虑，即市区和郊区分两个部分，火车站建立在郊区费用要小得多，郊区共用管线与非共用管线的费用相同所以可以用最短路径的方法来考虑，同时又要求费用最小，可以解出最低费用及对应的铺设线路。通过建立坐标系设两个点的坐标，同时也是表达管线的长度，然后再与各自的

6、费用之积确定总的费用，从而算出两点的坐标值。即确定了管线的路线。五、模型的建立与求解5.1关于问题1的模型建立与求解对于管线布置的分析，分为两种情况：1. 两厂分别在铁路的两侧如下图：那么连接两厂A、B与铁路的交点C即为火车站的位置。2. 当两厂位于铁路的同一侧时,此时要分有公用管线与没有公用管线两种情况。a.当没有公用管线时，此时找出两厂与铁路交点连线的最近路线即可，如图：过铁路作A厂的对称点A，连接AB与铁路交于一点C，该点C即为火车站的位置。b.当有共用管线时又要分为共线管线费用与非共线管线费用相同与不同两种情况：当共线管线与非共线管线相同时，费用为m万元/千米如图所示：假设共线管线的长

7、度为h，A厂到铁路的距离为a，B厂到铁路的距离为b，则总的管线长度为：则总费用：c.当共线管线与非共线管线不同时，共用管线费用为m万元/千米 ，非共用管线费用为n万元/千米，如图所示：总费用为：其中实际的费用可以根据已知道的常量a、b、l再结合x、y的取值范围可以得出最小费用。5.2关于问题2的模型建立与求解因为在城区和郊区铁路管线的费用相同，但城区要增加拆迁和工程补偿等费用，因此城区和郊区要分为两部分来考虑。我们考虑三家咨询公司给出的三个方案，我们考虑到甲级资质和乙级资质的评估准确性，首先排除掉公司二的预算，对于公司一和公司三的预算，我们将分别求出最小费用，考察两者的差别。1假设共用管线在郊

8、区把该模型看作是一束光从B点发射在分界处G点发生了折射，把左边的问题看作是最短路径问题，如图所示：设共用管线的长度为h1，G点到O2B的距离为h2。在区域中即BG段每千米的费用为：20+7.2=27.2万元。由以上分析数据可得如下关系式：总费用： W1（最小）= （式1）参数的取值范围： （式2） （式3）利用Matlab将式（1）（2）（3）联立关系式绘图：用Microsoft Visual C+ 6.0解:W1（最小）= 275.13404万元运行结果：在这种情况下采用公司一的预算，只需要在上式中将27.2增加为28.2即可，计算得到总费用：280.177831万元2假设共用管线在城区同理

9、，如图所示：W2（最小）=（式1）参数用Microsoft Visual C+ 6.0解得W2（最小）= 355.25587显然W1（最小）W2（最小）方案一费用少于方案二，因此舍掉这种方案。 最终求得的结果为，如果采用一咨询公司的估算价格，则最终费用为275.134304万元，如果采用三咨询公司的估算价格，则最终费用为280.177831万元，考虑到公司一具有高级资质，因此我们采用公司一的价格方案，将最终预算设为280.177831万元，但是实际铺设管道的价格有可能在两种估算价格之间。5.3关于问题3的模型建立与求解1、O点为B管线与分界线的交点，O1点为A管与B管的交点，如下图建立坐标轴，

10、采用公司三的估算费用，总费用等于各段路线的长度与各段费用的积为：坐标法解答，A01，OO1，OB，如图：O（x1，y1）,O1（x2，y2）A厂到管道交接点O1的长度：AO1=管道交点O1到B厂与城郊分界线交点O的长度：OO1= （式2）B厂到交点O的长度：OB=铁路站点O2到交叉管道O1的长度：O1O2= （式4） （式5） （式6） （式7）由以上式子利用Microsoft Visual C+ 6.0软件求得最小经费：W3（最小值）= 244.386494万元。在同种情况下，用公司一的预算费的总费用：当火车站建在市区费用太高同二题中的方案二故不选用那种方案。因此类似于问题二，我们采取公司一

11、的估算价格，最终预算为249.468791万元。但是公司三的价格也具有一定参考性，实际铺设管道价格应为244.386494万元到249.468791万元之间。六、模型的评价与应用从实际的生活出发输油管道是石油生产过程中的重要环节，是石油工业的动脉。在石油的生产过程中，至始至终都离不开输油管道。我们可以把石油的生产过程简单的表示为： dO_;vcg_vb 油 计量站 井联合站 转油站 矿场油库 炼油厂 用户M_.P（/FS9 从油井出来的油气通过管道输送到计量站，经过计量后又由管道输送往联合站，在联合站生产出合格的原油，合格原油通过管道和转油站输到矿场油库或外输到管道首站，通过长输原油管道输到炼

12、油厂加工精练，生产出各种产品，通过成品油管道或铁路、公路、水路将各种产品送往用户，其中成品油管道就需要用到管道的布置设计。qetP93_N_*优点：模型使问题由复杂变简单，方便运输，提高输油效率，规划线路。管线布置和规划及相应的费用减到最小，在不同的环境下用这种环境中的最优模型，方便快捷，节约开支，使实际问题更加精确。同时对于题目的三个问题都设计了合适的模型，并且当给出具体数值的时候能够给出足够精确的解，具有一定的普遍性。OOokh_Zd其四：其四 e # 5BPI_缺点：该模型在提出的时候将部分因素没有考虑进来，例如管线接头处的费用，以及工作工程中的一些意外情况等等，使得该模型在实际应用中会

13、缺少精确性。应用：模型在实际运用中，不仅仅可以用在成品油运输管布置，还可运用到原油输送和污水处理，电线电缆的布置还有公路铁路的修建等一些列的线路布置问题。b8rp8 _r_b8,I6 Rx7QlA_;七、参考文献【1】 赵静 但琦 数学建模与数学实验第三版 22-29页，178-194页 高等教育出版社 2008年1月【2】 曹戈 MATLAB教程及实训 37-60页 机械工业出版社 2008年5月【3】 邬学军 周凯数学建模竞赛辅导教程 73-96页 浙江大学出版社 2009年1月附录：问题2程序1、按照公司三的评估总费用为：#includemath.hvoid main（） double

14、h1,h2,w;double a,b;double min = 10000;for（h1=0;h1=8;h1+=0.001）for（h2=0;h28） continue; w=27.2*sqrt（25+h2*h2）+（sqrt（5-h1+8-h1-h2）*（5-h1+8-h1-h2）+225）+h1）*7.2; if（minw） min=w; a=h1;b=h2; printf（%f n,min）; printf（%f %f n,a,b）;2、按照公司一评估总费用为：问题3程序：1、按照公司三评估总费用为：#includestdio.hmath.hmain（） double x2,y2,y1,w,a,b,c,min=1000; for（x2=-15;x2=0;x2+=1）for（y1=0;y1y1+=1）for（y2=0;y2 min=w; a=x2;b=y2;c=y1;min=%f nx2=%f ny2=%f ny1=%fn,min,a,b,c）;将循环因子的步长降低之后，进一步精确求解值为： for（x2=-9;=-7;x2+=0.001）for（y1=6;y1+=0.001）=1;y2+=0.001）2、按照公司一的评估的总费用为：