CN111896027B - 一种考虑地形起伏的测距敏感器仿真建模方法 - Google Patents
一种考虑地形起伏的测距敏感器仿真建模方法 Download PDFInfo
- Publication number
- CN111896027B CN111896027B CN202010682507.5A CN202010682507A CN111896027B CN 111896027 B CN111896027 B CN 111896027B CN 202010682507 A CN202010682507 A CN 202010682507A CN 111896027 B CN111896027 B CN 111896027B
- Authority
- CN
- China
- Prior art keywords
- terrain
- coordinate system
- ranging
- distance measuring
- measuring sensor
- 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.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C25/00—Manufacturing, calibrating, cleaning, or repairing instruments or devices referred to in the other groups of this subclass
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Manufacturing & Machinery (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明涉及一种考虑地形起伏的测距敏感器仿真建模方法,通过下述方式实现:S1、以当前星下点地形高程所对应的待测星体表面作为基准平面,确定测距波束所在地形范围内地形最大高程;S2、根据确定的地形最大高程确定波束搜索起点,按照预设的步长进行粗搜索,搜索波束与地形交点,当波束某个位置对应高度小于地形高度时,采用二分法精细搜索波束与地形的交点,直到测距精度满足要求。
Description
技术领域
本发明涉及一种考虑地形起伏的测距敏感器仿真建模方法,属于航天器制导、导航与控制技术领域。本发明可以应用于采用测距敏感器导航的各类地外天体软着陆探测任务的数学仿真验证,具有广泛的应用价值和市场前景。
背景技术
测距敏感器作为保证火星进入舱安全着陆的关键敏感器之一,需要根据数字地形精确模拟其在着陆过程中的相对于地面的距离测量,以真实有效地验证基于测距敏感器的导航修正方案和算法。在嫦娥系列月球软着陆动力下降过程中测距敏感器模拟测量均未考虑地形起伏的影响。对于崎岖地形,不考虑地形起伏对测距敏感器的影响显然与实际情况不符,难以真实验证***导航方案和算法,因此,需要在继承嫦娥系列月球软着陆测距敏感器数学仿真模型的基础上,设计考虑地形起伏影响的高精度测距敏感器测量建模方法。
发明内容
本发明解决的技术问题是:针对火星进入、下降、着陆过程精确模拟测距敏感器测量的需求,提出了一种考虑地形起伏的测距敏感器仿真建模方法。
本发明解决技术的方案是:一种考虑地形起伏的测距敏感器仿真建模方法,通过下述方式实现:
S1、以当前星下点地形高程所对应的待测星体表面作为基准平面,确定测距波束所在地形范围内地形最大高程;
S2、根据确定的地形最大高程确定波束搜索起点,按照预设的步长进行粗搜索,搜索波束与地形交点,当波束某个位置对应高程小于地形高程时,采用二分法精细搜索波束与地形的交点,直到测距精度满足要求。
优选的,S1具体通过下述方式实现:
利用测距敏感器安装以及着陆器的位置和姿态信息确定测距波束与基准平面交点的经纬度;
以当前星下点经纬度和测距波束与基准面交点经纬度为矩形对角点确定测距波束所在地形范围;
搜索该范围内数字地形最大高程确定测距波束所在地形范围内地形的最大高程。
优选的,通过下述方式确定测距敏感器星下点经纬度lon0、lat0和测距敏
感器到星下点的高度h:
rSF=[rSFx rSFy rSFz]T=rF+CFBrS
lon0=atan2(rSFy,rSFx)
h=||rSF||-hTerrain(lon0,lat0)
其中,||·||表示向量求模运算,rF为着陆器在待测星体固连坐标系下位置矢量,CFB为从着陆器本体坐标系到待测星体固连坐标系的转换阵,rS为测距敏感器在着陆器本体坐标系下安装位置,rSF为测距敏感器在待测星体固连坐标系下的位置,[rSFx rSFy rSFz]T为rSF的三个分量,hTerrain(lon0,lat0)表示在经纬度lon0、lat0处地形高程。
优选的,通过下述方式确定第i个测距波束与基准面交点经纬度loni、lati
loni=arctan2(rSiFy,rSiFx)
其中,rSiF为测距敏感器第i个波束与基准平面交点在待测星体固连坐标系下的位置,[rSiFx rSiFy rSiFz]T为rSiF的三个分量,rSF为测距敏感器在待测星体固连坐标系下的位置,CFB为从着陆器本体坐标系到待测星体固连坐标系的转换阵,dSi为测距敏感器第i个波束在着陆器本体坐标系下的指向,h为测距敏感器到星下点的高度。
优选的,通过下述方式确定测距波束所在地形范围内地形的最大高程:
hTerrainMaxi=max(hTerrain(lon,lat))lon0≤lon≤loni且lat0≤lat≤lati
其中,lon0、lat0为测距敏感器星下点经纬度,loni、lati为测距敏感器第i个波束与基准平面交点处经纬度,max(·)表示求最大值运算。
优选的,通过下述方式确定波束搜索起点rSiF0:
其中,||·||表示向量求模运算,hTerrainMaxi为测距波束所在地形范围内的地形最大高程、CFB为从着陆器本体坐标系到待测星体固连坐标系的转换阵,dSi为测距敏感器第i个波束在本体系下的指向,rSF为测距敏感器在待测星体固连坐标系下的位置。
优选的,所述的粗搜索过程如下:
按照下式计算第j次搜索时波束i方向上位置rSiFj及其对应经纬度lonij、latij
rSiFj=rSiF0+(j·Δl)CFBdSi=[rSiFjx rSiFjy rSiFjz]T
lonij=arctan2(rSiFjy,rSiFjx)
若||rSiFj||-hTerrain(lonij,latij)≤0则停止搜索,粗搜索得到波束i与地形交点位置即为rSiFj,否则将j+1值赋给j继续进行下一次搜索,直至得到满足条件的rSiFj;
上述公式中||·||表示向量求模运算,rSiF0为波束搜索起点,CFB为从着陆器本体坐标系到待测星体固连坐标系的转换阵,dSi为测距敏感器第i个波束在本体系下的指向,Δl为粗搜索步长,[rSiFjx rSiFjy rSiFjz]T为rSiFj的三个分量,hTerrain(lonij,latij)表示在经纬度lonij、latij处地形高程。
优选的,粗搜索步长Δl按下述方式确定:
其中,||·||表示向量求模运算,hTerrainMaxi为测距波束所在地形范围内的地形最大高程、hTerrain(lon0,lat0)表示在经纬度lon0、lat0处地形高程,CFB为从着陆器本体坐标系到待测星体固连坐标系的转换阵,dSi为测距敏感器第i个波束在本体系下的指向,rSF为测距敏感器在待测星体固连坐标系下的位置,max(·)表示取最大值运算。
优选的,所述的精细搜索过程如下:
按照下式进行第k次精细搜索得到波束i方向上位置rSiFjk及其对应经纬度lonijk、latijk;
当k=0时
rSiFj0=rSiFj
lonij0=lonij
latij0=latij
当k>0时
lonijk=arctan2(rSiFjky,rSiFjkx)
其中
上述公式中||·||表示向量求模运算,rSiFjk-1为第k-1次搜索得到的波束i方向上位置、CFB为从着陆器本体坐标系到待测星体固连坐标系的转换阵,dSi为测距敏感器第i个波束在本体系下的指向,Δl为粗搜索步长,k为精搜索次数,hTerrain(lonijk,latijk)表示在经纬度lonijk、latijk处地形高程。
优选的,通过下述公式计算测距波束i的距离测量值li:
li=||rSiFjk-rSF||
式中,rSiFjk为精细搜索下波束i与地形交点在待测星体固连坐标系下的位置;rSF为测距敏感器在待测星体固连坐标系下的位置。
本发明与现有技术相比的有益效果是:
针对火星进入、下降、着陆过程精确模拟测距敏感器测量的需求,本发明设计了测距敏感器波束作用范围内地形最大高程确定策略,极大缩小了波束与地形交点的搜索范围,为快速搜索测距波束与地形交点确定测距测量奠定了基础;设计了测距敏感器波束与地形交点快速精确搜索策略,通过按一定步长进行粗搜索和采用二分法进行精确搜索,实现了按要求精度真实模拟沿测距波束方向上的距离测量,从而更加真实有效的通过数学仿真验证基于测距敏感器的导航方案和算法等关键技术。
附图说明
图1为本发明流程图。
具体实施方式
下面结合实施例及附图1对本发明作进一步阐述。本发明一种考虑地形起伏的测距敏感器仿真建模方法,步骤如下:
第一步,测距敏感器波束作用范围内地形最大高程确定
以当前星下点地形高程所对应的地形表面作为基准平面,利用测距敏感器安装以及着陆器的位置和姿态信息确定测距波束与基准面交点的经纬度;以当前星下点经纬度和测距波束与基准面交点经纬度为矩形对角点确定测距波束所在地形范围,搜索该范围内数字地形最大高程确定测距波束所在地形范围内地形的最大高程。
定义rF为着陆器在待测星体固连坐标系下位置矢量,CFB为从着陆器本体坐标系到待测星体固连坐标系的转换阵,rS为测距敏感器在着陆器本体坐标系下安装位置,dSi为测距敏感器第i个波束在本体系下的指向。
(1)根据着陆器位置和测距敏感器安装计算测距敏感器星下点对应的经纬度lon0、lat0和测距敏感器到星下点的高度h
其中||·||表示向量求模运算,[rSFx rSFy rSFz]T为rSF的三个分量,hTerrain(lon0,lat0)表示在经纬度lon0、lat0处地形高程。
(2)根据测距敏感器第i个波束在本体系下的指向确定第i个波束与基准平面交点处经纬度loni、lati;
其中,rSiF为测距敏感器第i个波束与基准平面交点在待测星体固连坐标系下的位置,[rSiFx rSiFy rSiFz]T为rSiF的三个分量。
(3)搜索确定经度lon0~loni和纬度lat0~lati范围内地形的最大高程hTerrainMaxi
hTerrainMaxi=max(hTerrain(lon,lat))lon0≤lon≤loni且lat0≤lat≤lati (4)
其中max(·)表示求最大值运算,hTerrain(lon,lat)表示在经纬度lon、lat处地形高程。
第二步,测距波束与地形交点确定
根据第一步确定的地形最大高程hTerrainMaxi计算搜索波束i与地形交点的起始点位置rSiF0,并由该点开始沿波束方向按一定步长Δl粗略搜索波束与地形交点。当波束某个位置对应高度小于地形高程时采用二分法继续精细搜索波束与地形交点,直到满足测距精度要求为止。
(1)确定测距波束搜索起始点rSiF0
(2)按步长Δl粗略搜索波束与地形交点
当进行第j次搜索时得到波束i方向上位置rSiFj(j=0,1,2,…)及其对应经纬度lonij、latij
其中,j为粗搜索次数,[rSiFjx rSiFjy rSiFjz]T为rSiFj的三个分量,Δl为粗搜索步长。
Δl可按下式范围选取:
其中,max(·)表示取最大值运算。
Δl取值并非一定按上式范围确定,可根据地形情况调整,当地形起伏较小时,为减小计算量粗搜索步长Δl可取较大值,当地形起伏变化较为剧烈时则粗搜索步长Δl取较小值以防止出现测距波束穿越地形情况出现。
若||rSiFj||-hTerrain(lonij,latij)≤0则停止搜索,粗搜索得到波束i与地形交点位置即为rSiFj,否则j=j+1继续进行下一次搜索。其中hTerrain(lonij,latij)表示在经纬度lonij、latij处地形高程。
(3)采用二分法精细搜索波束与地形交点
当由粗搜索得到波束i与地形交点位置rSiFj,令第k次精细搜索得到波束i方向上位置rSiFjk(k=0,1,2,…)及其对应经纬度lonijk、latijk,则
当k=0时
rSiFj0=rSiFj
lonij0=lonij
latij0=latij (8)
当k>0时
其中
rSiFjk-1为第k-1次搜索得到的波束i方向上位置、CFB为从着陆器本体坐标系到待测星体固连坐标系的转换阵,dSi为测距敏感器第i个波束在本体系下的指向,Δl为粗搜索步长,k为精搜索次数,hTerrain(lonijk,latijk)表示在经纬度lonijk、latijk处地形高程。
第三步,计算测距波束i的距离测量值li
li=||rSiFjk-rSF|| (11)
本发明说明书中未作详细描述的内容属本领域技术人员的公知技术。
Claims (7)
1.一种考虑地形起伏的测距敏感器仿真建模方法,其特征在于通过下述方式实现:
S1、以当前星下点地形高程所对应的待测星体表面作为基准平面,确定测距波束所在地形范围内地形最大高程,具体通过下述方式实现:
利用测距敏感器安装以及着陆器的位置和姿态信息确定测距波束与基准平面交点的经纬度;
以当前星下点经纬度和测距波束与基准面交点经纬度为矩形对角点确定测距波束所在地形范围;
搜索该范围内数字地形最大高程确定测距波束所在地形范围内地形的最大高程;
S2、根据确定的地形最大高程确定波束搜索起点,按照预设的步长进行粗搜索,搜索波束与地形交点,当波束某个位置对应高程小于地形高程时,采用二分法精细搜索波束与地形的交点,直到测距精度满足要求;
通过下述方式确定波束搜索起点rSiF0:
其中,||·||表示向量求模运算,hTerrainMaxi为测距波束所在地形范围内的地形最大高程、CFB为从着陆器本体坐标系到待测星体固连坐标系的转换阵,dSi为测距敏感器第i个波束在本体系下的指向,rSF为测距敏感器在待测星体固连坐标系下的位置;
通过下述方式确定测距波束所在地形范围内地形的最大高程:
hTerrainMaxi=max(hTerrain(lon,lat)) lon0≤lon≤loni且lat0≤lat≤lati
其中,lon0、lat0为测距敏感器星下点经纬度,loni、lati为测距敏感器第i个波束与基准平面交点处经纬度,max(·)表示求最大值运算,hTerrain(lon,lat)表示在经纬度lon、lat处地形高程。
2.根据权利要求1所述的方法,其特征在于:通过下述方式确定测距敏感器星下点经纬度lon0、lat0和测距敏感器到星下点的高度h:
rSF=[rSFx rSFy rSFz]T=rF+CFBrS
lon0=atan2(rSFy,rSFx)
h=||rSF||-hTerrain(lon0,lat0)
其中,||·||表示向量求模运算,rF为着陆器在待测星体固连坐标系下位置矢量,CFB为从着陆器本体坐标系到待测星体固连坐标系的转换阵,rS为测距敏感器在着陆器本体坐标系下安装位置,rSF为测距敏感器在待测星体固连坐标系下的位置,[rSFx rSFy rSFz]T为rSF的三个分量,hTerrain(lon0,lat0)表示在经纬度lon0、lat0处地形高程。
4.根据权利要求1所述的方法,其特征在于:所述的粗搜索过程如下:
按照下式计算第j次搜索时波束i方向上位置rSiFj及其对应经纬度lonij、latij
rSiFj=rSiF0+(j·Δl)CFBdSi=[rSiFjx rSiFjy rSiFjz]T
lonij=arctan2(rSiFjy,rSiFjx)
若||rSiFj||-hTerrain(lonij,latij)≤0则停止搜索,粗搜索得到波束i与地形交点位置即为rSiFj,否则将j+1值赋给j继续进行下一次搜索,直至得到满足条件的rSiFj;
上述公式中||·||表示向量求模运算,rSiF0为波束搜索起点,CFB为从着陆器本体坐标系到待测星体固连坐标系的转换阵,dSi为测距敏感器第i个波束在本体系下的指向,Δl为粗搜索步长,[rSiFjx rSiFjy rSiFjz]T为rSiFj的三个分量,hTerrain(lonij,latij)表示在经纬度lonij、latij处地形高程。
6.根据权利要求1所述的方法,其特征在于:所述的精细搜索过程如下:
按照下式进行第k次精细搜索得到波束i方向上位置rSiFjk及其对应经纬度lonijk、latijk;
当k=0时
rSiFj0=rSiFj
lonij0=lonij
latij0=latij
当k>0时
lonijk=arctan2(rSiFjky,rSiFjkx)
其中
上述公式中||·||表示向量求模运算,rSiFjk-1为第k-1次搜索得到的波束i方向上位置、CFB为从着陆器本体坐标系到待测星体固连坐标系的转换阵,dSi为测距敏感器第i个波束在本体系下的指向,Δl为粗搜索步长,k为精搜索次数,hTerrain(lonijk,latijk)表示在经纬度lonijk、latijk处地形高程。
7.根据权利要求1所述的方法,其特征在于:通过下述公式计算测距波束i的距离测量值li:
li=||rSiFjk-rSF||
式中,rSiFjk为精细搜索下波束i与地形交点在待测星体固连坐标系下的位置;rSF为测距敏感器在待测星体固连坐标系下的位置。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010682507.5A CN111896027B (zh) | 2020-07-15 | 2020-07-15 | 一种考虑地形起伏的测距敏感器仿真建模方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010682507.5A CN111896027B (zh) | 2020-07-15 | 2020-07-15 | 一种考虑地形起伏的测距敏感器仿真建模方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111896027A CN111896027A (zh) | 2020-11-06 |
CN111896027B true CN111896027B (zh) | 2022-07-29 |
Family
ID=73191293
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010682507.5A Active CN111896027B (zh) | 2020-07-15 | 2020-07-15 | 一种考虑地形起伏的测距敏感器仿真建模方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111896027B (zh) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103900576A (zh) * | 2014-03-31 | 2014-07-02 | 北京控制工程研究所 | 一种深空探测自主导航的信息融合方法 |
CN109633580A (zh) * | 2018-12-25 | 2019-04-16 | 西安空间无线电技术研究所 | 一种考虑地形高度的天基雷达杂波仿真方法及装置 |
CN110929402A (zh) * | 2019-11-22 | 2020-03-27 | 哈尔滨工业大学 | 一种基于不确定分析的概率地形估计方法 |
CN111319802A (zh) * | 2020-03-18 | 2020-06-23 | 北京理工大学 | 小天体附着轨迹自适应曲率匹配制导方法 |
Family Cites Families (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8321074B1 (en) * | 2012-08-03 | 2012-11-27 | Bell Douglas A | Altitude data generation system, device, and method |
US9404754B2 (en) * | 2013-03-25 | 2016-08-02 | Raytheon Company | Autonomous range-only terrain aided navigation |
CN103592632B (zh) * | 2013-10-16 | 2015-05-27 | 北京控制工程研究所 | 一种适用于月球着陆过程的测距测速波束指向确定方法 |
CN103674034B (zh) * | 2013-12-26 | 2015-12-30 | 北京控制工程研究所 | 多波束测速测距修正的鲁棒导航方法 |
CN108389255B (zh) * | 2018-02-12 | 2021-05-11 | 西安电子科技大学 | 基于分层高程云图的地形几何参数提取方法 |
CN109269512B (zh) * | 2018-12-06 | 2021-05-04 | 北京理工大学 | 行星着陆图像与测距融合的相对导航方法 |
CN110826174B (zh) * | 2019-09-25 | 2023-07-14 | 北京控制工程研究所 | 一种动力下降落月仿真过程中考虑三维地形的着陆器相对月面的距离确定方法 |
CN111076709A (zh) * | 2019-12-19 | 2020-04-28 | 中国农业大学 | 一种地形测量方法及*** |
-
2020
- 2020-07-15 CN CN202010682507.5A patent/CN111896027B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103900576A (zh) * | 2014-03-31 | 2014-07-02 | 北京控制工程研究所 | 一种深空探测自主导航的信息融合方法 |
CN109633580A (zh) * | 2018-12-25 | 2019-04-16 | 西安空间无线电技术研究所 | 一种考虑地形高度的天基雷达杂波仿真方法及装置 |
CN110929402A (zh) * | 2019-11-22 | 2020-03-27 | 哈尔滨工业大学 | 一种基于不确定分析的概率地形估计方法 |
CN111319802A (zh) * | 2020-03-18 | 2020-06-23 | 北京理工大学 | 小天体附着轨迹自适应曲率匹配制导方法 |
Non-Patent Citations (1)
Title |
---|
"Effect of DEM Interpolation Neighbourhood on Terrain Factors";Zhu, Y 等;《ISPRS INTERNATIONAL JOURNAL OF GEO-INFORMATION》;20191230;第8卷(第1期);正文第1-24页 * |
Also Published As
Publication number | Publication date |
---|---|
CN111896027A (zh) | 2020-11-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US11313698B2 (en) | Method for initial alignment of radar assisted airborne strapdown inertial navigation system | |
CN106052688B (zh) | 基于地形轮廓匹配的惯性导航***速度累积误差修正方法 | |
CN110779521A (zh) | 一种多源融合的高精度定位方法与装置 | |
US5345241A (en) | Self-contained method for correction of an inertial system over a body of water | |
CN103900576B (zh) | 一种深空探测自主导航的信息融合方法 | |
CN101246012B (zh) | 一种基于鲁棒耗散滤波的组合导航方法 | |
CN103674034B (zh) | 多波束测速测距修正的鲁棒导航方法 | |
CN112325886B (zh) | 一种基于重力梯度仪和陀螺仪组合的航天器自主定姿*** | |
Xu et al. | A reliable hybrid positioning methodology for land vehicles using low-cost sensors | |
Kim et al. | Error analysis of a low-low satellite-to-satellite tracking mission | |
CN113252038B (zh) | 基于粒子群算法的航迹规划地形辅助导航方法 | |
WO2016203744A1 (ja) | 測位装置 | |
CN111024074B (zh) | 一种基于递推最小二乘参数辨识的惯导速度误差确定方法 | |
CN110207721B (zh) | 一种考虑残差分布的无效地形匹配结果识别方法 | |
CN109855623A (zh) | 基于Legendre多项式和BP神经网络的地磁模型在线逼近方法 | |
CN106872962A (zh) | 一种用于星载激光测高仪检校的地面探测器的布设方法 | |
Dou et al. | A novel polarized skylight navigation model for bionic navigation with marginalized unscented Kalman filter | |
CN115265532A (zh) | 一种用于船用组合导航中的辅助滤波方法 | |
CN109855652B (zh) | 星载激光测高仪指向角误差为非常数时的在轨标定方法 | |
CN114111767B (zh) | 基于多信息融合对线路设计线型进行优化的方法 | |
CN109668562A (zh) | 一种考虑偏差时引入伪测量的重力梯度运动学导航方法 | |
CN111896027B (zh) | 一种考虑地形起伏的测距敏感器仿真建模方法 | |
Choe et al. | Lightweight marginalized particle filtering with enhanced consistency for terrain referenced navigation | |
Ye et al. | Deep Gaussian process regression for performance improvement of POS during GPS outages | |
CN115392540A (zh) | 一种用于月球轨道交会制导的快速预报方法 |
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 |