数学建模 飞机管理优化模型

飞行管理优化模型摘 要本文建立了关于飞行管理问题的简洁数学模型。

首先我们对新进入的飞机作出判断，通过模型Ⅰ给出了计算机模拟求解，看其是否会与其他六架飞机相碰，若会，则再次通过建型Ⅱ求出使各飞机安全通过区域应调整的方向角，对模型Ⅱ给出了非线性优化的具体算法。

在模型改进中，我们对风速、人的反应时间及飞机的实际速度等对方向角的调整的影响做了简单的分析与评价，使得模型更易在实际应用中推广。

模型Ⅰ：针对问题一，建立了碰撞检测模型。

首先，对已给数据进行分析，并利用VC 编程，模拟在6驾飞机都不改变飞行方向的条件下的飞行情况，结果是会相撞的。

模型Ⅱ：针对问题二，建立了多元非线性动态优化模型。

在确保互不相撞的前提下，要使得飞机调整的角度尽可能小，满足Min if i β261∑==，及最优解。

运用MATLAB 软件编程给出了具体算法。

第i 架飞机（i=1,2,……6）需调整角度分别为：0.0000 0.0000 2.0683 -0.4896 -0.0055 1.5611关键词：飞机碰撞 方向角调整 非线性优化 模拟仿真一、问题重述在约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轴方向的夹角。

试根据实际应用背景对你的模型进行评价与推广。

问题要求（1）判断飞机是否会碰撞；（2）若会碰撞，请给出方向角最优调整方案。

二、模型假设1、不碰撞的标准为任意两架飞机的距离大于8公里；2、飞机飞行方向角调整的幅度不超过30度；3、所有飞机在同一平面飞行，且速度均为每小时800公里；4、当飞机刚进入区域时就开始调整方向角，全过程只调整一次，且调整时间忽略不计；5、飞机在区域内飞行时不考虑其他一切外界因素的影响；6、相对正方形区域而言，将飞机作为质点考虑；7、最多需考虑6架飞机，不必考虑飞机离开此区域后的状况。

三、符号说明（i i y x ,）: 第i 架飞机的坐标（i=1，2，……6）单位：km ； （j x y x ,）：第j 架飞机的坐标（j=1，2，……6）单位：km ；v : 飞机飞行速度，其值为800 km/h ；i t : 第i 架飞机在区域内飞行时间 单位：h ； t : 飞机在区域内飞行的任一时刻 单位：h ；i α、j α : 第i 、j 架飞机的方向角 单位：°； i β、j β : 第i 、j 架飞机需调整的方向角 单位：°；ij d : 表示第i 架飞机与第j 架飞机之间的距离，单位：km 。

四、问题分析由上图可知：模型假设飞机在同一高度飞行，所以只需讨论飞机在二维平面的碰撞及调节管理问题。

以实现进入该区域的飞机都能安全通过该区域的目标。

当新飞入的飞机刚进入该正方形区域时，首先应通过模型Ⅰ对其作出判断，看其在区域内飞行时是否会与区域内的其他飞机发生碰撞。

并在二维坐标系中标出各飞机的位置坐标，结合其方向角借助计算机进行仿真模拟。

若通过模型Ⅰ验证不会碰撞，则使各飞机按原路线飞行即可安全通过正方形区域；否则，就需通过模型Ⅱ计算得出各飞机需要调整的角度，并且使各角度的调整方向角判断是否碰撞发生碰撞改变量尽可能小，即最优解。

在原方向角的基础之上调整计算得出的方向角，再用模型Ⅰ对其进行检验，提高结果的可信度，以确保各飞机能安全通过该正方形区域。

另外，为了使结果更加让人信服，在对时间的范围选取和最优解的处理问题上，我们进行了细致的分析。

（1）时间最小：i t 表示为第i 架飞机在区域内飞行的时间，j t 表示为第j 架飞机在区域内飞行的时间，取t=min （j i t t ,），即计算第i 架飞机与第j 架飞机间的距离时保证两架飞机都在区域内,且在t 时间内第i 架飞机不与第j 架飞机相撞，即保证所讨论范围在正方形区域内。

（2）平方和最小：由于所调角度i β∈（-30°，+30°），如果直接求和，正负将会抵消，所以目标函数定为：Min if i β261∑==。

五、模型的建立与求解模型一：碰撞检测模型在平面直角坐标系中设出各飞机（包括新进入的飞机）的坐标分别为()i i y x ,，方向角分别为i α()6,5,4,3,2,1=i ，则第i 架飞机t 时刻的位置为：()i i i i vt y vt x ααsin ,cos ++第i 架飞机在区域内时满足：160sin 0160cos 0≥≤+≤≤+≤t vt y vt x i i i i αα 求得i t t ≤≤0,i t 即为第i 架飞机在区域内的飞行时间。

所以，第6架飞机在区域内飞行的时间为6t 。

