CN107504972A - 一种基于鸽群算法的飞行器航迹规划方法及装置 - Google Patents
一种基于鸽群算法的飞行器航迹规划方法及装置 Download PDFInfo
- Publication number
- CN107504972A CN107504972A CN201710625878.8A CN201710625878A CN107504972A CN 107504972 A CN107504972 A CN 107504972A CN 201710625878 A CN201710625878 A CN 201710625878A CN 107504972 A CN107504972 A CN 107504972A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- path
- aircraft
- algorithm
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 21
- 238000005457 optimization Methods 0.000 claims abstract description 19
- 230000008859 change Effects 0.000 claims description 12
- 230000006870 function Effects 0.000 claims description 12
- 238000005516 engineering process Methods 0.000 abstract description 5
- 241000272201 Columbiformes Species 0.000 description 26
- 230000003044 adaptive effect Effects 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 230000005389 magnetism Effects 0.000 description 3
- 230000002068 genetic effect Effects 0.000 description 2
- 238000009434 installation Methods 0.000 description 2
- 239000011664 nicotinic acid Substances 0.000 description 2
- 230000006978 adaptation Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 210000004556 brain Anatomy 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000012804 iterative process Methods 0.000 description 1
- 230000004899 motility Effects 0.000 description 1
- 230000007935 neutral effect Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
- 239000013598 vector Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C21/00—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
- G01C21/20—Instruments for performing navigational calculations
-
- G—PHYSICS
- G08—SIGNALLING
- G08G—TRAFFIC CONTROL SYSTEMS
- G08G5/00—Traffic control systems for aircraft, e.g. air-traffic control [ATC]
- G08G5/003—Flight plan management
- G08G5/0034—Assembly of a flight plan
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05D—SYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
- G05D1/00—Control of position, course, altitude or attitude of land, water, air or space vehicles, e.g. using automatic pilots
- G05D1/10—Simultaneous control of position or course in three dimensions
- G05D1/101—Simultaneous control of position or course in three dimensions specially adapted for aircraft
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/004—Artificial life, i.e. computing arrangements simulating life
- G06N3/006—Artificial life, i.e. computing arrangements simulating life based on simulated virtual individual or collective life forms, e.g. social simulations or particle swarm optimisation [PSO]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
- G06Q10/047—Optimisation of routes or paths, e.g. travelling salesman problem
-
- G—PHYSICS
- G08—SIGNALLING
- G08G—TRAFFIC CONTROL SYSTEMS
- G08G5/00—Traffic control systems for aircraft, e.g. air-traffic control [ATC]
- G08G5/0004—Transmission of traffic-related information to or from an aircraft
- G08G5/0013—Transmission of traffic-related information to or from an aircraft with a ground station
-
- G—PHYSICS
- G08—SIGNALLING
- G08G—TRAFFIC CONTROL SYSTEMS
- G08G5/00—Traffic control systems for aircraft, e.g. air-traffic control [ATC]
- G08G5/0017—Arrangements for implementing traffic-related aircraft activities, e.g. arrangements for generating, displaying, acquiring or managing traffic information
- G08G5/0026—Arrangements for implementing traffic-related aircraft activities, e.g. arrangements for generating, displaying, acquiring or managing traffic information located on the ground
-
- G—PHYSICS
- G08—SIGNALLING
- G08G—TRAFFIC CONTROL SYSTEMS
- G08G5/00—Traffic control systems for aircraft, e.g. air-traffic control [ATC]
- G08G5/003—Flight plan management
- G08G5/0039—Modification of a flight plan
-
- G—PHYSICS
- G08—SIGNALLING
- G08G—TRAFFIC CONTROL SYSTEMS
- G08G5/00—Traffic control systems for aircraft, e.g. air-traffic control [ATC]
- G08G5/0047—Navigation or guidance aids for a single aircraft
- G08G5/0069—Navigation or guidance aids for a single aircraft specially adapted for an unmanned aircraft
-
- G—PHYSICS
- G08—SIGNALLING
- G08G—TRAFFIC CONTROL SYSTEMS
- G08G5/00—Traffic control systems for aircraft, e.g. air-traffic control [ATC]
- G08G5/04—Anti-collision systems
- G08G5/045—Navigation or guidance aids, e.g. determination of anti-collision manoeuvers
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Aviation & Aerospace Engineering (AREA)
- Business, Economics & Management (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Human Resources & Organizations (AREA)
- Strategic Management (AREA)
- Economics (AREA)
- Theoretical Computer Science (AREA)
- Quality & Reliability (AREA)
- Marketing (AREA)
- Entrepreneurship & Innovation (AREA)
- Operations Research (AREA)
- Game Theory and Decision Science (AREA)
- Tourism & Hospitality (AREA)
- Development Economics (AREA)
- General Business, Economics & Management (AREA)
- Automation & Control Theory (AREA)
- Life Sciences & Earth Sciences (AREA)
- Computational Linguistics (AREA)
- Mathematical Physics (AREA)
- General Engineering & Computer Science (AREA)
- Computing Systems (AREA)
- Molecular Biology (AREA)
- General Health & Medical Sciences (AREA)
- Evolutionary Computation (AREA)
- Data Mining & Analysis (AREA)
- Software Systems (AREA)
- Biophysics (AREA)
- Biomedical Technology (AREA)
- Artificial Intelligence (AREA)
- Health & Medical Sciences (AREA)
- Computer Networks & Wireless Communication (AREA)
- Traffic Control Systems (AREA)
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
- Navigation (AREA)
Abstract
一种基于鸽群算法的飞行器航迹规划方法及装置,属于无人飞行器航迹规划与多属性决策领域技术领域,本发明提供的基于鸽群算法的飞行器航迹规划方法,首先建立包含不确定性的轨迹预测模型,然后确定规定区域内的待优化路径,采用鸽群算法,通过地图和指南针操作和地标操作,迭代得到最优路径,最后将获得的最优路径的各个参数输出。对应地,本发明提供的基于鸽群算法的飞行器航迹规划装置,包括获取模块、构建模块、确定模块、优化模块和存储模块。本发明推导计算出轨迹预测模型,利用该模型获得的路径稳定性好,具有鲁棒性和可行性;并且采用鸽群智能优化方法,解决了复杂连续优化问题,计算搜索过程具有并行性、可行性、强鲁棒性的特点。
Description
技术领域
本发明属于无人飞行器航迹规划与多属性决策领域技术领域,具体是一种基于鸽群算法的飞行器航迹规划方法及装置。
背景技术
路径规划是按照一定的评价标准体系,寻找运动体从起始状态点到达目标状态点的满足特定约束的路径。航迹规划是路径规划的一种,比一般的路径规划更加困难复杂。航迹规划是在给定的规划空间内,寻找运动体从起始点到达目标点且满足某些约束条件和一定性能指标的最优或可行飞行航迹,使运动体安全地完成预定任务。由于飞行器的运动性能复杂、任务环境复杂,航迹规划***需要综合考虑飞行器的机动性能、任务时间、地形环境、敌控区域等因素。从数学意义上讲,航迹规划就是在众多约束条件中寻找最优解。
多目标优化问题(MOP)存在多个彼此冲突的目标。与单目标优化问题不同,多目标优化的本质在于,在大多数情况下,某目标的改善可能引起其他目标性能的降低,同时使多个目标达到最优是不可能的,只能在各目标间协调权衡折中处理,是所有目标函数尽可能达到最优。
航迹规划***主要包括航迹规划算法、轨迹跟踪控制、虚拟现实技术等内容,其中核心是航迹规划算法。目前应用于航迹规划的算法主要有A*算法、遗传算法、神经网络等,其中遗传算法常用来作为航迹规划算法。但是由于环境空间巨大、约束条件繁多且耦合性强等情况,之前的算法结果并没有很好地解决问题。
鸽群智能优化算法(Pigeon-Inspired Optimization,PIO)是段海滨教授于2014年提出的一种基于鸽群启发式的仿生智能优化算法。古代的飞鸽传书,鸽子主要借助地磁场、太阳和地标确定飞行方向,进而抵达目的地。鸽群算法主要包含两个数学模型:一个是基于地磁场和太阳的地图指南针模型,一个是基于地标的地标模型。鸽群智能优化较于其它仿生智能算法在搜索过程中具有并行性、可行性、强鲁棒性等特点,因此它可以用来解决复杂的连续优化问题。
目前,现有飞行器航迹规划的技术中,由于未考虑不确定性条件,例如风、航向角改变、开始和结束的操作位置点等,所规划得到的路径稳定性差,当稍有偏离时,所规划的路径不能适应,就需要再重新规划,耗费时间。
发明内容
本发明的目的是提供一种基于鸽群算法的飞行器航迹规划方法及装置,着重考虑不确定性条件下的航迹优化问题,相对于普通的航迹规划问题,本发明考虑不确定性对路径优化的显著影响,获得的路径稳定性好。
本发明提供的基于鸽群算法的飞行器航迹规划方法,首先建立包含不确定性的轨迹预测模型,然后确定规定区域内的待优化路径,采用鸽群算法对待优化路径进行优化获取最优路径。
所述的建立包含不确定性的轨迹预测模型,具体是:
设规定区域内起始点和目的地之间共有K个改变航向角的位置,所改变的航向角依次表示为θ1,θ2,...,θK,整个飞行器的航行轨迹共有K+1段路径组成,路径长度依次为d0,d1,...,dK,则建立飞行器的航行轨迹路径函数fL为:
K为正整数,k=0,1,…,K;
设规定区域内共有威胁中心m个,飞行器的航行轨迹点表示为p0,p1,...,pn,pn+1,p0,pn+1分别表示起始点与目的地,p1,...,pn为规划出的n个航迹点,n为正整数,为每个轨迹点考虑不确定性因素产生的椭圆凸包,设rij表示航迹点pi的椭圆凸包与第j个威胁中心的最短距离;则飞行器的航行轨迹威胁代价fTA为:
其中,j=1,2,...,m,i=1,2,...,n;m,n均为正整数;rsafe表示距离威胁中心的安全距离;
则包含不确定性的轨迹预测模型表示如下:
minfcost=wfL+(1-w)fTA;
约束条件包括:各航向角θ1,θ2,...,θK-1改变的值在设定范围内;各路径长度d0,d1,...,dK-1的最小值为最小步长L,最大值为设定的上限值;d0,d1,...,dK-1,θ1,θ2,...,θK-1不同时为0;rij≥rsafe;
其中,w表示权重系数。
根据所建立的轨迹预测模型,采用鸽群算法获取最优路径,输出d0,d1,...,dK-1,θ1,θ2,...,θK-1。
相应地,本发明提供的飞行器航迹规划装置,包括:
获取模块,用于获取规定区域内的路径信息;
构建模块,用于建立包含不确定性的轨迹预测模型;
确定模块,用于根据所述路径信息及轨迹预测模型来确定待优化路径;
优化模块,用于采用鸽群算法对所述待优化路径进行优化;
存储模块,用于存储最优路径的各个参数。
其中的构建模块所建立的轨迹预测模型为上面所述的minfcost=wfL+(1-w)fTA。
本发明提供的基于鸽群算法的飞行器航迹规划方法及装置,其优点和积极效果在于:考虑了不确定性条件下的航迹优化问题,本发明推导计算出轨迹预测模型,相对于现有方法,利用该模型获得的路径稳定性好,具有一定的鲁棒性和可行性。并且本发明采用鸽群智能优化方法,解决了所涉及的复杂连续优化问题,计算搜索过程中具有并行性、可行性、强鲁棒性等特点。
附图说明
图1是本发明实施例中改变航线角避开威胁的航迹操作示意图;
图2是本发明考虑不确定性因素后的航迹示意图;
图3是本发明的飞行器航迹规划实现步骤的流程框架图;
图4是本发明采用的鸽群算法中地图和指南针模型示意图;
图5是本发明采用的鸽群算法中地标模型示意图;
图6是本发明采用的鸽群算法模型示意图;
图7是本发明的路径优化装置示意图。
具体实施方式
下面将结合附图和实施例对本发明作进一步的详细说明。
本发明在飞行器的轨迹模型中,预测下一飞行位置时,考虑了不同的不确定性来源,如风、航向角改变、开始和结束的操作位置点等,然后运用鸽群算法进行路径优化,最终选择最优路径。
本发明实施例在轨迹模型中,为了限制搜索空间在一个“合理”的范围内,考虑从起始点到目的地存在三步操作,即三次改变航向角。如图1所示,在飞行距离d0进行第一步操作改变航向角α,对于距离d0存在不确定性参数距离误差ε0,这意味着飞机可能在(d0+ε0)或(d0-ε0)启动操作,对于航向角α存在不确定性参数εα,εα是指第一步操作航向角改变α时存在的误差;再飞行距离d1后进行第二步操作改变航向角β,同样,对于距离d1存在误差ε1,对于航向角β存在角度误差εβ;然后在飞行距离d2进行第三部操作改变航向角飞向目的地,对于距离d2存在距离误差ε2。
图1中,设O为起始点,D为目的地点,A、B、C分别是改变航向角的位置,起始点的坐标为(x0,y0),目的地点的坐标为(x4,y4),A点的坐标为(x1,y1),B点的坐标为(x2,y2),C点的坐标为(x3,y3)。A点航向角改变α,B点航向角改变β。
变量d0、d1、d2、α、β,相应的不确定性参数由ε0、ε1、ε2、εα、εβ表示,相应的上限值分别为d0max、d1max、d2max、αmax、βmax。dmin表示起始点到目的地的最短路径距离,航迹点用p0,p1,...,pn,pn+1表示,p0,pn+1分别表示起始点与目的地点,p1,...,pn为规划出的n个航迹点,n为正整数。考虑不确定性因素产生的椭圆凸包,如图2所示,设fTA表示椭圆凸包与威胁区域间所产生的威胁代价,设威胁中心有m个,表示为j=1,2,...,m,设rij表示航迹点pi的椭圆凸包与第j个威胁中心的最短距离,rsafe表示距离威胁中心的安全距离。
根据起始点O和目的地点D的二维坐标,可得:
起始点O和目的地点D直接连线的斜率
得到A点坐标(x1,y1)为:
得到A点与B点连线斜率k2=tan(α+γ);
得到B点坐标(x2,y2)为:
得到B点与C点的连线斜率k3=tan(α+β+γ);
得到C点坐标(x3,y3)为:
最终推导出C点到目的地的距离
则进一步,可获得飞行器的航行轨迹路径函数fL和威胁代价fTA;
本发明设置目标函数如下:
minfcost=wfL+(1-w)fTA;
约束条件如下:
L≤d0≤d0max,L≤d1≤d1max,L≤d2≤d2max;
d0、d1、d2、α、β不可同时为0;
rij≥rsafe。
其中,w表示权重系数,取值范围从0到1;L表示最小步长,最小步长是改变飞行航线进行拐弯前后都必须保持直飞的最短距离。
根据所建立的轨迹预测模型,来求取最优路径,即获得d0、d1、d2、α、β。
本发明的基于鸽群算法的飞行器航迹规划方法,整体可包括如图3所示的步骤,如下:
步骤1,建立包含不确定性的轨迹预测模型;
步骤2,根据规定区域内的路径信息,初始化鸽群算法待优化路径,以及对鸽群算法中搜索空间维度D、鸽群规模Np、迭代次数Ncmax、地磁因数R等参数进行初始化设置;每只鸽子代表一条待优化路径。
步骤3,随机设置每只鸽子的速度与位置,根据目标函数设置适应度函数,计算适应值,找到当前最优路径,并存储当前最优路径的各个参数d0、d1、d2、α、β。当前最优路径对应的适应值最大。
根据上面的轨迹预测模型,本发明是求解最小化问题,所使用的目标函数表示为
fmin(X)为目标函数minfcost=wfL+(1-w)fTA,X为某路径。分母不能为0,因为最小代价可能为0,所以ε表示一个较小常数。
步骤4,地图和指南针算子操作,更新每只鸽子的速度与位置;
步骤5:地标操作,根据适应值大小给所有鸽子排序,适应值较低的鸽子跟随适应值高的鸽子飞行,找到鸽群的中心位置(目的地),所有鸽子将直接飞向目的地。
计算各路径的适应值,更新当前最优路径的各个参数d0、d1、d2、α、β。
步骤6:判断是否达到最大迭达次数,若否,继续转步骤4执行,反复进行地图和指南针操作和地标操作,直到迭代次数大于地标算子的最大迭代次数时停止操作。
本发明在求解最优路径时,所采用的鸽群算法,参考文献1:Three-DimensionalPath Planning for Uninhabited Combat Aerial Vehicle Based on Predator-PreyPigeon-Inspired Optimization in Dynamic Environment;Bo Zhang,Haibin Duan;《IEEE/ACM Transactions on Computational Biology&Bioinformatics》,2017,PP(99):1-1。
步骤4中进行地图和指南针算子操作。鸽子可以使用磁性物体感知地磁场,然后在头脑中形成地图。它们把太阳高度作为指南针来调整飞行方向,当它们接近目的地的时候,它们对太阳和磁性物体的依赖性便减小。
在D维搜索空间中,设第i只鸽子的位置为Xi,速度为Vi,表示如下:
Xi=(Xi1,Xi2,...,XiD),Vi=(Vi1,Vi2,...,ViD);
第i只鸽子的位置和速度更新公式如下:
其中,R是地磁因数,取值范围为0到1;r1是0到1的随机数。Vi t、分别表示迭代至第t代时的第i只鸽子的速度与位置,Xg表示在t-1次迭代循环后,通过比较所有鸽子的位置而得到的全局最优位置。
如图4所示,地图指南针算子中粗箭头(速度矢量)指向适应值最优的鸽子,对应公式中的细箭头是原本鸽子飞行方向,对应公式中的Vi t-1e-Rt,两者矢量和表示鸽子下一时刻飞行方向。
步骤5中进行地标算子操作。地标算子是模仿地标对鸽子的影响。当鸽群接近目的地时鸽子会依靠熟悉的地标直接寻找目的地。而不熟悉地标的鸽子会跟随熟悉地标的鸽子飞行最终到达目的地。地标算子在工作过程中会在每次迭代循环中将鸽子的总数Np减半,然后找到所有鸽子的加权中心位置,这个位置就是目的地,如图5所示。
当鸽群接近目的地时,由于此时算法的解趋于收敛,且收敛速度较快,所以地标算子所起作用较小。因此本发明在每次迭代过程中,地图和指南针算子操作与地标算子操作同步进行。本发明利用地标算子对鸽子的位置进行更新,如下:
更新
其中,表示第t代鸽子的数量,Npmax表示最大鸽子总数,常数c是地标操作中鸽群数目因子,表示第t代鸽群的中心位置(目的地)。适应度函数f定义为鸽子个体的品质因数,即路径优化代价。表示为:本发明是求解最小化问题。参数q是地标操作的影响因子,表示为:s是0到1的常数,Ncmax表示最大迭代次数。r2是0到1的随机数。
接近中心位置的鸽群会依靠熟悉的地标直接寻找目的地,而离中心位置较远即不熟悉地标的鸽子会跟随熟悉地标的鸽子飞行最终到达目的地。
本发明所采用的鸽群算法模型如图6所示,通过地图和指南针操作和地标操作,迭代得到最优路径,最后将获得的最优路径的各个参数输出。
对应地,本发明提供的基于鸽群算法的飞行器航迹规划装置,如图7所示,包括获取模块、构建模块、确定模块、优化模块和存储模块。下面说明各模块。
获取模块用于获取规定区域内的路径信息,主要包括规定区域中起始点和目的地,障碍物信息等等。
构建模块用于建立包含不确定性的轨迹预测模型。建立模型的过程在此不在赘述。
确定模块用于根据所述路径信息及轨迹预测模型来确定待优化路径。
优化模块用于采用鸽群算法对所述待优化路径进行优化。采用如图4~6所示的鸽群算法对待优化路径进行优化,获取最优路径。
存储模块用于存储最优路径的参数,参数包括起始点和目的地之间改变航向角的位置以及改变的角度。
本发明考虑了不确定性影响,相对于现有方法,所获得的路径稳定性好,具有一定的鲁棒性和可行性。
Claims (4)
1.一种基于鸽群算法的飞行器航迹规划方法,在预测飞行位置时考虑不确定性来源,其特征在于,所述方法首先建立包含不确定性的轨迹预测模型,然后确定规定区域内的待优化路径,采用鸽群算法对待优化路径进行优化获取最优路径;
所述的建立包含不确定性的轨迹预测模型,具体是:
设规定区域内起始点和目的地之间共有K个改变航向角的位置,所改变的航向角依次表示为θ1,θ2,...,θK,整个飞行器的航行轨迹共有K+1段路径组成,路径长度依次为d0,d1,...,dK,则建立飞行器的航行轨迹路径函数fL为:
K为正整数,k=0,1,…,K;
设规定区域内共有威胁中心m个,飞行器的航行轨迹点表示为p0,p1,...,pn,pn+1,p0,pn+1分别表示起始点与目的地,p1,...,pn为规划出的n个航迹点,n为正整数,为每个轨迹点考虑不确定性因素产生的椭圆凸包,设rij表示航迹点pi的椭圆凸包与第j个威胁中心的最短距离;则飞行器的航行轨迹威胁代价fTA为:
<mrow>
<msub>
<mi>f</mi>
<mrow>
<mi>T</mi>
<mi>A</mi>
</mrow>
</msub>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>m</mi>
</munderover>
<mfrac>
<mn>1</mn>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>r</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>/</mo>
<msub>
<mi>r</mi>
<mrow>
<mi>s</mi>
<mi>a</mi>
<mi>f</mi>
<mi>e</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mfrac>
</mrow>
其中,j=1,2,...,m,i=1,2,...,n;m,n均为正整数;rsafe表示距离威胁中心的安全距离;
则包含不确定性的轨迹预测模型表示如下:
min fcost=wfL+(1-w)fTA;
约束条件包括:各航向角θ1,θ2,...,θK-1改变的值在设定范围内;各路径长度d0,d1,...,dK-1的最小值为最小步长L,最大值为设定的上限值;d0,d1,...,dK-1,θ1,θ2,...,θK-1不同时为0;rij≥rsafe;
其中,w表示权重系数;
根据所建立的轨迹预测模型,采用鸽群算法获取最优路径,输出d0,d1,...,dK-1,θ1,θ2,...,θK-1。
2.根据权利要求1所述的飞行器航迹规划方法,其特征在于,所述的建立包含不确定性的轨迹预测模型,设起始点和目的地之间有三个改变航向角的位置,依次表示为A点、B点和C点;在A点航向角改变α,B点航向角改变β;飞行器的航行轨迹包括4段路径d0,d1,d2,d3,则飞行器的航行轨迹路径函数fL=(d0+d1+d2+d3)2;
所建立的包含不确定性的轨迹预测模型为:
min fcost=wfL+(1-w)fTA;
约束条件如下:
<mrow>
<mo>-</mo>
<msub>
<mi>&alpha;</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msub>
<mo>&le;</mo>
<mi>&alpha;</mi>
<mo>&le;</mo>
<msub>
<mi>&alpha;</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msub>
<mo>,</mo>
<msub>
<mi>&alpha;</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mi>&pi;</mi>
<mn>6</mn>
</mfrac>
<mo>;</mo>
<mo>-</mo>
<msub>
<mi>&beta;</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msub>
<mo>&le;</mo>
<mi>&beta;</mi>
<mo>&le;</mo>
<msub>
<mi>&beta;</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msub>
<mo>,</mo>
<msub>
<mi>&beta;</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mi>&pi;</mi>
<mn>6</mn>
</mfrac>
<mo>;</mo>
</mrow>
L≤d0≤d0max,L≤d1≤d1max,L≤d2≤d2max;
d0、d1、d2、α、β不同时为0;
rij≥rsafe;
其中,d0max、d1max、d2max、αmax、βmax分别是变量d0、d1、d2、α、β的上限值。
3.一种基于鸽群算法的飞行器航迹规划装置,其特征在于,包括:
获取模块,用于获取规定区域内的路径信息;
构建模块,用于建立包含不确定性的轨迹预测模型;
确定模块,用于根据所述路径信息及轨迹预测模型来确定待优化路径;
优化模块,用于采用鸽群算法对所述待优化路径进行优化;
存储模块,用于存储最优路径的各个参数。
4.根据权利要求3所述的飞行器航迹规划装置,其特征在于,所述的构建模块
建立包含不确定性的轨迹预测模型,具体是:
设规定区域内起始点和目的地之间共有K个改变航向角的位置,所改变的航向角依次表示为θ1,θ2,...,θK,整个飞行器的航行轨迹共有K+1段路径组成,路径长度依次为d0,d1,...,dK,则建立飞行器的航行轨迹路径函数fL为:
K为正整数,k=0,1,…,K;
设规定区域内共有威胁中心m个,飞行器的航行轨迹点表示为p0,p1,...,pn,pn+1,p0,pn+1分别表示起始点与目的地,p1,...,pn为规划出的n个航迹点,n为正整数,为每个轨迹点考虑不确定性因素产生的椭圆凸包,设rij表示航迹点pi的椭圆凸包与第j个威胁中心的最短距离;则飞行器的航行轨迹威胁代价fTA为:
<mrow>
<msub>
<mi>f</mi>
<mrow>
<mi>T</mi>
<mi>A</mi>
</mrow>
</msub>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>m</mi>
</munderover>
<mfrac>
<mn>1</mn>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>r</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>/</mo>
<msub>
<mi>r</mi>
<mrow>
<mi>s</mi>
<mi>a</mi>
<mi>f</mi>
<mi>e</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mfrac>
</mrow>
其中,j=1,2,...,m,i=1,2,...,n;m,n均为正整数;rsafe表示距离威胁中心的安全距离;
则包含不确定性的轨迹预测模型表示如下:
min fcost=wfL+(1-w)fTA;
约束条件包括:各航向角θ1,θ2,...,θK-1改变的值在设定范围内;各路径长度d0,d1,...,dK-1的最小值为最小步长L,最大值为设定的上限值;d0,d1,...,dK-1,θ1,θ2,...,θK-1不同时为0;rij≥rsafe;
其中,w表示权重系数;
根据所建立的轨迹预测模型,采用鸽群算法获取最优路径,输出d0,d1,...,dK-1,θ1,θ2,...,θK-1。
Priority Applications (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710625878.8A CN107504972B (zh) | 2017-07-27 | 2017-07-27 | 一种基于鸽群算法的飞行器航迹规划方法及装置 |
US16/048,206 US20190035286A1 (en) | 2017-07-27 | 2018-07-27 | Airplane flight path planning method and device based on the pigeon-inspired optimization |
US17/033,785 US11727812B2 (en) | 2017-07-27 | 2020-09-26 | Airplane flight path planning method and device based on the pigeon-inspired optimization |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710625878.8A CN107504972B (zh) | 2017-07-27 | 2017-07-27 | 一种基于鸽群算法的飞行器航迹规划方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107504972A true CN107504972A (zh) | 2017-12-22 |
CN107504972B CN107504972B (zh) | 2018-08-07 |
Family
ID=60689550
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710625878.8A Active CN107504972B (zh) | 2017-07-27 | 2017-07-27 | 一种基于鸽群算法的飞行器航迹规划方法及装置 |
Country Status (2)
Country | Link |
---|---|
US (1) | US20190035286A1 (zh) |
CN (1) | CN107504972B (zh) |
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108592921A (zh) * | 2018-05-02 | 2018-09-28 | 山东理工大学 | 一种分段最速下降混合航线规划方法 |
CN109254588A (zh) * | 2018-10-17 | 2019-01-22 | 北京航空航天大学 | 一种基于交叉变异鸽群优化的无人机集群协同侦察方法 |
CN109917806A (zh) * | 2019-03-14 | 2019-06-21 | 北京航空航天大学 | 一种基于非劣解鸽群优化的无人机集群编队控制方法 |
CN110069075A (zh) * | 2019-04-19 | 2019-07-30 | 北京航空航天大学 | 一种仿鸽群应急避障机制的集群超机动避障方法 |
CN110081875A (zh) * | 2019-04-25 | 2019-08-02 | 北京航空航天大学 | 一种仿鸽子智能的无人机自主导航***及方法 |
CN110109477A (zh) * | 2019-05-10 | 2019-08-09 | 北京航空航天大学 | 基于鸽群智能反向学习的无人机集群多目标控制优化方法 |
CN110672088A (zh) * | 2019-09-09 | 2020-01-10 | 北京航空航天大学 | 一种仿信鸽地貌感知归巢机制的无人机自主导航方法 |
CN110749325A (zh) * | 2019-11-29 | 2020-02-04 | 北京京东乾石科技有限公司 | 航迹规划方法和装置 |
CN110989690A (zh) * | 2019-12-24 | 2020-04-10 | 北京航空航天大学 | 一种基于诱导信息的多无人机寻路方法 |
CN111089593A (zh) * | 2019-12-24 | 2020-05-01 | 山东科技大学 | 一种基于鸽群算法的机器人路径规划方法 |
CN111913482A (zh) * | 2020-07-14 | 2020-11-10 | 徐州工程学院 | 一种基于多Agent的鸽群机器人运动控制方法 |
CN112596574A (zh) * | 2020-12-21 | 2021-04-02 | 广东工业大学 | 基于分层鸽群算法的光伏最大功率跟踪方法及装置 |
CN113821054A (zh) * | 2021-09-29 | 2021-12-21 | 北京航空航天大学 | 基于鸽群智能优化动态逆控制的无人机轨迹跟踪制导方法 |
CN113885320A (zh) * | 2021-09-26 | 2022-01-04 | 北京航空航天大学 | 一种基于混合量子鸽群优化的飞行器随机鲁棒控制方法 |
Families Citing this family (29)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US11727812B2 (en) * | 2017-07-27 | 2023-08-15 | Beihang University | Airplane flight path planning method and device based on the pigeon-inspired optimization |
US11398158B2 (en) * | 2018-06-29 | 2022-07-26 | Satcom Direct, Inc. | System and method for forecasting availability of network services during flight |
US10860115B1 (en) | 2019-09-19 | 2020-12-08 | Bao Tran | Air transportation systems and methods |
US11675324B2 (en) | 2019-09-19 | 2023-06-13 | Bao Tran | Air transportation systems and methods |
CN111176324B (zh) * | 2019-12-31 | 2023-08-11 | 航天时代飞鸿技术有限公司 | 一种多无人机分布式协同编队规避动态障碍的方法 |
CN111207752B (zh) * | 2020-01-17 | 2022-10-04 | 西北工业大学 | 一种基于动态切点调整的无人机航迹规划方法 |
US11804052B2 (en) * | 2020-03-26 | 2023-10-31 | Seiko Epson Corporation | Method for setting target flight path of aircraft, target flight path setting system, and program for setting target flight path |
CN111562786B (zh) * | 2020-05-19 | 2023-03-28 | 智慧航海(青岛)科技有限公司 | 多级优化的无人船舶路径规划方法及装置 |
CN111580556B (zh) * | 2020-05-22 | 2023-04-07 | 中国人民解放军国防科技大学 | 一种时空约束下多无人机协同路径规划与制导方法 |
CN113095538B (zh) * | 2020-06-08 | 2024-03-19 | 华北电力大学 | 面向灵活性运行的热电联产机组宽负荷运行动态特性建模方法 |
CN112150059B (zh) * | 2020-06-23 | 2022-09-02 | 国网天津市电力公司电力科学研究院 | 一种基于乌鸦算法的计量器具智能仓库调度优化方法 |
CN112015199B (zh) * | 2020-07-17 | 2024-05-28 | 煤炭科学技术研究院有限公司 | 应用于煤矿井下智能巡检无人机的航迹规划方法及装置 |
CN114372603A (zh) * | 2020-11-13 | 2022-04-19 | 北京航空航天大学 | 一种仿鸽群多学习智能的无人靶机协同航路动态规划方法 |
CN112666981B (zh) * | 2021-01-04 | 2022-11-18 | 北京航空航天大学 | 基于原鸽群动态群组学习的无人机集群动态航路规划方法 |
CN113155880B (zh) * | 2021-05-08 | 2023-06-27 | 电子科技大学 | 一种采用无人机和xrf技术对土壤重金属污染的检测方法 |
CN113741508B (zh) * | 2021-06-29 | 2023-03-17 | 南京航空航天大学 | 基于改进狼群算法的无人机任务分配方法 |
CN113188562B (zh) * | 2021-07-01 | 2022-03-01 | 新石器慧通(北京)科技有限公司 | 可行驶区域的路径规划方法、装置、电子设备及存储介质 |
CN113885566B (zh) * | 2021-10-21 | 2024-01-23 | 重庆邮电大学 | 面向多无人机数据采集时间最小化的v型轨迹规划方法 |
CN113867369B (zh) * | 2021-12-03 | 2022-03-22 | 中国人民解放军陆军装甲兵学院 | 一种基于交流学习的海鸥算法的机器人路径规划方法 |
CN113867368B (zh) * | 2021-12-03 | 2022-02-25 | 中国人民解放军陆军装甲兵学院 | 一种基于改进海鸥算法的机器人路径规划方法 |
CN114545976A (zh) * | 2022-03-08 | 2022-05-27 | 中山大学 | 一种飞行器的避障飞行控制方法、装置及*** |
CN114440896B (zh) * | 2022-04-08 | 2022-06-21 | 中国人民解放军96901部队 | 一种基于威胁场景动态识别的高速飞行器飞行管道规划方法 |
CN115240475B (zh) * | 2022-09-23 | 2022-12-13 | 四川大学 | 融合飞行数据和雷达图像的航空器进近规划方法及装置 |
CN116614827B (zh) * | 2023-06-02 | 2024-05-17 | 昆明理工大学 | 一种面向无人机辅助通信的航迹优化方法 |
CN116518982B (zh) * | 2023-07-03 | 2023-09-19 | 江西啄木蜂科技有限公司 | 一种低空林业监测遥感无人机路径多目标规划方法 |
CN117371640B (zh) * | 2023-12-08 | 2024-04-12 | 山东省地质测绘院 | 基于无人机遥感的测绘路线优化方法及*** |
CN117373263B (zh) * | 2023-12-08 | 2024-03-08 | 深圳市永达电子信息股份有限公司 | 一种基于量子鸽群算法的交通流量预测方法及装置 |
CN117591794B (zh) * | 2024-01-18 | 2024-03-22 | 吉林省林业科学研究院(吉林省林业生物防治中心站) | 基于时间序列的鸟类迁徙轨迹预测方法 |
CN118102414B (zh) * | 2024-04-23 | 2024-07-09 | 南京理工大学 | 一种基于谱聚类-离散鸽群算法的无人机集群分簇方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102509476A (zh) * | 2011-11-01 | 2012-06-20 | 四川川大智胜软件股份有限公司 | 一种短期飞行冲突解脱方法 |
CN102854883A (zh) * | 2012-07-04 | 2013-01-02 | 北京航空航天大学 | 一种无人机动态紧急避撞区的建模方法 |
CN103971180A (zh) * | 2014-05-09 | 2014-08-06 | 北京航空航天大学 | 一种基于鸽群启发式优化的连续优化问题求解方法 |
CN105043400A (zh) * | 2015-06-30 | 2015-11-11 | 百度在线网络技术(北京)有限公司 | 路径规划方法及装置 |
CN105930307A (zh) * | 2016-04-22 | 2016-09-07 | 大连理工大学 | 一种新型群智能优化算法-鸽群算法 |
CN106441308A (zh) * | 2016-11-10 | 2017-02-22 | 沈阳航空航天大学 | 一种基于自适应权重鸽群算法的无人机航路规划方法 |
-
2017
- 2017-07-27 CN CN201710625878.8A patent/CN107504972B/zh active Active
-
2018
- 2018-07-27 US US16/048,206 patent/US20190035286A1/en not_active Abandoned
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102509476A (zh) * | 2011-11-01 | 2012-06-20 | 四川川大智胜软件股份有限公司 | 一种短期飞行冲突解脱方法 |
CN102854883A (zh) * | 2012-07-04 | 2013-01-02 | 北京航空航天大学 | 一种无人机动态紧急避撞区的建模方法 |
CN103971180A (zh) * | 2014-05-09 | 2014-08-06 | 北京航空航天大学 | 一种基于鸽群启发式优化的连续优化问题求解方法 |
CN105043400A (zh) * | 2015-06-30 | 2015-11-11 | 百度在线网络技术(北京)有限公司 | 路径规划方法及装置 |
CN105930307A (zh) * | 2016-04-22 | 2016-09-07 | 大连理工大学 | 一种新型群智能优化算法-鸽群算法 |
CN106441308A (zh) * | 2016-11-10 | 2017-02-22 | 沈阳航空航天大学 | 一种基于自适应权重鸽群算法的无人机航路规划方法 |
Non-Patent Citations (3)
Title |
---|
RAN HAO等: "Multiple UAVs Mission Assignment Based on Modified Pigeon-Inspired Optimization Algorithm", 《IEEE CHINESE GUIDANCE,NAVIGATION AND CONTROL CONFERENCE》 * |
ZHANG BO等: "Three-Dimensional Path Planning for Uninhabited Combat Aerial Vehicle Based on Predator-Prey Pigeon-Inspired Optimization in Dynamic Environment", 《IEEE/ACM TRANSACTIONS ON COMPUTATIONAL BIOLOGY AND BIOINFORMATICS》 * |
段海滨等: "基于捕食逃逸鸽群优化的无人机紧密编队协同控制", 《中国科学:技术科学》 * |
Cited By (23)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108592921A (zh) * | 2018-05-02 | 2018-09-28 | 山东理工大学 | 一种分段最速下降混合航线规划方法 |
CN108592921B (zh) * | 2018-05-02 | 2021-07-27 | 山东理工大学 | 一种分段最速下降混合航线规划方法 |
CN109254588A (zh) * | 2018-10-17 | 2019-01-22 | 北京航空航天大学 | 一种基于交叉变异鸽群优化的无人机集群协同侦察方法 |
CN109917806A (zh) * | 2019-03-14 | 2019-06-21 | 北京航空航天大学 | 一种基于非劣解鸽群优化的无人机集群编队控制方法 |
CN109917806B (zh) * | 2019-03-14 | 2020-08-11 | 北京航空航天大学 | 一种基于非劣解鸽群优化的无人机集群编队控制方法 |
CN110069075B (zh) * | 2019-04-19 | 2020-06-09 | 北京航空航天大学 | 一种仿鸽群应急避障机制的集群超机动避障方法 |
CN110069075A (zh) * | 2019-04-19 | 2019-07-30 | 北京航空航天大学 | 一种仿鸽群应急避障机制的集群超机动避障方法 |
CN110081875A (zh) * | 2019-04-25 | 2019-08-02 | 北京航空航天大学 | 一种仿鸽子智能的无人机自主导航***及方法 |
CN110081875B (zh) * | 2019-04-25 | 2020-10-09 | 北京航空航天大学 | 一种仿鸽子智能的无人机自主导航***及方法 |
CN110109477A (zh) * | 2019-05-10 | 2019-08-09 | 北京航空航天大学 | 基于鸽群智能反向学习的无人机集群多目标控制优化方法 |
CN110672088A (zh) * | 2019-09-09 | 2020-01-10 | 北京航空航天大学 | 一种仿信鸽地貌感知归巢机制的无人机自主导航方法 |
CN110749325A (zh) * | 2019-11-29 | 2020-02-04 | 北京京东乾石科技有限公司 | 航迹规划方法和装置 |
CN110749325B (zh) * | 2019-11-29 | 2021-11-02 | 北京京东乾石科技有限公司 | 航迹规划方法和装置 |
CN110989690A (zh) * | 2019-12-24 | 2020-04-10 | 北京航空航天大学 | 一种基于诱导信息的多无人机寻路方法 |
CN111089593A (zh) * | 2019-12-24 | 2020-05-01 | 山东科技大学 | 一种基于鸽群算法的机器人路径规划方法 |
CN111913482B (zh) * | 2020-07-14 | 2022-06-21 | 徐州工程学院 | 一种基于多Agent的鸽群机器人运动控制方法 |
CN111913482A (zh) * | 2020-07-14 | 2020-11-10 | 徐州工程学院 | 一种基于多Agent的鸽群机器人运动控制方法 |
CN112596574A (zh) * | 2020-12-21 | 2021-04-02 | 广东工业大学 | 基于分层鸽群算法的光伏最大功率跟踪方法及装置 |
CN112596574B (zh) * | 2020-12-21 | 2022-06-24 | 广东工业大学 | 基于分层鸽群算法的光伏最大功率跟踪方法及装置 |
CN113885320A (zh) * | 2021-09-26 | 2022-01-04 | 北京航空航天大学 | 一种基于混合量子鸽群优化的飞行器随机鲁棒控制方法 |
CN113885320B (zh) * | 2021-09-26 | 2024-02-23 | 北京航空航天大学 | 一种基于混合量子鸽群优化的飞行器随机鲁棒控制方法 |
CN113821054A (zh) * | 2021-09-29 | 2021-12-21 | 北京航空航天大学 | 基于鸽群智能优化动态逆控制的无人机轨迹跟踪制导方法 |
CN113821054B (zh) * | 2021-09-29 | 2024-04-02 | 北京航空航天大学 | 基于鸽群智能优化动态逆控制的无人机轨迹跟踪制导方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107504972B (zh) | 2018-08-07 |
US20190035286A1 (en) | 2019-01-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107504972A (zh) | 一种基于鸽群算法的飞行器航迹规划方法及装置 | |
CN106125764B (zh) | 基于a*搜索的无人机路径动态规划方法 | |
Li et al. | Autonomous underwater vehicle optimal path planning method for seabed terrain matching navigation | |
CN103557867B (zh) | 一种基于稀疏a*搜索的三维多uav协同航迹规划方法 | |
CN103528585B (zh) | 一种不等距分割可通行区域的路径规划方法 | |
CN104155998B (zh) | 一种基于势场法的航迹规划方法 | |
CN105700549A (zh) | 一种基于序列小生境粒子群算法的无人机多航迹规划方法 | |
CN104121903B (zh) | 一种基于边界值问题的滚动航路规划方法 | |
CN103697895A (zh) | 基于自适应a星算法的飞行器最优路径确定方法 | |
CN106444835A (zh) | 基于Lazy Theta星和粒子群混合算法的水下潜器三维路径规划方法 | |
CN108280866A (zh) | 道路点云数据处理方法及*** | |
CN107966153A (zh) | 水下航行器路径规划算法 | |
CN103631264A (zh) | 一种同时定位与地图创建方法及装置 | |
Xu et al. | Heuristic and random search algorithm in optimization of route planning for Robot’s geomagnetic navigation | |
US20210325891A1 (en) | Graph construction and execution ml techniques | |
CN108489491A (zh) | 一种水下自主航行器的三维航迹智能规划方法 | |
Sangeetha et al. | An intelligent gain-based ant colony optimisation method for path planning of unmanned ground vehicles | |
Liu et al. | Application of dijkstra algorithm in path planning for geomagnetic navigation | |
Chu et al. | Joint vehicular localization and reflective mapping based on team channel-SLAM | |
CN114088098A (zh) | 一种用于极区水下航行器数据库辅助导航路径规划方法 | |
CN117850471A (zh) | 三维环境下考虑雷达威胁的多智能体协同航迹规划方法及其规划*** | |
Cai et al. | Multiple paths planning for UAVs using particle swarm optimization with sequential niche technique | |
Song et al. | Search and tracking strategy of autonomous surface underwater vehicle in oceanic eddies based on deep reinforcement learning | |
CN101833702B (zh) | 一种基于行人可视范围的导航点动态更换方法 | |
Li et al. | An enhanced direction calibration based on reinforcement learning for indoor localization system |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |