CN104280025B - 基于无色卡尔曼滤波的深海机器人超短基线组合导航方法 - Google Patents
基于无色卡尔曼滤波的深海机器人超短基线组合导航方法 Download PDFInfo
- Publication number
- CN104280025B CN104280025B CN201310284641.XA CN201310284641A CN104280025B CN 104280025 B CN104280025 B CN 104280025B CN 201310284641 A CN201310284641 A CN 201310284641A CN 104280025 B CN104280025 B CN 104280025B
- Authority
- CN
- China
- Prior art keywords
- ultra
- short baseline
- kalman filter
- information
- robot
- 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
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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/38—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
- G01S19/39—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/42—Determining position
- G01S19/45—Determining position by combining measurements of signals from the satellite radio beacon positioning system with a supplementary measurement
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/38—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
- G01S19/39—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/42—Determining position
- G01S19/45—Determining position by combining measurements of signals from the satellite radio beacon positioning system with a supplementary measurement
- G01S19/47—Determining position by combining measurements of signals from the satellite radio beacon positioning system with a supplementary measurement the supplementary measurement being an inertial measurement, e.g. tightly coupled inertial
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Computer Networks & Wireless Communication (AREA)
- Navigation (AREA)
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
Abstract
本发明涉及一种基于无色卡尔曼滤波的深海机器人超短基线组合导航方法,利用全球定位系统获取深海机器人的初始绝对位置作为航迹推算的初始点,并采集初始信息;采用无色卡尔曼滤波方法对采集到的初始信息进行数据融合,得出融合后的信息。本方法可以提高使用超短基线定位系统的深海机器人的导航精度,同时能够平滑深海机器人控制系统所需的航向、深度以及载体坐标系下的速度信息。
Description
技术领域
本发明涉及海洋工程领域的一种深海机器人导航技术,具体地说是一种基于无色卡尔曼滤波的超短基线组合导航方法。
背景技术
精确的导航能力是水下机器人有效作业和安全回收的关键。导航系统必须提供远距离及长时间范围内的精确定位、速度和姿态信息。但是由于受到水下环境的复杂性、机器人自身体积、重量、能源以及隐蔽性等因素的影响,实现高精度的水下机器人导航仍然是一项艰巨的任务。
现有的水下导航方法主要是多传感器组合导航方法,与陆地组合导航方法无异,采用运动学方程对多传感器信息进行滤波。但是针对深海机器人的导航问题,机器人在水中的运动状态更加复杂多变,声学定位数据有较多野值而且有很大的延迟,采用水中动力学方程才能更精确地刻画机器人在水中的运动状态,对于有延迟的位置信息正确融合对提高导航精度也至关重要。
发明内容
针对现有技术的不足,本发明提供一种基于无色卡尔曼滤波的深海机器人组合导航方法,目的在于提高使用超短基线定位系统的深海机器人的导航精度,同时能够平滑深海机器人控制系统所需的航向、深度以及载体坐标系下的速度信息。
本发明为实现上述目的所采用的技术方案是:
一种基于无色卡尔曼滤波的深海机器人超短基线组合导航方法,利用全球定位系统获取深海机器人的初始绝对位置作为航迹推算的初始点,并采集深海机器人的初始信息;采用无色卡尔曼滤波方法对采集到的初始信息进行数据融合,得出融合后的信息。
所述初始信息包括多普勒计程仪采集的线速度信息,运动传感器采集的姿态信息,超短基线声学定位装置采集的位置信息和深度计采集的深度信息。
所述超短基线声学定位装置通过位移修正方法采集位置信息。
所述位移修正方法为用当前时刻接收到的超短基线声学定位装置提供的位置信息加上当前时刻距上次水声通讯时刻时间段内深海机器人走过的位移作为当前时刻超短基线定位装置提供的位置观测值。
在无色卡尔曼滤波过程中,取深海机器人的动力学方程作为系统的状态方程。
所述状态方程建立如下:
其中,M为6×6维的惯性矩阵,包含离心力和哥氏力的6×6维矩阵,为6×6维的水动力矩阵,G(q)包含重力和浮力的6×1维矩阵,q和u分别为6×1维的状态和电机转矩矢量。
在无色卡尔曼滤波过程中的测量方程选择超短基线数据有效时刻的切换观测方程。
所述有效时刻为超短基线数据有更新并且不是野值的时刻。
所述测量方程定义如下:
当超短基线定位装置提供的观测值有效时为:
yk=xk+Nn
x=[ux,vx,wx,rx,ξx,ηx,ζx,ψx]T为状态矢量,其中ux、vx、wx和rx分别是前向速度,侧向速度,垂向速度和转艏角速度的状态量,ξx、ηx、ζx和ψx分别是北向位置、东向位置、深度和艏向角的状态量;
y=[uy,vy,wy,ry,ξy,ηy,ζy,ψy]T为观测矢量,其中uy、vy、wy、ry分别是前向速度、侧向速度、垂向速度和转艏角速度的观测量,ξy、ηy、ζy、ψy是北向位置、东向位置、深度和艏向角的观测量;
Nn为观测噪声;
当超短基线定位装置提供的观测值无效时为:
其中,y=[uy,vy,wy,ry,ζy,ψy]T为观测矢量。
本发明具有以下有益效果及优点:
1.本发明采用无色卡尔曼滤波作为数据融合方法,计算量小,在过程噪声和测量噪声的先验信息已知的情况下,可以达到较高的估计精度;
2.本发明使用深海机器人动力学模型作为系统的状态方程,更精确地描述深海机器人在水下的运动状态,提高了导航精度;
3.本发明提出了根据超短基线定位数据的有效性切换观测方程,解决了超短基线定位数据野值较多以及超短基线定位装置与其他传感器采集数据更新率不一致的问题;
4.本发明利用上一次水声通讯时刻位置的新息信息更新当前时刻的位置信息,解决了超短基线定位系统的位置观测值有延迟的问题。
附图说明
图1是本发明的组合导航方法的工作原理图;
图2是本发明的无色卡尔曼滤波算法的流程图;
图3是基于无色卡尔曼滤波算法的机器人超短基线组合导航的实施例结果。
具体实施方式
下面结合附图及实施例对本发明做进一步的详细说明。
如图1所示为本发明的组合导航方法的工作原理图:利用多普勒计程仪采集线速度信息,利用超短基线声学定位装置采集平面位置信息,利用深度计采集深度信息,利用运动传感器采集姿态信息,利用数据融合方法对上述信息进行融合,得到深海机器人在海中的三维运动轨迹,同时能够平滑深海机器人控制系统所需的航向以及载体坐标系下的速度信息。
(1)取载体坐标系下的速度矢量和大地坐标系下的位置矢量和艏向角作为滤波器状态矢量,定义为:x=[ux,vx,wx,rx,ξx,ηx,ζx,ψx]T,其中ux,vx,wx,rx是前向速度,侧向速度,垂向速度和转艏角速度的状态量,ξx,ηx,ζx,ψx是北向位置,东向位置,深度和艏向角的状态量。
取深海机器人的动力学方程作为系统的状态方程,深海机器人的动力学方程如下:
式中,M为6×6维的惯性矩阵,包含离心力和哥氏力的6×6维矩阵,为6×6维的水动力矩阵,G(q)包含重力和浮力的6×1维矩阵,q和u分别为6×1维的状态和电机转矩矢量。
本发明中所采用的深海机器人动力学方程中的矢量为上述矢量中对应前向速度ux,侧向速度vx,垂向速度wx和转艏角速度rx的四维。
除ux,vx,wx,rx状态矢量x中的ξx,ηx,ζx,ψx的状态方程如下:
Nv为系统噪声,Δt为传感器采样间隔。
(2)所述根据超短基线定位数据的有效性切换观测方程步骤如下:
超短基线数据有效时刻是指超短基线数据更新并且不是野值的时刻,
1.超短基线定位装置提供的观测值有效时:
观测矢量定义为:y=[uy,vy,wy,ry,ξy,ηy,ζy,ψy]T,其中uy,vy,wy,ry是前向速度,侧向速度,垂向速度和转艏角速度的观测量,ξy,ηy,ζy,ψy是北向位置,东向位置,深度和艏向角的观测量。观测值来自于多普勒计程仪和运动传感器,东向和北向位置观测值来自于超短基线定位装置,深度观测值来自于深度计,艏向角观测值来自于运动传感器。
测量方程定义为:yk=xk+Nn,Nn为观测噪声。
2.超短基线定位装置提供的观测值无效时:
观测矢量定义为:y=[uy,vy,wy,ry,ζy,ψy]T
测量方程定义为:
(3)由于超短基线定位装置的工作原理基于水声通讯,超短基线定位装置提供的位置观测值是上一次水声通讯时刻机器人的位置,所述利用上一次水声通讯时刻位置的新息信息更新当前时刻的位置信息步骤如下:
用当前时刻接收到的超短基线定位装置提供的位置信息加上当前时刻距上次水声通讯时刻这个时间段内机器人走过的位移作为当前时刻超短基线定位装置提供的位置观测值,经过无色卡尔曼滤波算法得到当前时刻位置的最优估计。
(4)如图2所示为本发明的无色卡尔曼滤波算法的流程图,所述的无色卡尔曼滤波估计算法是以基本的卡尔曼滤波算法为框架,通过非线性的无色变换来进行非线性系统的状态和误差协方差的递推和更新。
无色卡尔曼滤波估计实现过程:
对上一时刻即k-1由被估计量“先验”的均值和方差Pk-1经非线性无色变换得到的一组离散Sigma点χk-1,χk-1经系统状态方程计算更新后得到离散的Sigma点对计算得到一步估计的扩展状态均值和方差Pk|k-1;对和Pk|k-1再进行无色变化得到一组离散的Sigma点χk|k-1,χk|k-1经测量方程计算更新后得到离散的Sigma点γk|k-1,对γk|k-1计算得到深海机器人可观测状态的均值和方差的估计值;所得到的一步估计值和得到的可观测状态的估计值与传感器采集观测值yk,经计算估计得到系统估算值和Pk。
如图3所示为基于无色卡尔曼滤波算法的机器人超短基线组合导航的实施例结果实施例:
外场实验中,机器人以不超过2m/s的速度在水下航行并上浮,由于航行期间机器人处于深水,所以没有GPS的定位信息。图3给出了利用本发明基于无色卡尔曼滤波的深海机器人超短基线组合导航方法进行机器人导航的结果,并同超短基线设备的定位结果和航迹推算法的导航结果作比较。其中圆圈表示航迹推算法的导航结果,星号表示超短基线设备的定位结果,点表示基于无色卡尔曼滤波的深海机器人超短基线组合导航方法的导航结果。从图中可以明显看出,采用基于无色卡尔曼滤波的深海机器人超短基线组合导航方法估计出的机器人的轨迹可以克服超短基线设备定位值少、野值多和航迹推算法随时间增长的累积误差。此实施例表明本发明基于无色卡尔曼滤波的深海机器人超短基线组合导航方法可显著提高组合导航系统精度。
Claims (4)
1.一种基于无色卡尔曼滤波的深海机器人超短基线组合导航方法,其特征在于:
利用全球定位系统获取深海机器人的初始绝对位置作为航迹推算的初始点,并采集深海机器人的初始信息;
采用无色卡尔曼滤波方法对采集到的初始信息进行数据融合,得出融合后的信息;
在无色卡尔曼滤波过程中,取深海机器人的动力学方程作为系统的状态方程;
所述初始信息包括多普勒计程仪采集的线速度信息,运动传感器采集的姿态信息,超短基线声学定位装置采集的位置信息和深度计采集的深度信息;
所述超短基线声学定位装置通过位移修正方法采集位置信息;
所述位移修正方法为用当前时刻接收到的超短基线声学定位装置提供的位置信息加上当前时刻距上次水声通讯时刻时间段内深海机器人走过的位移作为当前时刻超短基线定位装置提供的位置观测值;
所述状态方程建立如下:
其中,M为6×6维的惯性矩阵,包含离心力和哥氏力的6×6维矩阵,为6×6维的水动力矩阵,G(q)包含重力和浮力的6×1维矩阵,q和u分别为6×1维的状态和电机转矩矢量。
2.根据权利要求1所述的基于无色卡尔曼滤波的深海机器人超短基线组合导航方法,其特征在于:在无色卡尔曼滤波过程中的测量方程选择超短基线数据有效时刻的切换观测方程。
3.根据权利要求2所述的基于无色卡尔曼滤波的深海机器人超短基线组合导航方法,其特征在于:所述有效时刻为超短基线数据有更新并且不是野值的时刻。
4.根据权利要求2所述的基于无色卡尔曼滤波的深海机器人超短基线组合导航方法,其特征在于:所述测量方程定义如下:
当超短基线定位装置提供的观测值有效时为:
x=[ux,vx,wx,rx,ξx,ηx,ζx,ψx]T为状态矢量,其中ux、vx、wx和rx分别是前向速度,侧向速度,垂向速度和转艏角速度的状态量,ξx、ηx、ζx和ψx分别是北向位置、东向位置、深度和艏向角的状态量;
y=[uy,vy,wy,ry,ξy,ηy,ζy,ψy]T为8维观测矢量,其中uy、vy、wy、ry分别是前向速度、侧向速度、垂向速度和转艏角速度的观测量,ξy、ηy、ζy、ψy是北向位置、东向位置、深度和艏向角的观测量;
Nn为观测噪声;
当超短基线定位装置提供的观测值无效时为:
其中,y'=[uy,vy,wy,ry,ζy,ψy]T为6维观测矢量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310284641.XA CN104280025B (zh) | 2013-07-08 | 2013-07-08 | 基于无色卡尔曼滤波的深海机器人超短基线组合导航方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310284641.XA CN104280025B (zh) | 2013-07-08 | 2013-07-08 | 基于无色卡尔曼滤波的深海机器人超短基线组合导航方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104280025A CN104280025A (zh) | 2015-01-14 |
CN104280025B true CN104280025B (zh) | 2017-06-23 |
Family
ID=52255153
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310284641.XA Active CN104280025B (zh) | 2013-07-08 | 2013-07-08 | 基于无色卡尔曼滤波的深海机器人超短基线组合导航方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104280025B (zh) |
Families Citing this family (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104964689A (zh) * | 2015-07-06 | 2015-10-07 | 南京航空航天大学 | 基于声源定位、ins的爬壁机器人导航定位系统及方法 |
CN105547290B (zh) * | 2015-12-29 | 2018-02-13 | 北京理工大学 | 一种基于超短基线定位系统的从潜器导航方法 |
CN107024226B (zh) * | 2016-02-01 | 2021-03-16 | 北京自动化控制设备研究所 | 一种基于惯导/dvl/usbl组合的惯导误差估计方法 |
CN107505597A (zh) * | 2017-08-03 | 2017-12-22 | 浙江大学 | 一种基于深度信息的超短基线定位系统及方法 |
CN107839859A (zh) * | 2017-10-19 | 2018-03-27 | 丁建玲 | 一种海底光缆巡检水下航行器及巡检方法 |
CN108490396B (zh) * | 2018-03-16 | 2022-09-16 | 哈尔滨工业大学(威海) | 一种基于水声信号卡尔曼滤波的超短基线定位方法 |
CN109001829B (zh) * | 2018-07-12 | 2019-11-05 | 中国人民解放军国防科技大学 | 一种捷联式水下动态重力测量仪 |
CN109737958A (zh) * | 2019-03-08 | 2019-05-10 | 哈尔滨工程大学 | 一种声学测速辅助的极区格网惯性导航误差抑制方法 |
CN110132281B (zh) * | 2019-05-21 | 2023-10-20 | 哈尔滨工程大学 | 一种基于询问应答模式的水下高速目标高精度自主声学导航方法 |
CN110207694A (zh) * | 2019-05-27 | 2019-09-06 | 哈尔滨工程大学 | 一种基于相对位置信息的极区格网惯导/超短基线组合导航方法 |
CN111076728A (zh) * | 2020-01-13 | 2020-04-28 | 东南大学 | 基于dr/usbl的深潜载人潜水器组合导航方法 |
CN113074725B (zh) * | 2021-05-11 | 2022-07-22 | 哈尔滨工程大学 | 一种基于多源信息融合的小型水下多机器人协同定位方法及系统 |
CN114111840B (zh) * | 2021-11-12 | 2022-08-26 | 哈尔滨工业大学 | 一种基于组合导航的dvl误差参数在线标定方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1361431A (zh) * | 2000-12-23 | 2002-07-31 | 林清芳 | 完全整合式导航定位方法和系统 |
CN1779485A (zh) * | 2004-11-17 | 2006-05-31 | 中国科学院沈阳自动化研究所 | 载人潜水器的组合导航定位方法 |
CN102519450A (zh) * | 2011-12-12 | 2012-06-27 | 东南大学 | 一种用于水下滑翔器的组合导航装置及方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20080077325A1 (en) * | 2006-09-25 | 2008-03-27 | Honeywell International Inc. | Systems and methods for a hybrid transition matrix |
-
2013
- 2013-07-08 CN CN201310284641.XA patent/CN104280025B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1361431A (zh) * | 2000-12-23 | 2002-07-31 | 林清芳 | 完全整合式导航定位方法和系统 |
CN1779485A (zh) * | 2004-11-17 | 2006-05-31 | 中国科学院沈阳自动化研究所 | 载人潜水器的组合导航定位方法 |
CN102519450A (zh) * | 2011-12-12 | 2012-06-27 | 东南大学 | 一种用于水下滑翔器的组合导航装置及方法 |
Also Published As
Publication number | Publication date |
---|---|
CN104280025A (zh) | 2015-01-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104280025B (zh) | 基于无色卡尔曼滤波的深海机器人超短基线组合导航方法 | |
CN102323586B (zh) | 一种基于海流剖面的uuv辅助导航方法 | |
CN103376452B (zh) | 一种用单台声信标修正水下机器人位置误差的方法 | |
CN104075715A (zh) | 一种结合地形和环境特征的水下导航定位方法 | |
CN103968838B (zh) | 一种基于极坐标系的auv曲线运动状态下的协同定位方法 | |
CN111596333B (zh) | 一种水下定位导航方法及系统 | |
CN1325932C (zh) | 载人潜水器的组合导航定位方法 | |
CN109579850B (zh) | 基于对水速度辅助惯导的深水智能导航方法 | |
CN103217160A (zh) | 一种船舶动力定位卡尔曼滤波方法 | |
CN105549611B (zh) | 自主水下机器人的海洋环境自识别的航迹精确跟踪方法 | |
CN107990891A (zh) | 基于长基线和信标在线标定的水下机器人组合导航方法 | |
CN105628016A (zh) | 一种基于超短基线的导航定位方法 | |
CN111076728A (zh) | 基于dr/usbl的深潜载人潜水器组合导航方法 | |
CN102795323A (zh) | 一种基于ukf的水下机器人状态和参数联合估计方法 | |
CN106643723A (zh) | 一种无人艇安全航行船位推算方法 | |
CN103115625A (zh) | 一种浮体横纵荡及升沉运动的测量方法及系统 | |
Liu et al. | Research into the integrated navigation of a deep-sea towed vehicle with USBL/DVL and pressure gauge | |
CN104280024B (zh) | 一种深水机器人组合导航装置和方法 | |
CN104280026B (zh) | 基于aukf的深海机器人长基线组合导航方法 | |
CN109813316B (zh) | 一种基于地形辅助的水下载体紧组合导航方法 | |
CN103292813A (zh) | 一种提高水面艇编队导航精度的信息滤波方法 | |
CN103090864A (zh) | 一种基于通信延时补偿的多水面无人艇协同定位方法 | |
CN113190025B (zh) | 一种适用于变结构水下机器人的运动控制方法 | |
CN102830408B (zh) | 一种惯性信息辅助的卫星接收机基带信号处理方法 | |
CN105890589B (zh) | 一种水下机器人单目视觉定位方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |