飞行管理_数学建模.docx

B题：飞行管理问题摘要：飞行管理问题是一个既现实又重耍的课题，本文利用偏转角度尽可能的小建 立两个非线性规划模型模型一：时间模型考虑到各架飞机的偏转角有正有负，在此模型中，对于 各架E机调整选取各个偏转角的绝对值的和作为目标函数，要求任意两架0机任 意时刻的距离大于8公里，则可以求出任意两架飞机的距离如此吋，两架飞机距离是时间/与各个飞机偏转角A0的函数，编写程序时将/离散化，且/有 最大值0. 2828 5 （沿对角线飞过的时间）,这样可得到表1T的结果:表1T飞机序号123456总和偏转角（度）002.540629001.0745263.615155模型二：闭塞区域模型在两架飞机中，将其中一架看成“静止”，另一架 相对于它而运动而以“静止”飞机为圆心，8如7为半径的圆形区域构成该飞机 的闭塞区域，任意一架飞机的方向角均不能在此区域内，则为不相撞为此，木 文用复变函数的知识表示各架飞机的速度，从而算出相对速度，再求出相对位移, 以相对速度与相对位移的夹角大丁•每两架飞机的临界夹角來刻画不相撞冃标函 数为每架飞机偏转角的平方和利用计算机编程得到表1-2的结果：表1-2飞机序号123456总和偏转角（度）002.5566001.06843.6250对于上述两个非线性规划，在理论方面，木文利用SUTM内点法（障碍函数法） 进行算法描述，在操作方面，分别利用lingo语言与MATLAB语言宜接编写程序 进行计算关键词：非线性规划、复变函数、SUMT内点法、闭塞区域、禁飞角一、问题重述1. 背景知识与其他交通工具相比，飞机以其速度快、安全舒适等特点在交通领域占据了 绝对地位。

而近年来飞机事故的频繁发生也预示着飞机存在一定的安全隐患经 调杳造成飞机相撞事故的原因主要是人、飞机（设备）、环境，而人的因素是事 故中通常起主体作用的因素，直接影响事故的发生和结局飞机事故的发生难以 预测口死亡率极高，所以航空安全机制的健全，航空人员素质的提高已变得刻不 容缓2•问题重述在约10000米的高空某边长为160公里的止方形区域内，经常有若干架飞机 作水平飞行区域内每架飞机的位置和速度向量均由计算机记录其数据，以便进 行飞行管理当一驾欲进入该区域的飞机到达区域边缘时，记录其数据后，要立 即计算并判断是否会与区域内的飞机发生碰撞如果会碰撞，则应计算如何调整 各架（包括新进入的）飞机飞行的方向角，以避免碰撞现假定条件如下：1） 不碰撞的标准为任意两架飞机的距离大于8公里；2） 飞机飞行方向角调整的幅度不应超过30度；3） 所有飞机飞行速度均为每小时800公里；4） 进入该区域的飞机在到达区域边缘时，与区域内飞机的距离应在60公里以上;5） 最多需考虑6架飞机；6） 不必考虑飞机离开此区域后的状况请你对这个避免碰撞的飞机管理问题建立数学模型，列岀计算步骤，对以下 数据进行计算（方向角谋差不超过0・01度），要求飞机飞行方向角调整的幅度 尽量小。

设该区域内4个顶点的坐标为（0, 0）, （160, 0）, （160, 160）, （0, 160）记录数据为：飞机编号横坐标X纵坐标y方向角（度）115040243285852363150155220.54145501595130150230新进入0052 1二、 基本假设1•飞机在规定区域内飞行，不会发生任何自然灾害；2. 飞行人员能够探知飞行区域内所有飞机的飞行轨迹，并能够及吋作出相应调 整；3. 飞机在规定区域内最多只调整一次，且在第六架飞机刚进入区域时调整；4•飞机在调整过程中是瞬间完成，飞行轨迹可看成直线；5. 不碰撞的标准是任何两架飞机的距离大于8如2 ；6. 飞机飞行方向的调整幅度不应超过30 ；7. 所有飞机调整角度总和最小为最优方案；&暂时不考虑飞机离开此区域的情况；9. 0机在调整之后按固定速度飞行；10. 飞机在飞行过程中可看为质点三、 参数说明1. （心,兀）表示笫i架飞机的初始位置；2. 0表示第i架飞机初始飞行角度；3. 表示第i架飞机的调整角度；4. 5..表示第，架飞机和第丿架飞机的距离；5. 表示第i架飞机和第/架飞机距离的平方；5. 岭表示笫i架飞机的飞行速度；6. /表示飞机的飞行吋间；7. 厶表示第/架飞机飞出飞行区域的时间；10. 0•表示第/架飞机和对于第j架飞机的飞行角度；11. Z&表示第R架飞机的速度矢量；12. &表示第k架飞机指向第/飞机的飞行矢量；13. 叫/•表示第R架飞机对第/架飞机的禁飞角的大小。

四、问题分析■第？^ 三■/第二架 _飞机四廟1 1 1 1 114012010080架60 新进入的020 40 60 30 100 120 140 16020图4-1：飞机初始状态从距离约束条件来说，在初始状态下，我们可以通过调整飞行方向，口的 使飞机在飞行区域内满足任意两架飞机的距离都大于8公里调整角度我们可以 假设已知，将上述初始角替换成调整角，通过计算，我们知道存在调整方法能够 使飞机安全飞出飞行区域，月•调整方式多样若我们假设调整角度对于每架飞机 来说，其满意度是一样的，那在这一系列调整方式中，必定存在调整角度之和达 到最小的调整方式在飞行过程中，初始位置和飞行方向对飞机的飞行路径有决 定性作用，而在速度已知的情况下，时间的确定就可以直接决定飞机的具体位置 由此，我们想到，如果取适当的一段时间为步长，确定检测时间点，若这些时间 点各架飞机都能满足不相撞的条件，那我们可以近似认为，飞机在这些飞行方向 上是符合安全飞行条件的，之后只需找出调整最小的方式即可从角度约束条件来说，将飞机看做质点，在飞行区域内，以飞机为圆心，4 公里为半径的圆形区域内是不能存在其他飞机的即每架飞机都冇其自身的闭寒 区域，而且是动态的闭塞区域。

在边长为160公里的正方形二维平面内，这6 个闭塞区域在运动过程中是不允许存在重合部分的飞机在飞行区域内是按照同 一方向飞行，也就是说在该正方形区域内，闭塞区域是按直线行走的运动初期, 各圆的位置已知，飞行方向可调整（先不考虑调整角度的限制问题），而调整的 廿的是为了让这些圆在以相同速度且沿直线的行走过程中不存在重合部分将模 型简化，仅看两架飞机的运动形式，以一个I员I的I员I心出发，引发出两条相对于另 一个圆的切线，而这两条切线在包含另一个圆的一侧所构成的夹域是禁飞域，夹 角为禁飞角o在这种解法中不仅要考虑飞机之间的相对运动速度还要考虑飞机之 间初始相对位置的具体关系讨论初期，我们确定飞机之间相对运动的速度还是引用原冇的坐标系，飞机Z间相对位置关系也在原冇的坐标系中讨论，该模型的 约束条件就是相对运动方向与相对位置方向之间的相对夹角只要满足不在禁飞 夹角中即可首先我们选择通过几何方法计算相对夹角，而在进一步计算过程中, 我们发现相对夹角若以角度制表示要分不同种情况讨论，过程复杂难以计算之 后我们又将模型进一步改进，对于坐标系的选取，若以飞机Z间相对位置所在向 量确定相对坐标则使约束条件简化。

