CN105573109B - 基于引力频移的授时方法和系统 - Google Patents

基于引力频移的授时方法和系统 Download PDF

Info

Publication number
CN105573109B
CN105573109B CN201610064927.0A CN201610064927A CN105573109B CN 105573109 B CN105573109 B CN 105573109B CN 201610064927 A CN201610064927 A CN 201610064927A CN 105573109 B CN105573109 B CN 105573109B
Authority
CN
China
Prior art keywords
local
satellite system
beidou satellite
frequency
strange land
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.)
Expired - Fee Related
Application number
CN201610064927.0A
Other languages
English (en)
Other versions
CN105573109A (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.)
Jianghan University
Original Assignee
Jianghan University
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 Jianghan University filed Critical Jianghan University
Priority to CN201610064927.0A priority Critical patent/CN105573109B/zh
Publication of CN105573109A publication Critical patent/CN105573109A/zh
Application granted granted Critical
Publication of CN105573109B publication Critical patent/CN105573109B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G04HOROLOGY
    • G04RRADIO-CONTROLLED TIME-PIECES
    • G04R20/00Setting the time according to the time information carried or implied by the radio signal
    • G04R20/02Setting the time according to the time information carried or implied by the radio signal the radio signal being sent by a satellite, e.g. GPS
    • G04R20/04Tuning or receiving; Circuits therefor

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Position Fixing By Use Of Radio Waves (AREA)

Abstract

本发明公开了一种基于引力频移的授时方法和系统,属于时间频率比对领域。所述方法包括:获取本地时间频率标准T1;记录北斗卫星系统的参考时间频率T0;计算所述本地时间频率标准T1与本地记录的所述参考时间频率T0的差值ΔT1;获取异地计算出的异地时间频率标准T2与异地记录的北斗卫星系统的参考时间频率T0的差值ΔT2;计算两地时间频率比对结果ΔT,ΔT=ΔT1‑ΔT2‑ΔT0,ΔT0=(U1‑U2)/C,其中,所述ΔT0为频移量,所述U1为北斗卫星系统上测得的本地的引力势,所述U2为北斗卫星系统上测得的异地的引力势,所述C为光速,根据上述公式算得的ΔT中除去了引力频移,提高了时间频率的比对的精度。

Description

