CN105527642B - 一种单星定位装置与方法 - Google Patents

一种单星定位装置与方法 Download PDF

Info

Publication number
CN105527642B
CN105527642B CN201610044470.7A CN201610044470A CN105527642B CN 105527642 B CN105527642 B CN 105527642B CN 201610044470 A CN201610044470 A CN 201610044470A CN 105527642 B CN105527642 B CN 105527642B
Authority
CN
China
Prior art keywords
user terminal
mrow
satellite
msup
positioning
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
Application number
CN201610044470.7A
Other languages
English (en)
Other versions
CN105527642A (zh
Inventor
施浒立
刘成
程涛
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing Riyuejiutian Technology Co Ltd
Original Assignee
Beijing Riyuejiutian Technology Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Beijing Riyuejiutian Technology Co Ltd filed Critical Beijing Riyuejiutian Technology Co Ltd
Priority to CN201610044470.7A priority Critical patent/CN105527642B/zh
Publication of CN105527642A publication Critical patent/CN105527642A/zh
Application granted granted Critical
Publication of CN105527642B publication Critical patent/CN105527642B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/38Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
    • G01S19/39Determining 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/42Determining position
    • G01S19/45Determining position by combining measurements of signals from the satellite radio beacon positioning system with a supplementary measurement
    • G01S19/47Determining 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

本发明涉及一种单星定位装置与方法。所述装置包括用户终端和一颗具有测距功能的卫星或单基站或单信号源,用户终端安装了气压测高芯片、测角部件、计步器、囚禁原子钟CPT。所述方法利用气压测高芯片与一颗具有测距功能的卫星或单基站或单信号源分别获得高程值与距离测量值,通过两个曲面交会确定用户所在位置的环。然后基于航位外推原理,利用测角部件获得的方位角信息不断外推、修正,连续获得用户终端的三维位置坐标。本发明的定位原理不同于传统GPS类卫星***,摆脱了四星定位的限制,既能用作卫星数目较少时传统GPS类卫星定位***的补充,亦可作为一种独立的定位方法用于手机定位、基于位置服务、海洋浮标和轮船的定位与监测领域。

Description

一种单星定位装置与方法
技术领域
本发明属于导航定位领域,尤其涉及一种包括微型传感测量器件的单星定位导航装置与方法。
背景技术
以GPS(Global Positioning System,全球定位***)为代表的各种GNSS(GlobalNavigation Satellite System,全球卫星导航***)作为现代科学与技术的结晶,为人们的生活与工作带来了巨大的便利和变革。然而,GPS类卫星定位***也并非尽善尽美,其本身仍存在着一些明显的不足。从测量学原理讲,GNSS采用基于伪距测量的距离交会原理,需要观测4颗以上的卫星才能顺利进行三维定位和导航,所以GNSS依赖于完整和可靠的卫星星座。这使得GNSS的建设与维护需要花费大量的人力与物力,建设周期较长,***复杂度较高。而实际上,***本身又比较脆弱,信号容易受到干扰。在诸如城市、高楼和山区等遮挡物较多的环境下,观测卫星数目容易少于4颗,此时传统的导航定位模型将无法进行实时定位。若在战时等特殊时期,其可用性更会大大降低。一份来自美国海军研究GPS卫星***未来发展的研究报告曾明确指出:如果不对GPS***进行大规模的建设与改造,那么由于其脆弱性与易干扰性,卫星星座与信号将很容易受到干扰或破坏,因此也将无法适用于未来战争。
发明内容
为了解决现有技术中存在的上述问题,本发明提出一种包括微传感器件的单星或单基站或单信号源定位装置与方法,在仅有一颗测距卫星或单基站时也能实现具有一定精度水平的连续三维定位。
为达到上述目的,本发明采用如下技术方案:
一种单星定位装置,包括用户终端和一颗具有测距功能的卫星或单基站或单信号源,用户终端安装了气压测高芯片、测角部件、计步器、囚禁原子钟CPT(CoherentPopulation Trapping)。其中,
所述单卫星、单基站或单信号源为导航定位提供基准点,测量与用户终端的伪距值。
所述气压测高芯片用以测量用户终端所处位置的气压值,将所述气压值转换后得到大地高程值。
所述测角部件利用惯性或磁场测量用户终端运动的方位角度及方位角度的变化量。
所述囚禁原子钟CPT用于提供时间和频率值。
所述计步器用于测量用户终端行进的步距值。
进一步地,所述测角部件为一指北计。
进一步地,所述测角部件包括陀螺仪和加速度计。
进一步地,所述囚禁原子钟CPT提供的时间和频率值的时频稳定度优于10-10
一种应用所述单星定位装置进行定位的方法,包括以下步骤:
步骤1,用户终端接收来自具有所述测距功能的卫星或单基站或单信号源的导航信号,解调后得到导航电文和一个伪距测量值,消除伪距测量值中包括接收终端钟差的误差,得到用户终端与所述卫星或单基站之间的伪距测量值ρ1。用户终端位于以所述卫星或基站或信号源为球心、以ρ1为半径的球面上。
步骤2,由所述气压测高芯片测得用户终端所处位置的大地高程值h,得到高程值为h的等高程约束曲面,用户终端位于所述等高程约束曲面上。
步骤3,步骤1所述的球面与步骤2所述的等高程约束曲面交会得到一个圆S,用户终端位于所述圆S上。
步骤4,由测角部件测得用户终端运动的方位角及变化量,所述方位角为用户终端运动方向与正北的夹角。若上一观测历元用户终端位于U0(x0,y0,z0)处,上一观测历元测得的方位角为α1,过点U0、方向与用户终端运动方向相同的直线与步骤3所述圆S的交点即为当前观测历元的用户终端位置U1(x,y,z)。
步骤5,由所述计步器获得步距,采用广义融合方法修正终端位置U1(x,y,z)。
进一步地,步骤1所述的球面方程为:
其中,(x1,y1,z1)为卫星或基站的坐标。
进一步地,步骤2所述的等高程约束曲面的方程为:
其中,a、b分别为地球椭圆的长轴和短轴的长度。
进一步地,步骤4所述的直线的方程为:
y-y0=tanα1×(x-x0)
进一步地,用户终端处于海洋环境时,所述等高程约束曲面为海平面。
与现有技术相比,本发明具有以下有益效果:
本发明所述的定位装置采用气压测高芯片、陀螺仪、计步器等传感器器件,实现了一类基于单颗卫星测距的导航定位方法。本发明所述方法的定位原理不同于传统的GPS类卫星***,摆脱了四星定位的限制,方法简单易行,既能用作卫星数目较少时传统GPS类卫星定位***的补充和辅助,亦可作为一种独立的定位方法,能够广泛适用于手机定位、LBS(Location Based Service,基于位置的服务)、海洋浮标和轮船的定位与监测等领域。
附图说明
图1为曲面交会示意图;
图2为本发明单星定位示意图;
图3为定位累积误差消除示意图。
具体实施方式
下面结合附图和实施例对本发明做进一步说明。
近年来,随着导航定位、电子信息、LBS产业等技术领域的迅速发展,出现了许多崭新的定位设备、技术与方法。例如,目前廉价且微型的气压测高芯片、微陀螺仪芯片、计步器等传感器,已经可以测量得到较高精度的高程值、角度和步距等信息,并已被集成到手机等定位终端之中,参与导航与定位功能。CPT(Coherent Population Trapping,囚禁原子钟或高精度晶振)已经可以达到优于10-10的时频稳定度,且价格可以控制在千元量级。因此,在某些定位终端设备中已经可以配备CPT原子钟,以便直接获得接收机时钟与卫星(或基站或单信号源)时钟之间的差值。
本发明所述的单星定位装置,包括用户终端和一颗具有测距功能的卫星或单基站或单信号源(任意一个可利用的信号辐射源)。用户终端至少包括气压测高芯片、测角部件、计步器、囚禁原子钟CPT。其中,
气压测高芯片用以测量用户终端所处位置的气压值。气压值与大地高程相关联,将测量得到的气压值转换成大地高程值,解决大地高程的度量问题。
测角部件利用惯性原理或磁场测量得到方位角度及方位角度的变化量,解决行进时的方向问题。测角部件可以采用指北计,也可以由陀螺仪和加速度计完成测角。
计步器利用加速度计获得用户终端的速度和步距。根据速度和步距信息,采用广义融合方法修正终端位置U1(x,y,z)。
用户终端可以是一部普通的手机,因此气压测高芯片、测角部件和测速部件均采用微型芯片,它们的尺寸小到足以使它们能够安装到手机的线路板上。
囚禁原子钟CPT用以提供高稳定性的时间和频率值,能够提供稳定度优于10-10的时频信号。高稳定度的时频信号能够为用户终端提供高精度的伪距测量值。也可以采用高精度晶振代替囚禁原子钟CPT。
单卫星或单基站或单信号源为导航定位提供基准点,同时还用以测量伪距值,为导航定位提供最基本的条件。
本发明提出一种不同于传统GPS类卫星***原理的定位方法。利用气压测高芯片与一颗具有测距功能的卫星或单基站或单信号源,分别获得高程值与距离测量值,通过两个曲面交会的方法确定用户所在位置的环。然后,基于航位外推原理,利用指北计或陀螺仪及加速度计所获得的方位角信息,不断外推、修正,从而可以连续获得用户终端的三维位置坐标。所述方法包括以下步骤:
步骤1,用户终端接收来自一颗导航或通信卫星或单基站或单信号源的导航信号,解调后得到导航电文和一个伪距测量值。
由于定位终端中的CPT囚禁原子钟(或高精度晶振)具有较高的精度和稳定度,因此可以通过直接比对得到伪距测量值中的接收机的钟差值,从而获得终端与该卫星之间的高精度伪ρ1。以卫星或基站或单信号源为圆心、该终端与该卫星或基站或单信号源之间的距离ρ1为半径可以作球,用户位置即位于球面上。该球面的方程为:
其中,(x1,y1,z1)为卫星或基站或单信号源的坐标。
还可以采用以下两种方法获得高精度伪距测量值:(1)采用钟差模型外推加上卫星定期比对授时修正的方法;(2)利用已知点或者差分基准点/站的伪距测量差分信息进行修正。
步骤2,利用安装在定位终端上的气压测高芯片,测量得到用户终端的高程值h,因此为用户终端提供了一个等高程的约束曲面,用户终端位于该约束曲面上。该曲面的方程为:
其中,a、b分别为地球椭圆的长轴和短轴的长度。
联立以上两个曲面方程,可以在地面上相交得到一个圆S。用户位于该圆S上。
当本发明应用于海洋浮标或舰船的导航定位领域时,由于存在着海水面这个天然的高程约束面,可以不采用气压测高芯片进行测高从而得到等高程约束曲面。
步骤3,由测角部件测量用户终端运动的方位角度及方位角度的变化量。所述方位角为用户终端运动方向与正北的夹角。根据上一观测历元用户终端的位置和运动方向确定当前观测历元用户终端位置。具体方法如下:
若上一观测历元的用户位置在U0(x0,y0,z0)处,而在上一观测历元测得的方位角为α1,如图2所示,过点U0、方向与运动方向相同的直线与所述圆S的交点,即为当前观测历元的用户终端位置U1(x,y,z)。所述直线方程为:
y-y0=tanα1×(x-x0)
因此,单星定位方法的导航量测方程组模型可表述为:
通过本方法得到的U1的坐标在其法线方向n1上的误差较小,误差主要分布在点U1与圆S1切线方向l1附近,如附图2所示。由于用户终端与卫星的位置都在不停地变化与改变,因此在该法线方向上的定位误差一般不会持续累积。可以根据用户终端与卫星之间的运动趋势与定位点运动状态的改变对其进行修正。
如附图3所示,在某一观测历元i得到约束圆Si与用户坐标位置Ui,Ui点切线li的方位角为αi,则此时定位点误差将沿li方向上较大。而如果在另一观测历元j得到约束圆Sj与用户坐标位置Uj,Uj点法线nj的方位角为αj,则此时定位点误差将沿nj方向上较小。因此,当αi与αj相等或者近似相等时,即可有效消除该方向的累积误差。由于用户在运动过程中会不断改变其前进方向,一直沿着准li方向漂移的概率较小,所以误差会得到修正和消除。也可以利用计步器测量步距,进行融合修正,从而使定位精度得到改善。
本发明不限于上述实施方式,本领域技术人员所做出的对上述实施方式任何显而易见的改进或变更,都不会超出本发明的构思和所附权利要求的保护范围。

Claims (8)

1.一种单星定位装置,包括用户终端和一颗具有测距功能的卫星或单基站或单信号源,其特征在于,所述用户终端安装了气压测高芯片、测角部件、囚禁原子钟CPT、计步器;其中,
所述单卫星或单基站或单信号源为导航定位提供基准点,测量与用户终端的伪距值;
所述气压测高芯片用于测量用户终端所处位置的气压值,将所述气压值转换后得到大地高程值;
所述测角部件用于测量用户终端运动的方位角度及方位角度的变化量;
所述囚禁原子钟CPT用于提供时间和频率值;
所述计步器用于测量用户终端行进的步距值;所述囚禁原子钟CPT提供的时间和频率值的时频稳定度优于10-10
2.根据权利要求1所述单星定位装置,其特征在于,所述测角部件包括陀螺仪和加速度计。
3.根据权利要求1所述单星定位装置,其特征在于,所述测角部件为一指北计。
4.一种应用权利要求1所述单星定位装置进行定位的方法,包括以下步骤:
步骤1,用户终端接收来自具有所述测距功能的卫星或单基站或单信号源的导航信号,解调后得到导航电文和一个伪距测量值,消除伪距测量值中包括接收终端钟差的误差,得到用户终端与所述卫星或单基站之间的伪距测量值ρ1;用户终端位于以卫星或单基站为球心、以ρ1为半径的球面上;
步骤2,由所述气压测高芯片测得用户终端所处位置的大地高程值h,得到高程值为h的等高程约束曲面,用户终端位于所述等高程约束曲面上;
步骤3,步骤1所述的球面与步骤2所述的等高程约束曲面交会得到一个圆S,用户终端位于所述圆S上;
步骤4,由测角部件测得用户终端运动的方位角及变化量,所述方位角为用户终端运动方向与正北的夹角;若上一观测历元用户终端位于U0(x0,y0,z0)处,上一观测历元测得的方位角为α1,过点U0、方向与用户终端运动方向相同的直线与步骤3所述圆S的交点即为当前观测历元的用户终端位置U1(x,y,z);
步骤5,由所述计步器获得步距,采用广义融合方法修正终端位置U1(x,y,z)。
5.根据权利要求4所述的定位方法,其特征在于,步骤1所述的球面方程为:
<mrow> <msqrt> <mrow> <msup> <mrow> <mo>(</mo> <mi>x</mi> <mo>-</mo> <msub> <mi>x</mi> <mn>1</mn> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msup> <mrow> <mo>(</mo> <mi>y</mi> <mo>-</mo> <msub> <mi>y</mi> <mn>1</mn> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msup> <mrow> <mo>(</mo> <mi>z</mi> <mo>-</mo> <msub> <mi>z</mi> <mn>1</mn> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> </msqrt> <mo>=</mo> <msub> <mi>&amp;rho;</mi> <mn>1</mn> </msub> </mrow>
其中,(x1,y1,z1)为卫星或单基站的坐标。
6.根据权利要求4所述的定位方法,其特征在于,步骤2所述的等高程约束曲面的方程为:
<mrow> <mfrac> <mrow> <msup> <mi>x</mi> <mn>2</mn> </msup> <mo>+</mo> <msup> <mi>y</mi> <mn>2</mn> </msup> </mrow> <msup> <mrow> <mo>(</mo> <mi>a</mi> <mo>+</mo> <mi>h</mi> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mfrac> <mo>+</mo> <mfrac> <msup> <mi>z</mi> <mn>2</mn> </msup> <msup> <mrow> <mo>(</mo> <mi>b</mi> <mo>+</mo> <mi>h</mi> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mfrac> <mo>=</mo> <mn>1</mn> </mrow>
其中,a、b分别为地球椭圆的长轴和短轴的长度。
7.根据权利要求4所述的定位方法,其特征在于,步骤4所述的直线的方程为:
y-y0=tanα1×(x-x0)。
8.根据权利要求4所述的定位方法,其特征在于,所述用户终端处于海洋环境时,所述等高程约束曲面为海平面。
CN201610044470.7A 2016-01-22 2016-01-22 一种单星定位装置与方法 Active CN105527642B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610044470.7A CN105527642B (zh) 2016-01-22 2016-01-22 一种单星定位装置与方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610044470.7A CN105527642B (zh) 2016-01-22 2016-01-22 一种单星定位装置与方法

Publications (2)

Publication Number Publication Date
CN105527642A CN105527642A (zh) 2016-04-27
CN105527642B true CN105527642B (zh) 2018-03-30

Family

ID=55769978

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610044470.7A Active CN105527642B (zh) 2016-01-22 2016-01-22 一种单星定位装置与方法

Country Status (1)

Country Link
CN (1) CN105527642B (zh)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106054227B (zh) * 2016-07-04 2018-05-08 西北工业大学 惯导辅助下的伪距差值单星高动态定位方法
CN106019346B (zh) * 2016-07-04 2018-04-20 西北工业大学 基于双向通信协同定位辅助的单星定位方法
CN106842256B (zh) * 2016-12-12 2019-03-12 西安空间无线电技术研究所 一种利用gnss单星信号的导航定位方法
CN109429333B (zh) * 2017-08-29 2021-04-20 西安中兴新软件有限责任公司 一种定位方法及基站
CN110109163B (zh) * 2019-04-09 2022-01-07 天津大学 一种附有高程约束的精密单点定位方法
CN110109167B (zh) * 2019-04-18 2021-12-07 天津大学 一种基于高程约束的海上精密定位方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102445200A (zh) * 2011-09-30 2012-05-09 南京理工大学 微小型个人组合导航***及其导航定位方法
CN103383462A (zh) * 2012-05-03 2013-11-06 中国科学院国家天文台 一种卫星导航通信终端高程定位的方法
CN103968836A (zh) * 2014-05-16 2014-08-06 施浒立 一种基于时序伪距差分推算移动目标位置的方法及装置
CN104374368A (zh) * 2014-06-25 2015-02-25 中国科学院国家天文台 一种基于差分气压测高计及其测量方法
CN104459749A (zh) * 2014-12-02 2015-03-25 范子成 一种能对基站定位进行纠正的定位方法和装置
CN104931993A (zh) * 2015-05-25 2015-09-23 清华大学 微型定位导航授时***

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5901171A (en) * 1996-03-15 1999-05-04 Sirf Technology, Inc. Triple multiplexing spread spectrum receiver

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102445200A (zh) * 2011-09-30 2012-05-09 南京理工大学 微小型个人组合导航***及其导航定位方法
CN103383462A (zh) * 2012-05-03 2013-11-06 中国科学院国家天文台 一种卫星导航通信终端高程定位的方法
CN103968836A (zh) * 2014-05-16 2014-08-06 施浒立 一种基于时序伪距差分推算移动目标位置的方法及装置
CN104374368A (zh) * 2014-06-25 2015-02-25 中国科学院国家天文台 一种基于差分气压测高计及其测量方法
CN104459749A (zh) * 2014-12-02 2015-03-25 范子成 一种能对基站定位进行纠正的定位方法和装置
CN104931993A (zh) * 2015-05-25 2015-09-23 清华大学 微型定位导航授时***

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
一种气压测高辅助下的地面移动网定位方法;刘成等;《桂林理工大学学报》;20120229;第32卷(第1期);全文 *
星座通信***中单星定位研究;张计全等;《遥测遥控》;20090131;第30卷(第1期);全文 *
联合气压测高定位的高程转换问题;刘成等;《测绘科学》;20150630;第40卷(第6期);第6页左栏第2段 *

Also Published As

Publication number Publication date
CN105527642A (zh) 2016-04-27

Similar Documents

Publication Publication Date Title
CN105527642B (zh) 一种单星定位装置与方法
US10309786B2 (en) Navigational and location determination system
CN102176041B (zh) 一种基于gnss/sins组合的车辆导航监控***
CN103115623B (zh) 基于偏振光仿生导航的定位***的定位方法
CN104457754A (zh) 一种基于sins/lbl紧组合的auv水下导航定位方法
CN106908086A (zh) 一种多普勒计程仪测速误差的修正方法
CN113311436A (zh) 一种移动平台上激光测风雷达运动姿态测风订正方法
CN110174104A (zh) 一种组合导航方法、装置、电子设备及可读存储介质
CN110133700B (zh) 一种船载综合导航定位方法
CN106092094A (zh) 基于gnss/ins/激光测距组合定位的手持接收机及测量方法
CN103900580A (zh) 基于gis技术的北斗/gps与ins组合车载导航定位***
CN103389092A (zh) 一种系留飞艇姿态测量装置及测量方法
CN106093992A (zh) 一种基于cors的亚米级组合定位导航***及导航方法
CN112730883B (zh) 测量河流水位高度和坡度的gnss漂流浮毯及方法
CN103575297B (zh) 基于卫星导航接收机的gnss和mimu组合导航航向角估计方法
Wright et al. The effectiveness of global positioning system electronic navigation
JP2005241441A (ja) 移動式水上風観測システム
CN203364836U (zh) 一种系留飞艇姿态测量装置
KR19990017468A (ko) 상대 위치 측정장치
KR20050025037A (ko) 전파나침반
CN110646817A (zh) 计算定位误差的方法和高精度定位方法
Xiaoyu et al. GH-LPS/INS integration for precise UAV application
CN110988796B (zh) 一种基于频偏机会信号的定位方法及***
Moore et al. RiGHt: River level monitoring using GPS heighting
CN113933879B (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