飞行管理数学建模.doc

1、摘要近年来，随着现代航空运输不断发展，为了维护航空器的航空秩序，保证飞机飞行安全，对于同一区域的飞行管理问题提出了要求。本文讨论了在一定区域范围内飞机飞行管理的最优化问题，通过建立数学模型计算求解，对飞机是否发生碰撞冲突进行预测，根据计算机求解结果对如何解脱冲突给出了较好的解决方法。对于飞机是否发生碰撞冲突问题，本文提出了基于飞机位置速度矢量关系的碰撞冲突检测方案，证明了只有位置差与速度差矢量内积小于零，即这样的航迹才存在潜在碰撞冲突,并根据安全飞行间隔规定,采用线性预测方法对冲突进行有效性确认,解决了飞机碰撞冲突检测的同时也避免了碰撞虚警问题。在此基础上，对于存在潜在碰撞冲突的飞行问题，运用

2、航向调整的方法解脱冲突，建立非线性数学模型通过引入新的决策变量、，将原来的非线性模型转换成线性模型 其中，。再运用LINGO11编程求得该模型最优解为，第3架飞机的调整角为，第6架飞机（新进入的飞机）的调整角为，其余飞机不进行调整，从而给出了冲突解决方案。之后，本文对计算结果做出了分析和评价，同时还分析了滞后时间和转弯半径和限定在区域范围内对飞机航向调整的影响，使问题更符合实际情况。在对模型进行评价与分析的同时，本文又对模型进行了推广，对速度不同、飞行高度不同的情况下进行了分析，并给出了合理的解释；增强了模型的实际应用意义。关键词：飞行管理 碰撞冲突 线性规划一问题重述本题主要分析了在同一高度

3、，一定范围内的飞行管理问题。在10000米高空、边长为160公里的正方形区域内最多有6架飞机做水平飞行，其中飞机以每小时800公里的速度匀速飞行。为了便于飞行管理，在每一架飞机刚刚进入此飞行区域边界时，飞机的位置和速度均由计算机记录其数据，并立即进行计算和判断是否会与区域内的其它飞机发生潜在碰撞，如果存在碰撞危险，就应该计算如何调整各架飞机的飞行方向角度，以避免发生碰撞。因此，根据题目的条件和假设，对该避免碰撞的飞机管理问题建立数学模型，列出计算步骤并对已经给出的数据进行计算。在保证进入该区域的飞机到达区域边缘时, 与区域内飞机的距离在 60 公里以上的同时，要求满足飞机飞行的方向角度调整幅度

4、不能超过30度，而且要尽量小。在该区域内建立直角坐标系，4个顶点的坐标为 (0,0)， (160,0)， (160,160)， (0,160)。 记录数据如下：飞机编号横坐标x纵坐标y方向角(度)1150140243285852363150155220.54145501595130150230新进入0052最后，对建立的模型进行评价和推广。二问题分析2.1、问题分析针对问题1，一架新飞机飞入飞行区域边缘时，计算机进行记录其数据并进行判断该飞机是否和区域内的其它飞机存在潜在碰撞危险，由此作出分析，若存；在相互碰撞危险就要对各架飞机的方向角进行调整；否则维持原状，不做任何变化。 对两飞机间距离的理

5、解在本题中，根据题目假设所给的条件即在安全飞行过程中任意两架飞机间的距离要在8公里之上，这个距离远大于飞机本身的尺寸大小，所以，在以后的问题求解过程中我们将飞机当作质点看待。这样，根据计算机记录的数据可以知道每架飞机的坐标，根据两间距离公式即可得到两飞机之间的距离。对飞机方向角与调整角的理解根据题目得知，飞机的方向角即为飞机的飞行方向与x轴正向的夹角。对于两架飞机、而言，要求飞机在飞行过程中相对速度的方向不能在飞机与飞机的碰撞角（即两切线交角中指向圆心的那个角）范围内。由于调整角以两飞机的连线为对称轴，左右由正负之分，故以调整角的绝对值最为目标函数。2.2、问题2分析 针对问题2的论述，由于在

