CN102981174A - 一种gps天线相位中心变化对相对定位精度的修正方法 - Google Patents
一种gps天线相位中心变化对相对定位精度的修正方法 Download PDFInfo
- Publication number
- CN102981174A CN102981174A CN2012103342961A CN201210334296A CN102981174A CN 102981174 A CN102981174 A CN 102981174A CN 2012103342961 A CN2012103342961 A CN 2012103342961A CN 201210334296 A CN201210334296 A CN 201210334296A CN 102981174 A CN102981174 A CN 102981174A
- Authority
- CN
- China
- Prior art keywords
- centerdot
- delta
- phase center
- antenna phase
- gps
- 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
Images
Landscapes
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明公开了一种GPS天线相位中心变化对相对定位精度的修正方法,包括以下步骤:两个目标均载有GPS接收机及GPS天线,利用GPS卫星信号对两个目标进行相对定位;确定两个目标的GPS天线相位中心变化;根据所述天线相位中心变化确定与目标相对定位精度修正量的关系;利用所述修正量对目标相对定位精度进行修正。采用本发明实现了根据天线相位中心变化修正目标相对定位结果,提高了目标相对定位精度。
Description
技术领域
本发明涉及一种GPS相对定位精度修正方法,特别涉及一种利用GPS天线相位中心变化对相对定位精度的修正方法。
背景技术
GPS(Global Position System,全球定位系统)天线与GPS接收机共同完成对全球导航卫星(GPS卫星)信号的接收、处理及解算工作。GPS接收机通过天线接收GPS导航卫星信号及相关数据信息。由于GPS天线本身的特性,天线几何中心与相位中心不可能完全一致,两者在水平方向和高程方向均存在偏差,并且随着接收到的GPS卫星的方位和高度的变化而变化。天线相位中心的误差,必将对测量结果带来影响,尤其是在高精度的GPS测量中影响不可忽视。
天线相位中心是指天线的辐射场分量等效的由这点辐射出去,辐射波为一球面。理想天线存在唯一的相位中心,其等相面为球面,因此接收来自不同方向的卫星信号时,天线因自身产生额外的相位差而造成定位测量结果的偏差。然而在整个波束空间存在唯一相位中心的天线几乎是不存在的,只可能在主瓣某一范围内相位保持相对恒定,因此接收天线在接收不同方向的卫星信号时会引入额外的相位差异,而引入定位测量结果的误差。显然,天线类型不同,其相位特性不同,对接收机测量精度的影响也不同。天线相位中心与几何中心以及等相面关系示意图见图2,ARP为天线几何中心(天线参考点)一般定义为天线上表面的中心点,E为天线的平均相位中心。可以看出,正常情况下相位中心存在两个偏差。天线相位中心E与天线参考点ARP的偏差称为天线相位中心偏移PCO(Phase Center Offset),实际相位中心P点与平均相位中心E点的偏差称为天线相位中心变化PCV(Phase Center Variations),其数值等效于实际等相位面与拟合等相位球面的偏差。
发明内容
本发明的技术解决问题是:提供了一种利用GPS天线相位中心变化对相对定位精度的修正方法,实现了根据天线相位中心变化修正目标相对定位结果,提高了目标相对定位精度。
本发明的技术解决方案是:
一种GPS天线相位中心变化对相对定位精度的修正方法,包括以下步骤:
1)两个目标均载有GPS接收机及GPS天线,利用GPS卫星信号对两个目标进行相对定位;
2)确定2个目标的GPS天线相位中心变化;
3)根据GPS天线相位中心变化确定目标相对定位精度修正量:
单差方式进行相对定位时,根据所述天线相位中心变化通过下式确定与目标相对定位精度修正量的关系:
其中,σδX为用于表示目标相对定位精度修正量的标准差;PDOP为空间位置精度因子;σ为用于表示GPS天线相位中心变化在均值为0时的测量误差的标准差;σpcv为用于表示GPS天线相位中心变化在均值为0时的标准差。
双差方式进行相对定位时,根据所述天线相位中心变化通过下式确定与目标相对定位精度修正量的关系:
其中,σδX为用于表示目标相对定位精度修正量的标准差;DOP为精度因子;σ为用于表示GPS天线相位中心变化在均值为0时的测量误差的标准差;σpcv为用于表示GPS天线相位中心变化在均值为0时的标准差。
4)利用目标相对定位精度修正量对目标相对定位精度进行修正。
单差方式进行相对定位时,所述GPS天线相位中心变化与所述目标相对定位修正量的关系通过下式在GPS天线到各卫星产生的天线相位变化呈正态分布,且不同卫星到接收机之间的相位中心变化互不相关时获得:
其中,δXP为只考虑天线相位中心变化影响的δX2,
δt12=δt2-δt1;
δt1、δt2分别为两个接收机的钟差,c为光速;
A为δX2的系数矩阵,
其中[xp(t1),yp(t1),zp(t1)]为t1时刻卫星p的位置坐标,p=1,2,…n。D12为只含有天线相位中心变化信息,
双差方式进行相对定位时,所述GPS天线相位中心变化与所述目标相对定位修正量的关系通过下式在GPS天线到各GPS卫星产生的天线相位变化呈正态分布,且不同卫星到GPS接收机之间的相位中心变化互不相关时获得:
δXP=-(ATA)-1AT·ΔD
其中,δXP为只考虑天线相位中心变化影响的δX2,
A为δX2的系数矩阵,
其中[xp(t1),yp(t1),zp(t1)]为t1时刻卫星p的位置坐标,p=1,2,…n。
ΔD为只含有天线相位中心变化信息,
本发明与现有技术相比的优点在于:
1)采用本发明获得了利用GPS天线相位中心变化对GPS相对定位精度的影响,实现了对目标相对定位精度的修正。
2)本发明给出了星载GPS天线相位中心变化与GPS相对定位精度之间的定量关系,给出了精确的数学模型。现有的技术虽然分析到GPS天线相位中心变化的存在,对GPS测量的水平位置和高程位置精度都会产生影响。但从来没有提出过明确的定量关系和数学模型。
附图说明
图1为本发明流程图。
图2为本发明天线相位中心与几何中心以及等相面关系示意图
图3为本发明卫星信号的方位角和仰角示意图
图4为本发明天线相位中心变化导致的径向距离变化
图5为GPS单差相对定位示意图
图6为GPS双差相对定位示意图
图7为本发明天线相位中心变化对单差相对定位影响总误差实例
图8为本发明天线相位中心变化对单差相对定位各个方向上误差实例
图9为本发明天线相位中心变化对双差相对定位影响总误差实例
图10为本发明天线相位中心变化对双差相对定位各个方向上误差实例
图11为不同PDOP值下相对定位精度统计曲线(单差方式)
图12为不同DOP值下相对定位精度统计曲线(双差方式)
具体实施方式
下面就结合附图对本发明做进一步介绍。
由于星载GPS天线几何中心与相位中心不一致,随着接收到的GPS卫星的高度和方位的变化而变化,天线相位中心的误差直接对测量结果(目标相对定位精度)带来影响,尤其是在高精度的GPS测量中需要分析并消除部分影响。
天线相位中心偏移PCO是一个固定不变的矢量,而天线相位中心变化PCV则随着卫星入射信号高度角和方位角的变化而变化,因此会给目标相对定位精度引入额外的相位偏差。如图1所示为本发明流程图,下面就针对流程图对本发明进行介绍。
天线平均相位中心E点,PCO、PCV、天线参考点APR的关系如图2所示。对于每个天线,向量表示天线相位中心偏移PCO,其长度为一个固定的值。向量表示卫星信号方向上的单位向量,可以用卫星信号相对于大地局部坐标系的方位角θ、仰角表示,如图3所示。在图3中,S为卫星位置,O为接收机位置,向量ON为大地正北方向,OE为大地正东方向,OU为大地天顶方向。向量OS在大地平面EON上的投影为向量OA,那么方位角θ为向量ON到OA的角度,范围为0~360度;仰角为向量OA到向量OS的角度,范围为0~90度。
GPS相对定位的基本原理是:用2台接收机分别安置在基线的2个端点,同步观测相同的4颗以上GPS卫星,确定基线2个端点在地球坐标系中的相对位置。
在2个观测站或多个观测站同步观测相同卫星的情况下,卫星的轨道误差、卫星钟差、接收机钟差以及电离层和对流层的折射误差等对观测值的影响具有一定的相关性。利用这些观测量的不同组合(求差)进行相对定位,可有效削弱相关误差的影响,从而提高相对定位的精度。对GPS载波相位观测值的求差,可以分别在卫星之间、接收机之间和不同历元之间进行,根据求差次数的不同又可分为单差,双差等方式。实际工程应用中,可以采用单差或者双差方式。
单差(Single Different,SD)是指不同观测站同步观测相同卫星所得观测量之差。
在图5中,若在t1时刻在观测站1、观测站2同时对卫星p进行了载波相位测量,得观测方程:
公式(3)和公式(4)中,表示t1时刻接收机i到卫星j的相位观测量,表示t1时刻接收机i到卫星j的几何距离,分别表示t1时刻接收机i到卫星j的电离层延时距离误差和对流层延时距离误差,δti表示接收机i的钟差,δtj表示卫星j的钟差,表示观测初始时刻t0接收机i到卫星j的整周未知数,f为载波频率,c为光速。
将公式(3)和公式(4)代入公式(2)中,可处理得:
其中,
其中[xp(t1),yp(t1),zp(t1)]为t1时刻卫星p的位置坐标,p=1,2,…n。
表示观测初始时刻t0接收机i到卫星p的整周未知数,i=1,2。
对于观测站1和观测站2,对其相位观测值进行修正,并经过一系列的处理,可得
L=Ln-D12 (6)
则有,
其中,只有D12矩阵包含天线相位中心变化因素,因此,天线相位中心变化引起的相对定位单差观测方式的误差为:
双差(Double Different,DD)是指不同观测站同步观测同一组卫星,所得单差观测值之差。
如图6所示,观测站1、观测站2在t1时刻同时观测了2颗卫星p、q,那么对p、q卫星分别有如公式(2)所示的单差方程。若忽略大气折射残差,双差观测方程为
经过一系列处理,可得
δX2=λ(ATA)-1AT(ΔL-CΔN) (10)
对于观测站1和观测站2,对其相位观测值进行修正,并经过一系列的处理,可得
λΔL=λΔLn-ΔD (11)
参考双差观测的解算公式(9),得到
其中,只有ΔD矩阵包含天线相位中心变化因素,因此,天线相位中心变化引起的双差观测方式的相对定位精度为
δXP=-(ATA)-1AT·ΔD (13)
利用载波相位测量的星载GPS天线相位中心变化与GPS接收机相对定位精度的数学关系公式(8)和(13),分别建立星载GPS天线相位中心变化对相对定位精度影响的单差、双差数学模型,并设计仿真软件,定量的分析天线相位中心变化对相对定位精度的影响。从而实现接收机天线相位中心变化对相对定位精度的影响数学模型的验证。该软件要求输入天线相位中心变化数据(自定义格式或者标准格式),卫星位置数据(自定义格式或者标准格式),以及模型验证数据。用于数学模型的仿真验证。
GPS天线相位中心变化引起的单差相对定位精度公式(8)中,假设GPS接收机天线到各个卫星产生的相位中心变化均呈相同的正态分布,其均值为0,并且不同卫星到接收机之间的相位中心变化互不相关。则三维空间单差相对定位精度的标准差:
其中,PDOP定义为空间位置精度因子,它的值为
如果用hii代表权系数矩阵H的对角元素,其中i=1,2,...,4,H矩阵为H=(CTC)-1,H通常称为权系数阵,他是一个4*4的对称矩阵。
因此,在样本数足够多的情况下,对于单差相对定位方式,三维空间相对定位精度的标准差σδX等于天线相位中心变化产生的测量误差的标准差σ乘以空间位置精度因子(PDOP),理论上σ为倍σpcv,σpcv为天线相位中心变化标准差。
在GPS星座中,已知在遮蔽角为5度时,PDOP值小于3的概率约为99%。进一步分析相对定位精度的分布区间,从图11可以看出,在单差相对定位时,定位误差的标准差小于3倍天线相位中心变化产生的测量误差的标准差的置信度为99%。因此,可以认为:在不考虑其他误差影响的情况下,对于单差相对定位方式,相对定位精度的标准差小于3倍天线相位中心变化引起的测量误差的标准差。
GPS天线相位中心变化引起的双差相对定位精度公式(13)中,假设GPS天线到各个卫星产生的相位中心变化均呈相同的正态分布,其均值为0,并且不同卫星到接收机之间的相位中心变化互不相关。则三维空间双差相对定位精度的标准差:
其中,DOP定义为精度因子,它的值为
如果用hii代表权系数矩阵H的对角元素,其中i=1,2,3,H矩阵为H=(ATA)-1,H通常称为权系数阵,它是一个3*3的对称矩阵。
因此,在样本数足够多的情况下,对于双差相对定位方式,三维空间相对定位精度的标准差σδX等于天线相位中心变化产生的测量误差的标准差σ乘以精度因子(DOP),理论上σ为倍σpcv,σpcv为天线相位中心变化标准差。
通常双差定位的精度因子DOP小于单差定位的精度因子PDOP。由于在遮蔽角为5度时,单差定位的精度因子PDOP小于3的概率为99%。我们可以得到DOP通常小于3。
进一步分析相对定位精度的分布区间,从图12可以看出,在双差相对定位时,相对定位精度的标准差小于3倍天线相位中心变化产生的测量误差的标准差的置信度为99%。因此,可以认为:在不考虑其他误差影响的情况下,对于双差相对定位方式,相对定位精度的标准差小于3倍天线相位中心变化引起的测量误差的标准差。
对于该统计结果,现有技术虽然分析了GPS天线几何中心与其电学相位中心不一致,不仅在水平、高程方向存在一个固定的偏差,而且随着接收到的GPS卫星的高度和方位的变化而变化。但没有做出过详细的统计特性估计。
实施例
观测站、参考站接收机均静止,其中观测站的空间坐标为[-113402.1865,-5504362.8394,3209404.3787](单位:米),参考站的空间坐标为[-112409.5336,-5504317.7352,3209516.6545](单位:米)。接收机天线类型为:XX-XX,PCV的三个分量分别表示North、East、Up方向上的长度。选择的观测历元区间为2011-7-1 0:0:00至2011-7-1 1:0:00,历元采样间隔为2分钟。考查天线相位中心变化PCV对相对定位精度的影响,PCV的数据采取标准格式输入,单差方式的仿真结果如图7、图8所示,双差方式的仿真结果如图9、图10所示。
图7和图9是相对定位精度总的长度,表示定位位置与实际位置之间的距离误差;图8和图10中,“定位精度(东方向)”表示定位精度在大地正东方向上的分量,“定位精度(北方向)”表示定位精度在大地正北方向上的分量,“定位精度(水平方向)”表示定位精度在大地水平面上的分量,是东方向和北方向分量的合成,“定位精度(高程方向)”表示定位精度在大地垂直方向上的分量,即天顶方向上的分量。
在图7、图8、图9和图10中,星号表示数学模型分析得到的误差,圆圈表示实际定位程序分析得到的误差。从仿真结果中可以看出,在仅讨论天线相位中心变化对GPS相对定位精度的影响时,单差方式和双差方式,其相对定位精度很相近,差别不大。采用单差方式或者双差方式,由实际工程应用的具体情况决定,这里不作讨论。可以看出,相对定位精度修正量很小,约为0.0147毫米,基本可以忽略,因此相对定位可以大大提高定位精度。
本发明说明书中未作详细描述的内容属本领域技术人员的公知技术。
Claims (2)
1.一种GPS天线相位中心变化对相对定位精度的修正方法,其特征在于,包括以下步骤:
1)两个目标均载有GPS接收机及GPS天线,利用GPS卫星信号对两个目标进行相对定位;
2)确定两个目标的GPS天线相位中心变化;
3)根据GPS天线相位中心变化确定目标相对定位精度修正量:
单差方式进行相对定位时,根据所述天线相位中心变化通过下式确定与目标相对定位精度修正量的关系:
其中,σδX为用于表示目标相对定位精度修正量的标准差;PDOP为空间位置精度因子;σ为用于表示GPS天线相位中心变化在均值为0时的测量误差的标准差;σpcv为用于表示GPS天线相位中心变化在均值为0时的标准差。
双差方式进行相对定位时,根据所述天线相位中心变化通过下式确定与目标相对定位精度修正量的关系:
其中,σδX为用于表示目标相对定位精度修正量的标准差;DOP为精度因子;σ为用于表示GPS天线相位中心变化在均值为0时的测量误差的标准差;σpcv为用于表示GPS天线相位中心变化在均值为0时的标准差。
4)利用目标相对定位精度修正量对目标相对定位精度进行修正。
2.如权利要求1所述的一种GPS天线相位中心变化对相对定位精度的修正方法,其特征在于:
单差方式进行相对定位时,所述GPS天线相位中心变化与所述目标相对定位修正量的关系通过下式在GPS天线到各卫星产生的天线相位变化呈正态分布,且不同卫星到接收机之间的相位中心变化互不相关时获得:
其中,δXP为只考虑天线相位中心变化影响的δX2,
δt12=δt2-δt1;
δt1、δt2分别为两个接收机的钟差,c为光速;
A为δX2的系数矩阵,
其中[xp(t1),yp(t1),zp(t1)]为t1时刻卫星p的位置坐标,p=1,2,…n。
D12为只含有天线相位中心变化信息,
双差方式进行相对定位时,所述GPS天线相位中心变化与所述目标相对定位修正量的关系通过下式在GPS天线到各卫星产生的天线相位变化呈正态分布,且不同卫星到接收机之间的相位中心变化互不相关时获得:
δXP=-(ATA)-1AT·ΔD
其中,δXP为只考虑天线相位中心变化影响的δX2,
A为δX2的系数矩阵,
其中[xp(t1),yp(t1),zp(t1)]为t1时刻卫星p的位置坐标,p=1,2,…n。
ΔD为只含有天线相位中心变化信息,
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210334296.1A CN102981174B (zh) | 2012-09-12 | 2012-09-12 | 一种gps天线相位中心变化对相对定位精度的修正方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210334296.1A CN102981174B (zh) | 2012-09-12 | 2012-09-12 | 一种gps天线相位中心变化对相对定位精度的修正方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102981174A true CN102981174A (zh) | 2013-03-20 |
CN102981174B CN102981174B (zh) | 2014-10-08 |
Family
ID=47855378
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210334296.1A Active CN102981174B (zh) | 2012-09-12 | 2012-09-12 | 一种gps天线相位中心变化对相对定位精度的修正方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102981174B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105740517A (zh) * | 2016-01-25 | 2016-07-06 | 西安电子科技大学 | 一种考虑在轨热环境的星载网状可展开天线结构优化方法 |
CN106443722A (zh) * | 2016-09-09 | 2017-02-22 | 中国人民解放军国防科学技术大学 | 一种天线相位中心偏差检测方法 |
CN112327340A (zh) * | 2021-01-06 | 2021-02-05 | 腾讯科技(深圳)有限公司 | 终端定位精度评估方法、装置、设备以及介质 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7382313B1 (en) * | 2004-11-03 | 2008-06-03 | Topcon Gps, Llc | Method for absolute calibration of global navigation satellite system antennas |
CN102323489A (zh) * | 2011-08-18 | 2012-01-18 | 北京华力创通科技股份有限公司 | 对接收机天线的相位中心进行校准的方法和装置 |
-
2012
- 2012-09-12 CN CN201210334296.1A patent/CN102981174B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7382313B1 (en) * | 2004-11-03 | 2008-06-03 | Topcon Gps, Llc | Method for absolute calibration of global navigation satellite system antennas |
CN102323489A (zh) * | 2011-08-18 | 2012-01-18 | 北京华力创通科技股份有限公司 | 对接收机天线的相位中心进行校准的方法和装置 |
Non-Patent Citations (3)
Title |
---|
S.Y.ZHU ET AL.: "Satellite antenna phase center offsets and scale errors in GPS solutions", 《JOURNAL OF GEODESY》 * |
李东皤 等: "天线相位中心偏差对精密定位精度影响的研究", 《第三届中国卫星导航学术年会电子文集》 * |
黄繁荣 等: "星载GPS天线相位中心偏差对定位精度影响的试验方法研究", 《遥测遥控》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105740517A (zh) * | 2016-01-25 | 2016-07-06 | 西安电子科技大学 | 一种考虑在轨热环境的星载网状可展开天线结构优化方法 |
CN105740517B (zh) * | 2016-01-25 | 2019-02-12 | 西安电子科技大学 | 一种考虑在轨热环境的星载网状可展开天线结构优化方法 |
CN106443722A (zh) * | 2016-09-09 | 2017-02-22 | 中国人民解放军国防科学技术大学 | 一种天线相位中心偏差检测方法 |
CN106443722B (zh) * | 2016-09-09 | 2019-05-17 | 中国人民解放军国防科学技术大学 | 一种天线相位中心偏差检测方法 |
CN112327340A (zh) * | 2021-01-06 | 2021-02-05 | 腾讯科技(深圳)有限公司 | 终端定位精度评估方法、装置、设备以及介质 |
Also Published As
Publication number | Publication date |
---|---|
CN102981174B (zh) | 2014-10-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108445518B (zh) | 一种基于双差模糊度固定解约束的gnss精密时间传递方法 | |
CN104714244B (zh) | 一种基于抗差自适应Kalman滤波的多系统动态PPP解算方法 | |
CN105589087B (zh) | 一种高精度实时卫星定位装置及其方法 | |
CN102753991B (zh) | 短基线和超短基线相位图 | |
CN109613585A (zh) | 一种对基站天线超短基线gnss双天线实时测向的方法 | |
CN103814307A (zh) | 包括抗干扰天线并且利用经相位中心校正的载波的gnss定位系统 | |
CN102998690A (zh) | 一种基于gps载波双差方程的姿态角直接求解方法 | |
CN111399020A (zh) | 一种定向测姿系统及方法 | |
Dong et al. | Multi-antenna synchronized global navigation satellite system receiver and its advantages in high-precision positioning applications | |
CN102830406B (zh) | 一种gps天线相位中心变化对绝对定位精度的修正方法 | |
CN105549048B (zh) | 基于三角形接收机阵列的高精度卫星定位装置及其方法 | |
Seepersad | Reduction of initial convergence period in GPS PPP data processing | |
ZhG | Entanglement difference of GNSS carrier phase for vehicle attitude determination | |
Zhang et al. | Initial performance evaluation of precise point positioning with triple-frequency observations from BDS-2 and BDS-3 satellites | |
Kersten et al. | Comparison concept and quality metrics for GNSS antenna calibrations: Cause and effect on regional GNSS networks | |
CN102981174A (zh) | 一种gps天线相位中心变化对相对定位精度的修正方法 | |
US10935669B2 (en) | Using SDP relaxation for optimization of the satellites set chosen for positioning | |
Lee et al. | Analysis of pseudolite augmentation for GPS airborne applications | |
CN109254270A (zh) | 一种星载x波段合成孔径雷达干涉定标方法 | |
Kim et al. | Precise phase calibration of a controlled reception pattern GPS antenna for JPALS | |
Bisnath | Relative Positioning and Real‐Time Kinematic (RTK) | |
Peng et al. | Precise orbit determination for Jason-1 satellite using on-board GPS data with cm-level accuracy | |
Grgac et al. | Quality assessment of Locata positioning system | |
Milbert | Influence of pseudorange accuracy on phase ambiguity resolution in various GPS modernization scenarios | |
Zhao et al. | Error analysis for the baseline estimation and calibration of distributed InSAR satellites |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |