一个飞行管理问题数模竞赛

一个飞行管理问题

摘要

在某一空域里对飞机的飞行合理管理事关重大，比如乘客及机上工作人员生命财产安全和航空公司的运作效益等。本文通过对飞机飞行管理问题的研究，得到了调整飞机架数较少同时调整幅度均最小（平和最小）的飞行管理最优安排的非线性模型，这样既使得乘客所受影响达到最少，也便于飞机调整，还有利于飞机回到原来的航线，同时还在决策时间上对模型进行了优化和调整。

本文不仅一般性地将不相撞的问题转化为欧式距离控制，而且很巧妙的将不碰撞条件转化成简单的二次函数标准形式进行含参讨论，建立一个只含有转向角变量的模型。并且再次很妙的具体化区域受控时间形成矩阵，大大得简化运算，节约了大量运算的时间，便于管理人员控制操作，从而确保飞机的安全。更重要的是最后结合实际缩短了搜索区间，并优化算法，使得决策更加高效。最后的延时检验也充分体现了模型的可靠性。

关键字：欧氏距离约束转化缩短搜索区间时间矩阵延时检验

一、问题重述

在约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初步分析

根据问题容易知道，这显然是一个优化问题，当两架飞机可能发生碰撞时，即在规定区域某一时刻两架飞机之间的距离小于8公里，因此要调整飞行向一定角度，保证任意两架飞机在区域任意时刻，两者的距离均不小于8公里，避免相撞。考虑到调整角度应尽量小，可以简化飞行向调整策略，降低调整难度，同时减轻机乘客及工作人员的不适。此外由此初步确定了调整目标：所有六架飞机的飞行向调整角度均尽量小。

2解决案

由于所有飞机均处于1000米得高空作水平飞行，可将飞机飞行的空域视为二维平面xoy中的一个正形区域,顶点为(0,0),(160,0),(160,160),(0,160)。于是可以引入时间变量后，确定每架飞机在任意时刻的坐标，列出任意两点的欧氏距离，令其恒大于8公里，则得出一个重要约束条件。再结合变化角度应小于30度，即可得出约束条件，然后运用LINGO软件编辑程序进行求解。

为提高决策效率，在反复试验中又可对约束条件进行调整。

三、条件假设

1.不碰撞的标准为任意两架飞机的距离在以后任一个时间里大于8公里；

2.飞机飞行向角调整的幅度不应超过30度；

3.所有飞机飞行速度均为每小时800公里；

4.进入该区域的飞机在到达区域边缘时，与区域飞机的距离应在60公里以上。即在

计算如最优地调整各架（包括新进入的）飞机飞行的向角时，飞行管理中心得出合理的最优调整措施，；

5.最多需考虑6架飞机。

6.此处忽略飞机在执行过程中所需耗费的时间，即假设从飞机管理中心发出的调整信

息飞机马上可以接收并执行，不存在滞后或延迟；

7.飞行管理中心在计算飞行调整信号和发出信号所需时间，忽略各架飞机（包括刚进

入的飞机）调整航向前飞行数据的变化；

8.假定飞机在该区域完全依赖飞行管理中心调度；

9.假设飞机在飞出区域之后，飞行员可以自觉调整飞行策略，回归原始航线，即飞行

管理中心不必考虑飞机离开此区域后的状况。

四、符号说明

五、模型建立与求解

1问题简化

首先，如果对六架飞机在区域做实时监控，再做多次调整，则每作一次航向调整都要进行一次决策，这将使问题复杂化，总体计算量较大，同时实际问题中计算也要耗费

时间，效率大大降低，飞机控制的安全性必然会降低。并且对问题所给原始数据利用MATLAB软件（程序见附录1）作出原始航线图，如图1，可以粗略验证一次调整可行，既可以避免相撞，又简单易行。

图1

结论一：我们认为只做一次调整是优于多次调整的。

其次，分析飞机飞行轨迹，作如图2，

图2

某架飞机在点A的航向是由A到C，此时次飞机在点C可能遭碰撞，D点是一个安全位置，如果飞机在A点早做航向的调整，向角变化角为α；如果待飞机飞到B点时再做航向的调整，则向的变化角为β，很显然，β> α。由此我们得出结论，早调整航向优于晚调整航向，因此我们可将问题简化。

结论二：在新进入飞机到达区域边缘时便对所有飞机做航向调整，避免碰撞。

2模型的导出过程

1）模型推导

首先引入时间变量t,其中t的取值区间为（0，

T）；另外，为了提高程序运行效率，简

ij

化计算，我们引入一个矩阵来控制搜索循环，较合理的设定了变量的上下界，从而大大提高决策效率。

如图3所示，

T表示第i架飞机在规定区域可能飞行的最长时间，即该架飞机飞至区域

i

边界四个角上点所用时间的最大值，令

T=min{i T，j T}，则所有的ij T可以构成一个N*N

ij

的矩阵。此矩阵可以运用MATLAB软件求解，程序见附录2，计算结果见表1。