基于引力频移的授时方法和系统
技术领域
本发明涉及时间频率比对领域,特别涉及一种基于引力频移的授时方法和系统。
背景技术
随着原子钟和通信技术的发展,时间频率比对技术得到不断进步。时间频率比对技术主要包括时间频率的传递和同步。
按照时间频率传递技术出现的先后顺序,远距离的时间频率传递包括短波无线电播时、LORAN-C长波导航系统、以GPS(Global Positioning System,全球定位系统)为代表的星基导航、及网络传递时间频率。其中,由于高精度和高可靠性,GPS在远距离的时间频率传递中占主导地位。通过GPS进行时间频率传递的方法为,处于地球上任何两个地点的原子钟利用同一时刻收到的同一颗GPS导航卫星的时间信号进行时间频率比对。
在实现本发明的过程中,发明人发现现有技术至少存在以下问题:
现有技术中时间频率比对结果精度低,不能满足研究需求。
发明内容
为了解决现有技术的问题,本发明实施例提供了一种基于引力频移的授时方法和系统。所述技术方案如下:
第一方面,本发明实施例提供了一种基于引力频移的授时方法,所述方法包括:
获取本地时间频率标准T1;
记录北斗卫星系统的参考时间频率T0;
计算所述本地时间频率标准T1与本地记录的所述参考时间频率T0的差值ΔT1;
获取异地计算出的异地时间频率标准T2与异地记录的北斗卫星系统的参考时间频率T0的差值ΔT2;
计算两地时间频率比对结果ΔT,ΔT=ΔT1-ΔT2-ΔT0,ΔT0=(U1-U2)/C,其中,所述ΔT0为频移量,所述U1为北斗卫星系统上测得的本地的引力势,所述U2为北斗卫星系统上测得的异地的引力势,所述C为光速;
所述方法还包括:
获取所述北斗卫星系统通过重力梯度仪测得的U1和U2,所述U1是北斗卫星系统中处在地心和本地连线的延长线上的卫星测得的,所述U2是北斗卫星系统中处在地心和异地连线的延长线上的卫星测得的。
在本发明实施例的一种实现方式中,所述获取本地时间频率标准T1,包括:
将处于本地的本地接收机同步到本地参考源上,使所述本地接收机获取所述本地参考源提供的所述本地时间频率标准T1。
在本发明实施例的另一种实现方式中,所述获取异地计算出的异地时间频率标准T2与异地记录的北斗卫星系统的参考时间频率T0的差值ΔT2,包括:
接收处于异地的异地接收机发送的所述ΔT2。
第二方面,本发明实施例还提供了一种基于引力频移的授时系统,所述系统包括:
本地参考源,用于提供本地时间频率标准T1;
本地接收机,用于记录北斗卫星系统的参考时间频率T0;计算所述本地时间频率标准T1与本地记录的所述北斗卫星系统的参考时间频率T0的差值ΔT1;
异地参考源,用于提供异地时间频率标准T2;
异地接收机,用于记录的北斗卫星系统的参考时间频率T0;计算所述异地时间频率标准T2与异地记录的所述北斗卫星系统的参考时间频率T0的差值ΔT2;
所述本地接收机和所述异地接收机,还用于将所述ΔT1或所述ΔT2传递给对端的接收机;
所述本地接收机和所述异地接收机,还用于计算两地时间频率比对结果ΔT,ΔT=ΔT1-ΔT2-ΔT0,ΔT0=(U1-U2)/C,其中,所述ΔT0为频移量,所述U1为北斗卫星系统上测得的本地的引力势,所述U2为北斗卫星系统上测得的异地的引力势,所述C为光速;所述U1是北斗卫星系统中处在地心和本地连线的延长线上的卫星测得的,所述U2是北斗卫星系统中处在地心和异地连线的延长线上的卫星测得的;
所述本地接收机和异地接收机分别用于与所述北斗卫星系统无线连接,所述本地接收机和异地接收机之间无线连接,所述本地接收机和所述本地参考源电连接,所述异地接收机和所述异地参考源电连接。
在本发明实施例的一种实现方式中,所述本地接收机和异地接收机之间通过卫星通信网络相互连接。
在本发明实施例的另一种实现方式中,所述本地接收机和所述异地接收机均为全球定位系统接收机。
在本发明实施例的另一种实现方式中,所述本地参考源和所述异地参考源均为原子钟。
在本发明实施例的另一种实现方式中,所述本地接收机和异地接收机均包括:
相连的第一接收单元和第二接收单元,所述第一接收单元通过所述北斗卫星系统的通信及处理装置与所述北斗卫星系统上的高精度空间冷原子钟连接,所述第二接收单元通过所述北斗卫星系统的通信及处理装置与所述北斗卫星系统上的重力梯度仪连接。
本发明实施例提供的技术方案带来的有益效果是:
通过计算本地时间频率标准T1与本地记录的所述参考时间频率T0的差值ΔT1,然后接收异地计算出的异地时间频率标准T2与异地记录的北斗卫星系统的参考时间频率T0的差值ΔT2,然后计算两地时间频率比对结果ΔT,ΔT=ΔT1-ΔT2-ΔT0,ΔT0=(U1-U2)/C,其中,所述ΔT0为频移量,所述U1为北斗卫星系统上测得的本地的引力势,所述U2为北斗卫星系统上测得的异地的引力势,所述C为光速;根据爱因斯坦的等效原理,处在两个不同的引力势U1和U2中的辐射源相对于观测者的辐射频率会有一定的引力频移,频移量即为前述ΔT0,因此,上述公式算得的ΔT中除去了引力频移,提高了时间频率的比对的精度。
附图说明
为了更清楚地说明本发明实施例中的技术方案,下面将对实施例描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1是本发明实施例提供的一种基于引力频移的授时方法的流程图;
图2是本发明实施例提供的北斗卫星系统测量引力势的位置示意图;
图3是本发明实施例提供的一种基于引力频移的授时系统的结构示意图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,下面将结合附图对本发明实施方式作进一步地详细描述。
图1是本发明实施例提供的一种基于引力频移的授时方法的流程图,参见图1,方法包括:
步骤101:获取本地时间频率标准T1。
在本发明实施例中,获取本地时间频率标准T1,包括:
将处于本地的本地接收机同步到本地参考源上,使本地接收机获取本地参考源提供的本地时间频率标准T1。
步骤102:记录北斗卫星系统的参考时间频率T0。
即接收北斗卫星系统下发的参考时间频率T0,该参考时间频率T0是处在地心和本地连线的延长线上的卫星提供的。
步骤103:计算本地时间频率标准T1与本地记录的参考时间频率T0的差值ΔT1。
步骤104:获取异地计算出的异地时间频率标准T2与异地记录的北斗卫星系统的参考时间频率T0的差值ΔT2。
在本发明实施例中,获取异地计算出的异地时间频率标准T2与异地记录的北斗卫星系统的参考时间频率T0的差值ΔT2,包括:
接收处于异地的异地接收机发送的ΔT2。
在本发明实施例中,异地记录的北斗卫星系统的参考时间频率T0与本地记录的北斗卫星系统的参考时间频率T0为同一时刻北斗卫星系统下发的参考时间频率,异地的参考时间频率T0是处在地心和异地连线的延长线上的卫星提供的。但是本地记录的是地心和本地连线的延长线上的卫星下发的T0,异地记录的是地心和异地连线的延长线上的卫星下发的T0,由于两个位置受到不同的引力,造成引力频移等原因,使得两地获得的T0并不相等,存在一定差值。
步骤105:计算两地时间频率比对结果ΔT,ΔT=ΔT1-ΔT2-ΔT0,ΔT0=(U1-U2)/C,其中,ΔT0为频移量,U1为北斗卫星系统上测得的本地的引力势,U2为北斗卫星系统上测得的异地的引力势,C为光速。
在本发明实施例中,本地和异地用于表示两个不同的地点。
进一步地,该方法还包括:
获取北斗卫星系统通过重力梯度仪测得的U1和U2,U1是北斗卫星系统中处在地心和本地连线的延长线上的卫星测得的,U2是北斗卫星系统中处在地心和异地连线的延长线上的卫星测得的。
如图2所示,北斗卫星系统30中一卫星处在地心0和本地1连线D1的延长线d1上时,可以测得本地的引力势U1;北斗卫星系统30中另一卫星处在地心0和异地2连线D2的延长线d2上时,可以测得异地的引力势U2。
本发明通过计算本地时间频率标准T1与本地记录的参考时间频率T0的差值ΔT1,然后接收异地计算出的异地时间频率标准T2与异地记录的北斗卫星系统的参考时间频率T0的差值ΔT2,然后计算两地时间频率比对结果ΔT,ΔT=ΔT1-ΔT2-ΔT0,ΔT0=(U1-U2)/C,其中,ΔT0为频移量,U1为北斗卫星系统上测得的本地的引力势,U2为北斗卫星系统上测得的异地的引力势,C为光速;根据爱因斯坦的等效原理,处在两个不同的引力势U1和U2中的辐射源相对于观测者的辐射频率会有一定的引力频移,频移量即为前述ΔT0,因此,上述公式算得的ΔT中除去了引力频移,提高了时间频率的比对的精度。
图3是本发明实施例提供的一种基于引力频移的授时系统的结构示意图,参见图3,系统包括:
本地参考源10,用于提供本地时间频率标准T1;
本地接收机11,用于记录北斗卫星系统30的参考时间频率T0,即接收北斗卫星系统下发的参考时间频率T0;计算本地时间频率标准T1与本地记录的北斗卫星系统30的参考时间频率T0的差值ΔT1;
异地参考源20,用于提供异地时间频率标准T2;
异地接收机21,用于记录的北斗卫星系统30的参考时间频率T0;计算异地时间频率标准T2与异地记录的北斗卫星系统30的参考时间频率T0的差值ΔT2;
本地接收机11和异地接收机21,还用于将ΔT1或ΔT2传递给对端的接收机,本地接收机11的对端为异地接收机21,异地接收机21的对端为本地接收机11;
本地接收机11和异地接收机21,还用于计算两地时间频率比对结果ΔT,ΔT=ΔT1-ΔT2-ΔT0,ΔT0=(U1-U2)/C,其中,ΔT0为频移量,U1为北斗卫星系统30上测得的本地的引力势,U2为北斗卫星系统30上测得的异地的引力势,C为光速;
本地接收机11和异地接收机21分别用于与北斗卫星系统30无线连接,本地接收机11和异地接收机21之间无线连接,本地接收机11和本地参考源10电连接,异地接收机21和异地参考源20电连接。
在本发明实施例中,异地记录的北斗卫星系统的参考时间频率T0与本地记录的北斗卫星系统的参考时间频率T0为同一时刻北斗卫星系统下发的参考时间频率。但由于引力频移等原因,两地获得的T0并不相等,存在一定差值。
在本发明实施例中,本地接收机11和异地接收机21之间通过卫星通信网络相互连接。
在本发明实施例中,本地接收机11和异地接收机21均为全球定位系统接收机。
在本发明实施例中,本地参考源10和异地参考源20均为原子钟。
在本发明实施例中,本地接收机11和异地接收机21均包括:
相连的第一接收单元和第二接收单元,第一接收单元通过北斗卫星系统30的通信及处理装置与北斗卫星系统30上的高精度空间冷原子钟连接,第二接收单元通过北斗卫星系统30的通信及处理装置与北斗卫星系统30上的重力梯度仪连接。
其中,通信及处理装置用于与接收机通信,以及对高精度空间冷原子钟和重力梯度仪的输出信号进行处理。
在本发明实施例中,U1是北斗卫星系统30中处在地心和本地连线的延长线上的卫星测得的,U2是北斗卫星系统30中处在地心和异地连线的延长线上的卫星测得的。
如图2所示,北斗卫星系统30中一卫星处在地心0和本地1连线D1的延长线d1上时,可以测得本地的引力势U1;北斗卫星系统30中另一卫星处在地心0和异地2连线D2的延长线d2上时,可以测得异地的引力势U2。
在本发明实施例中,本地和异地用于表示两个不同的地点。具体地,本地参考源10和本地接收机11设置在同一地理位置,异地参考源20和异地接收机21设置在同一地理位置,本地参考源10和异地参考源20设置在不同地理位置。
进一步地,由于硬件环境和工作机制,接收机与参考源间的时间频率传递过程中会产生误差(参考源的理论时间频率与接收机实际获得的时间频率之间的差值)。将本地和异地的误差分别记为L1和L2,则ΔT=ΔT1-ΔT2-ΔT0+ΔL,ΔL=L1-L2。
本发明通过计算本地时间频率标准T1与本地记录的参考时间频率T0的差值ΔT1,然后接收异地计算出的异地时间频率标准T2与异地记录的北斗卫星系统的参考时间频率T0的差值ΔT2,然后计算两地时间频率比对结果ΔT,ΔT=ΔT1-ΔT2-ΔT0,ΔT0=(U1-U2)/C,其中,ΔT0为频移量,U1为北斗卫星系统上测得的本地的引力势,U2为北斗卫星系统上测得的异地的引力势,C为光速;根据爱因斯坦的等效原理,处在两个不同的引力势U1和U2中的辐射源相对于观测者的辐射频率会有一定的引力频移,频移量即为前述ΔT0,因此,上述公式算得的ΔT中除去了引力频移,提高了时间频率的比对的精度。
本领域普通技术人员可以理解实现上述实施例的全部或部分步骤可以通过硬件来完成,也可以通过程序来指令相关的硬件完成,的程序可以存储于一种计算机可读存储介质中,上述提到的存储介质可以是只读存储器,磁盘或光盘等。
以上仅为本发明的较佳实施例,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (8)

