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

飞行管理优化模型

摘 要

本文建立了关于飞行管理问题的简洁数学模型。首先我们对新进入的飞机作出判断，通过模型Ⅰ给出了计算机模拟求解，看其是否会与其他六架飞机相碰，若会，则再次通过建型Ⅱ求出使各飞机安全通过区域应调整的方向角，对模型Ⅱ给出了非线性优化的具体算法。在模型改进中，我们对风速、人的反应时间及飞机的实际速度等对方向角的调整的影响做了简单的分析与评价，使得模型更易在实际应用中推广。

模型Ⅰ：针对问题一，建立了碰撞检测模型。首先，对已给数据进行分析，并利用VC 编程，模拟在6驾飞机都不改变飞行方向的条件下的飞行情况，结果是会相撞的。

模型Ⅱ：针对问题二，建立了多元非线性动态优化模型。在确保互不相撞的前提下，要使得飞机调整的角度尽可能小，满足Min i

f i β2

6

1∑

==，及最优解。运

用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 i

f i β2

6

1∑

==。

五、模型的建立与求解

模型一：碰撞检测模型

在平面直角坐标系中设出各飞机（包括新进入的飞机）的坐标分别为()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 0160

cos 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 架飞机都在正方形区域内， 第六架飞机与其他五架飞机的距离为

()[]()[]2

662666sin 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 ≤≤）