CN103033828B - 一种高灵敏度北斗辅助授时装置和授时接收机及授时方法 - Google Patents
一种高灵敏度北斗辅助授时装置和授时接收机及授时方法 Download PDFInfo
- Publication number
- CN103033828B CN103033828B CN201210543822.5A CN201210543822A CN103033828B CN 103033828 B CN103033828 B CN 103033828B CN 201210543822 A CN201210543822 A CN 201210543822A CN 103033828 B CN103033828 B CN 103033828B
- Authority
- CN
- China
- Prior art keywords
- signal
- module
- time
- time service
- week
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 37
- 238000009825 accumulation Methods 0.000 claims abstract description 125
- 238000000605 extraction Methods 0.000 claims abstract description 18
- 230000035945 sensitivity Effects 0.000 claims description 67
- DMBHHRLKUKUOEG-UHFFFAOYSA-N diphenylamine Chemical compound C=1C=CC=CC=1NC1=CC=CC=C1 DMBHHRLKUKUOEG-UHFFFAOYSA-N 0.000 claims description 57
- 230000001360 synchronised effect Effects 0.000 claims description 17
- 238000004364 calculation method Methods 0.000 claims description 8
- 239000012634 fragment Substances 0.000 claims description 8
- 101100268668 Caenorhabditis elegans acc-2 gene Proteins 0.000 claims description 3
- 230000001427 coherent effect Effects 0.000 abstract 2
- 238000012937 correction Methods 0.000 description 12
- 238000001514 detection method Methods 0.000 description 6
- 230000008569 process Effects 0.000 description 6
- 230000008901 benefit Effects 0.000 description 4
- 230000005540 biological transmission Effects 0.000 description 4
- 230000000694 effects Effects 0.000 description 3
- 230000007257 malfunction Effects 0.000 description 3
- 239000000203 mixture Substances 0.000 description 3
- 241001269238 Data Species 0.000 description 2
- 238000004891 communication Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 230000036541 health Effects 0.000 description 2
- 230000013011 mating Effects 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 230000002123 temporal effect Effects 0.000 description 2
- 238000012546 transfer Methods 0.000 description 2
- 108010076504 Protein Sorting Signals Proteins 0.000 description 1
- 230000000903 blocking effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- YGANSGVIUGARFR-UHFFFAOYSA-N dipotassium dioxosilane oxo(oxoalumanyloxy)alumane oxygen(2-) Chemical compound [O--].[K+].[K+].O=[Si]=O.O=[Al]O[Al]=O YGANSGVIUGARFR-UHFFFAOYSA-N 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 229910052627 muscovite Inorganic materials 0.000 description 1
- 238000004064 recycling Methods 0.000 description 1
- 238000010561 standard procedure Methods 0.000 description 1
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/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/13—Receivers
- G01S19/24—Acquisition or tracking or demodulation of signals transmitted by the system
- G01S19/246—Acquisition or tracking or demodulation of signals transmitted by the system involving long acquisition integration times, extended snapshots of signals or methods specifically directed towards weak signal acquisition
-
- G—PHYSICS
- G04—HOROLOGY
- G04R—RADIO-CONTROLLED TIME-PIECES
- G04R20/00—Setting the time according to the time information carried or implied by the radio signal
- G04R20/02—Setting 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
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)
Abstract
本发明公开了一种高灵敏度北斗辅助授时装置和授时接收机及授时方法,在恶劣信号条件下,即便只跟踪到一颗卫星,在不能完成定位解算的情况下,依靠从接收到的导航信号中提取的周内时信息,仍然能够完成对北斗***时间的同步,提供了高灵敏度的授时服务。高灵敏度北斗辅助授时装置包括顺次连接的帧头搜索相关器、帧头搜索门限判决模块、周内时提取相关器以及周内时门限判决模块,帧头搜索门限判决模块的输出端和周内时门限判决模块的输出端还和帧头搜索相关器的输入端连接;帧头搜索相关器包括相互连接的帧头相参积累模块和帧头非相参积累模块,周内时提取相关器包括相互连接的周内时相参积累模块和周内时非相参积累模块。
Description
技术领域
本发明涉及卫星导航授时技术领域,特别是涉及一种高灵敏度北斗辅助授时装置和内嵌有所述高灵敏度辅助授时装置的高灵敏度北斗导航授时接收机及授时方法。
背景技术
授时的基本含义就是通过标准或者定制的接口和协议,为其他设备或***提供时间信息。最初的授时基本都依靠对太阳和恒星等天体规则运动的观测来实现,称为天文授时。天文测时所依赖的是地球自转,而地球自转的不均匀性使得天文方法所得到的时间(世界时)精度无法满足二十世纪中叶社会经济各方面的需求。一种更为精确和稳定的时间标准应运而生,这就是“原子钟”。目前世界各国都采用原子钟来产生和保持标准时间,这就是“时间基准”,然后,通过各种手段和媒介将时间信号送达用户,这些手段包括:短波、长波、电话网、互联网、卫星等。这一整个工序,就称为“授时***”。
短波授时的基本方法是由无线电台发播时间信号(简称时号),用户用无线电接收机接收时号,然后进行本地对时。长波授时利用长波(低频)进行时间频率传递与校准,是一种覆盖能力比短波强,校准的准确度更高的授时方法。网络授时和电话授时,采用用户询问方式向用户提供标准时间信号。卫星授时可以实现发播信号大面积的覆盖,而且比起前两种授时方法,它的精度更高。根据卫星在授时中所起的作用卫星授时分为主动式和中转式。主动式卫星有精密时钟,可发播标准时间信号;中转式仅转发由地面时间基准通过卫星地面站送来的标准时间信号。
卫星导航定位的基本原理是测量出已知位置的卫星到用户接收机之间的距离,然后综合多颗卫星的数据就可知道接收机的具***置。要达到这一目的,卫星的位置可以根据星载时钟所记录的时间在卫星星历中查出,而用户到卫星的距离则通过记录卫星信号传播到用户所经历的时间,再将其乘以光速得到。由于卫星轨道误差、大气层电离层的干扰,星钟与本地时钟的不同步等原因,这一距离并不是用户与卫星之间的真实距离,而是伪距。当卫星正常工作时,会不断地用1和0二进制码元组成的伪随机码发射导航电文。导航电文包括卫星星历、工作状况、时钟改正、电离层时延修正、大气折射修正等信息。用户接受到导航电文,提取出卫星时间并将其与自己的时钟做对比便可得知卫星与用户的距离,再利用导航电文中的卫星星历数据推算出卫星发射电文时所处位置,用户在大地坐标系中的位置速度等信息便可得知。
导航***卫星部分的作用就是不断地发射导航电文,卫星导航电文包含有***时间、时钟改正参数、电离层延迟模型参数、卫星星历及卫星健康状况信息等,这是为了给用户提供时间、位置坐标、速度等结果数据,而用于导航卫星信号处理的有关信息。然而,由于用户接收机使用的时钟与卫星星载时钟存在时间差,称为钟差,所以除了用户的三维坐标x、y、z外,还要引进一个Δt即卫星与接收机之间的时间差作为未知数,然后用4个方程将这4个未知数解出来。所以,如果想知道接收机所处的位置,至少要能接收到4颗卫星的信号。利用解算出的本地接收机时钟与卫星导航***时的时间差Δt,可以实现卫星导航接收机本地时钟与以原子钟为基础的导航***时间的同步。因此,每一个处于定位解算成功状态的卫星导航接收机都是一个高精度的时钟。
当前,卫星导航授时应用已发展成为全球性的高新技术产业,并且随着北斗卫星导航***日趋成熟,基于该***的授时应用也越来越得到人们的关注。北斗卫星导航***是中国自行研制的全球卫星定位与通信***,是继美国的全球定位***(GPS)和俄罗斯的GLONASS之后第三个成熟的卫星导航***,可在全球范围内全天候、全天时为各类用户提供高精度、高可靠的定位、导航、授时服务,并兼具短报文通信能力,截止2012年5月在轨卫星12颗,已经初步具备区域导航、定位和授时能力。北斗卫星导航***将在2020年形成全球覆盖能力。随着中国北斗卫星***的发展,覆盖范围的扩大,移动终端内部的时钟可以得到北斗卫星信号的有力辅助。但是由于目前大量移动终端用户都居住在城市环境中,移动终端在严重遮挡的环境下,可能无法保证对卫星连续的跟踪,任何一颗卫星的信号都可能时断时续,以致无法保证接收到完整的4颗卫星信号。在卫星信号强度较弱和遮挡的情况下,对北斗卫星导航***的授时服务灵敏度提出了更高的要求。
发明内容
本发明公开了一种高灵敏度北斗辅助授时装置和内嵌有所述高灵敏度辅助授时装置的高灵敏度北斗导航授时接收机及授时方法,在恶劣的信号遮挡和衰减条件下,即便只跟踪到一颗卫星,在不能完成位置、速度和时间的定位解算情况下,依靠从接收到的导航信号中提取的周内时信息,仍然能够完成对北斗***时间的同步,提供高灵敏度的授时服务。
本发明的技术方案是:
一种高灵敏度北斗辅助授时装置,其特征在于,包括顺次连接的帧头搜索相关器、帧头搜索门限判决模块、周内时提取相关器以及周内时门限判决模块,所述帧头搜索门限判决模块的输出端和周内时门限判决模块的输出端还和帧头搜索相关器的输入端连接;所述帧头搜索相关器包括相互连接的帧头相参积累模块和帧头非相参积累模块,将截取接收到的基带信号片段Sr(n+m)和本地复现信号Sl(m)先进行相参积累再进行非相参积累,将得到的非相参积累值通过帧头搜索门限判决模块进行计算和比较,确定北斗导航电文数据帧的子帧帧头所处位置;所述周内时提取相关器包括相互连接的周内时相参积累模块和周内时非相参积累模块,当接收到确定的帧头位置信号后,将截取的固定位置上的一段基带信号Sr(m)和本地复现信号Sk(m)先进行相参积累再进行非相参积累,将得到的非相参积累值通过周内时门限判决模块进行计算和比较,提取导航电文的周内时信息,输出1PPS秒脉冲信号,完成对北斗***时间的同步。
一种低成本型高灵敏度北斗导航授时接收机,其特征在于,包括与天线连接的射频模块、信号捕获模块、信号跟踪模块以及所述高灵敏度北斗辅助授时装置,所述射频模块分别连接信号捕获模块和信号跟踪模块的输入端,信号捕获模块的输出端连接信号跟踪模块的输入端,所述信号跟踪模块的输出端连接所述高灵敏度北斗辅助授时装置;当信号捕获模块捕获少于四颗卫星信号时,由所述高灵敏度北斗辅助授时装置从接收到的卫星信号中提取导航电文的周内时信息,输出1PPS秒脉冲信号,完成对北斗***时间的同步。
当信号捕获模块只捕获一颗卫星信号时,由所述高灵敏度北斗辅助授时装置从接收到的卫星信号中提取导航电文的周内时信息,输出1PPS秒脉冲信号,完成对北斗***时间的同步。
一种高灵敏度北斗导航授时接收机,其特征在于,包括与天线连接的射频模块、信号捕获模块、信号跟踪模块、定位授时模块以及所述高灵敏度北斗辅助授时装置和授时模式选择模块,所述射频模块分别连接信号捕获模块和信号跟踪模块的输入端,信号捕获模块的输出端连接信号跟踪模块的输入端,所述信号跟踪模块的输出端分别连接定位授时模块和所述高灵敏度北斗辅助授时装置,定位授时模块和高灵敏度北斗辅助授时装置的输出端分别连接授时模式选择模块;当信号捕获模块捕获四颗以上卫星信号时,由定位授时模块完成定位解算,通过授时模式选择模块输出定位结果及1PPS秒脉冲信号;当信号捕获模块捕获少于四颗卫星信号时,由所述高灵敏度北斗辅助授时装置从接收到的卫星信号中提取导航电文的周内时信息,输出1PPS秒脉冲信号,完成对北斗***时间的同步。
当信号捕获模块只捕获一颗卫星信号时,由所述高灵敏度北斗辅助授时装置从接收到的卫星信号中提取导航电文的周内时信息,输出1PPS秒脉冲信号,完成对北斗***时间的同步。
一种基于上述高灵敏度北斗辅助授时装置的授时方法,其特征在于,包括以下步骤:
1)在完成对北斗卫星信号的捕获和跟踪的基础上,根据捕获和跟踪到的至少一颗的北斗卫星信号,得到PRN码相位和多普勒频率信息;截取接收到的基带信号片段Sr(n+m),与本地复现信号Sl(m)进行点积运算,进行相参积累,再将所述点积运算结果取绝对值,进行非相参积累,所述非相参积累值
搜索非相参积累最大值所在的位置,找到北斗导航电文数据帧的子帧帧头所处位置;其中,L表示帧头的长度,n是搜索的起始位置;
2)根据帧头位置所确定的周内时字段的位置,截取固定位置上的一段基带信号Sr(m),与本地复现信号Sk(m)先进行相参积累再进行非相参积累,所述非相参积累值
搜索非相参积累最大值,提取导航电文的周内时信息,获取北斗***时间;其中,M表示存贮周内时间占用的比特数,Sk对应于第k个可能的周内时间序列。
所述步骤2)中,对于没有关于时间先验知识的首次提取,求取所有可能的本地复现信号Sk(m)对应的非相参积累值,搜索非相参积累最大值,完成首次周内时的提取;对于获得了先验知识的非首次提取,先估算非相参积累值的搜索范围,在估算的非相参积累值的搜索范围内,求取所有可能的本地复现信号Sk(m)对应的非相参积累值,搜索非相参积累最大值,完成周内时的提取。
所述步骤1)中,所述非相参积累值按下式计算:
所述步骤1)中,所述非相参积累值为多个子帧所对应的非相参积累值进行二次累加得到的最终统计量,Acc2nd(n)=Acc(n)+Acc(n+Psf),Psf是子帧长度或周期。
所述步骤2)中,所述非相参积累值按下式计算:
所述步骤2)中,所述非相参积累值为多个子帧所对应的非相参积累值进行二次累加得到的最终统计量, Psf是子帧长度或周期。
本发明的技术效果:
本发明提供了一种高灵敏度北斗辅助授时装置,信号捕获通道可以只捕获一颗或少于四颗导航卫星,在恶劣的信号遮挡或衰减条件下,即便只跟踪到一颗卫星,在不能完成位置、速度和时间的定位解算情况下,依靠从接收到的导航信号中提取的周内时信息,仍然能够完成对北斗***时间的同步,能够适应标准授时方法不能正常工作的卫星信号严重衰减、可见卫星数量不足等恶劣信号环境,提供高灵敏度的授时服务。
本发明的高灵敏度北斗辅助授时装置的授时方法,根据导航电文采用数据帧的结构传送,数据帧由多个子帧构成,每一子帧的开始部分是一个固定图案的帧头,随后的数据中包含了周内时信息的特点,在完成对北斗卫星信号的捕获和跟踪的基础上,根据捕获和跟踪到的至少一颗的北斗卫星信号,首先搜索帧头位置,帧头位置确定后,即可提取北斗周内时,输出1PPS秒脉冲信号,完成对北斗***时间的同步。不必进行PVT解算,仅仅提取周内时的方法突出的优点是与标准PVT解算相比大为缩减的运算负荷。
本发明还提供了内嵌有上述高灵敏度北斗辅助授时装置的高灵敏度北斗导航授时接收机,包括低成本型高灵敏度北斗导航授时接收机和标准型高灵敏度北斗导航授时接收机两种方案。低成本型高灵敏度北斗导航授时接收机包括与天线连接的射频模块、信号捕获模块、信号跟踪模块以及所述高灵敏度辅助授时装置,其优点是信号捕获通道可以只捕获一颗或少于四颗导航卫星,不必实现PVT解算,只要不是所有的卫星信号都被阻断,则授时服务就不会中断,卫星信号中的周内时仍然提供了相对精确的时间,这使得卫星导航接收机的时钟在无法定位的情况下依旧可以被视为一个可取的精确参考时间。标准型高灵敏度北斗导航授时接收机包括与天线连接的射频模块、信号捕获模块、信号跟踪模块、定位授时模块以及所述高灵敏度辅助授时装置和授时模式选择模块,其优点是能够在提供基于PVT解算的标准授时方法的基础上,还具有提供高灵敏度授时服务的功能,在卫星信号严重衰减、可见卫星数量不足等恶劣信号环境下,标准授时方法不能正常工作时,仍然能够获取北斗***时间,提供授时服务。
附图说明
图1为本发明高灵敏度北斗辅助授时装置的结构示意图。
图2为本发明低成本型高灵敏度北斗导航授时接收机的结构示意图。
图3为本发明标准型高灵敏度北斗导航授时接收机的结构示意图。
图4为本发明基于高灵敏度北斗辅助授时装置的授时方法流程图。
图5为北斗导航电文结构。
具体实施方式
以下结合附图对本发明的实施例做进一步的详细说明。
如图1所示,一种高灵敏度北斗辅助授时装置,包括顺次连接的帧头搜索相关器、帧头搜索门限判决模块、周内时提取相关器以及周内时门限判决模块,帧头搜索门限判决模块的输出端和周内时门限判决模块的输出端还和帧头搜索相关器的输入端连接;帧头搜索相关器包括相互连接的帧头相参积累模块和帧头非相参积累模块,帧头相参积累模块的输入端输入截取接收到的基带信号片段Sr(n+m)和本地复现信号Sl(m),先进行相参积累再进行非相参积累,将得到的非相参积累值输入帧头搜索门限判决模块,通过帧头搜索门限判决模块的计算和比较,确定是否找到北斗导航电文的子帧帧头的位置,若未找到帧头,重复上述过程,若找到帧头,将帧头位置信号输入给周内时提取相关器;周内时提取相关器包括相互连接的周内时相参积累模块和周内时非相参积累模块,当接收到确定的帧头位置信号后,将截取的固定位置上的一段基带信号Sr(m)和本地复现信号Sk(m)输入给周内时提取相关器,先进行相参积累再进行非相参积累,将得到的非相参积累值通过周内时门限判决模块的计算和比较,提取导航电文的周内时信息,输出1PPS秒脉冲信号,完成对北斗***时间的同步。其中,Sr(n+m)中的r表示received,即接收信号的意思,Sl(m)中的l表示local,即本地信号的意思。
如图2所示,描述了本发明内嵌有上述高灵敏度北斗辅助授时装置的低成本型高灵敏度北斗导航授时接收机的基本组成、结构和信号流向。包括与天线连接的射频模块、信号捕获模块、信号跟踪模块以及所述高灵敏度北斗辅助授时装置,射频模块分别连接信号捕获模块和信号跟踪模块的输入端,信号捕获模块的输出端连接信号跟踪模块的输入端,信号跟踪模块的输出端连接所述高灵敏度北斗辅助授时装置;其中,天线主要负责信号的接收;射频模块负责信号的下变频,在当前的数字化接收机中还包括A/D转换,它是所有后端处理的基础,其信号的品质直接影响接收机的性能;信号捕获模块和信号跟踪模块即为接收机的基带信号处理部分,主要负责信号的捕获、跟踪与锁定,包括各种原始数据和测量数据的输出;当信号捕获模块只捕获一颗或少于四颗卫星信号时,由所述高灵敏度北斗辅助授时装置从接收到的卫星信号中提取导航电文的周内时信息,输出1PPS秒脉冲信号,完成对北斗***时间的同步。
在低成本实现方案中,可以选用较低质量的时钟源。其缺点是可能延长卫星信号捕获任务完成所需的时间。假如对于导航卫星辅助时钟***的初次定位时间,即TTFF指标不做较高要求,则可以考虑低成本时钟方案。从获取周内时的角度出发,信号捕获通道可以只捕获一颗(或少于四颗)导航卫星,这是因为不必实现PVT解算,也就无需捕获多颗卫星的信号。因此,仅仅提取周内时的设计突出的优点是与标准PVT解算相比大为缩减的运算负荷。
如图3所示,描述了本发明内嵌有上述高灵敏度北斗辅助授时装置的标准型高灵敏度北斗导航授时接收机的基本组成、结构和信号流向。包括与天线连接的射频模块、信号捕获模块、信号跟踪模块、定位授时模块以及所述高灵敏度北斗辅助授时装置和授时模式选择模块,射频模块分别连接信号捕获模块和信号跟踪模块的输入端,信号捕获模块的输出端连接信号跟踪模块的输入端,信号跟踪模块的输出端分别连接定位授时模块和所述高灵敏度北斗辅助授时装置,定位授时模块和高灵敏度北斗辅助授时装置的输出端分别连接授时模式选择模块;当信号捕获模块捕获四颗或四颗以上卫星信号时,由定位授时模块完成定位解算,通过授时模式选择模块输出定位结果及1PPS信号;当信号捕获模块只捕获一颗或少于四颗卫星信号时,由所述高灵敏度北斗辅助授时装置从接收到的卫星信号中提取导航电文的周内时信息,输出1PPS秒脉冲信号,完成对北斗***时间的同步。
一种基于所述高灵敏度辅助授时装置的授时方法,包括以下步骤:
1)在完成对北斗卫星信号的捕获和跟踪的基础上,根据捕获和跟踪到的至少一颗的北斗卫星信号,得到PRN码相位和多普勒频率信息;截取接收到的基带信号片段Sr(n+m),与本地复现信号Sl(m)进行点积运算,进行相参积累,再将所述点积运算结果取绝对值,进行非相参积累,所述非相参积累值
搜索非相参积累最大值所在的位置,找到北斗导航电文数据帧的子帧帧头所处位置;其中,L表示帧头的长度,n是搜索的起始位置;
2)根据帧头位置所确定的周内时字段的位置,截取的固定位置上的一段基带信号Sr(m),与本地复现信号Sk(m)先进行相参积累再进行非相参积累,所述非相参积累值
搜索非相参积累最大值,提取导航电文的周内时信息,获取北斗***时间,其中,M表示存贮周内时间占用的比特数,Sk对应于第k个可能的周内时间序列。
其中,步骤2)中,对于没有关于时间先验知识的首次提取,求取所有可能的本地复现信号Sk(m)对应的非相参积累值,搜索非相参积累最大值,完成首次周内时的提取;对于获得了先验知识的非首次提取,先估算非相参积累值的搜索范围,在估算的非相参积累值的搜索范围内,求取所有可能的本地复现信号Sk(m)对应的非相参积累值,搜索非相参积累最大值,完成周内时的提取。
如图4所示,为本发明基于高灵敏度北斗辅助授时装置的授时方法流程图。在完成对北斗卫星信号的捕获和跟踪的基础上,得到PRN码相位和多普勒频率信息,然后进行相参积累和非相参积累,并搜索非相参积累最大值所在的位置,确定帧头位置;帧头位置确定之后,即可以提取北斗周内时,提取周内时之前,须首先明确有无关于时间的先验知识,如果无时间先验知识,进行周内时全组合非相参积累值的计算及搜素;如果有时间先验知识,则先估算非相参积累值的搜索范围,尽可能减小运算量,在估算的非相参积累值的搜索范围内,进行周内时非相参积累值的计算及搜素,搜索非相参积累最大值,完成北斗导航电文周内时的提取。
卫星导航电文包含有***时间、时钟改正参数、电离层延迟模型参数、卫星星历及卫星健康状况信息等。这是为了给用户提供时间、位置坐标、速度等结果数据,而用于导航卫星信号处理的有关信息。导航电文同样以二进制码的形式播送给用户,因此又叫数据码,或称D码。电文采用数据帧的结构传送,数据帧由多个子帧构成。每一子帧的开始部分是一个固定图案的帧头,随后的数据中包含了周内时间,即下一子帧的发射时刻。
固定图案的帧头会重复出现在各子帧的开始部分。在一个实现中,高灵敏度北斗辅助授时装置持续搜索每个子帧,使用非相参积累确定帧头的位置。更明确地说,各子帧的帧头将以非相参积累加以检测。在一个实现中,非相参积累将确定接收到的卫星信号(Sr(m))与本地序列(Sl(m)),如帧头,之间的相关函数。
接收到的卫星信号经相位同步和相参积累后以比特流的形式保存在接收机基带信号处理部分存贮器单元中,每一个比特可以表示为一个复数ArD(m)eiθ。其中,Ar是一个未知的幅度,θ是未知的相位,D(m)是电文序列。设定一个宽度为L的滑动窗截取接收到的卫星信号片断Sr(n+m),计算该片段与本地序列Sl(m)点积(相乘再求和)的绝对值,即可得到非相参积累值。具体的计算方法见方程(1):
L是信号片断或目标序列长度。假如待检测电文序列为0b0011110,则本地目标电文序列是复序列:
接收到的卫星信号也是一个复序列
Sr(m)=[A,A,-A,-A,-A,-A,A]eiθ+N ()
其中,N是复噪声序列。各序列中对应元素相乘(复数乘法),然后将乘积的实部和虚部分别相加,最后求取绝对值。显然,当目标序列与滑动窗截取到的输入信号序列匹配时,非相参积累取得最大值;而序列不匹配时得到的非相参积累值较小。在另一个实现中非相参积累也可以按下式计算:
两种计算方法的效果是类似的,但从实现的角度出发,第二种算法可以避免运算量较大的开方计算,更节省软/硬件开销。改变滑动窗位置/n的取值,计算和比较非相参积累值可以确定接收信号中帧头所处的位置。搜索范围,即n的变化范围,取决于北斗电文子帧的长度。在一个实现中,计算接收序列一个子帧内每一个元素处开始的子序列的非相参积累值,最大值所在的位置即为帧头。在弱信号或强噪声条件下,只对单个帧头的位置进行非相参检测可能无法得到正确的结果。解决的办法是在与第一个相关运算相隔固定长度(子帧长度)位置处计算额外的非相参积累值。图5展示了导航电文的结构以及前述多个非相参积累值的关系。在第一子帧进行相关运算得到第一非相参积累值,在第二子帧进行相关运算得到第二非相参积累值,第三子帧进行相关运算得到第三非相参积累值。图中只画出3个子帧,事实上用同样的思路可以得到第四非相参积累值、第五非相参积累值,以此类推。因为帧头会不断重复出现在后续子帧的相同位置上,将这些在对应位置处取得的非相参积累值相加得到最终的检测统计量。所述非相参积累值为多个子帧所对应的非相参积累值进行二次累加得到的最终统计量,Acc2nd(n)=Acc(n)+Acc(n+Psf),Psf是子帧长度或周期。形成最终检测统计量的对应非相参积累值越多,越有助于平滑检测统计量中的噪声,提高检测灵敏度。
周内时的提取遵循类似的思路,同样使用非相参积累来实现。但是具体的步骤有所不同,原因在于根据导航电文结构,通过前述帧头检测确定了帧头位置,周内时字段的位置也就明确了,因而不再使用滑动窗而是截取固定位置上的一段基带信号。不同于帧头,周内秒字段一直变化/递增,在相邻子帧不重复且一开始是未知的。
周内时提取分无先验知识/首次提取和有先验知识提取。首次提取由于没有先验知识,必须求取所有可能的二进制组合Sk(m)对应的非相参积累值,搜索最大值,认为最大非相参积累值出现在与真实时间匹配的二进制组合处。授时应用可能只关心以24小时为模的时分秒,从而减小搜索范围。搜索范围,即k的变化范围,取决于北斗电文周内时间编码组合的总数。在弱信号情况下,可以运用与帧头检测类似的非相参积累值再累积的办法抑制噪声,增大检测概率。所述非相参积累值为多个子帧所对应的非相参积累值进行二次累加得到的最终统计量, Psf是子帧长度或周期。考虑到相邻子帧周内时字段递增的特性,产生参与再累积的来自相邻周内时的非相参积累值的本地序列应该具有确定的对应关系。
一旦有一颗卫星成功完成了首次周内时提取,终端即取得了关于时间的先验知识。这一先验知识的辅助将大大减少待搜索二进制组合的数量,减轻对该卫星和其它卫星周内时提取的工作量。移动终端在严重遮挡的环境下,可能无法保证对卫星连续的跟踪,任何一颗卫星的信号都可能时断时续。但是只要不是所有的卫星信号都被阻断,则授时服务不会中断,因为终端可能交替提取不同卫星的周内时间。
与基于PVT解算的标准授时方法相比,本方法的钟差修正频度较低。标准方法的修正频度取决于PVT解算的频率,通常高于1Hz,某些北斗接收机可以达到20Hz或更高。本方法的钟差校正与子帧的发射频率相等,一个北斗信号子帧的持续时间通常在秒级,故而本方法的钟差修正频度小于1Hz。由于时钟随时间逐渐漂移,时间越长,钟漂越大。钟差修正频度越高,时钟的钟漂越小。因此,本方法的授时精度低于以PVT解算为基础的标准授时方法,但是能够适应标准授时方法不能正常工作的卫星信号严重衰减、可见卫星数量不足等恶劣信号环境。
应当指出,以上所述具体实施方式可以使本领域的技术人员更全面地理解本发明创造,但不以任何方式限制本发明创造。因此,尽管本说明书和实施例对本发明创造已进行了详细的说明,但是,本领域技术人员应当理解,仍然可以对本发明创造进行修改或者等同替换;而一切不脱离本发明创造的精神和范围的技术方案及其改进,其均涵盖在本发明创造专利的保护范围当中。
Claims (11)
1.一种高灵敏度北斗辅助授时装置,其特征在于,包括顺次连接的帧头搜索相关器、帧头搜索门限判决模块、周内时提取相关器以及周内时门限判决模块,所述帧头搜索门限判决模块的输出端和周内时门限判决模块的输出端还和帧头搜索相关器的输入端连接;所述帧头搜索相关器包括相互连接的帧头相参积累模块和帧头非相参积累模块,将截取接收到的基带信号片段Sr(n+m)和本地复现信号Sl(m)先进行相参积累再进行非相参积累,将得到的非相参积累值通过帧头搜索门限判决模块进行计算和比较,确定北斗导航电文数据帧的子帧帧头所处位置;所述周内时提取相关器包括相互连接的周内时相参积累模块和周内时非相参积累模块,当接收到确定的帧头位置信号后,将截取的固定位置上的一段基带信号Sr(m)和本地复现信号Sk(m)先进行相参积累再进行非相参积累,将得到的非相参积累值通过周内时门限判决模块进行计算和比较,提取导航电文的周内时信息,输出1PPS秒脉冲信号,完成对北斗***时间的同步。
2.一种低成本型高灵敏度北斗导航授时接收机,其特征在于,包括与天线连接的射频模块、信号捕获模块、信号跟踪模块以及根据权利要求1所述的高灵敏度北斗辅助授时装置,所述射频模块分别连接信号捕获模块和信号跟踪模块的输入端,信号捕获模块的输出端连接信号跟踪模块的输入端,所述信号跟踪模块的输出端连接所述高灵敏度北斗辅助授时装置;当信号捕获模块捕获少于四颗卫星信号时,由所述高灵敏度北斗辅助授时装置从接收到的卫星信号中提取导航电文的周内时信息,输出1PPS秒脉冲信号,完成对北斗***时间的同步。
3.根据权利要求2所述的低成本型高灵敏度北斗导航授时接收机,其特征在于,当信号捕获模块只捕获一颗卫星信号时,由所述高灵敏度北斗辅助授时装置从接收到的卫星信号中提取导航电文的周内时信息,输出1PPS秒脉冲信号,完成对北斗***时间的同步。
4.一种高灵敏度北斗导航授时接收机,其特征在于,包括与天线连接的射频模块、信号捕获模块、信号跟踪模块、定位授时模块以及根据权利要求1所述的高灵敏度北斗辅助授时装置和授时模式选择模块,所述射频模块分别连接信号捕获模块和信号跟踪模块的输入端,信号捕获模块的输出端连接信号跟踪模块的输入端,所述信号跟踪模块的输出端分别连接定位授时模块和所述高灵敏度北斗辅助授时装置,定位授时模块和高灵敏度北斗辅助授时装置的输出端分别连接授时模式选择模块;当信号捕获模块捕获四颗以上卫星信号时,由定位授时模块完成定位解算,通过授时模式选择模块输出定位结果及1PPS秒脉冲信号;当信号捕获模块捕获少于四颗卫星信号时,由所述高灵敏度北斗辅助授时装置从接收到的卫星信号中提取导航电文的周内时信息,输出1PPS秒脉冲信号,完成对北斗***时间的同步。
5.根据权利要求4所述的高灵敏度北斗导航授时接收机,其特征在于,当信号捕获模块只捕获一颗卫星信号时,由所述高灵敏度北斗辅助授时装置从接收到的卫星信号中提取导航电文的周内时信息,输出1PPS秒脉冲信号,完成对北斗***时间的同步。
6.一种基于权利要求1所述的高灵敏度北斗辅助授时装置的授时方法,其特征在于,包括以下步骤:
1)在完成对北斗卫星信号的捕获和跟踪的基础上,根据捕获和跟踪到的至少一颗的北斗卫星信号,得到PRN码相位和多普勒频率信息;截取接收到的基带信号片段Sr(n+m),与本地复现信号Sl(m)进行点积运算,进行相参积累,再将所述点积运算结果取绝对值,进行非相参积累,所述非相参积累值
搜索非相参积累最大值所在的位置,找到北斗导航电文数据帧的子帧帧头所处位置;其中,L表示帧头的长度,n是搜索的起始位置;
2)根据帧头位置所确定的周内时字段的位置,截取固定位置上的一段基带信号Sr(m),与本地复现信号Sk(m)先进行相参积累再进行非相参积累,所述非相参积累值
搜索非相参积累最大值,提取导航电文的周内时信息,获取北斗***时间;其中,M表示存贮周内时间占用的比特数,Sk对应于第k个可能的周内时间序列。
7.根据权利要求6所述的授时方法,其特征在于,所述步骤2)中,对于没有关于时间先验知识的首次提取,求取所有可能的本地复现信号Sk(m)对应的非相参积累值,搜索非相参积累最大值,完成首次周内时的提取;对于获得了先验知识的非首次提取,先估算非相参积累值的搜索范围,在估算的非相参积累值的搜索范围内,求取所有可能的本地复现信号Sk(m)对应的非相参积累值,搜索非相参积累最大值,完成周内时的提取。
8.根据权利要求6所述的授时方法,其特征在于,所述步骤1)中,用 代替 计算所述非相参积累值。
9.根据权利要求6或8所述的授时方法,其特征在于,在弱信号或强噪声条件下,将所述步骤1)中得到的多个子帧所对应的非相参积累值进行二次累加,得到非相参积累值的最终统计量,Acc2nd(n)=Acc(n)+Acc(n+Psf),Psf是子帧长度或周期。
10.根据权利要求6所述的授时方法,其特征在于,所述步骤2)中,用 代替 计算所述非相参积累值。
11.根据权利要求10所述的授时方法,其特征在于,在弱信号条件下,将所述步骤2)中得到的多个子帧所对应的非相参积累值进行二次累加,得到非相参积累值的最终统计量, Psf是子帧长度或周期。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210543822.5A CN103033828B (zh) | 2012-12-14 | 2012-12-14 | 一种高灵敏度北斗辅助授时装置和授时接收机及授时方法 |
PCT/CN2013/081844 WO2014089988A1 (zh) | 2012-12-14 | 2013-08-20 | 一种高灵敏度北斗辅助授时装置和授时接收机及授时方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210543822.5A CN103033828B (zh) | 2012-12-14 | 2012-12-14 | 一种高灵敏度北斗辅助授时装置和授时接收机及授时方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103033828A CN103033828A (zh) | 2013-04-10 |
CN103033828B true CN103033828B (zh) | 2015-02-04 |
Family
ID=48020898
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210543822.5A Expired - Fee Related CN103033828B (zh) | 2012-12-14 | 2012-12-14 | 一种高灵敏度北斗辅助授时装置和授时接收机及授时方法 |
Country Status (2)
Country | Link |
---|---|
CN (1) | CN103033828B (zh) |
WO (1) | WO2014089988A1 (zh) |
Families Citing this family (33)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103033828B (zh) * | 2012-12-14 | 2015-02-04 | 北京东方联星科技有限公司 | 一种高灵敏度北斗辅助授时装置和授时接收机及授时方法 |
CN105388488A (zh) * | 2015-11-27 | 2016-03-09 | 江苏艾倍科科技股份有限公司 | 一种北斗自行车 |
CN106338739A (zh) * | 2016-08-31 | 2017-01-18 | 江苏艾倍科科技股份有限公司 | 一种北斗飞行器 |
CN107171761B (zh) * | 2017-04-19 | 2018-11-16 | 西安电子科技大学 | 利用短波授时信号校准本地时钟的方法 |
CN109001970B (zh) * | 2017-06-07 | 2021-09-24 | 精工爱普生株式会社 | 计时装置、电子设备以及移动体 |
CN107396300A (zh) * | 2017-06-09 | 2017-11-24 | 武汉米风通信技术有限公司 | 基于北斗地基增强***和物联网的基站定位方法 |
CN107357160A (zh) * | 2017-07-13 | 2017-11-17 | 星展测控科技股份有限公司 | 一种卫星智能时钟 |
CN107515407A (zh) * | 2017-09-14 | 2017-12-26 | 北京航宇创通技术有限公司 | 定位授时卫星通信终端及其控制方法 |
CN107450083A (zh) * | 2017-09-20 | 2017-12-08 | 西安康帕斯电子科技有限公司 | 北斗中心型指挥机 |
CN109143283A (zh) * | 2018-02-24 | 2019-01-04 | 上海华测导航技术股份有限公司 | 一种为精密授时提供高精度pps的卫星接收机的设计方法 |
CN110662174A (zh) * | 2018-06-30 | 2020-01-07 | 北京华泰安信科技有限公司 | 基于北斗短信的指挥代码传输实现方法 |
CN109639382A (zh) * | 2018-12-20 | 2019-04-16 | 郑州轻工业学院 | 一种获取高精度定时信号的*** |
CN110308466A (zh) * | 2019-06-21 | 2019-10-08 | 南京理工大学 | 基于Zynq-7020的微型双模接收机及其导航方法 |
CN110244324A (zh) * | 2019-06-30 | 2019-09-17 | 西安寰宇赛特信息科技有限公司 | 一种智能超低功耗北斗信标示位一体化终端 |
CN110515109B (zh) * | 2019-08-16 | 2023-05-23 | 中国航空工业集团公司西安飞行自动控制研究所 | 一种基于多信息源融合自主pnt时间的方法及装置 |
CN111078393B (zh) * | 2019-11-07 | 2022-11-15 | 河北德冠隆电子科技有限公司 | 一种交互型边缘计算服务器及其控制方法 |
CN111176107B (zh) * | 2020-01-19 | 2021-10-22 | 成都天奥电子股份有限公司 | 一种基于北斗三号***的快速授时手表 |
CN112511254B (zh) * | 2020-11-18 | 2022-09-06 | 西安航光卫星测控技术有限公司 | 基于北斗定位授时和lora通信的移动母钟及授时方法 |
CN112600608B (zh) * | 2020-11-25 | 2022-07-22 | 四川安迪科技实业有限公司 | 一种多帧叠加的卫星信号捕获方法及装置 |
US11683771B2 (en) * | 2020-11-30 | 2023-06-20 | Viettel Group | Method and apparatus for data frame synchronization of 5G base station |
CN112925033B (zh) * | 2021-01-23 | 2022-12-13 | 中国科学院国家授时中心 | 一种长波授时等效大地电导率数据的差分测算方法 |
CN113189626B (zh) * | 2021-03-05 | 2023-06-27 | 长沙学院 | 一种载体旋转运动下北斗卫星导航信号的快速位同步方法 |
CN113365232B (zh) * | 2021-05-17 | 2022-03-08 | 福建吉星智能科技股份有限公司 | 一种北斗通讯应用***短报文通讯的终端控制方法 |
CN113541914B (zh) * | 2021-07-05 | 2022-05-24 | 中国人民解放军国防科技大学 | 基于关联规则的cnav-2电文帧同步方法 |
CN113960919A (zh) * | 2021-10-15 | 2022-01-21 | 合肥移瑞通信技术有限公司 | Gnss授时精度测试***及方法 |
CN114095070B (zh) * | 2021-10-29 | 2023-12-29 | 中国电子科技集团公司第五十四研究所 | 一种基于北斗卫星导航的箭体信息回传装置 |
CN114268582B (zh) * | 2021-11-24 | 2024-04-05 | 上海机电工程研究所 | 弹动转发***及装置 |
CN114637033B (zh) * | 2022-03-11 | 2024-04-05 | 重庆市计量质量检测研究院 | 一种基于北斗的远程实时时间校准方法 |
CN114706105B (zh) * | 2022-06-07 | 2022-08-23 | 尚禹河北电子科技股份有限公司 | 一种长波引导北斗授权信号捕获的方法、装置及*** |
CN115616890A (zh) * | 2022-09-28 | 2023-01-17 | 广西自贸区见炬科技有限公司 | 一种热电***钟差预报误差修正*** |
CN115865301B (zh) * | 2022-10-08 | 2023-05-26 | 中国科学院空天信息创新研究院 | 一种卫星通信中抗大码偏、低功耗的定时同步方法 |
CN115932921B (zh) * | 2022-12-22 | 2023-07-04 | 中国人民解放军军事科学院***工程研究院 | 一种多卫星多载波的非相参联合累积处理方法 |
CN118112624A (zh) * | 2024-04-25 | 2024-05-31 | 江西省军民融合研究院 | 一种基于北斗定位的无人机位置通信方法及*** |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN201540463U (zh) * | 2009-11-13 | 2010-08-04 | 湖南创越电子科技有限公司 | 北斗一号无源定位授时终端 |
CN101969690A (zh) * | 2010-10-09 | 2011-02-09 | 东莞市泰斗微电子科技有限公司 | 一种利用导航卫星进行授时的方法及装置 |
CN102073056A (zh) * | 2009-11-20 | 2011-05-25 | 郑州威科姆科技股份有限公司 | 北斗/gps双***授时接收机 |
CN102129218A (zh) * | 2010-12-14 | 2011-07-20 | 中电普瑞科技有限公司 | 一种使用北斗卫星的精确授时方法 |
CN102354101A (zh) * | 2011-08-17 | 2012-02-15 | 东莞市泰斗微电子科技有限公司 | 一种利用导航卫星的授时方法及装置 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102062862B (zh) * | 2010-12-14 | 2012-11-14 | 东莞市泰斗微电子科技有限公司 | 一种导航卫星混合定位或授时的模式选择方法 |
CN202394056U (zh) * | 2011-10-19 | 2012-08-22 | 成都金本华科技有限公司 | 一种北斗一号卫星单向授时*** |
CN103033828B (zh) * | 2012-12-14 | 2015-02-04 | 北京东方联星科技有限公司 | 一种高灵敏度北斗辅助授时装置和授时接收机及授时方法 |
-
2012
- 2012-12-14 CN CN201210543822.5A patent/CN103033828B/zh not_active Expired - Fee Related
-
2013
- 2013-08-20 WO PCT/CN2013/081844 patent/WO2014089988A1/zh active Application Filing
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN201540463U (zh) * | 2009-11-13 | 2010-08-04 | 湖南创越电子科技有限公司 | 北斗一号无源定位授时终端 |
CN102073056A (zh) * | 2009-11-20 | 2011-05-25 | 郑州威科姆科技股份有限公司 | 北斗/gps双***授时接收机 |
CN101969690A (zh) * | 2010-10-09 | 2011-02-09 | 东莞市泰斗微电子科技有限公司 | 一种利用导航卫星进行授时的方法及装置 |
CN102129218A (zh) * | 2010-12-14 | 2011-07-20 | 中电普瑞科技有限公司 | 一种使用北斗卫星的精确授时方法 |
CN102354101A (zh) * | 2011-08-17 | 2012-02-15 | 东莞市泰斗微电子科技有限公司 | 一种利用导航卫星的授时方法及装置 |
Non-Patent Citations (2)
Title |
---|
北斗卫星导航***授时应用;陈洪卿 等;《数字通信世界》;20110630(第6期);54-58 * |
基于FPGA和DSP的北斗1PPS授时***的设计与实现;陈芬芬;《黑龙江科技信息》;20081231(第36期);77 * |
Also Published As
Publication number | Publication date |
---|---|
WO2014089988A1 (zh) | 2014-06-19 |
CN103033828A (zh) | 2013-04-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103033828B (zh) | 一种高灵敏度北斗辅助授时装置和授时接收机及授时方法 | |
CN101646955B (zh) | 时延测量 | |
CN100397094C (zh) | 卫星定位***接收机中的时间确定及其方法 | |
CN101821643B (zh) | 全球导航卫星***的信号捕获方法及接收器 | |
US8391341B2 (en) | Process for sub-microsecond time transfer using weak GPS/GNSS signals | |
CN101755223B (zh) | 无定时信息的导航定位 | |
CN101458321B (zh) | 基于pc机的gps/caps双模软件接收机 | |
CN102435999B (zh) | Gps接收机基带模块和gps信号捕获及跟踪方法 | |
EP2056121B1 (en) | Mobile unit positioning device | |
CN1987515B (zh) | 实现在微弱信号环境下gps定位的方法和*** | |
RU2625819C1 (ru) | Способ вспомогательного захвата в приемнике объединенной навигационной системы | |
CN113009524B (zh) | 长时间相干积分捕获的导航电文比特翻转估计方法及*** | |
CN102183770A (zh) | 一种抗多径干扰的gps伪随机码跟踪环路及其抗多径干扰方法 | |
CN105492927A (zh) | 云卸载的全球卫星定位 | |
US8472504B2 (en) | Method for optimizing an acquisition of a spread-spectrum signal from a satellite by a mobile receiver | |
CN103176189A (zh) | 高灵敏度卫星导航接收机的远近效应抑制器及其方法 | |
CN101937089B (zh) | 用于在定位***中弱数据位同步的方法和装置 | |
CN111458730B (zh) | 基于多普勒残差估计的gnss载波跟踪方法及接收机 | |
CN2854622Y (zh) | 实现在微弱信号环境下gps定位的*** | |
CN102680985B (zh) | 在全球导航卫星***中抑制多径信号的鉴相方法 | |
US10838072B2 (en) | Fast fix using a perfect reference in a satellite-based positioning system | |
CN1690726B (zh) | 按比例缩减装置及方法,全球定位***同步方法及接收器 | |
CN104237912A (zh) | 导航比特同步方法及检查导航比特同步的方法 | |
CN109541651A (zh) | 一种gpsl1l2双频信号联合捕获方法 | |
Akopian et al. | A network aided iterated LS method for GPS positioning and time recovery without navigation message decoding |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20150204 |