数学建模报告-飞行问题

《数学建模》课程设计报告课题名称：___飞行管理问题系（院）：理学院专业：数学与应用数学班级：10122111学生姓名：邵仁和学号：1012211122指导教师：陈宏宇开课时间：2011-2012学年二学期1012211122邵仁和2飞行管理问题的优化模型摘要为了避免较多飞机在区域内会发生碰撞，让飞机在某正方形区域内安全飞行，便于进行飞行管理，所以在飞机飞行过程中，要适当调整各架飞机的方向角（调整幅度尽量小），所以这是个优化问题。本文我们根据题目所给的数据，利用matlab软件绘制出飞机的位置图标及飞行路径，并利用lingo软件找出了碰撞发生的飞机、碰撞发生的点和时间。同时再寻找判断两架飞机是否会相撞的方法，我们发现可以在飞机飞出区域之前每隔一段较短的时间对飞机进行监控，看是否与别的飞机相撞。然后，我们根据问题讨论了飞行方向角的调整时间和次数对最优解的影响，发现调整时间越早，调整角度就越小，所以我们决定在第六架飞机刚飞到区域边缘的时候就进行飞行角度的调整，并且达到了优化目标：61|)(|miniia。由题意，我们找到约束条件，然后把这些约束条件在lingo中用语言描述出来，再针对运算方面进行改进，得到我们的lingo程序，运行后我们得到了飞机调整的飞行方向角和方案。关键词：简化，最小调整幅度，最优1012211122邵仁和3一、问题重述6.飞行管理问题(优化模型)在约10000米高空的某边长160km的正方形区域内,经常有若干架飞机作水平飞行.区域内飞行的每架飞机的位置和速度向量均由计算机记录其数据,以便进行飞行管理.当一架欲进入该区域的飞机到达区域的边界时,记录其数据后,须立即判断是否将与区域内的飞机相碰撞.若可能发生碰撞,则应计算如何调整各架飞机的飞行的方向角，以避免碰撞。作如下假设:(1)任意两架飞机的安全飞行距离为8公里;(2)所有飞机的飞行速度为800公里/小时;(3)进入该区域的飞机在到达区域边界时,与区域内的飞机的距离应在60公里以上;(4)最多考虑6架飞机;(5)不必考虑飞机离开此区域后的情况.请你对这个避免碰撞的飞行管理问题建立数学模型，列出计算步骤，对以下数据进行计算（方向角误差不超过0.01），要求飞机飞行方向角调整的幅度尽量小。设该区域四个顶点的坐标为：(0,0)，(160,0)，(160,160)，(0,160)记录数据为：飞机编号横坐标x纵坐标y方向角（度）1150140243285852363150155220.54145501595130150230新进入0052注：方向角指飞行方向与x轴正向的夹角。1012211122邵仁和4二、问题分析根据问题容易知道，这显然是一个优化问题，当两架飞机可能发生碰撞时，即在规定区域内某一时刻两架飞机之间的距离小于8公里，因此要调整飞行方向一定角度，保证任意两架飞机在区域内任意时刻，两者的距离均不小于8公里，避免相撞。考虑到调整角度应尽量小，可以简化飞行方向调整策略，降低调整难度，同时减轻机内乘客及工作人员的不适。此外由此初步确定了调整目标：所有六架飞机的飞行方向调整角度均尽量小。三、模型假设与符号约定（1）假设飞机进入控制区域后完全服从地面控制台的调度，其他任何因素或人都不能改变飞机的飞行方向角。（2）假设从飞机管理处发出的信息飞机马上可以接收并执行（此处忽略飞机在执行过程中所需耗费的时间），不存在滞后或延迟，即可以实现实时控制。（3）不考虑本组设计的程序在实时控制中运行的时间。(4)不碰撞的标准为任意两架飞机的距离大于8km。（5）所有飞机飞行速度均为每小时800km。（6）飞机在区域外靠雷达自动与其他飞机保持距离大于60km，进入区域后由地面控制台进行统一控制，保证飞机距离大于8km。（7）假设飞机在区域内改变方向，在飞出区域后驾驶人员会自觉调整方向回归原航线继续飞行。（8）为了表达清晰，我们对符号作出以下说明：符号含义备注pi区域内的第i架飞机ttyxii,第i架飞机在时刻t时的坐标此处记第6架飞机进入区域时为t=0v飞机飞行速度按题中所给条件为800hkm/ia各飞机未调整方向前的的方向角飞机头与x轴正向间的夹角ia第i架飞机方向角的调整值1012211122邵仁和5四、模型的建立和求解1、问题简化明显我们不可能每时每刻都对每架飞机进行调整，当然希望调整的次数越少越好，调整的飞机数越少越好，角度越小越好。第一，我们需要知道哪些飞机会在飞行过程中与其他的飞机相撞（距离小于8km）。第二，我们需要调整哪些飞机和多少角度来使61|)(|iia最小。假使飞机要从A点向C飞行，途中想要偏离C点，由图αβ（外角大于内角），明显是从一开始就改变α角度使得61|)(|iia更小，所以越早调整越好，即在第六架飞机进入时即可调整角度。如图：2、模型建立由问题一，我们首先判断在第6架飞机进入正方形区域后会否发生飞机碰撞。我们依照数据，用matlab画出大致的航线图形（程序见附录1）。其中灰线代表飞机向上飞行，黑线代表飞机向下飞行。（如图1）ABαβCD1012211122邵仁和6图1：把飞机速度hkmv/800转化成skmv/...2222.0，由此我们可以知道飞机在5秒内仅飞行1km多，根据这个我们将飞机飞行时间分成5秒一段，由这种分段来达到近乎实时监控的目的，知晓每一架飞机每一个5秒的飞行情况，尽可能精确地求出有可能发生碰撞的飞机及碰撞发生的时间段，方便调整。我们对这个思路整合后用lingo编出程序（程序见附录2），运行得出结果：如果飞机在区域内按照原飞行方向角沿原航线继续飞行，那么是会发生一次飞机碰撞的。是在第88个5秒时，第5架和第6架飞机会发生碰撞。所以需要调整飞机的飞行方向角来避过这次相撞。下面进行模型求解：根据题意，目标函数应该为飞行方向角的调整幅度取最小值，所以我们可以得到目标函数：1012211122邵仁和761|)(|miniia题目告诉我们，不碰撞的标准为任意两架飞机的距离都大于8km。我们将飞机看做一个质点，求出两个质点之间的欧式距离：22)()()()(ttttyyxxjiji即为两架飞机间的距离（如图2）：图2可得到：6422ttttyyxxjiji，ji6|)(|iani，，，21为了方便计算并且更加直观，我们对本题理解为在第6架飞入时飞机管理处便进行计算，将得出的最优调整方案马上提供给飞机执行，即每架飞机最多在t=0时改变一次躲过碰撞，之后在区域内便不再改变方向角。根据这种理解，我们又得到一个约束条件：iaiavttxxiicos0ni，，，21iaiavttyyiisin0ni，，，21综上分析，可以得到第一个模型：yx0ttyxii,ttyxjj,ttyyjittxxji1012211122邵仁和861|)(|miniia6422ttttyyxxjiji，ji..ts)()(cos)0()(iaiavttxxiini，，，21iaiavttyyiisin0ni，，，216|)(|iani，，，21由此编写lingo程序（见附录3），可知最优解为624.023°，4-0.4954°，61.567°（其他调整角度为0）使得61|)(|iia最小。五、模型改进与推广模型优点：在调整各飞机偏转的方向角时，我们只调整少数飞机，能把问题的计算简单化，而且借助各架飞机初始位点的运动路线，更直观地初步判断各飞机于某时刻的基本碰撞情况。模型缺点：1这个模型是平面的，尚未考虑空间上的问题。2在这个模型中我们没有考虑飞机接受命令到执行命令之间的时间，事实上这段时间是存在且不可忽略的。3在实际的飞机航行中，改变飞行角度后，飞机便离开了原航线，在以后的飞行中是要矫正过来的，但是在我们的模型中这个问题并没有列入考虑范围。4在飞机的飞行过程中，它会受到风速、仪器、天气的影响，进而导致飞行员在调整角度的过程中可能出现偏差。并且在模型求解过程中用lingo求解最优解时，程序运行时间有些长，这在现实中是没有那么多的等待时间，因为飞机也在飞行，坐标还在改变，所以模型求解时间应越短越好。模型推广：此类模型可以用于航海中轮船及各种相遇问题等。总体来说这个模型比较简单易懂，符合一般的调整需要，不足的是模型需要强大的实时数据支持。六、参考文献[1]朱道元，《数学建模精品案例》，南京：东南大学出版社，1999[2]程极泰，《最优设计的数学方法》，北京：国防工业出版社，19941012211122邵仁和9附录（1）x1=150;y1=140;a1=243;b1=y1-x1*tan(a1/180*pi);xa=1:5:160;ya=xa*tan(a1/180*pi)+b1;plot(xa,ya,'r',x1,y1,'+')holdonx2=85;y2=85;a2=236;b2=y2-x2*tan(a2/180*pi);xb=1:5:160;yb=xb*tan(a2/180*pi)+b2;plot(xb,yb,'r',x2,y2,'+')holdonx3=150;y3=155;a3=220.5;b3=y3-x3*tan(a3/180*pi);xc=1:5:160;yc=xc*tan(a3/180*pi)+b3;plot(xc,yc,'r',x3,y3,'+')holdonx4=145;y4=50;a4=159;b4=y4-x4*tan(a4/180*pi);xd=1:5:160;yd=xd*tan(a4/180*pi)+b4;plot(xd,yd,'g',x4,y4,'+')holdonx5=130;y5=150;a5=230;b5=y5-x5*tan(a5/180*pi);xe=1:5:160;ye=xe*tan(a5/180*pi)+b5;plot(xe,ye,'r',x5,y5,'+')1012211122邵仁和10holdonx6=0;y6=0;a6=52;b6=y6-x6*tan(a6/180*pi);xf=1:5:160;yf=xf*tan(a6/180*pi)+b6;plot(xf,yf,'g',x6,y6,'+')axis([01600160])holdon（2）model:sets:jiaodu/1..6/:a,a1,b,k,s,c,x,y;bianhua/1..204/:t;zuobiao(jiaodu,bianhua):x1,y1;chongdie(jiaodu,jiaodu):cu;juli(bianhua,jiaodu,jiaodu):ju;jiao/1..3/:ja;endsets@for(bianhua(i):t(i)=5*i;);@for(jiaodu(i):@free(k(i));@free(s(i));@free(c(i)););@for(jiaodu(i):@free(a1(i));@free(b(i)););@for(jiaodu(i):@for(bianhua(j):@free(x1(i,j));@free(y1(i,j));););min=(@sum(jiaodu(i):a1(i)));@for(jiaodu(i):@bnd(-30,a1(i),30););@for(jiao(e):a1(6)=-3+e;@for(jiao(z):a1(3)=-2+z;@for(jiaodu(i):a(i)=a(i)+a1(i);1012211122邵仁和11k(i)=@tan(a(i)/180*3.14159);c(i)=@cos(a(i)/180*3.14159);b(i)=y(i)-x(i)*k(i););@for(jiaodu(i):@for(bianhua(j):x1(i,j)=x(i)+c(i)*