CN104199280A - 一种基于差分gps的时间同步误差测量方法 - Google Patents
一种基于差分gps的时间同步误差测量方法 Download PDFInfo
- Publication number
- CN104199280A CN104199280A CN201410488224.1A CN201410488224A CN104199280A CN 104199280 A CN104199280 A CN 104199280A CN 201410488224 A CN201410488224 A CN 201410488224A CN 104199280 A CN104199280 A CN 104199280A
- Authority
- CN
- China
- Prior art keywords
- mobile node
- time
- differential gps
- synchronization error
- coordinate
- 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
- 238000000691 measurement method Methods 0.000 title abstract 2
- 238000000034 method Methods 0.000 claims abstract description 25
- 230000005540 biological transmission Effects 0.000 claims abstract description 18
- 238000004364 calculation method Methods 0.000 claims description 12
- 238000004891 communication Methods 0.000 claims description 8
- 238000006243 chemical reaction Methods 0.000 claims description 5
- 238000005259 measurement Methods 0.000 abstract description 10
- 230000001360 synchronised effect Effects 0.000 description 10
- 238000010586 diagram Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 238000012360 testing method Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000004069 differentiation Effects 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 238000010998 test method Methods 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
Landscapes
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明提供了一种基于差分GPS的时间同步误差测量方法。具体方法为:定义各个移动节点,并测量移动节点之间的通道延时补偿值;计算出各个移动节点的差分GPS坐标值;依据传输时间计算出各个移动节点之间的传输时间伪距,依据差分GPS坐标值计算三维坐标伪距;根据所述两个伪距求出移动节点A和B之间的时间同步误差。解决了处于高速机动状态,距离达到百公里量级节点间纳秒级时间同步误差实时测量的问题,保证了时间同步测量精度在2ns之内,可以用于10ns左右同步误差的测量。
Description
技术领域
本发明涉及一种通信领域中基于差分GPS的时间同步误差测量方法,特别是涉及一种适用于测量机动节点的时间同步误差测量方法。
背景技术
在通信领域,组网通信的各节点之间对时间同步具有较高要求,有时甚至达到几纳秒量级。针对处于高速机动状态,通信距离达到百公里量级的节点,如何实时、高精度测量多节点之间纳秒级的同步误差,是该领域的一个难点。时间同步误差作为时间统一系统的关键指标,其测试方法有多种。
目前,精度较高的时间同步误差测量基本采用GPS辅助技术,利用GPS,并采用复杂可编程逻辑器件实现的时间同步测量系统。该方法主要利用GPS秒脉冲作为基准信号,测试时间统一系统秒脉冲与基准秒脉冲之间的时间差。
该方法的步骤主要包括:
(1)将GPS秒脉冲和时间统一系统秒脉冲做抗干扰处理,使得小于20ns的干扰信号被滤掉;
(2)通过一系列逻辑门操作,得到GPS秒脉冲和时间统一系统秒脉冲前沿之间的误差计数值;
(3)采用D触发器判断GPS秒脉冲与时间统一系统秒脉冲之间的先后关系。
最终得到GPS秒脉冲和时间统一系统秒脉冲时间的误差值。
该方法的测试误差主要为GPS秒脉冲同步误差,当节点之间距离较远(如达到百公里量级),且处于相对运动时,GPS秒脉冲同步误差仅能达到几十纳秒量级。
差分GPS技术可以获得很高的定位精度,从而得到精度较高的节点间距离,但如何利用这一优势测量机动节点时间同步误差,目前尚无文献给出详细的方法。
发明内容
本发明要解决的技术问题是提供一种基于差分GPS的机动点时间同步误差测量方法。
本发明采用的技术方案如下:一种基于差分GPS的时间同步误差测量方法,其方法步骤为:
步骤一、定义各个移动节点,并测量移动节点之间的通道延时补偿值;
步骤二、计算出各个移动节点的差分GPS坐标值;
步骤三、计算出各个移动节点之间的传输时间伪距和三维坐标伪距;
步骤四、根据所述两个伪距求出移动节点A和B之间的时间同步误差。
作为优选,所述各个移动节点的差分GPS坐标值的计算方法为:基准站与各个移动节点通过各自的GPS接收机同步观测GPS卫星,并得到各自的载波相位;基准站将测算到的载波相位实时发送给各个移动节点,各个移动节点分别将此载波相位与各自测算的载波相位求差分,进而实时解算出各自的差分GPS坐标值。
作为优选,所述传输时间伪距的计算方法为:在所述各个移动节点获得差分GPS坐标值的时间段,按照一定周期重复下述操作:移动节点发送信息到另外一个移动节点,发送信息包含发送时间,所述另外一个移动节点记录到达的本地时间,根据到达的本地时间和发送时间之间的时间间隔计算出传输时间伪距。
作为优选,所述传输时间伪距的计算方法还包括:在进行伪距计算时,对光速进行修正,根据射频通道延迟补偿值对传输时间进行修正。
作为优选,所述三维坐标伪距的计算方法为:将所述各个移动节点的差分GPS坐标转换为大地直角坐标,利用插值算法对各个移动节点的坐标进行插值,并依据大地直角坐标进行距离计算;同时,移动节点实时将获得的差分GPS坐标值发送到另外一个移动节点,实时计算出三维坐标伪距。
作为优选,所述插值算法为均匀插值算法。
与现有技术相比,本发明的有益效果是:解决了处于高速机动状态,距离达到百公里量级节点间纳秒级时间同步误差实时测量的问题,保证了时间同步测量精度在2ns之内,可以用于10ns左右同步误差的测量。
附图说明
图1为本发明其中一具体实施例的基准站O原理性结构框图。
图2为图1所示实施例的移动节点原理性结构框图。
图3为本发明其中一具体实施例的误差测量流程图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅用以解释本发明,并不用于限定本发明。
本说明书(包括任何附加权利要求、摘要和附图)中公开的任一特征,除非特别叙述,均可被其他等效或者具有类似目的的替代特征加以替换。即,除非特别叙述,每个特征只是一系列等效或类似特征中的一个例子而已。
在本具体实施例中,以移动节点A和移动节点B为例进行具体说明。
如图1和图2所示,在本具体实施例中,将一台GPS接收机设置在位置固定的基准站内,该基准站在WGS-84统一坐标系中拥有已知的精密坐标。同时每个移动节点配备一台GPS差分接收机。基准站和移动节点配置有收发装置,可以实现相互通信。
如图3所示,定义移动节点A、移动节点B,对节点A和B进行时间同步误差测量,其方法步骤为:
步骤一、测量移动节点之间的通道延时补偿值。安装地面辅助测试基准站O,其具有精确的GPS坐标(B,L,H),节点A和B通道延时补偿值为Δt'。
步骤二、计算出两个移动节点的差分GPS坐标值。在本具体实施例中,所述两个移动节点的GPS坐标值的计算方法为:基准站O与两个移动节点A和B通过各自的GPS接收机同步观测GPS卫星,并从接收的信号中得到各自的载波相位;基准站O将测算到的载波相位实时发送给移动节点A和B,移动节点A和B分别将此载波相位与各自测算的载波相位求差分,进而实时解算出各自的差分GPS坐标值(Ba,La,Ha)和(Bb,Lb,Hb)。
步骤三、计算出移动节点A和B之间的传输时间伪距L2和三维坐标伪距L1。
在本具体实施例中,所述传输时间伪距L2的计算方法为:移动节点A和B获得差分GPS坐标值的时间段,按照一定周期重复下述操作::移动节点A发送信息到移动节点B,发送信息包含发送时间Nt,移动节点B记录到达的本地时间Nr,根据到达的本地时间和发送时间之间的时间间隔计算出传输时间伪距L2;同时,移动节点A实时将获得的差分GPS坐标值发送到移动节点B,实时计算出三维坐标伪距。
在本具体实施例中,所述三维坐标伪距的计算方法为:将移动节点A和B的差分GPS坐标转换为大地直角坐标,并利用插值算法(在本具体实施例中为均匀插值算法)对移动节点A和B的坐标插值。
在本具体实施例中,计算三维坐标伪距时,使用地心空间大地直角坐标系,因此首先需要进行坐标转换,将节点A的差分GPS坐标(Bb,Lb,Hb)转换为大地直角坐标(Xb,Yb,Zb),节点B的差分GPS坐标(Ba,La,Ha)转换为大地直角坐标(Xa,Ya,Za),GPS坐标向大地直角坐标的转换公式为:
其中,e为地球扁率,取值1/298.257223563,N为卯酉面内曲率半径,a为地球长半轴,取值为6378137米。另外由于GPS信息刷新较慢,节点A相对于节点B处于机动状态,为获得更精确的实时坐标值样本,需要利用插值算法对节点A的坐标插值,即从(Xa,Ya,Za)插值得到(Xa j,Ya j,Za j)。两个伪距的数学表达式为:
L2=c0(Δt+δt+Δt')=c0((Nr-Nt)T-Δt') (4)
其中c0为光在大气中的传播速率,Δt为信号在节点A与节点B之间的传输时间,δt为节点A与B之间的时间同步误差。
所述传输时间伪距的计算方法还包括:在进行伪距计算时,对光速进行修正,使得大气折射带来的测量误差小于0.2ns。在实际计算时,需要根据当地气温、气压等环境因素进行修正,修正的公式为:c0=c/n,其中c为光在真空中的传输速率,取值为2.99792458×108m/s,n为传输环境下的大气折射率。
步骤四、根据所述两个伪距求出移动节点A和B之间的时间同步误差。
由于采用了载波相位求差法的差分GPS技术,多数共模误差都可以消除或者减弱,使得移动节点的定位精度高,伪距L1误差很小,一般可以控制在几个厘米之内,因此可以使用L1作为真实距离的参考值,求取两个节点之间的时间同步误差δt:
本具体实施例保证了时间同步测量精度在2ns之内,可以用于10ns左右同步误差的测量。
Claims (6)
1.一种基于差分GPS的时间同步误差测量方法,其方法步骤为:
步骤一、定义各个移动节点,并测量移动节点之间的通道延时补偿值;
步骤二、计算出各个移动节点的差分GPS坐标值;
步骤三、计算出各个移动节点之间的传输时间伪距和三维坐标伪距;
步骤四、根据所述两个伪距求出移动节点A和B之间的时间同步误差。
2.根据权利要求1所述的时间同步误差测量方法,所述各个移动节点的差分GPS坐标值的计算方法为:基准站与各个移动节点通过各自的GPS接收机同步观测GPS卫星,并得到各自的载波相位;基准站将测算到的载波相位发送给各个移动节点,各个移动节点分别将此载波相位与各自测算的载波相位求差分,进而解算出各自的差分GPS坐标值。
3.根据权利要求1所述的时间同步误差测量方法,所述传输时间伪距的计算方法为:在所述各个移动节点获得差分GPS坐标值的时间段,按照一定周期重复下述操作:移动节点发送信息到另外一个移动节点,发送信息包含发送时间,所述另外一个移动节点记录到达的本地时间,根据到达的本地时间和发送时间之间的时间间隔计算出传输时间伪距。
4.根据权利要求3所述的时间同步误差测量方法,所述传输时间伪距的计算方法还包括:在进行伪距计算时,对光速进行修正,根据射频通道延迟补偿值对传输时间进行修正。
5.根据权利要求1到4之一所述的时间同步误差测量方法,所述三维坐标伪距的计算方法为:将所述各个移动节点的差分GPS坐标转换为大地直角坐标,利用插值算法对各个移动节点的坐标进行插值,并依据大地直角坐标进行距离计算;同时,移动节点实时将获得的差分GPS坐标值发送到另外一个移动节点,实时计算出三维坐标伪距。
6.根据权利要求5所述的时间同步误差测量方法,所述插值算法为均匀插值算法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410488224.1A CN104199280B (zh) | 2014-09-23 | 2014-09-23 | 一种基于差分gps的时间同步误差测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410488224.1A CN104199280B (zh) | 2014-09-23 | 2014-09-23 | 一种基于差分gps的时间同步误差测量方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104199280A true CN104199280A (zh) | 2014-12-10 |
CN104199280B CN104199280B (zh) | 2017-12-15 |
Family
ID=52084583
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410488224.1A Expired - Fee Related CN104199280B (zh) | 2014-09-23 | 2014-09-23 | 一种基于差分gps的时间同步误差测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104199280B (zh) |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105911858A (zh) * | 2016-04-11 | 2016-08-31 | 中国科学院国家授时中心 | 一种亚纳秒级远程时间比对/传递系统误差测试方法 |
CN108540207A (zh) * | 2018-04-20 | 2018-09-14 | 中国电子科技集团公司第二十九研究所 | 一种卫星载荷分布式时间同步方法及系统 |
CN108632979A (zh) * | 2017-03-20 | 2018-10-09 | 中兴通讯股份有限公司 | 一种优化时间同步误差的方法、装置和设备 |
CN108761434A (zh) * | 2018-03-27 | 2018-11-06 | 中国电子科技集团公司第二十九研究所 | 一种基于cdma下行信号的伪距测量方法及系统 |
WO2018205811A1 (zh) * | 2017-05-10 | 2018-11-15 | 中兴通讯股份有限公司 | 时钟同步的方法、时间基准源设备和时钟复现设备 |
CN109655846A (zh) * | 2019-01-30 | 2019-04-19 | 清华大学 | 一种多站差分后处理高精度时间同步方法及系统 |
CN110850367A (zh) * | 2019-11-14 | 2020-02-28 | 南京邮电大学 | 一种在定位系统中测量空气折射率及晶振实际频率并应用于距离测量的方法 |
CN110988947A (zh) * | 2019-02-20 | 2020-04-10 | 以见科技(上海)有限公司 | 一种基于实时动态载波相位差分技术的增强现实定位方法 |
CN111459009A (zh) * | 2020-04-21 | 2020-07-28 | 哈尔滨工业大学 | 一种多数字电子设备同步的随机误差估计系统及其估计方法 |
CN111605734A (zh) * | 2020-04-28 | 2020-09-01 | 北京控制工程研究所 | 一种航天器三超控制星时准确同步的滤波校时系统及方法 |
CN112229401A (zh) * | 2020-09-15 | 2021-01-15 | 北京菲斯罗克光电技术有限公司 | 适用于ins-gps伪距融合的量测信息同步外推方法及系统 |
CN112260790A (zh) * | 2020-10-28 | 2021-01-22 | 中标慧安信息技术股份有限公司 | 一种基于ntp协议的时间同步方法及系统 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1425601A (zh) * | 2001-12-12 | 2003-06-25 | 上海振华港口机械股份有限公司 | 用于轮胎式龙门集装箱起重机的卫星定位系统移动站 |
CN101750620A (zh) * | 2009-12-25 | 2010-06-23 | 三一重工股份有限公司 | 臂架系统的定位方法、定位装置及混凝土泵车 |
CN201567797U (zh) * | 2009-12-25 | 2010-09-01 | 三一重工股份有限公司 | 臂架系统的定位装置、臂架系统、混凝土泵车及布料架 |
CN102137488A (zh) * | 2011-02-16 | 2011-07-27 | 北京邮电大学 | 定位方法及装置 |
CN102196559A (zh) * | 2011-05-04 | 2011-09-21 | 西安电子科技大学 | 基于tdoa定位的通道时延误差消除方法 |
CN103472459A (zh) * | 2013-08-29 | 2013-12-25 | 镇江青思网络科技有限公司 | 一种基于gps伪距差分的车辆协作定位方法 |
CN103210690B (zh) * | 2011-09-09 | 2014-12-31 | 华为技术有限公司 | 一种时间同步的方法和系统及节点设备 |
-
2014
- 2014-09-23 CN CN201410488224.1A patent/CN104199280B/zh not_active Expired - Fee Related
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1425601A (zh) * | 2001-12-12 | 2003-06-25 | 上海振华港口机械股份有限公司 | 用于轮胎式龙门集装箱起重机的卫星定位系统移动站 |
CN101750620A (zh) * | 2009-12-25 | 2010-06-23 | 三一重工股份有限公司 | 臂架系统的定位方法、定位装置及混凝土泵车 |
CN201567797U (zh) * | 2009-12-25 | 2010-09-01 | 三一重工股份有限公司 | 臂架系统的定位装置、臂架系统、混凝土泵车及布料架 |
CN102137488A (zh) * | 2011-02-16 | 2011-07-27 | 北京邮电大学 | 定位方法及装置 |
CN102196559A (zh) * | 2011-05-04 | 2011-09-21 | 西安电子科技大学 | 基于tdoa定位的通道时延误差消除方法 |
CN103210690B (zh) * | 2011-09-09 | 2014-12-31 | 华为技术有限公司 | 一种时间同步的方法和系统及节点设备 |
CN103472459A (zh) * | 2013-08-29 | 2013-12-25 | 镇江青思网络科技有限公司 | 一种基于gps伪距差分的车辆协作定位方法 |
Cited By (20)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105911858B (zh) * | 2016-04-11 | 2019-05-10 | 中国科学院国家授时中心 | 一种亚纳秒级远程时间比对/传递系统误差测试方法 |
CN105911858A (zh) * | 2016-04-11 | 2016-08-31 | 中国科学院国家授时中心 | 一种亚纳秒级远程时间比对/传递系统误差测试方法 |
CN108632979A (zh) * | 2017-03-20 | 2018-10-09 | 中兴通讯股份有限公司 | 一种优化时间同步误差的方法、装置和设备 |
WO2018205811A1 (zh) * | 2017-05-10 | 2018-11-15 | 中兴通讯股份有限公司 | 时钟同步的方法、时间基准源设备和时钟复现设备 |
CN108882356A (zh) * | 2017-05-10 | 2018-11-23 | 深圳市中兴软件有限责任公司 | 时钟同步的方法、时间基准源设备和时钟复现设备 |
CN108761434A (zh) * | 2018-03-27 | 2018-11-06 | 中国电子科技集团公司第二十九研究所 | 一种基于cdma下行信号的伪距测量方法及系统 |
CN108761434B (zh) * | 2018-03-27 | 2022-06-28 | 中国电子科技集团公司第二十九研究所 | 一种基于cdma下行信号的伪距测量方法及系统 |
CN108540207A (zh) * | 2018-04-20 | 2018-09-14 | 中国电子科技集团公司第二十九研究所 | 一种卫星载荷分布式时间同步方法及系统 |
CN109655846A (zh) * | 2019-01-30 | 2019-04-19 | 清华大学 | 一种多站差分后处理高精度时间同步方法及系统 |
CN110988947B (zh) * | 2019-02-20 | 2020-09-22 | 以见科技(上海)有限公司 | 一种基于实时动态载波相位差分技术的增强现实定位方法 |
CN110988947A (zh) * | 2019-02-20 | 2020-04-10 | 以见科技(上海)有限公司 | 一种基于实时动态载波相位差分技术的增强现实定位方法 |
CN110850367A (zh) * | 2019-11-14 | 2020-02-28 | 南京邮电大学 | 一种在定位系统中测量空气折射率及晶振实际频率并应用于距离测量的方法 |
CN111459009A (zh) * | 2020-04-21 | 2020-07-28 | 哈尔滨工业大学 | 一种多数字电子设备同步的随机误差估计系统及其估计方法 |
CN111459009B (zh) * | 2020-04-21 | 2021-08-17 | 哈尔滨工业大学 | 一种多数字电子设备同步的随机误差估计系统及其估计方法 |
CN111605734B (zh) * | 2020-04-28 | 2021-11-16 | 北京控制工程研究所 | 一种航天器三超控制星时准确同步的滤波校时系统及方法 |
CN111605734A (zh) * | 2020-04-28 | 2020-09-01 | 北京控制工程研究所 | 一种航天器三超控制星时准确同步的滤波校时系统及方法 |
CN112229401A (zh) * | 2020-09-15 | 2021-01-15 | 北京菲斯罗克光电技术有限公司 | 适用于ins-gps伪距融合的量测信息同步外推方法及系统 |
CN112229401B (zh) * | 2020-09-15 | 2022-08-05 | 北京菲斯罗克光电技术有限公司 | 适用于ins-gps伪距融合的量测信息同步外推方法及系统 |
CN112260790A (zh) * | 2020-10-28 | 2021-01-22 | 中标慧安信息技术股份有限公司 | 一种基于ntp协议的时间同步方法及系统 |
CN112260790B (zh) * | 2020-10-28 | 2021-08-31 | 中标慧安信息技术股份有限公司 | 一种基于ntp协议的时间同步方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN104199280B (zh) | 2017-12-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104199280A (zh) | 一种基于差分gps的时间同步误差测量方法 | |
CN103399326B (zh) | 一种gnss动态测量精度检验系统及方法 | |
US8077089B2 (en) | Precision geolocation of moving or fixed transmitters using multiple observers | |
EP2634593B1 (en) | Positioning using a local wave-propagation model | |
US20150369924A1 (en) | Method and system for high-accuracy differential tracking of global positioning system (gps) receivers | |
CN103344942B (zh) | 控制节点、异步定位方法与系统 | |
JP5078082B2 (ja) | 測位装置、測位システム、コンピュータプログラム及び測位方法 | |
CN104316947B (zh) | Gnss/ins超紧组合导航装置及相对导航系统 | |
CN108842544B (zh) | 一种利用光纤惯导的高铁轨道静态参数检测系统和方法 | |
CN102253399B (zh) | 一种利用载波相位中心值的多普勒差分补偿测速方法 | |
EP2689268A1 (en) | Method, apparatus and system for determining a position of an object having a global navigation satellite system receiver by processing undifferenced data like carrier phase measurements and external products like ionosphere data | |
CN103576175A (zh) | 一种双频多星座gnss整周模糊度otf解算方法 | |
CN104536027A (zh) | 一种实时北斗精密相对定位的方法 | |
CN103293947A (zh) | 一种星地激光时间比对系统 | |
CN104237920A (zh) | 高反射环境中的交通工具定位 | |
CN110716217A (zh) | 一种阵列伪卫星室内定位方法及系统 | |
CN103645485A (zh) | 一种基于双星时差频差无源定位的伪距差分方法 | |
CN111970002B (zh) | 基于北斗geo卫星的原子钟远程频率传递与比对方法 | |
CN104880717A (zh) | 一种卫星测量值多路径误差检测装置及其算法 | |
CN108345018A (zh) | 一种自适应卫星导航多模式相对定位方法 | |
CN104730551B (zh) | 一种星地双基地差分干涉基线坐标以及形变量测量方法 | |
KR101851853B1 (ko) | 항체 사이의 상대 위치 결정시스템 및 이의 방법 | |
Gogoi et al. | On the cooperative ranging between android smartphones sharing raw GNSS measurements | |
CN102087363B (zh) | 一种中轨道卫星搜救系统定位方法 | |
CN109870672A (zh) | 一种基于锚节点差分时间同步和Taylor协同的定位算法 |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20171215 |