6、2.1.2中提出的绝对值目标函数为非线性函数，给求解带来了麻烦，因此应用变量代换的数学方法将其转换成非线性函数，运用MATLAB、LINGO软件编程，调用求解函数得出最优调整方向角度。三模型的基本假设1.新飞机进入边缘时，立即做出计算，每架飞机按照计算机计算后的指示立即作出方向角改变（有的飞机方向角可不变）。2.每架飞机在新飞机刚进入到下一架飞机要进入前整个过程中最多只改变一次方向角。3.忽略飞机转向角度，即认为飞机在接收到指令后立即对方向进行调整，且忽略调整时间。4.新飞机进入该区域前,在区域中飞行的飞机方向已调合适，不会相撞。5.对方向角的相同调整量的满意程度是一样的，且方向角调整越少，满

7、意程度越高。6.飞机在安全飞行是两两之间的最短距离为8公里，此距离远大于机身尺寸，因此将飞机当作质点看待，即假设两飞机之间的距离为两点之间的距离。7.最多考虑6架飞机。8.飞机飞行方向角调整的幅度不应超过30度。9.所有飞机匀速飞行，速度为。说明： 假设6中假设飞机为质点是为了将问题简化，忽略次要问题以便得出合理的结论。假设7中6架飞机的假设是足够多的。以世界上最繁忙的国际航空港之一希思罗机场邻近区域为例，因假设飞机在区域160160作水平飞行，即知该区域内无机场。设在希思罗机场起降的飞机有一半穿过该区域，希思罗机场起驾总架次为22.5万次，则平均每小时有15架飞机穿过该区域。而一架飞机穿过该

8、区域最多需小时，则任一时刻该区域上空飞机架数的期望值不超过架。另外事实上不同飞机的飞行高度是不同的，这就进一步减少了该区域同一水平面上飞机的数目。以上讨论虽然略显粗略，但是足以证明架飞机的假设是合理的。 假设9中假定所有飞机飞行速度均为，是出于对问题的简化。我们将在模型的推广中给出飞机速度各不相同的方案。四.符号说明：第 架飞机与第架飞机的碰撞角，即两圆公切线交角中指向圆的那个角，规定；：第架飞机相对于第架飞机的相对飞行速度；：第架飞机与第架飞机间距；：第架飞机相对于第架飞机的位置矢量与x轴的夹角：第架飞机相对于第架飞机的相对飞行速度与x轴的夹角。：第架飞机相对于第架飞机的相对飞行速度与两架飞

9、机圆心连线的交角，规定：以第架飞机为原点，连线从指向为正方向，逆时针旋转正，顺时针旋转为负；：第架飞机相对于直角坐标系旋转的角（即方向角改变量），为代数量；：第架飞机相对于第架飞机的改变量；：第架飞机的当前位置矢量；：第架飞机的当前速度矢量；：时间参数;: :五.模型的建立5.1 预测是否存在碰撞危险 潜在碰撞危险检测方案当相对速度矢量与相对位置矢量的矢量内积小于时，即从到夹角余弦值小于时，两飞机存在潜在碰撞的危险，否则，无碰撞危险。具体证明如下：对任意两架飞机，用位置矢量和速度矢量分别表示为：、，则经时间后两飞机的位置差为： （1）则 （2）为了求得两飞机相撞时刻，令：=得到 （3）当0时，

10、说明两飞机具有潜在碰撞的危险，而当0，即时，两飞机存在潜在碰撞冲突；反之，0 if X(j)-X(i)0 beita(i,j)=atan(Y(j)-Y(i)/(X(j)-X(i); end if X(j)-X(i)0 beita(i,j)=atan(Y(j)-Y(i)/(X(j)-X(i); beita(i,j)=beita(i,j)+pi; end if X(j)-X(i)=0 beita(i,j)=pi/2; end end if Y(j)-Y(i)0 beita(i,j)=atan(Y(j)-Y(i)/(X(j)-X(i); beita(i,j)=beita(i,j)+2*pi; end

11、 if X(j)-X(i)theta(j) beita(i,j)=beita(i,j)*180/pi; A(i,j)=90+(theta(i)+theta(j)/2-beita(i,j); end if theta(i)arfa(i,j);(ming(i)-ning(i)+ming(j)-ning(j)/2+beita(i,j)0);for(xx:ning0); for(xx:ning3.14/6);for(xx:ming3.14/6);!数据;data:data:beita = 0 -1.2346 -2.2375 -2.7196 -0.1210 0.2526 0 0 -1.5502 -0.7373 -1.6110 0.1571 0 0 0 -2.9230 -1.0251 0.0063 0 0 0 0 0.1042 -0.0615 0 0 0 0 0 0.0334