发明名称 一种适用于广域空中交通流量调控的方法
摘要 本发明公开了一种适用于广域空中交通流量调控的方法,针对广域空中交通流量的特点,根据协同进化的思想先将大问题划分成多个小问题,然后采用并行进化的方法同时对多个子种群进行进化,向不同的方向搜索,避免了陷入局部最优,而且改进了子种群之间信息交互的策略,实验证明能够找到比传统方法更优的解。本发明提出了一种并行进化中子种群之间的交换间隔动态变化的策略;在并行进化子种群环形拓扑结构交换的基础上提出了一种左右概率交换的策略,针对本问题采用了混合交换策略,能够得到较好的结果;本发明能够解决广域范围内的空中交通流量调控,比传统的遗传算法能够得到更为满意的结果。
申请公布号 CN103489336B 申请公布日期 2015.10.14
申请号 CN201310446753.0 申请日期 2013.09.26
申请人 北京航空航天大学 发明人 张学军;雷佳兴;管祥民
分类号 G08G5/00(2006.01)I 主分类号 G08G5/00(2006.01)I
代理机构 北京永创新实专利事务所 11121 代理人 官汉增
主权项 一种适用于广域空中交通流量调控的方法,其特征在于:假设,(1)所有的航班的飞行速度都一样,而且在飞行过程中保持不变;(2)所有航班的起飞时刻都是一个含有有限个元素的集合;(3)所有航班的飞行路径在一定的范围内可选;(4)同一个起止点之间的航班,可选路径集合是一样的;基于上述假设,所述方法包括如下步骤:第一步,数学模型建立,包括根据中国航路网模型建立约束条件和目标函数;第二步,采用合作型协同进化算法进行优化;所述的约束条件为,每个航班都包括两个变量(δ<sub>i</sub>,r<sub>i</sub>),δ<sub>i</sub>表示飞机起飞提前或延后的时间,r<sub>i</sub>表示重新选择的飞行路径,δ<sub>i</sub>和r<sub>i</sub>的可选集合表示为:Δ<sub>i</sub>=‑δ<sub>m</sub>,‑δ<sub>m</sub>+1,...,‑1,0,1,...,δ<sub>p</sub>‑1,δ<sub>p</sub>R<sub>i</sub>=r<sub>0</sub>,r<sub>1</sub>,r<sub>2</sub>,...,r<sub>max</sub>其中,i=1,2,…,N,N为航班数量,δ<sub>m</sub>,δ<sub>p</sub>分别表示航班提前和延后的最大时间,r<sub>0</sub>表示最优路径,r<sub>max</sub>表示最差路径;所述的目标函数包括两个,第一个目标函数表示为:<img file="FDA0000728771720000011.GIF" wi="889" he="209" />其中,<img file="FDA0000728771720000012.GIF" wi="160" he="133" />表示扇区S<sub>k</sub>在时间T内的总负荷,<img file="FDA0000728771720000013.GIF" wi="206" he="107" />表示扇区S<sub>k</sub>在时间T内的最大负荷,P表示扇区的数量,φ和<img file="FDA0000728771720000014.GIF" wi="53" he="54" />是权重系数;第二个目标函数表示为:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>y</mi><mn>2</mn></msub><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><mi>&delta;</mi><msup><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mn>2</mn></msup></mrow>]]></math><img file="FDA0000728771720000015.GIF" wi="325" he="167" /></maths>其中,如果航班i计划在t<sub>k</sub>时刻起飞,而实际上它在t<sub>n</sub>时刻起飞,这样航班的地面延误表示为:δ<sub>s</sub>(i)=t<sub>n</sub>‑t<sub>k</sub>,为了保证δ<sub>s</sub>(i)是正数,这里取绝对值,将δ<sub>s</sub>(i)表示为:δ<sub>s</sub>(i)=|t<sub>n</sub>‑t<sub>k</sub>|;航班的空中延误表示为:δ<sub>r</sub>(i)=3*(T<sub>r</sub>‑T<sub>0</sub>),T<sub>r</sub>表示实际飞行路径需要的时间,T<sub>0</sub>表示最短路径需要的时间;总的延误δ(i)表示为:δ(i)=δ<sub>s</sub>(i)+δ<sub>r</sub>(i)。
地址 100191 北京市海淀区学院路37号