




免费预览已结束,剩余12页可下载查看
下载本文档
版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
飞 行 管 理 问 题的非线性优化模型摘要本文研究了飞行管理问题的非线性优化模型。在约10000米高空的某边长160公里的正方形区域内,经常有若干架飞机作水平飞行。区域内每架飞机的位置和速度向量均由计算机记录其数据,以便进行飞行管理。当一架欲进入该区域的飞机到达区域边缘时,记录其数据后,要立即计算并判断是否会与区域内的其它飞机发生相撞。如果发生相撞,则应计算如何调整各架(包括新进入的)飞机的飞行方向角,以避免碰撞。借鉴物理学中质点的概念,我们在空域平面建立了二维直角坐标系;我们讨论飞机不相撞的的条件时,将每一个飞机看作半径4km的圆,任意两架飞机不相撞的条件为两圆外切;进一步讨论约束条件,我们利用点的速度合成定理,任意两架飞机不相撞的条件化为两架飞机的相对速度和相对位移的夹角不小于两飞机相撞的临界角;模型求解过程中,我们结合多种数学软件(AutoCAD、MATLAB、SPSS、lingo)的特点,简化了算法。用Lingo数学软件,我们求得了最优解,各个飞机的调整角均不大于17,调整幅度平方和为9.46 ,检验后比较符合实际情况,每架飞机的初始方向角和调整角及调整幅度如下表:表1 每架飞机的初始方向角和调整角及调整幅度初始方向角4.2410004.1189003.8483002.7750004.0141000.9075000调整后的方向角4.2409994.1188993.6191952.8583343.7226261.051778调整幅度-0.0000-0.0000-0.22910.0833-0.29150.1443调整的角度00-13.134.77-16.78.27飞行管理问题的非线性优化模型一、问题的重述 在约10000米高空的某边长160公里的正方形区域内,经常有若干架飞机作水平飞行。区域内每架飞机的位置和速度向量均由计算机记录其数据,以便进行飞行管理。当一架欲进入该区域的飞机到达区域边缘时,记录其数据后,要立即计算并判断是否会与区域内的其它飞机发生相撞。如果发生相撞,则应计算如何调整各架(包括新进入的)飞机的飞行方向角,以避免碰撞。现假设条件如下:1 不相撞的标准为任意两架飞机的距离大于8公里;2 飞机飞行方向角调整的幅度不应超过30度;3 所有飞机的飞行速度均为每小时800公里;4 进入该区域的飞机在到达区域边缘时,与区域内飞机的距离应在60公里以上;5 最多需考虑6架飞机;6 不必考虑飞机离开此区域后的情况。请你对这个避免碰撞的飞行管理问题建立数学模型,列出计算步骤,对以下数据进行计算(方向角误差不超过0.01度),要求飞机飞行方向角调整的幅度尽量小。二、符号说明符号含义第i架飞机的初始方向角第i架飞机调整后的方向角第i架飞机与第j架飞机初始位置的距离(,)第i架飞机的初始位置(,)调整后第i架飞机相对于第j架飞机的速度(,)调整后第i架飞机相对于第j架飞机的位移调整后第i架飞机相对于第j架飞机的速度与位移之间的夹角调整后第i架飞机与第j架飞机不相撞的临界角备注:(1)i=1,2,,6 (2)j=1,2,,6 (3)ij (4)角度均采用弧度值三、模型的假设(1)不相撞的标准为任意两架飞机的距离大于8公里;(2)飞机飞行方向角调整的幅度不应超过30度;(3)所有飞机的飞行速度均为每小时800公里;(4)进入该区域的飞机在到达区域边缘时,与区域内飞机的距离应在60公里以上;(5)最多需考虑6架飞机;(6)不必考虑飞机离开此区域后的情况;(7)计算机从记录新进入飞机数据到给飞机发指令的间隔时间忽略不计;(8)新飞机进入时,每架飞机(包括新飞机)立即改变方向角,且至多改变一次;(9)新飞机进入空域时,在空域中飞行的飞机(包括新飞机)方向已调合适,之后不会相撞。四、模型的建立1、模型分析该模型是基于假设条件的非线性规划模型(参考文献一)。飞机飞行的整个空域可以看成是一个二维平面,我们可以建立直角坐标系,顶点是(0,0),(160,0),(160,160),(0,160),飞机飞行的方向角是飞行方向与x轴正方向的夹角。我们可以用向量、分别表示横纵坐标,则V可以表示成V=,S可以表示成S=。每架飞机的位置都是时间T的函数,直接研究比较困难,为此我们引入相对运动。我们根据点的速度合成定理(参考文献二):动点在某瞬时的绝对速度等于它在该瞬时的牵连速度与相对速度的矢量和,即动点的相对速度可以由绝对速度与牵连速度的矢量差来确定。 (是绝对速度,是相对速度,是相对速度)。任意两架飞机不相撞的条件是其距离大于8千米,因此我们考虑将每架飞机视为半径4千米的圆心,在相对运动的情形下,两架飞机不相撞的条件就为两圆相离,不相撞的临界条件是在它们之间的距离最小的时候,两圆相切。如图1: yB C A x 图1飞机相对飞机运动,图中虚线所示的位置既是与不相撞的临界位置。这里我用到了解析几何中的向量夹角定理(参考文献三)求解两个向量之间的夹角:COS(A,B)=(x1*x2+y1*y2)/(*) (其中A=x1,y1,B=x2,y2)要使两架飞机不相撞,我们只需满足=。由飞机飞行角度调整幅度不超过30度(即/6),我们得到约束:|-|= (i=1,2,,6 j=1,2,,6 ij) (1) |-|=/6 (i=1,2,,6) (2)0=,同时,调整后的方向角与初始方向角的差值不超过/6,即|-|=/6。(2)lingo编程源代码见附录2(3)由lingo编程即可得最优解(见附录3),运行结果整理后如下表6 调整前后的方向角及调整幅度(弧度制)飞机的编号123456初始方向角4.2410004.1189003.8483002.7750004.0141000.9075000调整后的方向角4.2409994.1188993.6191952.8583343.7226261.051778调整幅度-0.0000-0.0000-0.22910.0833-0.29150.1443调整的角度00-13.134.77-16.78.27总的平方代价为:0.1652067(弧度制)。六、模型的检验及分析1、模型的检验我们对调整前的方向角度与调整的角度进行比较:初始方向角24323522.515923052调整角度00-13.1394.77-16.78.27用SPSS做出调整角度的变化曲线,见图2图2从结果中我们就可以看出,第一架和第二架调整的幅度为0,其它几架飞机调整的幅度=G(i,j);!约束相对飞行角度大于临界角(即两圆恰好相切的角度); for(row(i):for(col(j)|i#EQ#j:H(i,j)=0); for(row(i):for(col(j)|i#EQ#j:V0(i,j)=0); for(row(i):for(col(j)|i#EQ#j:V1(i,j)=0); for(row(i):for(col(j):free(V1(i,j);!取消非零约束; for(row(i):for(col(j):free(V0(i,j);!取消非零约束; for(row(i):for(col(j):free(X(i,j);!取消非零约束; for(row(i):for(col(j):free(Y(i,j);!取消非零约束; for(col:free(V00);!取消非零约束; for(col:free(V11);!取消非零约束; data: A=4.2410 4.1189 3.8483 2.7750 4.0141 0.9075;!初始方向角; X = 0 65 0 5 20 150 -65 0 -65 -60 -45 85 0 65 0 5 20 150 -5 60 -5 0 15 145 -20 45 -20 -15 0 130 -150 -85 -150 -145 -130 0; Y = 0 55 -15 90 -10 140 -55 0 -70 35 -65 85 15 70 0 105 5 155 -90 -35 -105 0 -100 50 10 65 -5 100 0 150 -140 -85 -155 -50 -150 0; G = 0 0.0941 0.5625 0.0889 0.3659 0.0390 0.0941 0 0.0838 0.1154 0.1014 0.0666 0.5625 0.0838 0 0.0762 0.3985 0.0371 0.0889 0.1154 0.0762 0 0.0792 0.0522 0.3659 0.1014 0.3985 0.0792 0 0.0403 0.0390 0.0666 0.0371 0.0522 0.0403 0; enddataend附录3、lingo结果 Local optimal solution found at iteration: 915 Objective value: 0.1652067 Variable Value Reduced Cost A( 1) 4.241000 0.000000 A( 2) 4.118900 0.000000 A( 3) 3.848300 0.000000 A( 4) 2.775000 0.000000 A( 5) 4.014100 0.000000 A( 6) 0.9075000 0.000000 B( 1) 4.240999 0.000000 B( 2) 4.118899 0.000000 B( 3) 3.619195 0.000000 B( 4) 2.858334 0.000000 B( 5) 3.722626 0.000000 B( 6) 1.051778 0.000000 V00( 1) -363.3001 0.000000 V00( 2) -447.4061 0.000000 V00( 3) -710.4796 0.000000 V00( 4) -768.1198 0.000000 V00( 5) -668.7167 0.000000 V00( 6) 396.8225 0.000000 V11( 1) -712.7503 0.000000 V11( 2) -663.1951 0.000000 V11( 3) -367.7210 0.000000 V11( 4) 223.5889 0.000000 V11( 5) -439.1105 0.000000 V11( 6) 694.6451 0.000000 G( 1, 1) 0.000000 0.000000 G( 1, 2) 0.9410000E-01 0.000000 G( 1, 3) 0.5625000 0.000000 G( 1, 4) 0.8890000E-01 0.000000 G( 1, 5) 0.3659000 0.000000 G( 1, 6) 0.3900000E-01 0.000000 G( 2, 1) 0.9410000E-01 0.000000 G( 2, 2) 0.000000 0.000000 G( 2, 3) 0.8380000E-01 0.000000 G( 2, 4) 0.1154000 0.000000 G( 2, 5) 0.1014000 0.000000 G( 2, 6) 0.6660000E-01 0.000000 G( 3, 1) 0.5625000 0.000000 G( 3, 2) 0.8380000E-01 0.000000 G( 3, 3) 0.000000 0.000000 G( 3, 4) 0.7620000E-01 0.000000 G( 3, 5) 0.3985000 0.000000 G( 3, 6) 0.3710000E-01 0.000000 G( 4, 1) 0.8890000E-01 0.000000 G( 4, 2) 0.1154000 0.000000 G( 4, 3) 0.7620000E-01 0.000000 G( 4, 4) 0.000000 0.000000 G( 4, 5) 0.7920000E-01 0.000000 G( 4, 6) 0.5220000E-01 0.000000 G( 5, 1) 0.3659000 0.000000 G( 5, 2) 0.1014000 0.000000 G( 5, 3) 0.3985000 0.000000 G( 5, 4) 0.7920000E-01 0.000000 G( 5, 5) 0.000000 0.000000 G( 5, 6) 0.4030000E-01 0.000000 G( 6, 1) 0.3900000E-01 0.000000 G( 6, 2) 0.6660000E-01 0.000000 G( 6, 3) 0.3710000E-01 0.000000 G( 6, 4) 0.5220000E-01 0.000000 G( 6, 5) 0.4030000E-01 0.000000 G( 6, 6) 0.000000 0.000000 H( 1, 1) 0.000000 0.000000 H( 1, 2) 1.906896 0.000000 H( 1, 3) 2.353088 0.000000 H( 1, 4) 2.678020 0.000000 H( 1, 5) 2.874664 0.000000 H( 1, 6) 0.3246630 0.000000 H( 2, 1) 1.906896 0.000000 H( 2, 2) 0.000000 0.000000 H( 2, 3) 1.665760 0.000000 H( 2, 4) 2.445895 0.000000 H( 2, 5) 1.756878 0.000000 H( 2, 6) 0.2291439 0.000000 H( 3, 1) 2.353088 0.000000 H( 3, 2) 1.665760 0.000000 H( 3, 3) 0.000000 0.000000 H( 3, 4) 2.996838 0.000000 H( 3, 5) 1.286457 0.000000 H( 3, 6) 0.3710000E-01 0.000000 H( 4, 1) 2.678020 0.000000 H( 4, 2) 2.445895 0.000000 H( 4, 3) 2.996838 0.000000 H( 4, 4) 0.000000 0.000000 H( 4, 5) 3.141590 -1.485291 H( 4, 6) 0.5220000E-01 0.000000 H( 5, 1) 2.874664 0.000000 H( 5, 2) 1.756878 0.000000 H( 5, 3) 1.286457 0.000000 H( 5, 4) 3.141590 0.000000 H( 5, 5) 0.000000 0.000000 H( 5, 6) 0.4030000E-01 0.000000 H( 6, 1) 0.3246630 0.000000 H( 6, 2) 0.2291439 0.000000 H( 6, 3) 0.3710000E-01 0.000000 H( 6, 4) 0.5220000E-01 0.000000 H( 6, 5) 0.4030000E-01 0.000000 H( 6, 6) 0.000000 0.000000 V0( 1, 1) 0.000000 0.000000 V0( 1, 2) 84.10597 0.000000 V0( 1, 3) 347.1795 0.000000 V0( 1, 4) 404.8197 0.000000 V0( 1, 5) 305.4165 0.000000 V0( 1, 6) -760.1227 0.000000 V0( 2, 1) -84.10597 0.000000 V0( 2, 2) 0.000000 0.000000 V0( 2, 3) 263.0735 0.000000 V0( 2, 4) 320.7137 0.000000 V0( 2, 5) 221.3106 0.000000 V0( 2, 6) -844.2286 0.000000 V0( 3, 1) -347.1795 0.000000 V0( 3, 2) -263.0735 0.000000 V0( 3, 3) 0.000000 0.000000 V0( 3, 4) 57.64018 0.000000 V0( 3, 5) -41.76292 0.000000 V0( 3, 6) -1107.302 0.000000 V0( 4, 1) -404.8197 0.000000 V0( 4, 2) -320.7137 0.000000 V0( 4, 3) -57.64018 0.000000 V0( 4, 4) 0.000000 0.000000 V0( 4, 5) -99.40311 0.000000 V0( 4, 6) -1164.942 0.000000 V0( 5, 1) -305.4165 0.000000 V0( 5, 2) -221.3106 0.000000 V0( 5, 3) 41.76292 0.000000 V0( 5, 4) 99.40311 0.000000 V0( 5, 5) 0.000000 0.000000 V0( 5, 6) -1065.539 0.000000 V0( 6, 1) 760.1227 0.000000 V0( 6, 2)
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 公司协议书范本格式字体
- 精细化管理班组施工安全责任合同范本
- 烧烤店试用期配菜员烧烤技艺聘用合同
- 留学行前准备协议:出国留学全方位辅导服务
- 科技园区场地租赁与研发创新合作合同
- 场地监管与廉政责任落实合同
- 超声回弹综合法测强技术
- 溶洞区桩基注浆压力调控
- 大数据管理能力评估和数据工程资料分享
- 医院重点部位消防安全培训
- 2025年吉林省国资委出资企业招聘笔试参考题库含答案解析
- 2025年全国安全生产月安全生产知识竞赛考试题库及答案(共四套)
- 基于MATLABsimulink同步发电机突然三相短路仿真
- 《标准的制定》课件
- 国土空间规划环评培训
- 风动凿岩机操作规程(4篇)
- 四川省成都市九县区2023-2024学年高一下学期期末调研考试化学试题(解析版)
- (完整版)python学习课件
- 联塑管材检验报告模板
- 高钠血症护理查房
- 小学数学练习设计的有效性研究结题报告
评论
0/150
提交评论