输油管道问题.docxVIP

数学模型课程设计 报 告 题目：输油管的布局 院 系： 班 级： 学 号： 姓 名： 华 南 理 工 大 学 2013.06 输油管的布局 问题描述 登峰油田计划在铁路(X轴)北侧建造两家炼油厂和(单位：km)，同时在铁路边上增建一个火车站来运送成品油。的区域是郊区，的区域是城区。所有输油管的铺设费用均为C1万元/km。 铺设在城区的输油管还需增加拆迁费用C2万元/km。 请你设计出从炼油厂到火车站的输油管最优布局方案使总费用最低。 M M N A(0,yA) B(20,yB) H(15, yH) 已知：yA，yB，C1，C2，求：yH，xN，yN，S(总费用)和铺管图。 算法原理简介。A( A(0,YA)) B(20,YB b) E M N H(15, yH) 1、 设点M是车站，点H是输油管在两个区域的分界点，点N是三角形AHM的费尔马点，三角形AHE是等边三角形，则点E的坐标 ， 而且，，则总费用 S=C1*Ye+(C1+C2)*HB 即S=C1*((Ya+Yh+15*sqrt(3))/2)+(C1+C2)*sqrt(25+(Yh-Yb)^2) 令Yh=t(i), 由题目要求易知Yh取值必定在0和Yb之间。通过一个for循环，使t(i)以0.01的步长增加，从而求得最小的总费用s，以及此时的Yh。 2、通过点A、点H和点M可以确定N的坐标（原理如下） 若ΔABC的三个内角都小于120o，分别在ΔABC之外作等边ΔABE和等边ΔADC，则直线BD与CE的交点F就是ΔABC 的费尔马点，且BD与CE就是费尔马距离，即BD=CE=AF+BF+CF. 当我们知道ΔABC三顶点的坐标时，借助平面几何可如下求点D和点E的坐标以及费尔马点： A A B T E C D M N y x 图1 设T是AB的中点，AN与EM平行与y轴,MT与BN平行与x轴. 以下我们用点的字母作坐标的下标，例如等等。则 (1) 可见，点E的坐标是很容易通过点A与点B的坐标计算出来的。同理得 (2) 然后求线段CE与BD的交点就得费尔马点F的坐标。若分别用与表示线段CE与BD的斜率，即，则易求得 (3) 一般地，A、B、C三点的位置未必与图1一致，从而影响结果的计算公式。因此，我们还要对此作进一步讨论。我们不妨设，即把横坐标最小者称为点A, 最大者称为点C, 中间者称为点B.这样仍按上述推导方法，仍设的AB边外与点E构成等边三角形，AC边外与点D构成等边三角形，可得如下结论。 当点B在线段AC之上时 (4) 其他情况 (5) 不论哪种情况，费尔马点的坐标计算，仍用(3)式。 如果三角形其中一个角大于或者等于120o则费尔马点与此顶点重合。 求得Xn= (15-(Yh-Ya)*3^(1/2))/2 Yn=(Ya+Yh)/2-(15*3^(1/2))/6 主要变量说明。 Ya、Yb、Yn、Ye、Yh分别表示点A、B、N、E和H的纵坐标，Xn、Xe分别表示点N和点E的横坐标，S表示总费用 程序开始算法流程图。 程序开始 输入 输入Ya、Yb、C1、C2 通过循环求得最小总费用 通过循环求得最小总费用s 以及此时对应Yh 求出 求出Xn、Yn Xn Xn0 ？ 是否 是 否 Xn15 Xn15 ？ Xn=0，Yn=Ya 是 是 否 否 Xn Xn=15，Yn=Yh 输出 输出s、Yh、Xn、Yn并绘制相应的管道布局图 四．操作说明： 1．打开程序即进入如下页面。 2按下问题描述，显示如下（再点击一次问题描述则隐藏问题描述）。 3根据提示，在文本框中输入YA、YA、C1、C2，输入完毕后按确认输入并且求解按钮即可在右下方得到所求的数据。左下方是铺管图。（以下是我求解时的多种情况）。 3． 重新算时，必须重新输入YA、YB、C1、C2，再按求解。求解完毕，需要退出时直接按退出键就好了。 总结 此次课程设计是我第一次如此近距离与matlab亲密接触，从而花了我不少时间，也遇到不少困难。 首先，通过网上学习、翻阅资料、请教同学学习matlab语言程序编写。在了解了matlab语言程序基本编写及运用之后，针对课程设计题目要求，经反复操作得到一相对较精确的解。 在此次课程设计过程中，学习了matlab gui及一些重要的函数和基本操作。将平时课堂上学习到的理论知识运用到实际的操作当中，自然受益颇多。在以后的学习当中还需不断将理论转化为实际操作。 附录 求解总费用s、Yh、Xn、Yn程序源代码。 输入Ya、Yb、C1、C2 t=0:0.001:Yb; s=C1*((Ya+15*sqrt(3))/2)+(C1+C2)*sqrt(25+Yb^2); for i=1:Yb*1000+1 W