1.一种基于引力频移的授时方法,其特征在于,所述方法包括:
获取本地时间频率标准T1;
记录北斗卫星系统的参考时间频率T0;
计算所述本地时间频率标准T1与本地记录的所述参考时间频率T0的差值ΔT1;
获取异地计算出的异地时间频率标准T2与异地记录的北斗卫星系统的参考时间频率T0的差值ΔT2;
计算两地时间频率比对结果ΔT,ΔT=ΔT1-ΔT2-ΔT0,ΔT0=(U1-U2)/C,其中,所述ΔT0为频移量,所述U1为北斗卫星系统上测得的本地的引力势,所述U2为北斗卫星系统上测得的异地的引力势,所述C为光速;
所述方法还包括:
获取所述北斗卫星系统通过重力梯度仪测得的U1和U2,所述U1是北斗卫星系统中处在地心和本地连线的延长线上的卫星测得的,所述U2是北斗卫星系统中处在地心和异地连线的延长线上的卫星测得的。
2.根据权利要求1所述的方法,其特征在于,所述获取本地时间频率标准T1,包括:
将处于本地的本地接收机同步到本地参考源上,使所述本地接收机获取所述本地参考源提供的所述本地时间频率标准T1。
3.根据权利要求1所述的方法,其特征在于,所述获取异地计算出的异地时间频率标准T2与异地记录的北斗卫星系统的参考时间频率T0的差值ΔT2,包括:
接收处于异地的异地接收机发送的所述ΔT2。
4.一种基于引力频移的授时系统,其特征在于,所述系统包括:
本地参考源,用于提供本地时间频率标准T1;
本地接收机,用于记录北斗卫星系统的参考时间频率T0;计算所述本地时间频率标准T1与本地记录的所述北斗卫星系统的参考时间频率T0的差值ΔT1;
异地参考源,用于提供异地时间频率标准T2;
异地接收机,用于记录的北斗卫星系统的参考时间频率T0;计算所述异地时间频率标准T2与异地记录的所述北斗卫星系统的参考时间频率T0的差值ΔT2;
所述本地接收机和所述异地接收机,还用于将所述ΔT1或所述ΔT2传递给对端的接收机;
所述本地接收机和所述异地接收机,还用于计算两地时间频率比对结果ΔT,ΔT=ΔT1-ΔT2-ΔT0,ΔT0=(U1-U2)/C,其中,所述ΔT0为频移量,所述U1为北斗卫星系统上测得的本地的引力势,所述U2为北斗卫星系统上测得的异地的引力势,所述C为光速;所述U1是北斗卫星系统中处在地心和本地连线的延长线上的卫星测得的,所述U2是北斗卫星系统中处在地心和异地连线的延长线上的卫星测得的;
所述本地接收机和异地接收机分别用于与所述北斗卫星系统无线连接,所述本地接收机和异地接收机之间无线连接,所述本地接收机和所述本地参考源电连接,所述异地接收机和所述异地参考源电连接。
5.根据权利要求4所述的系统,其特征在于,所述本地接收机和异地接收机之间通过卫星通信网络相互连接。
6.根据权利要求4所述的系统,其特征在于,所述本地接收机和所述异地接收机均为全球定位系统接收机。
7.根据权利要求4所述的系统,其特征在于,所述本地参考源和所述异地参考源均为原子钟。
8.根据权利要求4所述的系统,其特征在于,所述本地接收机和异地接收机均包括:
相连的第一接收单元和第二接收单元,所述第一接收单元通过所述北斗卫星系统的通信及处理装置与所述北斗卫星系统上的高精度空间冷原子钟连接,所述第二接收单元通过所述北斗卫星系统的通信及处理装置与所述北斗卫星系统上的重力梯度仪连接。
CN201610064927.0A 2016-01-29 2016-01-29 基于引力频移的授时方法和系统 Expired - Fee Related CN105573109B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610064927.0A CN105573109B (zh) 2016-01-29 2016-01-29 基于引力频移的授时方法和系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610064927.0A CN105573109B (zh) 2016-01-29 2016-01-29 基于引力频移的授时方法和系统

