CN111272177A - 基于时间对准的间接滤波相对导航方法及*** - Google Patents
基于时间对准的间接滤波相对导航方法及*** Download PDFInfo
- Publication number
- CN111272177A CN111272177A CN202010129670.9A CN202010129670A CN111272177A CN 111272177 A CN111272177 A CN 111272177A CN 202010129670 A CN202010129670 A CN 202010129670A CN 111272177 A CN111272177 A CN 111272177A
- Authority
- CN
- China
- Prior art keywords
- relative
- navigation
- time
- position information
- satellite
- 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
- 238000001914 filtration Methods 0.000 title claims abstract description 61
- 238000000034 method Methods 0.000 title claims abstract description 32
- 238000005259 measurement Methods 0.000 claims abstract description 43
- 239000011159 matrix material Substances 0.000 claims description 14
- 238000007781 pre-processing Methods 0.000 claims description 10
- 238000009434 installation Methods 0.000 claims description 6
- 238000003032 molecular docking Methods 0.000 claims description 6
- 238000012958 reprocessing Methods 0.000 claims description 5
- 230000008569 process Effects 0.000 abstract description 7
- 238000012545 processing Methods 0.000 abstract description 6
- 238000004088 simulation Methods 0.000 description 7
- 238000013459 approach Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- JXASPPWQHFOWPL-UHFFFAOYSA-N Tamarixin Natural products C1=C(O)C(OC)=CC=C1C1=C(OC2C(C(O)C(O)C(CO)O2)O)C(=O)C2=C(O)C=C(O)C=C2O1 JXASPPWQHFOWPL-UHFFFAOYSA-N 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 239000000446 fuel Substances 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Images
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
- 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/10—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration
- G01C21/12—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration executed aboard the object being navigated; Dead reckoning
- G01C21/16—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration executed aboard the object being navigated; Dead reckoning by integrating acceleration or speed, i.e. inertial navigation
- G01C21/165—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration executed aboard the object being navigated; Dead reckoning by integrating acceleration or speed, i.e. inertial navigation combined with non-inertial navigation instruments
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Automation & Control Theory (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Navigation (AREA)
Abstract
本发明公开一种基于时间对准的间接滤波相对导航方法及***,利用测量单机对目标星进行相对测量,对相对测量的结果进行滤波处理,便可以得到需要的相对导航信息。测量单机的输出频率与导航滤波算法周期的不一致,导致数据之间存延时,通过对两者的时间进行对准,降低延迟对导航滤波算法精度的影响;受限于测量范围的限制,通常追踪星在超近距离段会配备多台测量单机,每台测量单机输出的相对测量信息的定于不一致,为了在切换不同测量单机进行相对导航滤波的过程中,导航算法不需要重新收敛,对不同测量单机的输出量进行预先处理,使用间接滤波的方法进行相对导航的滤波输出。
Description
技术领域
本发明涉及基于时间对准的间接滤波相对导航方法及***,属于空间导航技术领域。
背景技术
造价昂贵的航天器受恶劣空间环境影响,导致部分设备老化损坏,失去一些功能,造成巨额经济损失,且失效的航天器长期占用宝贵的轨道资源,对其他航天器造成潜在的安全威胁。面对这样的情形,如何对这些失效航天器进行设备维修、轨道修正以及燃料加注等以使其能够继续正常运行工作,已经成为各国航天领域重点研究的问题之一,在轨服务技术随之产生,具有广泛的应用价值。
对目标航天器进行在轨服务的前提是完成对它的相对测量和相对导航。远距离和近距离阶段可以利用微波雷达对目标星进行点目标的探测,以此来完成相对导航。但是超近距离段,需要进行高精度的体目标相对导航,利用微波雷达已经不能满足使用需求,因此可以采用激光雷达和视觉相机获得目标星与追踪星之间的相对信息,以此来完成高精度的相对导航,方法简单可靠性高。
目前点目标的相对导航以及体目标的合作目标相对导航方法日趋成熟。目前超近距离非合作目标的相对导航算法处理过程中,测量单机输出频率低于***导航算法的处理频率,限制了导航算法的精度;不同测量单机测量输出的相对信息定义不一致,导致使用不同单机测量数据时需要设计多个滤波器,增加了算法的复杂程度。
发明内容
本发明解决的技术问题是:克服现有技术的不足,提供了基于时间对准的间接滤波相对导航方法及***,降低了使用多种测量单机进行非合作目标超近距离相对导航算法的复杂度的同时,提高了算法的精度。
本发明的技术解决方案是:基于时间对准的间接滤波相对导航方法,包括如下步骤:
利用追踪星上的测量单机获取目标星与追踪星之间相对位置信息和相对姿态信息;
对获取的相对位置信息和相对姿态信息进行预处理,使相对位置信息和相对姿态信息在时间维度上递推至追踪星导航滤波算法的时刻,得到时间对准后的相对位置信息和相对姿态信息;
对时间对准后的相对位置信息和相对姿态信息进行再处理,统一不同源测量单机输出的相对位置信息和相对姿态信息,并将统一后的相对位置信息和相对姿态信息作为基于扩展卡尔曼滤波的间接滤波相对导航算法的输入,得到的输出用于进行后续控制。
进一步地,所述对获取的相对位置信息和相对姿态信息进行预处理的方法为:利用导航滤波算法输出的追踪星和目标星的相对速度乘以测量单机输出数据的时间与追踪星导航滤波算法时刻的时间差,加上测量单机输出的相对位置信息,实现将测量单机输出的相对位置信息递推至追踪星导航滤波算法的时刻,完成时间维度上的对准。
进一步地,所述测量单机包括激光雷达和视觉相机。
进一步地,当测量单机为激光雷达时,目标星轨道下两星的相对位置为:
进一步地,当测量单机为视觉相机时,目标星轨道系下两星的相对位置为:
其中,为视觉相机在追踪星本体系下的安装位置,为经过时间对准后的视觉相机输出的相对位置信息,为目标星质心在目标星形心坐标系下的坐标,为目标星对接环中心在目标星形心坐标系下的坐标,为目标星形心坐标系到追踪星本体坐标系的旋转矩阵。
基于时间对准的间接滤波相对导航***,包括
第一模块,利用追踪星上的测量单机获取目标星与追踪星之间相对位置信息和相对姿态信息;
第二模块,对获取的相对位置信息和相对姿态信息进行预处理,使相对位置信息和相对姿态信息在时间维度上递推至追踪星导航滤波算法的时刻,得到时间对准后的相对位置信息和相对姿态信息;
第三模块,对时间对准后的相对位置信息和相对姿态信息进行再处理,统一不同源测量单机输出的相对位置信息和相对姿态信息,并将统一后的相对位置信息和相对姿态信息作为基于扩展卡尔曼滤波的间接滤波相对导航算法的输入,得到的输出用于进行后续控制。
进一步地,所述对获取的相对位置信息和相对姿态信息进行预处理,具体的方法为:利用导航滤波算法输出的追踪星和目标星的相对速度乘以测量单机输出数据的时间与追踪星导航滤波算法时刻的时间差,加上测量单机输出的相对位置信息,实现将测量单机输出的相对位置信息递推至追踪星导航滤波算法的时刻,完成时间维度上的对准。
进一步地,所述测量单机包括激光雷达和视觉相机。
进一步地,当测量单机为激光雷达时,目标星轨道下两星的相对位置为:
进一步地,当测量单机为视觉相机时,目标星轨道系下两星的相对位置为:
其中,为视觉相机在追踪星本体系下的安装位置,为经过时间对准后的视觉相机输出的相对位置信息,为目标星质心在目标星形心坐标系下的坐标,为目标星对接环中心在目标星形心坐标系下的坐标,为目标星形心坐标系到追踪星本体坐标系的旋转矩阵。
本发明与现有技术相比的优点在于:
(1)本发明算法处理的过程中考虑了测量单机输出频率与***导航算法处理频率不一致的问题,在时间维度上进行统一,提高了算法的精度;
(2)本发明设计的算法,将不同测量单机输出的相对信息进行预处理,设计基于扩展卡尔曼的间接滤波器,降低了算法的复杂度。
附图说明
图1为本发明方法流程示意图;
图2为本发明激光雷达测量输出示意图;
图3为本发明视觉相机测量输出示意图;
图4为本发明实施例导航输出相对位置误差曲线;
图5为本发明实施例导航输出相对位置误差曲线;
图6为本发明实施例单机切换处导航输出相对位置误差曲线。
具体实施方式
本发明为基于时间对准的间接滤波相对导航方法,利用时间对准的方法解决数据延时对***导航精度的影响,同时利用间接滤波的方法解决不同测量单机切换使用的过程中,导航滤波算法需要重新收敛导致***性能下降的问题。利用测量单机输出的相对信息,通过预处理将测量单机输出的数据递推到***导航滤波算法的时间,完成时间维度上的对准。对时间对准完成后的数据进一步进行预处理,统一不同测量单机输出的相对测量信息,相对导航滤波算法以此做为观测量,以目标星轨道系下两星的相对位置速度为状态量,最终可以获得相对轨道控制所需要的目标星轨道下的两星相对位置和相对速度。
如图1,基于时间对准的间接滤波相对导航方法,本发明包括:
步骤一,利用追踪星上的测量单机获取目标星与追踪星之间相对位置信息和相对姿态信息;
步骤二,对获取的相对位置信息和相对姿态信息进行预处理,使相对位置信息和相对姿态信息在时间维度上递推至追踪星导航滤波算法的时刻,得到时间对准后的相对位置信息和相对姿态信息;
步骤三,对时间对准后的相对位置信息和相对姿态信息进行再处理,统一不同源测量单机输出的相对位置信息和相对姿态信息,并将统一后的相对位置信息和相对姿态信息作为基于扩展卡尔曼滤波的间接滤波相对导航算法的输入,得到的输出用于进行后续控制。
上述步骤二,对测量单机输出的数据进行预处理,利用导航滤波算法上一个周期输出的目标星与追踪星之间的相对速度,即卡尔曼滤波上一周期的输出值的后三维量测量单机输出数据的时间tld或txj与***导航滤波算法时间tGNC的差记为Δtld=tGNC-tld或者Δtxj=tGNC-txj,时间对准完成后送给导航滤波算法的输入值为或者
上述步骤三,如图2所示,可知当测量单机为激光雷达时,目标星与追踪星质心连线在追踪星本体系下的投影:
如图3所示,可知当测量单机为视觉相机时,目标星与追踪星质心连线在追踪星本体系下的投影:
其中,为视觉相机在追踪星本体系下的安装位置,为经过时间对准后的视觉相机输出的相对位置信息,为目标星质心在目标星形心坐标系下的坐标,为目标星对接环中心在目标星形心坐标系下的坐标,为目标星形心坐标系到追踪星本体坐标系的旋转矩阵。
由以上分析可知,当测量单机为激光雷达时,目标星轨道下两星的相对位置为:
当测量单机为视觉相机时,目标星轨道系下两星的相对位置为:
以目标星轨道系下目标星与追踪星的相对位置为观测量,建立扩展卡尔曼滤波器,观测方程为:
其中,H=[I3×3 03×3],X=[x y z vx vy vz]为卡尔曼滤波算法中的状态量,本方法中为目标星轨道系下两星的相对位置和相对速度,
结合两星的相对轨道动力学方程采用滤波技术实时进行输出。
实施例
下面以某高轨卫星GNC***为对象,采用激光雷达与视觉相机作为测量单机,通过仿真实例验证本发明在超近距离非合作目标相对导航方面的应用。
(1)状态方程建立
将相对轨道的坐标系建立在目标星轨道系,以两星的相对位置、相对速度、目标质心坐标作为状态量,其状态方程为:
其中:X=[x y z vx vy vz]T;
x,y,z为两星相对位置在目标星轨道系下的投影;
vx,vy,vz为两星相对速度在目标星轨道系下的投影;
进行离散化处理后,求取一步预测值:
(2)观测方程建立
以单机输出的相对位置建立观测方程,由单机测量信息定义,可得
当测量单机为激光雷达时,目标星轨道下两星的相对位置为:
当测量单机为视觉相机时,目标星轨道系下两星的相对位置为:
以目标星轨道系下目标星与追踪星的相对位置为观测量,建立扩展卡尔曼滤波器,观测方程为:
其中H=[I3×3 03×3],X=[x y z vx vy vz]为卡尔曼滤波算法中的状态量,此算例中为目标星轨道系下两星的相对位置和相对速度。
(3)仿真条件
激光雷达测量精度:
相对位置测量精度:≤0.1m;
相对角度测量误差:3°。
视觉相机测量精度:
相对位置测量精度:≤0.05m;
相对角度测量误差:2°。
仿真轨迹:初始追踪星在目标后方60m→目标后方60m逼近到后方30m→目标后方30m逼近到后方10m→目标后方10m逼近到后方2m。
10m处单机由激光雷达切换成视觉相机(仿真时间1000s)
(3)仿真结果
如图4、5和6。
从仿真结果曲线可以看出:整个过程中相对位置误差优于0.05m;相对速度误差优于0.01m/s,单机切换处导航算法无需重新收敛,可以满足超近距离非合作目标相对导航的需求。以上仿真说明了基于时间对准的间接滤波相对导航方法的有效性。
本发明说明书中未作详细描述的内容属本领域技术人员的公知技术。
Claims (10)
1.基于时间对准的间接滤波相对导航方法,其特征在于,包括如下步骤:
利用追踪星上的测量单机获取目标星与追踪星之间相对位置信息和相对姿态信息;
对获取的相对位置信息和相对姿态信息进行预处理,使相对位置信息和相对姿态信息在时间维度上递推至追踪星导航滤波算法的时刻,得到时间对准后的相对位置信息和相对姿态信息;
对时间对准后的相对位置信息和相对姿态信息进行再处理,统一不同源测量单机输出的相对位置信息和相对姿态信息,并将统一后的相对位置信息和相对姿态信息作为基于扩展卡尔曼滤波的间接滤波相对导航算法的输入,得到的输出用于进行后续控制。
2.根据权利要求1所述的基于时间对准的间接滤波相对导航方法,其特征在于,所述对获取的相对位置信息和相对姿态信息进行预处理的方法为:利用导航滤波算法输出的追踪星和目标星的相对速度乘以测量单机输出数据的时间与追踪星导航滤波算法时刻的时间差,加上测量单机输出的相对位置信息,实现将测量单机输出的相对位置信息递推至追踪星导航滤波算法的时刻,完成时间维度上的对准。
3.根据权利要求1所述的基于时间对准的间接滤波相对导航方法,其特征在于,所述测量单机包括激光雷达和视觉相机。
6.基于时间对准的间接滤波相对导航***,其特征在于,包括
第一模块,利用追踪星上的测量单机获取目标星与追踪星之间相对位置信息和相对姿态信息;
第二模块,对获取的相对位置信息和相对姿态信息进行预处理,使相对位置信息和相对姿态信息在时间维度上递推至追踪星导航滤波算法的时刻,得到时间对准后的相对位置信息和相对姿态信息;
第三模块,对时间对准后的相对位置信息和相对姿态信息进行再处理,统一不同源测量单机输出的相对位置信息和相对姿态信息,并将统一后的相对位置信息和相对姿态信息作为基于扩展卡尔曼滤波的间接滤波相对导航算法的输入,得到的输出用于进行后续控制。
7.根据权利要求6所述的基于时间对准的间接滤波相对导航***,其特征在于,所述对获取的相对位置信息和相对姿态信息进行预处理,具体的方法为:利用导航滤波算法输出的追踪星和目标星的相对速度乘以测量单机输出数据的时间与追踪星导航滤波算法时刻的时间差,加上测量单机输出的相对位置信息,实现将测量单机输出的相对位置信息递推至追踪星导航滤波算法的时刻,完成时间维度上的对准。
8.根据权利要求6所述的基于时间对准的间接滤波相对导航***,其特征在于,所述测量单机包括激光雷达和视觉相机。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010129670.9A CN111272177B (zh) | 2020-02-28 | 2020-02-28 | 基于时间对准的间接滤波相对导航方法及*** |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010129670.9A CN111272177B (zh) | 2020-02-28 | 2020-02-28 | 基于时间对准的间接滤波相对导航方法及*** |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111272177A true CN111272177A (zh) | 2020-06-12 |
CN111272177B CN111272177B (zh) | 2022-03-04 |
Family
ID=70997422
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010129670.9A Active CN111272177B (zh) | 2020-02-28 | 2020-02-28 | 基于时间对准的间接滤波相对导航方法及*** |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111272177B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113607176A (zh) * | 2021-10-11 | 2021-11-05 | 智道网联科技(北京)有限公司 | 组合导航***轨迹输出方法及装置 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050060092A1 (en) * | 2003-08-05 | 2005-03-17 | The Boeing Company | Laser range finder closed-loop pointing technology of relative navigation, attitude determination, pointing and tracking for spacecraft rendezvous |
CN107529376B (zh) * | 2013-08-01 | 2015-12-30 | 上海新跃仪表厂 | 多模融合的微小卫星非合作目标相对导航的方法 |
CN105486305A (zh) * | 2014-09-17 | 2016-04-13 | 上海新跃仪表厂 | 一种估计加速度计漂移的近程相对导航滤波方法 |
CN110388941A (zh) * | 2019-08-28 | 2019-10-29 | 北京机械设备研究所 | 一种基于自适应滤波的车辆姿态对准方法 |
-
2020
- 2020-02-28 CN CN202010129670.9A patent/CN111272177B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050060092A1 (en) * | 2003-08-05 | 2005-03-17 | The Boeing Company | Laser range finder closed-loop pointing technology of relative navigation, attitude determination, pointing and tracking for spacecraft rendezvous |
CN107529376B (zh) * | 2013-08-01 | 2015-12-30 | 上海新跃仪表厂 | 多模融合的微小卫星非合作目标相对导航的方法 |
CN105486305A (zh) * | 2014-09-17 | 2016-04-13 | 上海新跃仪表厂 | 一种估计加速度计漂移的近程相对导航滤波方法 |
CN110388941A (zh) * | 2019-08-28 | 2019-10-29 | 北京机械设备研究所 | 一种基于自适应滤波的车辆姿态对准方法 |
Non-Patent Citations (6)
Title |
---|
QIYANG HU等: "A decoupled Method for Lidar-based Close-range Relative Navigation of Tumbling Non-cooperative Target", 《2019 CHINESE CONTROL AND DECISION CONFERENCE (CCDC)》 * |
刘翠翠等: "非等间隔滤波及其在导航***中的应用综述", 《航天控制》 * |
张晓彤等: "空间非合作机动目标相对导航方法概述", 《第二届中国空天安全会议论文集》 * |
彭杨等: "非合作目标超近距离光学相对导航方法及半物理仿真***研究", 《上海航天》 * |
王大轶: "非合作航天器自主相对导航研究综述", 《控制理论与应用》 * |
郭碧波等: "航天器最终逼近段的相对导航研究与半实物仿真验证", 《宇航学报》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113607176A (zh) * | 2021-10-11 | 2021-11-05 | 智道网联科技(北京)有限公司 | 组合导航***轨迹输出方法及装置 |
CN113607176B (zh) * | 2021-10-11 | 2021-12-10 | 智道网联科技(北京)有限公司 | 组合导航***轨迹输出方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN111272177B (zh) | 2022-03-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Peng et al. | Pose measurement and motion estimation of space non-cooperative targets based on laser radar and stereo-vision fusion | |
CN103852082B (zh) | 一种星间测量和陀螺的姿态轨道一体化滤波估计方法 | |
Felicetti et al. | Image-based attitude maneuvers for space debris tracking | |
CN107529376B (zh) | 多模融合的微小卫星非合作目标相对导航的方法 | |
Xu et al. | A pose measurement method of a non-cooperative GEO spacecraft based on stereo vision | |
CN113443173B (zh) | 运载火箭的导航***、方法、装置、电子设备及存储介质 | |
CN110412868B (zh) | 一种使用星间光学图像的非合作航天器轨道确定方法 | |
Wolf et al. | Toward improved landing precision on Mars | |
CN116540285B (zh) | 惯性辅助的gnss双天线定向方法、装置与电子设备 | |
Ivanov et al. | Relative pose and inertia determination of unknown satellite using monocular vision | |
CN111272177B (zh) | 基于时间对准的间接滤波相对导航方法及*** | |
Allerton et al. | Distributed data fusion algorithms for inertial network systems | |
Cropp | Pose estimation and relative orbit determination of a nearby target microsatellite using passive imagery | |
CN113359159B (zh) | 一种天基空间目标监视的多星定位定轨精度快速分析方法 | |
Capuano et al. | VINAG: A highly integrated system for autonomous on-board absolute and relative spacecraft navigation | |
CN105486305B (zh) | 一种估计加速度计漂移的近程相对导航滤波方法 | |
CN115855063A (zh) | 基于绝对姿态递推修正的交会对接敏感器数据预处理方法 | |
Xu et al. | Rotational motion estimation of non-cooperative target in space based on the 3D point cloud sequence | |
CN113608452A (zh) | 重复使用运载火箭导航***的半实物仿真平台 | |
Feetham et al. | Single camera absolute motion based digital elevation mapping for a next generation planetary lander | |
Du et al. | Pose determination of large non-cooperative satellite in close range using coordinated cameras | |
CN111928858B (zh) | 一种可回收航天运载器的导航方法 | |
Goecks et al. | Virtual reality for enhanced 3D astronaut situational awareness during robotic operations in space | |
CN112082551B (zh) | 一种可回收航天运载器的导航*** | |
CN116502399B (zh) | 基于stk和matlab联合仿真的卫星轨道生成方法和生成器 |
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 |