在计算方而，为了避免角度的多情况彩响， 我们选取以复数形式表示角度，这样就将相对夹角的形式统一起来，计算也有了 进一步突破在计算出具体结果后，就需要对结果的有效性进行检验，而我们选取的检验 条件就是飞机以调整后的角度飞行，在离开飞行区域时，是否与飞行区域内的其 他飞机的距离大于60公里，若都满足大于60公里，则认为结呆可行性较强，若 有岀入，在对模型进行相应调整五、模型建立与求解模型一设飞机初始状态下的坐标为(“」)易知任意吋刻飞机的坐标为Xi = xi + vt cos(^- + Ag)Yi = yi + wsin(Q +△&,)则两架飞机的距离：d{j = 5%- =( x 厂 x2+(y厂岭y讨论飞机在规定区威内的飞行状态，则飞机冇飞行时间的限制：若飞机沿对 角线飞行，则飞行距离最远，飞行时间最长，此时/唤二型返，即t e(0,0.2828), max 800在飞行时间内，如杲任意两架飞机的距离都超过8公里，则认为不会发生碰撞事 件，即：Vr g(0,0.2828) di} >64 (i = l・・.6J = l・・・6) (1)首先讨论如果不改变飞机偏转角是否有两架飞机相撞，即满足(1)式的当v = 800公里/小时，(兀,兀)和&心知的情况下，经计算：5362 =(150+ 800/cos220.5“ —cos52丁 + (155 + 800fsin220.5“ —sin52『<=645562 = (130+ 800fcos230—cos52F+(150+ 800/sin230—sin52$ <= 64在t e(0,0.2828)内有实解，贝lj如果按照初始方向飞行会发生碰撞事件，需要做出相应调整。

为了确保飞行时不会发生碰撞事件，则任意叨架飞机的飞行距离不得小于 8公里，在此条件下寻求6架飞机能够安全通过飞行区域的最小调整方式而对 于“最小的调整方式”，本文采取下述两种理解方式1 •偏转总和最小因为每架飞机偏转角可止可负，所以，本模型取偏转角的 平方和作为目标函数结合上述分析与条件，可建立如下非线性规划模型6min工(AQ尸1=1dq >64 (z=b--6,j = l-・6)皿|旳<30 (心1・・6)虫(0,0.282*其中:2.毎架飞机偏转都最小min = Tdj> 64t g (0,0.2828)编程时，我们将/离散化，分为200等分,Hn . 0.2828八八八八八彳即 dt = = 0.00094300误差 ds = vdt = 0.7541如？（1）求得的最终结果为：△q=0； △&2=0； △03 =2.540629； A^4 = 0； △比=0； \06 = 1.074526 o(2) 求得的最终结果为：= 2.2906 ； 込=0.2500 ； = 2.5406 ； A% = 0 ； Kd, = 0.2500 ；△佻=1・0745对比上述结果，结合题中误差在0.01。

之内，可只取第一种方案，这样既满 足飞机偏转总和最小，又满足每架飞机偏转最小下面检验对于上述结杲是否满足“进入该区域的飞机在到达区域边缘时，与 区域内飞机的距离应在60公里以上；”飞机的飞行路径方程为:（y _ X）/（兀~xi）= tan（Q + AQ）；求得每架飞机飞出飞行区域的坐标分别为（129.6190,0）； （27.666&0）； （0,15.1227）；（0,105.6603）；（4」351,0）； （120」314,160）飞机飞岀飞行区域的时间为：_（兀•一兀J+（”•一儿尸— V求得每架飞机的飞行时间为：斤=0.0561 ； t2 = 0.1282 ； = 0.2564 ；t4 =0.1941； t5 = 0.2448 ； t6 = 0.2501飞机飞出飞行区域的先后顺序为：1, 2, 4, 5, 6, 3o按照此飞行时间顺序 由少到多重新排序而飞机飞出区域要保证和飞行区域内的其他飞机距离都人于60公里；所以：U1 -认））2 + （儿-儿S）2）> 3600 。