2022年飞行管理数学建模 .pdf

《2022年飞行管理数学建模 .pdf》由会员分享，可在线阅读，更多相关《2022年飞行管理数学建模 .pdf（16页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

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

2、冲突的飞行问题，运用航向调整的方法解脱冲突，建立非线性数学模型61minii通过引入新的决策变量im、in，将原来的非线性模型转换成线性模型61mininmiiijijjjiinmnm2ijijjjiinmnm226/0pimi6/0pimj其中2iiim，2iiin。再运用 LINGO11编程求得该模型最优解为3.6326，第3架飞机的调整角为2.8419，第6架飞机（新进入的飞机）的调整角为0.7907，其余飞机不进行调整，从而给出了冲突解决方案。之后，本文对计算结果做出了分析和评价，同时还分析了滞后时间和转弯半径和限定在区域范围内对飞机航向调整的影响，使问题更符合实际情况。 在对模型进行

3、评价与分析的同时， 本文又对模型进行了推广， 对速度不同、 飞行高度不同的情况下进行了分析，并给出了合理的解释；增强了模型的实际应用意义。关键词： 飞行管理 碰撞冲突线性规划名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 1 页，共 16 页 - - - - - - - - - 一问题重述本题主要分析了在同一高度， 一定范围内的飞行管理问题。 在 10000米高空、边长为 160 公里的正方形区域内最多有6 架飞机做水平飞行， 其中飞机以每小时800 公里的速度匀速飞行。为了便于飞行

4、管理，在每一架飞机刚刚进入此飞行区域边界时，飞机的位置和速度均由计算机记录其数据，并立即进行计算和判断是否会与区域内的其它飞机发生潜在碰撞，如果存在碰撞危险， 就应该计算如何调整各架飞机的飞行方向角度，以避免发生碰撞。因此，根据题目的条件和假设，对该避免碰撞的飞机管理问题建立数学模型，列出计算步骤并对已经给出的数据进行计算。在保证进入该区域的飞机到达区域边缘时, 与区域内飞机的距离在60 公里以上的同时， 要求满足飞机飞行的方向角度调整幅度不能超过30 度，而且要尽量小。在该区域内建立直角坐标系， 4 个顶点的坐标为(0,0)， (160,0)， (160,160)，(0,160)。 记录数据

5、如下：飞机编号横坐标 x 纵坐标 y 方向角 (度) 1 150 140 243 2 85 85 236 3 150 155 220.5 4 145 50 159 5 130 150 230 新进入0 0 52 最后，对建立的模型进行评价和推广。二问题分析2.1、问题分析针对问题 1，一架新飞机飞入飞行区域边缘时，计算机进行记录其数据并进行判断该飞机是否和区域内的其它飞机存在潜在碰撞危险，由此作出分析，若存；在相互碰撞危险就要对各架飞机的方向角进行调整；否则维持原状， 不做任何变化。名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - -

6、 名师精心整理 - - - - - - - 第 2 页，共 16 页 - - - - - - - - - 对两飞机间距离的理解在本题中，根据题目假设所给的条件即在安全飞行过程中任意两架飞机间的距离要在 8 公里之上，这个距离远大于飞机本身的尺寸大小，所以，在以后的问题求解过程中我们将飞机当作质点看待。这样，根据计算机记录的数据可以知道每架飞机的坐标，根据两间距离公式即可得到两飞机之间的距离。对飞机方向角与调整角的理解根据题目得知， 飞机的方向角即为飞机的飞行方向与x 轴正向的夹角。 对于两架飞机 i 、j而言，要求飞机在飞行过程中相对速度的方向不能在飞机i 与飞机j的碰撞角（即两切线交角中指向

7、圆心的那个角）范围内。由于调整角以两飞机的连线为对称轴，左右由正负之分，故以调整角的绝对值最为目标函数。2.2 、问题 2 分析针对问题 2 的论述，由于在 2.1.2 中提出的绝对值目标函数为非线性函数，给求解带来了麻烦， 因此应用变量代换的数学方法将其转换成非线性函数，运用MATLAB 、LINGO 软件编程，调用求解函数得出最优调整方向角度。三模型的基本假设1. 新飞机进入边缘时， 立即做出计算， 每架飞机按照计算机计算后的指示立即作出方向角改变（有的飞机方向角可不变）。2. 每架飞机在新飞机刚进入到下一架飞机要进入前整个过程中最多只改变一次方向角。3. 忽略飞机转向角度， 即认为飞机在