在t=min （6t ，i t ）时间内，第6架飞机与第i 架飞机都在正方形区域内， 第六架飞机与其他五架飞机的距离为()[]()[]2662666sin sin cos cos ααααvt y vt y vt x vt x d i i i i i +-+++-+=如果满足 86>i d （i=1,2……5，),min(0j i t t t ≤≤）即认为不会碰撞。

根据题目所给数据，标出各飞机在二维坐标中的位置（如图 5.1）。

其中左上角为坐标原点，横轴水平向右为x轴正方向，纵轴竖直向下为y轴正方向。

1、2、3、4、5、6分别为飞机编号，红色圆圈表示两架飞机互不相撞应满足的范围，即圆的半径为4km。

图5.1 飞机6进入该区域时各飞机位置在VC模拟环境下对其进行判断，根据程序（见附录I）模拟结果的截图可直观的看出新进入的编号为6的飞机会与编号为5的飞机（见图5.2）和编号为3的飞机（见图5.3）碰撞：图5.2 飞机5和飞机6碰撞图表5.1：飞机5和飞机6碰撞时各架飞机间距（km）图5.3飞机3和飞机6碰撞图表5.2：飞机3和飞机6碰撞时各架飞机间距（km）根据以上判断，新进入的飞机会与区域内的飞机发生碰撞，所以需对各飞机的方向角进行调整。

模型二：多元非线性优化模型经过模型Ⅰ的判断，如果会碰撞，则需计算如何调整各架飞机的方向角，以避免碰撞，并且使各架飞机调整的角度尽可能小，即求一个最优解。

设各飞机飞行方向的调整角度为约束变量，建立模型，由于i β∈（-30°，+30°），如果直接求和，正负将会抵消，所以目标函数定为：Min if i β261∑==约束条件：（1）按调整后的角度飞行，任意两架飞机在区域内的距离大于8公里，即 )([]()[]{}()[]()[]{}2sin sin 2cos cos j j vt j y i i vt i y j j vt j x i i vt i x ij dβαβαβαβα++-+++++-++=ij d 64> ()6,5,4,3,2,1,,=<j i j i ；（2）根据模型Ⅰ知，第i 、j 架飞机在区域内飞行的时间分别为：i t 、j t ；在t=min （i t ，j t ）时间内，第i 架飞机与第j 架飞机都在正方形区域内。

（3）飞机飞行最大调整角度i β、j β满足：30303030≤≤-≤≤-j i ββ故飞机调整角度最优解模型如下：Min if i β261∑==64 ij d （ 6,5,4,3,2,1,,=j i j i ，()j i t t t ,min 0≤≤）s.t3030≤≤-i β 3030≤≤-j β这是一个非线性优化问题，借助MATLAB 的fimincon()函数编程（见附录Ⅱ）计算得出最优解。

由所求结果，利用模型一对结果进行检验，用VC 模拟调整角度后的飞机的动态飞行路线 ，检验结果的截图见图5.4。

图5.4 调整角度后模拟飞行路线截图通过图5.4与图5.3比较，各飞机均可安全飞出正方形区域，完全符合安全飞行标准。

六、模型的灵敏度分析由于飞机的内部性能和飞行员的技术水平等一系列内部因素及不断变化的外部条件都可能对结果产生影响，因此讨论这些因素对最优解的影响是十分必要的。

（1）t ∆表示新飞如的飞机进入区域t ∆后各飞机开始调整方向。

则)([]()()[]()2sin sin sin sin 2cos cos cos cos ⎭⎬⎫⎩⎨⎧⎥⎦⎤⎢⎣⎡⎭⎬⎫⎩⎨⎧⎥⎦⎤⎢⎣⎡++∆+-++∆++++∆+-++∆+=j j vt jt jyi i vt it i y jj vt jt jx i i vt it i x ijdβααβααβααβαα64 ij d （ 6,5,4,3,2,1,,=j i j i ，()j i t t t ,min 0≤≤）s.t 3030≤≤-i β 3030≤≤-j β由以上条件求得最优解。

显然，t ∆越大，模型Ⅱ无解的可能性越大，即各飞机在新飞入的飞机进入区域后越晚调整，飞机在正方形区域内发生碰撞的概率越大。

（2）0v 表示区域内平均风速，其方向角为0α。

所以此时飞机飞行的水平速度为00cos cos ααv v i +，竖直速度为00sin sin ααv v i +，即飞机此时沿曲线飞行。

风速与飞机飞行方向的夹角越大，飞机的操纵性越差，即飞机在侧风中飞行时，发生安全事故的可能性越大。

七、模型的评价1、优点（1）本模型成功解决了飞行管理问题，建立了较为满意的模型，并且应用了仿真模拟系统很好地检验了结果的正确性，可信度较高；（2）模型建立时在时间t 的处理问题上，精确地将t 限定在区域内，即只考虑飞机在区域内的碰撞情况，所以能够保证取得最优解；（3）本模型的约束条件较少，计算本应很复杂，但使用了Matlab 中的fmincon 函数使得计算大大简化，且计算结果够精确。

2、缺点（1）模型是在理想条件下进行的，不考虑天气、风向等一切因素的影响，且设定飞行员接受命令到执行的时间为零，飞行调整的角度为绝对精确，因此可能与实际情况有一定差距；（2）考虑飞机是在同一平面内飞行，未考虑飞机间的垂直距离；（3）由于在计算过程中题目给的飞机飞行速度800公里/小时，而现实中不可能每架飞机在每一时刻的速度都为800km/h ，所以使得模型计算的结果与实际情况有所差别。