运筹学课程设计-一个飞行管理问题

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

一个飞行管理问题模型2建立

在约10000m高空的某边长160km的正方形区域内,经常有若干架飞机作水平飞行,区域内每架飞机的位置和速度向量均由计算机记录其数据,以便进行飞行管理.当一架欲进入该区域的飞机到达边界区域边缘时, 记录其数据后,要立即计算并判断是否会与其区域内的飞机发生碰撞.如果会碰撞,则应计算如何调整各架(包括新进入的)飞机飞行的方向角,以避免碰撞.现假设条件如下:

1) 不碰撞的标准为任意两架飞机的距离大于8km;

2)飞机飞行方向角调整的幅度不应超过30度;

3)所有飞机飞行速度均为每小时为800km;

4)进入该区域的飞机在到达区域边缘时,与区域内飞机的距离应

在60km以上;

5)最多考虑6架飞机;

6)不必考虑飞机离开此区域后的状况;

对以下数据进行计算(方向角误差不超过0.01度),要求飞机飞行方向角调整的幅度尽量小.

飞机编号横坐标纵坐标方向角(度)

1 150 140 243

2 85 85 236

3 150 155 220.5

4 14

5 50 159

5 130 150 230

新进入0 0 52

注: 方向角指飞行方向与x轴正向的夹角

一、分析

在模型1的非线形规划模型是比较困难的,输入之后也很可能无权找不到好的解甚至出现错误,此外,演示软件还会受到求解规模的限制,尤其可能无权使用全局求解程序。因此,如果能把这个问题简单化成简单的规划模型,将是非常有价值的。下面就建立模型二。

二、模型的建立 如图1,把两架飞机i,j 分别看成半径为4km 的圆(图中i,j 为圆心),AB ,CD 为公切线,将AB 和CD 的夹角的一半称为碰撞角。在调整时刻,设第i 架飞机与第j 架飞机的碰撞角为,则易知

其中r ij 为当前这两架飞机连线的长度(距离)。

因为飞机的距离大于8km 就不会相撞,所以如果这两个圆随着时间推移不相交就可以了。为此,考虑第i 架飞机相对于第j 架飞机的相对速度(矢量,图中记为Vij )是比较方便的,因为相对速度的大小和方向在飞机飞行中会始终保持不变(除非调整飞行角度)。

社Bij 为调整前的相对速度Vij 与这架飞机连线(从i 指向j 的矢量)的夹角(以连线矢量为基准,逆时针方向为正,顺时针为负),则Bij=-Bji 。具体来说,应该计算如下:

注意:标准的反正切函数的符号是arctan,返回主值;我们这里使用arctan2表示一个特殊的返回四象限辐角的反正切函数,即arctan2(b,a)返回向量

,

8

sin 1

ij

ij r -=α

(a,b )的-π到π之间的辐角(或者是返回0到 2π之间也是可以的)。即使这样,也还能完全满足要求,因为这样的到ij β

图1:第i 架飞机与第j 架飞机的碰撞角

取值位于-2π到2π之间的,还需要将它转换到-π到π之间才行(超过π时就剪去2π,小于-π就加上2π)。

从图中可以看出(注意图中的两条辅助线in//CD 、im//AB ),两架飞机i ,j 不相撞的冲要条件是(实际上不只是在所考虑的区域内不相交,而是永远不相交)

|ij β|≥ij α

如果调整前这个关系式成立,则不需要调整。否则,仍用i θ∆表示第i 架飞机飞

行方向角的调整量,并记由此引起的ij β的改变量∆ij β。现在,问题的关键是如何弄清∆ij β如何和

i θ∆变化。可以证明

下面利用复数的知识证明上式

证明

由题知| i v |=800km=A.设改变前的速度分别为'

i

v = i i e θA ,

'

j i j e θv =A ,改变方向后速度分别为

2

i

v =A ()

i i i e

θθ+∆ ,

2

j

v =A ()

j j i e

θθ+∆

改变前相对速度

ij v ='i v -'j v =A (i i e θ-j

j e θ

),

改变后相对速度

'

ij v =2

i v -2

i v =A [()i

i

i e θθ+∆- ()j

j

i e θθ+∆],

'

ij

ij v v =

()

()

[]

[]

j j i i j i

i i i i A e

e

A e e θθθθθθ+∆+∆--

= (

)

2

sin

2sin 2

i j

i i j j

i i j

e

θθθθθθθθ∆+∆+∆--∆-

即'ij v 和ij v 辐角相差2

i j

θθ∆+∆.

因此,可以得到如下的数学规划模型:

.

6,,1,

30,,6,,1,,

)(..216

1

2

=≤∆≠=≥∆+∆+∑=i j i j i t s Min i ij j i ij i i

θαθθβθ

参考文献:

[1] 谢金星薛毅,《优化建模与LINDO_LINGO软件》,北京:清华大学出版社,2005年

[2] 姜启源谢金星叶俊,《数学模型》(第三版),北京:高等教育出版社,2003年

[3] 朱道元,《数学建模案例精选》,北京:科学出版社,2003年

相关文档
最新文档