Publications (2)

Publication Number Publication Date
CN105573109A CN105573109A (zh) 2016-05-11
CN105573109B true CN105573109B (zh) 2019-01-29

Family

ID=55883375

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610064927.0A Expired - Fee Related CN105573109B (zh) 2016-01-29 2016-01-29 基于引力频移的授时方法和系统

Country Status (1)

Country Link
CN (1) CN105573109B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113253307B (zh) * 2021-04-14 2024-06-18 长沙学院 旋转天线场景下的gnss多星粗校时方法、系统及存储介质

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN202661628U (zh) * 2012-05-28 2013-01-09 江汉大学 一种时间频率传递装置
CN103162844A (zh) * 2013-02-04 2013-06-19 江汉大学 一种光速各向异性测量方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN202661628U (zh) * 2012-05-28 2013-01-09 江汉大学 一种时间频率传递装置
CN103162844A (zh) * 2013-02-04 2013-06-19 江汉大学 一种光速各向异性测量方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
北斗导航卫星钟相对论效应误差的理论研究;达道安;《2005年全国时间频率学术交流会文集》;20051231;379-393 *

Also Published As

Publication number Publication date
CN105573109A (zh) 2016-05-11

Similar Documents

Publication Publication Date Title
US6492945B2 (en) Instantaneous radiopositioning using signals of opportunity
CN103293512B (zh) 使用本地波传播模型定位
JP4920079B2 (ja) Gps受信機を備えた地震測定システム
CN101867868B (zh) 一种组合导航装置和实施方法
US10520605B2 (en) Satellite signal reception characteristic estimation apparatus, method thereof, and program thereof
CN102033236B (zh) 一种卫星导航位置速度联合估计方法
BR112013017960A2 (pt) método e sistema para determinar correções de relógio de satélite
CN104181550A (zh) 一种基于北斗卫星导航系统的共视时间频率传递方法
CN103797727A (zh) 使用交联测距和准确时间源的用于卫星星座的先进定时和时间传递
Mahmoud et al. Integrated positioning for connected vehicles
CN109100746A (zh) 一种基于转发节点的隧道定位系统及方法
CN104122789A (zh) 高精度分布式同步时钟系统及方法
CN111970002B (zh) 基于北斗geo卫星的原子钟远程频率传递与比对方法
CN103148845A (zh) 一种卫星组合罗盘及方位角和俯仰角测量计算方法
CN105573109B (zh) 基于引力频移的授时方法和系统
CN105158781A (zh) 一种地球gnss卫星实现外层空间用户定位的方法
CN210742507U (zh) 一种基于全球导航卫星系统的标准时间频率源装置
KR20130135635A (ko) 글로벌 위성 항법 시스템을 이용한 대기층의 수증기량 측정시스템 및 측정방법
CN205620685U (zh) 基于引力频移的授时装置
CN104392108B (zh) 一种采用迭代差分算法的远程定位系统及方法
Baselga et al. GBDM+: an improved methodology for a GNSS-based distance meter
CN204065693U (zh) 高精度分布式同步时钟系统
CN110865200A (zh) 地球自转角速率测量方法
Hancock et al. Satellite mapping in cities and below cities: how good is it now?
CN106657419B (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
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20190129

CF01 Termination of patent right due to non-payment of annual fee