CN103926604B - 基于重叠差分循环相干积分的弱信号捕获方法 - Google Patents
基于重叠差分循环相干积分的弱信号捕获方法 Download PDFInfo
- Publication number
- CN103926604B CN103926604B CN201410153065.XA CN201410153065A CN103926604B CN 103926604 B CN103926604 B CN 103926604B CN 201410153065 A CN201410153065 A CN 201410153065A CN 103926604 B CN103926604 B CN 103926604B
- Authority
- CN
- China
- Prior art keywords
- data
- signal
- carrier
- carry out
- data block
- 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
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
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
本发明提供的是一种基于重叠差分循环相干积分的弱信号捕获方法。步骤1:取一段20Mms+10ms长度的中频信号数据,进行去载波处理,设定多普勒频移范围、搜索频率步进、检测阈值;步骤2:将去除载波之后的数据以L=10ms的长度进行数据分块处理,将每一个数据块进行自身叠加处理,并进行FFT变换;步骤3:将本地产生的C/A码进行FFT变换并取共轭,然后与步骤2中最终获得的每个结果分别相乘并进行IFFT处理;步骤4:将步骤3处理后的每个数据块与和它有数据重叠的相邻数据块进行复数共轭相乘,然后再进行分组相关求和,判断信号是否捕获成功。本发明的应用范围是全球定位***中单频民用软件接收机对于弱信号的捕获。
Description
技术领域
本发明涉及的是全球卫星导航***(GNSS——GlobalNavigationSatelliteSystem)中接收机的捕获处理方法,具体的应用范围是全球定位***(GPS——GlobalPositioningSystem)中单频民用软件接收机对于弱信号的捕获方法。
背景技术
全球卫星定位***简称GPS,能够在全球范围内为用户提供精确的位置、速度、时间等信息,并且已经得到了广泛的应用。GPS***是发展最早,也是目前发展最为成熟的导航定位***。由于GPS的优越性,使其受到了人们的普遍关注,它的应用范围非常广泛,包括可以为地球表面、近地表以及地球外太空任意地点的用户提供导航服务,但还是存在着局限性:比如在室内、森林、城市等复杂环境中,信号将会被衰减至很弱的情况下,GPS就不能很好的进行工作。如果使GPS接收机能够在信号较微弱的条件下进行定位工作,那么,就需要对GPS接收机最先启动的工作部分,即捕获功能模块的性能指标提出更高的要求,需要其有很好的对弱信号的捕获特性,即能够捕获到功率足够低的信号。当前GPS接收机所采用的捕获技术已经是难以满足用户的需求。如果不能成功的捕获到GPS信号,那么接收机后续的跟踪和解算等功能环节的进行就无从谈起,所以,为保证接收机的灵敏度,对微弱信号的捕获算法就成为了GPS接收机信号处理的关键部分。
GPS信号的捕获过程一般是通过对卫星信号的载波频率和码相位这两个参量值进行扫描式二维搜索来完成的。具体的过程就是在输入信号的多普勒频率和C/A码码相位的二维平面上寻找最大相关峰峰值,从而寻找到输入信号的多普勒频率和C/A码码相位所在位置,为后续信号的跟踪提供初始信息。目前典型的捕获方法有串行捕获算法和基于快速傅里叶变换(FFT)的循环相关算法。串行捕获算法的优点是结构比较简单且易于实现,但它需要对多组的码相位和多普勒频率的组合进行运算,运算量较大,处理速度缓慢。另一种比较典型的信号捕获方法是基于快速傅里叶变换的循环相关算法,该方法易于实现且效率高,通过按照多普勒频率范围内的每一种多普勒频移,对所有不确定的码时延进行快速傅里叶变换,计算相干积分值。这样二维搜索空间被简化成一维搜索空间,即只剩下对多普勒频移这一维的搜索,从而减少了计算量,节省了捕获所花费的时间,提高了捕获速度。但是该方法也是有缺陷的,那就是捕获的精度与接收数据所选取的长度密切相关,选取的数据长度越长,精度越高,但是随着数据长度的不断增长,处理过程的运算量和导航数据比特翻转的概率都会随之增大。
通常对于GPS弱信号的捕获处理,为了能够获得一个较高的载噪比,传统的解决方法是对接收数据进行较长时间的相干积分或者是进行多次的非相干积分累加。这两种方法都是存在着一定的缺陷的:对于相干积分来说,由于数据比特存在着跳变的限制,相干积分的长度一般是不能超过一个数据位即20ms,否则就有可能因为数据位的跳变,使得相干积分值变小甚至是出现抵消现象,所以并不能够无限制的增加相干积分的时间;由于非相干积分法虽然可以利用平方运算克服数据位跳变的影响,但是却也不可避免的引起了平方损耗,使信噪比的增益受到削弱,所以也不能无限制的进行非相干积分。但对于弱信号的捕获已经成为GPS接收机的关键技术,因此,GPS弱信号捕获成为当前的研究热点。
发明内容
本发明发明的目的在于提供一种可以有效捕获到弱信号,并且计算量不大的基于重叠差分循环相干积分的弱信号捕获方法。
本发明的目的是这样实现的:
步骤1:取一段20Mms+10ms长度的中频信号数据,M=1,2,3…,进行去载波处理,使其与本地复现的载波复正弦信号进行混频处理,设定合适的多普勒频移范围(ωDmin,ωDmax),以及搜索频率步进ΔωD,同时设定检测阈值γ;
步骤2:将去除载波之后的数据以L=10ms的长度进行数据分块处理,其中,第一个数据块与第二个数据块之间有5ms的数据重叠,第二个数据块与第三个数据块之间也有5ms的数据重叠,以此类推,共分为4M+1块的数据块,将每一个数据块进行自身叠加处理,并进行FFT变换;
步骤3:将本地产生的C/A码进行FFT变换并取共轭,然后与步骤2中最终获得的每个结果分别相乘并进行IFFT处理;
步骤4:将步骤3处理后的每个数据块与和它有数据重叠的相邻数据块进行复数共轭相乘,即进行重叠差分循环相关,然后再进行分组相关求和,对得到的4个结果进行判断,判断其大小,选取没有数据跳变的结果,获得峰值,将峰值与设定的门限值γ作比较,判断信号是否捕获成功,如果结果超出门限值,则捕获到信号,否则,没有捕获到信号,然后重新调制码相位和多普勒频率,返回步骤2,直到信号捕获到。
本发明的优点:
(1)本发明采用的是差分循环相干积分法。该方法将相邻的2个信号块进行共轭相乘,然后判断其中不含数据跳变的数据段,优点是能够较好的捕获到弱信号,使得接收机对于弱信号的适应性得到提高。
(2)本发明中对接收数据首先进行了叠加处理,减小了相干积分的运算量,从而提高了捕获速度。
(3)本发明采用了重叠差分的理念,可以实现较长时间的相干积分与差分相干运算的结合。
附图说明
图1是重叠差分分块示意图。
图2是基于傅里叶变换的循环相关算法示意图。
图3是传统弱信号捕获原理结构图。
图4是差分循环相干算法原理图。
具体实施方式
本发明描述的方法是一种适用于GPS弱信号的捕获方法,该发明采用了差分循环相干积分法,与一般的弱信号捕获方法相比,该方法的创新之处是首次采用了重叠差分循环相关积分,实现了较长时间的相干积分与差分相干运算的结合,本发明将相邻的2个数据块进行重叠共轭相乘,然后在4个结果中选取出不含数据跳变的数据段,使其与设定门限作比较,而且该发明对噪声的放大相对来说比较小,从而改善了信号的信噪比,并且与一般的方法相比计算量也不大。本发明方法如图4所示,步骤如下:
步骤1:取一段20Mms+10ms长度的中频信号数据作为接收数据,M=1,2,3…,选取的数据越长,捕获效果越好,将该数据信号进行去载波处理,使其与本地复现的载波复正弦信号进行混频处理,设定合适的多普勒频移范围(ωDmin,ωDmax),以及搜索频率步进ΔωD,同时设定检测阈值γ。
具体过程如下:
设接收信号的信号模型为:
式(1)中:s(k)是在采样时间tk时的接收信号,c(tk-ts)为初始码相位为ts的C/A码序列;为初始载波相位;频率ωIF是中频,ωD是载波信号的多普勒频移;n(k)是高斯白噪声;k=0,1,…N-1,为处理数据段的采样点数。
将接收信号s(k)与本地载波复正弦信号进行混频处理。滤除高频分量后得到:
式(2)中:为第i次载波频率搜索点与接收信号间的频率差值;是第i次的本地载波的多普勒频移;为对应的相位差值;当频率差值Δωd=0时,即本地复现的载波频率与接收信号的载波频率相等时,便可实现载波的完全剥离,最终会得到包含常数因子的接收到的伪码序列为:
步骤2:将去除载波之后的数据以L=10ms的长度进行数据分块处理,其中,第一个数据块与第二个数据块之间有5ms的数据重叠,第二个数据块与第三个数据块之间也有5ms的数据重叠,以此类推,该段信号共可分为4M+1块的数据块,将每一个数据块进行自身叠加处理,并进行FFT变换。
具体过程如下:
将去除载波之后的数据分块,如图1所示,每个块与每个块之间有5ms数据的重叠,依次记为其中m=1,2,3,…,M。例如,第一个20ms的数据长度(并不一定刚好是一个数据位)被分割成4块,依次记为需要用到下一个20ms数据中的前5ms,第二个数据块分割后的数据依次记为依次类推,最后一个块的数据表示为
将这些数据块进行自身叠加处理,即将每一个10ms长度的数据块叠加成1ms的长度,分别记为m=1,2…M。这样再对该1ms数据进行处理时实际上是对10ms的数据进行了处理,将数据进行FFT处理。
步骤3:将本地产生的C/A码进行FFT变换并取共轭,然后与步骤2中的每一个结果分别相乘并进行IFFT处理。
C=(FFT(c))*(4)
本地产生的C/A码进行FFT变换并取共轭,如公式(4)所示,其中c表示C/A码序列,C表示C/A码序列进行FFT变换之后取共轭,*表示共轭。然后将步骤2中的数据与该变换相乘并进行IFFT处理。
步骤4:将步骤3处理后的每个数据块与和它有数据重叠的相邻数据块进行复数共轭相乘,即进行重叠差分循环相关,然后再进行分组相关求和,对得到的4个结果进行判断,判断其大小,选取没有数据跳变的结果,获得峰值,将峰值与设定的门限值γ作比较,判断信号是否捕获成功,如果结果超出门限值,则捕获到信号,否则,没有捕获到信号,然后重新调制码相位和多普勒频率,返回步骤2,重复过程,直到信号捕获到。
具体过程如下:
假设在20Mms+6ms位置处发生数据跳变,M=1,2,3…,由于数据是进行10ms分块,相邻两个数据块之间有5ms的数据重合,所以含有数据跳变的块为Z1m和Z2m这两项。
通过观察上述4个等式,可以看出,等式(6)、(7)、(8)、(9)中不包含Z1m和Z2m这两项的只有等式(8),由于导航数据比特跳变出现在这些分块中时,那么该分块的前后差分相关将会出现较大的削弱,导致最终的累积结果出现较大的衰减,而不存在跳变的数据块,其差分相关累积结果将有比较大的累计值。通过比较这4个式子峰值的大小,就可以确定数据跳变位于哪组分块中。将不存在跳变的分块组所对应的差分相干累积结果的峰值与预先设定的检测阈值γ作比较,判断信号是否捕获成功。将得到的结果与检测阈值γ作比较,如果结果Y大于检测阈值γ则表明捕获成功;若小于检测阈值γ,则表明当前的卫星信号并没有捕获成功,调节本振和本地C/A码,重复以上的搜索过程,直到卫星信号捕获成功。
Claims (4)
1.一种基于重叠差分循环相干积分的弱信号捕获方法,其特征是:
步骤1:取一段20Mms+10ms长度的中频信号数据,M=1,2,3…,进行去载波处理,使其与本地复现的载波复正弦信号进行混频处理,设定多普勒频移范围(ωDmin,ωDmax),以及搜索频率步进ΔωD,同时设定检测阈值γ;
步骤2:将去除载波之后的数据以L=10ms的长度进行数据分块处理,其中,第一个数据块与第二个数据块之间有5ms的数据重叠,第二个数据块与第三个数据块之间也有5ms的数据重叠,以此类推,共分为4M+1块的数据块,将每一个数据块进行自身叠加处理,并进行FFT变换;
步骤3:将本地产生的C/A码进行FFT变换并取共轭,然后与步骤2中最终获得的每个结果分别相乘并进行IFFT处理;
步骤4:将步骤3处理后的每个数据块与和它有数据重叠的相邻数据块进行复数共轭相乘,即进行重叠差分循环相关,然后再进行分组相关求和,对得到的4个结果进行判断,判断其大小,选取没有数据跳变的结果,获得峰值,将峰值与设定的检测阈值γ作比较,判断信号是否捕获成功,如果结果超出检测阈值γ,则捕获到信号,否则,没有捕获到信号,然后重新调制码相位和多普勒频率,返回步骤2,直到信号捕获到;
所述步骤1具体包括:
设接收信号的信号模型为:
其中:s(k)是在采样时间tk时的接收信号,c(tk-ts)为初始码相位为ts的C/A码序列;为初始载波相位;频率ωIF是中频,ωD是载波信号的多普勒频移;n(k)是高斯白噪声;k=0,1,…N-1,为处理数据段的采样点数;
将接收信号s(k)与本地载波复正弦信号进行混频处理,滤除高频分量后得到:
其中:为第i次的本地载波的多普勒频移与接收信号间的载波信号的多普勒频移差值;是第i次的本地载波的多普勒频移;为对应的相位差值;当Δωd=0时,即本地复现的载波频率与接收信号的载波频率相等时,便实现载波的完全剥离,最终会得到包含常数因子的接收到的伪码序列为:
2.根据权利要求1所述的基于重叠差分循环相干积分的弱信号捕获方法,其特征是所述步骤2具体包括:
将去除载波之后的数据分块,每个块与每个块之间有5ms数据的重叠,依次记为其中m=1,2,3,…,M,第一个20ms的数据长度被分割成4块,依次记为 需要用到下一个20ms数据中的前5ms,第二个数据块分割后的数据依次记为 依次类推,最后一个块的数据表示为
将这些数据块进行自身叠加处理,即将每一个10ms长度的数据块叠加成1ms的长度,分别记为m=1,2…M,再对该1ms数据进行处理时实际上是对10ms的数据进行了处理,将数据进行FFT处理,
3.根据权利要求2所述的基于重叠差分循环相干积分的弱信号捕获方法,其特征是所述步骤3具体包括:
本地产生的C/A码进行FFT变换并取共轭,C=(FFT(c))*,其中c表示C/A码序列,C表示C/A码序列进行FFT变换之后取共轭,*表示共轭,然后将FFT处理的数据与该变换共轭后结果相乘并进行IFFT处理,
4.根据权利要求3所述的基于重叠差分循环相干积分的弱信号捕获方法,其特征是所述步骤4具体包括:
设在20Mms+6ms位置处发生数据跳变,M=1,2,3…,由于数据是进行10ms分块,相邻两个数据块之间有5ms的数据重合,所以含有数据跳变的块为和这两项,
比较这4个式子峰值的大小,确定数据跳变位于哪组分块中,将不存在跳变的分块组所对应的差分相干累积结果的峰值与预先设定的检测阈值γ作比较,判断信号是否捕获成功,将得到的结果与检测阈值γ作比较,如果结果Y大于检测阈值γ则表明捕获成功;若小于检测阈值γ,则表明当前的卫星信号并没有捕获成功,调制码相位与多普勒频率,重复以上的搜索过程,直到卫星信号捕获成功。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410153065.XA CN103926604B (zh) | 2014-04-17 | 2014-04-17 | 基于重叠差分循环相干积分的弱信号捕获方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410153065.XA CN103926604B (zh) | 2014-04-17 | 2014-04-17 | 基于重叠差分循环相干积分的弱信号捕获方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103926604A CN103926604A (zh) | 2014-07-16 |
CN103926604B true CN103926604B (zh) | 2016-05-04 |
Family
ID=51144892
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410153065.XA Expired - Fee Related CN103926604B (zh) | 2014-04-17 | 2014-04-17 | 基于重叠差分循环相干积分的弱信号捕获方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103926604B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105974448A (zh) * | 2016-06-20 | 2016-09-28 | 奇瑞汽车股份有限公司 | 捕获卫星信号的方法和装置 |
CN107783077B (zh) * | 2016-08-25 | 2020-11-17 | 大连楼兰科技股份有限公司 | 过门限峰值点的处理的方法 |
CN109917429B (zh) * | 2019-03-07 | 2023-09-12 | 西安开阳微电子有限公司 | 一种b1c弱信号的捕获方法、装置及计算机存储介质 |
CN110456393A (zh) * | 2019-08-21 | 2019-11-15 | 四川航天***工程研究所 | 北斗弱信号快速捕获方法 |
CN114280641A (zh) * | 2021-12-28 | 2022-04-05 | 上海恒岳智能交通科技有限公司 | 一种北斗卫星接收信号的数据块预处理方法及*** |
CN114614860B (zh) * | 2022-02-17 | 2023-06-23 | 中国电子科技集团公司第十研究所 | 一种高动态非相干直扩信号差分捕获*** |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5886749A (en) * | 1996-12-13 | 1999-03-23 | Cable Television Laboratories, Inc. | Demodulation using a time domain guard interval with an overlapped transform |
CN101699774A (zh) * | 2009-11-10 | 2010-04-28 | 京信通信***(中国)有限公司 | 用于CDMA2000 1x***的下行同步PN序列捕获方法 |
US7860476B1 (en) * | 2006-12-14 | 2010-12-28 | Alfred E. Mann Foundation For Scientific Research | Narrowband interference excision in the external controller of an implanted microstimulator network |
CN102323601A (zh) * | 2011-05-27 | 2012-01-18 | 哈尔滨工程大学 | 一种gnss-boc调制信号的捕获方法 |
CN102680986A (zh) * | 2012-05-15 | 2012-09-19 | 哈尔滨工程大学 | 一种基于频率可调的gps接收机信号捕获反馈回路及捕获方法 |
CN103645483A (zh) * | 2013-12-09 | 2014-03-19 | 西安电子科技大学昆山创新研究院 | 一种弱信号环境下北斗信号捕获方法 |
-
2014
- 2014-04-17 CN CN201410153065.XA patent/CN103926604B/zh not_active Expired - Fee Related
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5886749A (en) * | 1996-12-13 | 1999-03-23 | Cable Television Laboratories, Inc. | Demodulation using a time domain guard interval with an overlapped transform |
US7860476B1 (en) * | 2006-12-14 | 2010-12-28 | Alfred E. Mann Foundation For Scientific Research | Narrowband interference excision in the external controller of an implanted microstimulator network |
CN101699774A (zh) * | 2009-11-10 | 2010-04-28 | 京信通信***(中国)有限公司 | 用于CDMA2000 1x***的下行同步PN序列捕获方法 |
CN102323601A (zh) * | 2011-05-27 | 2012-01-18 | 哈尔滨工程大学 | 一种gnss-boc调制信号的捕获方法 |
CN102680986A (zh) * | 2012-05-15 | 2012-09-19 | 哈尔滨工程大学 | 一种基于频率可调的gps接收机信号捕获反馈回路及捕获方法 |
CN103645483A (zh) * | 2013-12-09 | 2014-03-19 | 西安电子科技大学昆山创新研究院 | 一种弱信号环境下北斗信号捕获方法 |
Non-Patent Citations (1)
Title |
---|
基于FFT_iFFT的GPS信号捕获算法研究等;任记达;《现代导航》;20131231(第6期);论文第409页第1段及附图1 * |
Also Published As
Publication number | Publication date |
---|---|
CN103926604A (zh) | 2014-07-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103926604B (zh) | 基于重叠差分循环相干积分的弱信号捕获方法 | |
CN103645483B (zh) | 一种弱信号环境下北斗信号捕获方法 | |
CN103926605B (zh) | 基于差分循环相干积分的gps弱信号捕获方法 | |
CN106646546A (zh) | 一种卫星信号多维快速捕获方法及*** | |
CN102662183B (zh) | Gps信号捕获方法与*** | |
CN104181556B (zh) | 一种基于重叠差分循环相干积分的boc调制信号捕获方法 | |
CN104035109B (zh) | 基于重叠1/5比特差分循环相干积分的弱信号捕获方法 | |
CN102282773A (zh) | 用于捕捉扩频信号的方法和*** | |
CN105759288B (zh) | 基于特征序列的北斗b1i弱信号捕获方法 | |
CN107085222A (zh) | 一种 boc 信号捕获方法和卫星导航接收机 | |
CN105607096A (zh) | 一种双星时差频差定位方法和定位装置 | |
CN105717522A (zh) | “北斗二代”b1频段弱信号捕获方法 | |
CN105553506B (zh) | 一种长码扩频信号的快速捕获方法及装置 | |
CN104597435B (zh) | 基于修正频域补偿和分数阶傅里叶变换的多帧相参tbd方法 | |
CN104765050A (zh) | 一种新的北斗信号二次捕获算法 | |
CN105911569A (zh) | 基于码元序列估计的北斗卫星导航***b1频带弱d1信号捕获方法 | |
CN104898136A (zh) | 一种北斗二代卫星b1信号的捕获方法及其*** | |
CN103969664B (zh) | 一种基于差分循环相干积分的gps弱信号捕获方法 | |
CN106291610A (zh) | 一种用于gnss信号压缩捕获处理装置的压缩并行相关模块及其实现方法 | |
CN103941269A (zh) | 用于卫星导航***的pn码捕获方法 | |
CN102680988B (zh) | 高灵敏度导航卫星信号非线性捕获方法及装置 | |
CN103454651B (zh) | 一种gps l2c cm码快速捕获方法 | |
JP4805930B2 (ja) | 共通の復調テンプレートを使用する信号収集のための装置、方法およびコンピュータプログラム | |
CN104360360A (zh) | 一种gps卫星捕获判决的高效算法 | |
CN105656511A (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 | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20160504 |