8、接收到指令后立即对方向进行调整，且忽略调整时间。4. 新飞机进入该区域前 , 在区域中飞行的飞机方向已调合适，不会相撞。5. 对方向角的相同调整量的满意程度是一样的，且方向角调整越少， 满意程度越高。6. 飞机在安全飞行是两两之间的最短距离为8 公里，此距离远大于机身尺寸， 因此将飞机当作质点看待，即假设两飞机之间的距离为两点之间的距离。7. 最多考虑 6 架飞机。8. 飞机飞行方向角调整的幅度不应超过30 度。9. 所有飞机匀速飞行，速度为800hkm/。说明： 假设6中假设飞机为质点是为了将问题简化， 忽略次要问题以便得出合理的结论。假设7中6架飞机的假设是足够多的。 以世界上最繁忙的国际

9、航空港之一希思罗机场邻近区域为例， 因假设飞机在区域 160160作水平飞行，即知该区域内无机场。设在希思罗机场起降的飞机有一半穿过该区域，希思罗机场起驾总架次为22.5 万次，则平均每小时有 15架飞机穿过该区域。 而一架飞机穿过该区域最多需名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 3 页，共 16 页 - - - - - - - - - 216080028.0小时，则任一时刻该区域上空飞机架数的期望值不超过5.4架。另外事实上不同飞机的飞行高度是不同的，这就进一步减少了该区

10、域同一水平面上飞机的数目。以上讨论虽然略显粗略，但是足以证明6架飞机的假设是合理的。 假设9中假定所有飞机飞行速度均为800hkm/，是出于对问题的简化。我们将在模型的推广中给出飞机速度各不相同的方案。四. 符号说明ij：第 i 架飞机与第j架飞机的碰撞角，即两圆公切线交角中指向圆的那个角，规定jiij；ijv ：第 i 架飞机相对于第j架飞机的相对飞行速度；jir：第i架飞机与第j架飞机间距；ji：第i架飞机相对于第j架飞机的位置矢量与 x轴的夹角0ji：第i架飞机相对于第j架飞机的相对飞行速度与 x轴的夹角。ji：第i架飞机相对于第j架飞机的相对飞行速度与两架飞机圆心连线的交角，规定：以第

11、 i 架飞机为原点，ji连线从 i 指向j为正方向，逆时针旋转正，顺时针旋转为负；i：第 i 架飞机相对于直角坐标系旋转的角（即方向角改变量），为代数量；i：第 i 架飞机相对于第j架飞机ij的改变量；ip：第 i 架飞机的当前位置矢量；iv ：第i架飞机的当前速度矢量；t ：时间参数 ; i : 654321j:654321五. 模型的建立5.1 预测是否存在碰撞危险 潜在碰撞危险检测方案当相对速度矢量v与相对位置矢量p的矢量内积小于0时，即从v到p名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - -

12、- - - 第 4 页，共 16 页 - - - - - - - - - 夹角余弦值小于0时，两飞机存在潜在碰撞的危险，否则，无碰撞危险。具体证明如下：对任意两架飞机，用位置矢量和速度矢量分别表示为：11,vp、2,2vp，则经时间 t后两飞机的位置差为：vtptd（1）则)()min(min2vtpvtptdT（2）为了求得两飞机相撞时刻t，令：ttd=0得到 v v vpt（3）当t 0时，说明两飞机具有潜在碰撞的危险，而当t 0，即0 vp时，两飞机存在潜在碰撞冲突；反之， t 0if X(j)-X(i)0 beita(i,j)=atan(Y(j)-Y(i)/(X(j)-X(i);end

13、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;endif X(j)-X(i)=0 beita(i,j)=pi/2;endendif 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;endif 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名师资料总

14、结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 14 页，共 16 页 - - - - - - - - - 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

15、 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 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 15 页，共 16 页 - - - - - - - - - 0 0 0 0 0 0 ; arfa = 0 0.0941 0.5625 0.0889 0.3659 0.0390 0 0 0.0838 0.1154 0.1014 0.0666 0 0 0 0.0762 0.3985 0.0371 0 0 0 0 0.0792 0.0522 0 0 0 0 0 0.0403 0 0 0 0 0 0 ; enddataend名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 16 页，共 16 页 - - - - - - - - -