CN108267765B - 使用不完整gnss信息的mems组合导航建模方法 - Google Patents
使用不完整gnss信息的mems组合导航建模方法 Download PDFInfo
- Publication number
- CN108267765B CN108267765B CN201810205371.1A CN201810205371A CN108267765B CN 108267765 B CN108267765 B CN 108267765B CN 201810205371 A CN201810205371 A CN 201810205371A CN 108267765 B CN108267765 B CN 108267765B
- Authority
- CN
- China
- Prior art keywords
- information
- navigation
- gnss
- mathematical model
- dimensional
- 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
- 238000000034 method Methods 0.000 title claims abstract description 19
- 238000013178 mathematical model Methods 0.000 claims abstract description 33
- 238000012417 linear regression Methods 0.000 claims abstract description 20
- LIWAQLJGPBVORC-UHFFFAOYSA-N ethylmethylamine Chemical compound CCNC LIWAQLJGPBVORC-UHFFFAOYSA-N 0.000 claims abstract 5
- 238000004458 analytical method Methods 0.000 claims description 6
- 238000000605 extraction Methods 0.000 claims description 3
- 238000005259 measurement Methods 0.000 description 4
- XUIMIQQOPSSXEZ-UHFFFAOYSA-N Silicon Chemical compound [Si] XUIMIQQOPSSXEZ-UHFFFAOYSA-N 0.000 description 2
- 239000000969 carrier Substances 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 229910052710 silicon Inorganic materials 0.000 description 2
- 239000010703 silicon Substances 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000003754 machining Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
Images
Classifications
-
- 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)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
- Navigation (AREA)
Abstract
本发明公开了使用不完整GNSS信息的MEMS组合导航建模方法,具体按照以下步骤实施,第一步:确定实现组合导航***所需的外部GNSSNMEA输入报文类型,其接收顺序可向组合***发送配置命令;第二步:从外部GNSSNMEA信息的GGA报文中提取三维位置信息;第三步:从外部GNSS NMEA信息的RMC报文中提取三维速度信息;第四步:建立三维位置误差的数学模型;第五步:线性回归算法确定三维位置误差的数学模型参数;第六步:建立三维速度误差的数学模型;第七步:线性回归算法确定速度误差模型参数;第八步:运行组合导航算法。本方法突破了常规MEMS组合导航设备的实现限制,将组合导航***的实现难度大大降低。
Description
技术领域
本发明属于惯性导航与卫星导航技术领域,具体涉及一种使用不完整GNSS信息的MEMS组合导航建模方法。
背景技术
北斗卫星导航***是与美国GPS、俄罗斯GLONASS以及欧洲GALILEO并驾齐驱的覆盖全球的星基导航***(GNSS)。卫星导航***的显著特点是直接提供绝对定位坐标,且误差不随时间积累,但存在无线电信号易受遮挡和干扰的不足。惯性导航***(InertialNavigation System,INS)能够自主地进行连续的三维空间导航与测姿,不存在信号的电磁干扰或遮挡问题,但惯性导航的不足在于其导航误差会随时间积累,当外界不给予其修正时,误差增长无界。
微机电***MEMS(Micro Electronic Mechanical System)是随着半导体集成电路微细加工技术和超精密机械加工技术的发展而发展起来的集微型传感器、执行器、信号处理与控制电路、接口电路、通信和电源一体的微型机电***。MEMS惯性传感器主要指硅微加速度计和硅微陀螺仪这两种惯性传感器。MEMS惯导是基于MEMS惯性传感器的惯性导航***,导航***通过加速度计测量载体的三维的线性运动,陀螺测量载体的三维角运动,进而在导航坐标系上进行积分来获取速度、位置和姿态。MEMS惯导继承了传统惯导***的完全自主性、保密性强、不存在信号的电磁干扰等特点,同时又具有尺寸小、重量轻、成本低、功耗小、可靠性高等传统惯导无法比拟的优点。但受制造工艺所限,其用作导航定位时,误差随时间快速积累。
GNSS与MEMS惯导所形成的组合导航***充分发挥了二者的优势,形成有机互补。优势在于:MEMS惯导在短时间内能够保持较高的精度,且其不受外界工作环境影响;同时,GNSS提供的绝对定位测速信息可以补偿惯导随时间累积的误差,而保证长距离运行中的测量精度。
在常规的GNSS与MEMS惯导的组合***中,卫星导航需向组合***提供三维位置、三维速度以及对应的不确定(即误差),这样就对组合在一起的***提出了较高的结合要求,组合***要求GNSS提供完整的三维定位测速信息。然而,在卫星导航接收机的输出信息/报文中,目前各个厂家普遍遵循美国国家海洋电子协会(NMEA)所指定的行业标准。NMEA报文所定义的卫星导航信息并非为完整的三维信息,缺少直接的三维速度信息、位置误差、速度误差等内容,不能满足用于组合导航***或辅助MEMS惯导***的需求。因此,GNSS与MEMS惯性导航的组合***无法具有通用性,即所实现的组合***不能任意更换GNSS接收机。
发明内容
本发明的目的是提供一种使用不完整GNSS信息的MEMS组合导航建模方法,该方法仅使用标准通用的NMEA信息报文作为输入,并对NMEA输入的不完整GNSS信息进行数学建模,进而将此模型应用于组合导航***的导航算法中。
本发明所采用的技术方案是:使用不完整GNSS信息的MEMS组合导航建模方法,其特征在于,具体按照以下步骤实施,
第一步:确定实现组合导航***所需的外部GNSSNMEA输入报文类型,所述外部GNSSNMEA输入报文类型包括GGA、RMC、GSA、GSV及GLL,其接收顺序可向组合***发送配置命令;
第二步:从外部GNSSNMEA信息的GGA报文中直接提取三维位置信息,包括经度、维度及高度;
第三步:从外部GNSS NMEA信息的RMC报文中直接提取三维速度信息;
第四步:建立三维位置误差的数学模型;
第五步:线性回归算法确定三维位置误差的数学模型参数;
第六步:建立三维速度误差的数学模型;
第七步:线性回归算法确定速度误差模型参数;
第八步:将第四步与第六步所建立的三维位置与三维速度模型应用于组合导航算法。
本发明的特点还在于,
第三步中,三维速度信息的提取算法为:
VGNSS,N=VGNSScos(A)·cos(p)
VGNSS,E=VGNSS·sin(A)·cos(p)
VGNSS,D=-VGNSS·sin(p)
其中,VGNSS为GNSS在NMEARMC信息中提供的地速,其沿导航坐标系-北、东、地的三个分量分别为VGNSS,N,VGNSS,E,VGNSS,D;A为NMEA信息中提供的航向角;p为载体的俯仰角,该俯仰角由组合导航***的姿态决定。
第四步中,三维位置误差的数学模型为:
其中,σp,i为三维位置误差在i方向上的投影,i由导航坐标系定义,即北、东、地导航坐标系;j=1,2,…,M标识卫星星座***,包括GPS、北斗、格洛纳斯、伽利略;fj,hj分别表示对于第j个卫星***的位置数学模型非线性方程;E,A,VGNSS,PDOP分别为NMEA信息所提供的水平位置误差、航向、地速以及位置几何精度因子;Cnj,ej,Nj为对应的参与定位的第j个卫星***的所涉及的卫星载噪比,仰角以及卫星数。
第五步中,线性回归算法确定三维位置误差的数学模型参数的具体算法如下:
fj=α1,j·E+α2,jti(A,VGNSS)+a3,j·hj(Cnj,ej)+α4,j·PDOP+α5,j·VGNSS·Nj+a6,j
其中,ti(A,VGNSS)为地速在导航***沿第i方向的投影,即由第三步所计算的三维速度信息;α1,j…α8,j分别表示模型中的系数,系数由批量数据的线性回归分析确定。
第六步中,三维速度误差的数学模型为:
其中,σv,i为速度误差在i方向上的投影,i由导航坐标系定义,即北、东、地导航坐标系;j=1,2,…,M标识卫星星座***,包括GPS、北斗、格洛纳斯、伽利略;fj,hj分别表示对于第j个卫星***的位置数学模型非线性方程;gj,qj分别表示对于第j个卫星***的速度数学模型非线性方程;E,A,VGNSS,PDOP分别为NMEA信息所提供的水平位置误差、航向、地速以及位置几何精度因子;Cnj,ej,Nj为对应的参与定位的第j个卫星***的所涉及的卫星载噪比,仰角以及卫星数。
第七步中,线性回归算法确定速度误差的数学模型参数的具体算法如下:
gj=b1,j·E+b2,jti(A,VGNSS)+b3,j·qj(Cnj,ej)+b4,j·PDOP+b5,j·VGNSS·Nj+b6,j
其中,ti(A,VGNSS)为地速在导航***沿第i方向的投影,b1,j…b8,j分别表示模型中的系数,系数由批量数据的线性回归分析确定。
本发明的有益效果是:该方法中仅使用标准通用的NMEA信息报文作为输入,并对NMEA输入的不完整GNSS信息进行数学建模,进而将此模型应用于组合导航***的导航算法中。该实现方法突破了常规MEMS组合导航设备的实现限制,将组合导航***的实现难度大大降低。对于实现GNSS/MEMS组合***在广泛载体上的应用有着重要意义。
附图说明
图1是本发明使用不完整GNSS信息的MEMS组合导航建模方法的***架构图。
具体实施方式
下面结合附图和具体实施方式对本发明进行详细说明。
本发明提供使用不完整GNSS信息的MEMS组合导航建模方法,如图1所示,图1为卫星导航NMEA信息输入的MEMS组合导航***架构图,具体按照以下步骤实施,
第一步:确定实现组合导航***所需的外部GNSS NMEA输入报文类型,所述外部GNSS NMEA输入报文类型包括GGA、RMC、GSA、GSV及GLL,其接收顺序可向组合***发送配置命令;
第二步:从外部GNSS NMEA信息的GGA报文中直接提取三维位置信息,包括经度、维度及高度;
第三步:从外部GNSS NMEA信息的RMC报文中直接提取三维速度信息;
第三步中,三维速度信息的提取算法为:
VGNSS,N=VGNSS·cos(A)·cos(p)
VGNSS,E==VGNSS·sin(A)·cos(p)
VGNSS,D=-VGNSS·sin(p)
其中,VGNSS为GNSS在NMEA RMC信息中提供的地速,其沿导航坐标系-北、东、地的三个分量分别为VGNSS,N,VGNSS,E,VGNSS,D;A为NMEA信息中提供的航向角;p为载体的俯仰角,该俯仰角由组合导航***的姿态决定。
第四步:建立三维位置误差的数学模型;
第四步中,所述三维位置误差的数学模型为:
其中,σp,i为三维位置误差在i方向上的投影,i由导航坐标系定义,即北、东、地导航坐标系;j=1,2,…,M标识卫星星座***,包括GPS、北斗、格洛纳斯、伽利略;fj,hj分别表示对于第j个卫星***的位置数学模型非线性方程;E,A,VGNSS,PDOP分别为NMEA信息所提供的水平位置误差、航向、地速以及位置几何精度因子;Cnj,ej,Nj为对应的参与定位的第j个卫星***的所涉及的卫星载噪比,仰角以及卫星数。
第五步:线性回归算法确定三维位置误差的数学模型参数;
第五步中,线性回归算法确定三维位置误差的数学模型参数的具体算法如下:
fj=a1,j·E+α2,j·ti(A,VGNSS)+α3,j·hj(Cnj,ej)+α4,j·PDOP+α5,j·VGNSS·Nj+α6,j
其中,ti(A,VGNSS)为地速在导航***沿第i方向的投影,即由第三步所计算的三维速度信息;a1,j…a8,j分别表示模型中的系数,系数由批量数据(多于20组)的线性回归分析确定。
第六步:建立三维速度误差的数学模型;
第六步中,所述三维速度误差的数学模型为:
其中,σv,i为速度误差在i方向上的投影,i由导航坐标系定义,即北、东、地导航坐标系;j=1,2,…,M标识卫星星座***,包括GPS、北斗、格洛纳斯、伽利略;fj,hj分别表示对于第j个卫星***的位置数学模型非线性方程;gj,qj分别表示对于第j个卫星***的速度数学模型非线性方程;E,A,VGNSS,PDOP分别为NMEA信息所提供的水平位置误差、航向、地速以及位置几何精度因子;Cnj,ej,Nj为对应的参与定位的第j个卫星***的所涉及的卫星载噪比,仰角以及卫星数。
第七步:线性回归算法确定速度误差模型参数;
第七步中,线性回归算法确定速度误差的数学模型参数的具体算法如下:
gj=b1,j·E+b2,j·ti(A,VGNSS)+b3,j·qj(Cnj,ej)+b4,jPDOP+b5,j·VGNSS·Nj+b6,j
其中,ti(A,VGNSS)为地速在导航***沿第i方向的投影,b1,j…b8,j分别表示模型中的系数,系数由批量数据(多于20组)的线性回归分析确定。
第八步:将第四步与第六步所建立的三维位置与三维速度模型应用于组合导航算法。
本发明的优点是:该方法中仅使用标准通用的NMEA信息报文作为输入,并对NMEA输入的不完整GNSS信息进行数学建模,进而将此模型应用于组合导航***的导航算法中。该实现方法突破了常规MEMS组合导航设备的实现限制,将组合导航***的实现难度大大降低。对于实现GNSS/MEMS组合***在广泛载体上的应用有着重要意义。
Claims (4)
1.使用不完整GNSS信息的MEMS组合导航建模方法,其特征在于,具体按照以下步骤实施,
第一步:确定实现组合导航***所需的外部GNSSNMEA输入报文类型,所述外部GNSSNMEA输入报文类型包括GGA、RMC、GSA、GSV及GLL,其接收顺序可向组合***发送配置命令;
第二步:从外部GNSSNMEA信息的GGA报文中直接提取三维位置信息,包括经度、维度及高度;
第三步:从外部GNSS NMEA信息的RMC报文中直接提取三维速度信息;
第四步:建立三维位置误差的数学模型;
第五步:线性回归算法确定三维位置误差的数学模型参数;
第六步:建立三维速度误差的数学模型;
第七步:线性回归算法确定速度误差模型参数;
第八步:将第四步与第六步所建立的三维位置与三维速度模型应用于组合导航算法;
第三步中,所述三维速度信息的提取算法为:
第四步中,所述三维位置误差的数学模型为:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810205371.1A CN108267765B (zh) | 2018-03-13 | 2018-03-13 | 使用不完整gnss信息的mems组合导航建模方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810205371.1A CN108267765B (zh) | 2018-03-13 | 2018-03-13 | 使用不完整gnss信息的mems组合导航建模方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108267765A CN108267765A (zh) | 2018-07-10 |
CN108267765B true CN108267765B (zh) | 2021-08-03 |
Family
ID=62774798
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810205371.1A Active CN108267765B (zh) | 2018-03-13 | 2018-03-13 | 使用不完整gnss信息的mems组合导航建模方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108267765B (zh) |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101000245A (zh) * | 2007-01-10 | 2007-07-18 | 北京航空航天大学 | 一种sins/gps/磁罗盘组合导航***的数据融合方法 |
CN103048671A (zh) * | 2011-11-16 | 2013-04-17 | 郭杭 | 全球导航卫星***RTD、MEMS-IMU、ZigBee组合高精度无缝导航新方法 |
CN103235328A (zh) * | 2013-04-19 | 2013-08-07 | 黎湧 | 一种gnss与mems组合导航的方法 |
CN103235327A (zh) * | 2013-04-07 | 2013-08-07 | 清华大学 | 一种gnss/mins超深组合导航方法、***及装置 |
CN103454652A (zh) * | 2012-05-31 | 2013-12-18 | 长沙威佳电子科技有限公司 | 采用多或双gnss接收***的高精度gnss定位*** |
CN105510935A (zh) * | 2015-12-15 | 2016-04-20 | 北京光华纺织集团有限公司 | 一种用于导航的装置 |
CN105607093A (zh) * | 2015-12-20 | 2016-05-25 | 上海华测导航技术股份有限公司 | 一种组合导航***及获取导航坐标的方法 |
CN105607104A (zh) * | 2016-01-28 | 2016-05-25 | 成都佰纳瑞信息技术有限公司 | 一种基于gnss与ins的自适应导航定位***及方法 |
CN105741665A (zh) * | 2016-01-29 | 2016-07-06 | 北京华云智联科技有限公司 | 卫星导航教学演示平台 |
CN106767787A (zh) * | 2016-12-29 | 2017-05-31 | 北京时代民芯科技有限公司 | 一种紧耦合gnss/ins组合导航装置 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP5398120B2 (ja) * | 2007-03-22 | 2014-01-29 | 古野電気株式会社 | Gps複合航法装置 |
US9978285B2 (en) * | 2015-06-10 | 2018-05-22 | Ecole Polytechnique Federale De Lausanne (Epfl) | Autonomous and non-autonomous dynamic model based navigation system for unmanned vehicles |
-
2018
- 2018-03-13 CN CN201810205371.1A patent/CN108267765B/zh active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101000245A (zh) * | 2007-01-10 | 2007-07-18 | 北京航空航天大学 | 一种sins/gps/磁罗盘组合导航***的数据融合方法 |
CN103048671A (zh) * | 2011-11-16 | 2013-04-17 | 郭杭 | 全球导航卫星***RTD、MEMS-IMU、ZigBee组合高精度无缝导航新方法 |
CN103454652A (zh) * | 2012-05-31 | 2013-12-18 | 长沙威佳电子科技有限公司 | 采用多或双gnss接收***的高精度gnss定位*** |
CN103235327A (zh) * | 2013-04-07 | 2013-08-07 | 清华大学 | 一种gnss/mins超深组合导航方法、***及装置 |
CN103235328A (zh) * | 2013-04-19 | 2013-08-07 | 黎湧 | 一种gnss与mems组合导航的方法 |
CN105510935A (zh) * | 2015-12-15 | 2016-04-20 | 北京光华纺织集团有限公司 | 一种用于导航的装置 |
CN105607093A (zh) * | 2015-12-20 | 2016-05-25 | 上海华测导航技术股份有限公司 | 一种组合导航***及获取导航坐标的方法 |
CN105607104A (zh) * | 2016-01-28 | 2016-05-25 | 成都佰纳瑞信息技术有限公司 | 一种基于gnss与ins的自适应导航定位***及方法 |
CN105741665A (zh) * | 2016-01-29 | 2016-07-06 | 北京华云智联科技有限公司 | 卫星导航教学演示平台 |
CN106767787A (zh) * | 2016-12-29 | 2017-05-31 | 北京时代民芯科技有限公司 | 一种紧耦合gnss/ins组合导航装置 |
Non-Patent Citations (3)
Title |
---|
INS/GPS组合导航下机载PD雷达DPCA性能分析;杨勇等;《雷达科学与技术》;20110430;第9卷(第2期);第166-171,179页 * |
典型弹道下的火箭弹MEMS-INS/GNSS组合导航姿态误差可观性分析;董进龙等;《兵工学报》;20140630;第35卷(第6期);第850-856页 * |
陆地车辆GNSS/MEMS惯性组合导航机体系约束算法研究;刘华等;《北京理工大学学报》;20130531;第33卷(第5期);第510-515页 * |
Also Published As
Publication number | Publication date |
---|---|
CN108267765A (zh) | 2018-07-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Cao et al. | GVINS: Tightly coupled GNSS–visual–inertial fusion for smooth and consistent state estimation | |
CN1322311C (zh) | 车载快速定位定向*** | |
De Agostino et al. | Performances comparison of different MEMS-based IMUs | |
CN108594283B (zh) | Gnss/mems惯性组合导航***的自由安装方法 | |
CN106767787A (zh) | 一种紧耦合gnss/ins组合导航装置 | |
CN112505737B (zh) | 一种gnss/ins组合导航方法 | |
CN103994763A (zh) | 一种火星车的sins/cns深组合导航***及其实现方法 | |
Li et al. | P $^{3} $-VINS: Tightly-coupled PPP/INS/visual SLAM based on optimization approach | |
CN111637888B (zh) | 基于惯导和激光雷达单点测距的掘进机定位方法及*** | |
Sun et al. | An adaptive weighting strategy for multisensor integrated navigation in urban areas | |
CN115327588A (zh) | 一种基于网络rtk的无人自动化作业特种车高精度定位方法 | |
CN115435779A (zh) | 一种基于gnss/imu/光流信息融合的智能***姿估计方法 | |
CN115388884A (zh) | 一种智能***姿估计器联合初始化方法 | |
Wen et al. | 3D LiDAR aided GNSS real-time kinematic positioning | |
CN105466423A (zh) | 一种无人机导航***及其运行方法 | |
CN204255368U (zh) | 一种适用于火星车的sins/cns深组合导航*** | |
CN116576849A (zh) | 一种基于gmm辅助的车辆融合定位方法及*** | |
CN108267765B (zh) | 使用不完整gnss信息的mems组合导航建模方法 | |
CN116105729A (zh) | 一种针对野外洞穴森林环境侦察的多传感器融合定位方法 | |
CN115523920A (zh) | 一种基于视觉惯性gnss紧耦合的无缝定位方法 | |
WO2002046699A1 (en) | Vehicle positioning and data integrating method and system thereof | |
CN112649001B (zh) | 一种小型无人机姿态与位置解算方法 | |
CN114646993A (zh) | 一种基于gnss、视觉以及imu进行精确定位的数据融合算法 | |
CN115047506A (zh) | 一种基于gnss-ins和视觉的农机组合导航*** | |
CN110764126B (zh) | 一种gps信息缺失情况下的无人车导航方法及*** |
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 |