资源描述：
《数学建模竞赛案例选讲.ppt》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、数学建模竞赛案例选讲飞行管理问题1995年A题飞行管理问题（1995年全国大学生数学建模竞赛试题A）问题：在约10000米高空的正方形区域内，有若干架飞机作水平飞行。区域内每架飞机的位置和速度向量均由计算机记录数据，以便进行飞行管理，当一架欲进入该区域的飞机到达区域边缘时，记录其数据后，要立即计算并判别是否会与区域内的飞机发生碰撞。如果会碰撞，则应计算如何调整各架（包括新进入的）飞机飞行的方向角，以避免碰撞。假定条件:1.不碰撞的标准是任意两架飞机的距离大于8km；2.飞机飞行方向角调整幅度不应超过30度，而要尽可能小；3.所有飞机的飞行速度为800km/h，
2、不受其他因素影响；4.进入该区域的飞机在到达边缘时，与该区域内的飞机的距离应在60km以上；5.不考虑飞机离开区域后的情况；6.建模时暂考虑6架飞机；问题的提出：请你对这个避免碰撞的飞行管理问题建立数学模型，列出计算步骤，对以下数据进行计算（方向角误差不超过0.01度），要求飞机的方向角调整的幅度尽量小。设该区域4个顶点的坐标为(0，0),(160，0),(160，160),(0，160），记录数据如下表（其中方向角指飞行方向与x轴正向夹角)飞机编号横坐标x纵坐标y方向角度1150140243285852363150155220.541455015951301
3、50230新进入0052问题分析根据题目的条件，可将飞机飞行的空域视为二维平面xoy中的一个正方形，顶点在(0，0),(160，0),(160，160),(0，160）。各架飞机的飞行方向角为飞行方向与x轴正向夹角（转角）。根据两飞机不碰撞的标准为二者距离大于8km，可将每架飞机视为一个以飞机为圆心、以4为半径的圆状物体（每架飞机在空域中的状态由圆心的位置矢量和飞行速度矢量确定）。这样两架飞机是否碰撞就化为两圆在运动中是否相交的问题。两圆是否相交只要讨论它们的相对运动即可。CABDnimilijαijαijαijβijvijij建模时补充假定条件:1.飞机在所
4、定区域内作直线飞行，不偏离航向；2.飞机管理系统内不发生意外，如发动机失灵，或其他意外原因迫使飞机改变航向；4.飞机管理系统发出的指令应被飞机立即执行，即认为转向是瞬间完成的（忽略飞机转向的影响，即转弯半径和转弯时间的影响）；3.飞机进入区域边缘时，立即作出计算，每架飞机按照计算后的指示立即作方向角改变；5.每架飞机在在整个过程中指点改变一次方向6.新飞机进入区域时，已在区域内部的飞机的飞行方向已调整合适，不会碰撞；7.对每架飞机方向角的相同调整量的满意程度是一样的。模型的建立（1）圆状模型ABDnimilijαijαijαijβijvijij采用相对速度作为
5、研究对象符号说明：i,j——第i，第j架飞机的圆心；αij——第i，第j架飞机的碰撞角，αij=αji;vij——第i架飞机相对第j架飞机的相对飞行速度;lij——第i，第j架飞机的圆心距θi——第i架飞机的飞行方向与x轴正向夹角（逆时针为正）xi——第i架飞机的位置矢量vi——第i架飞机的的速度矢量βij——第i飞机对第j架飞机的相对速度与两架飞机圆心连线的夹角（逆时针为正）不碰撞
6、βij
7、>αij(2)由圆状模型导出的方程讨论βij的改变量与第i第j两架飞机飞行方向角改变量Δθi,Δθj的关系由题目条件知
8、vi
9、=A=800，可用复数表示速度设第i，j飞机
10、飞行方向改变前的速度分别为改变后的速度分别为改变前后相对速度分别为两者之商的幅角就是Δβij定理：对第i，第j两架飞机，其相对速度方向βij的改变量Δβij等于两飞机飞行方向角改变量之和的一半，即模型目标函数：Minθ其中θ为各飞机方向角调整量的最大值或为约束条件：调整方向角时不能超过300：调整飞行方向后飞机不能碰撞：模型为或为化为线性规划模型由于Δθi可正可负，为使各变量均非负，引入新变量：模型化为模型求解αij的计算model:sets:plane/1..6/:x0,y0;link(plane,plane):alpha,sin2;endsets@for(
11、link(i,j)
12、i#ne#j:sin2(i,j)=64/((x0(i)-x0(j))^2+(y0(i)-y0(j))^2););@for(link(i,j)
13、i#ne#j:(@sin(alpha*3.14159265/180.0))^2=sin2;);data:x0=150,85,150,145,130,0;y0=140,85,155,50,150,0;enddataendALPHA(1,1)1.234568ALPHA(1,2)5.391190ALPHA(1,3)752.2310ALPHA(1,4)5.091816ALPHA(1,5)2000.963ALP
14、HA(1,6)2.234507ALPH
