CN104375139B - 一种基于一维集方法的脉冲多普勒雷达改进测距方法 - Google Patents
一种基于一维集方法的脉冲多普勒雷达改进测距方法 Download PDFInfo
- Publication number
- CN104375139B CN104375139B CN201410513115.0A CN201410513115A CN104375139B CN 104375139 B CN104375139 B CN 104375139B CN 201410513115 A CN201410513115 A CN 201410513115A CN 104375139 B CN104375139 B CN 104375139B
- Authority
- CN
- China
- Prior art keywords
- doppler radar
- pulse
- pulse doppler
- range
- transmission signal
- 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
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/02—Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
- G01S13/06—Systems determining position data of a target
- G01S13/08—Systems for measuring distance only
- G01S13/10—Systems for measuring distance only using transmission of interrupted, pulse modulated waves
- G01S13/18—Systems for measuring distance only using transmission of interrupted, pulse modulated waves wherein range gates are used
-
- 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
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/02—Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
- G01S13/06—Systems determining position data of a target
- G01S13/08—Systems for measuring distance only
- G01S13/10—Systems for measuring distance only using transmission of interrupted, pulse modulated waves
- G01S13/22—Systems for measuring distance only using transmission of interrupted, pulse modulated waves using irregular pulse repetition frequency
- G01S13/222—Systems for measuring distance only using transmission of interrupted, pulse modulated waves using irregular pulse repetition frequency using random or pseudorandom pulse repetition frequency
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明涉及一种基于一维集方法的脉冲多普勒雷达改进测距方法,其步骤为:建立查找表,令查找表中的每个数值取0;得出雷达发射信号的第i个脉冲重复频率Pi对应的视在距离门限si和最大模糊度fi;根据目标回波信号得出雷达发射信号的第i个脉冲重复频率Pi对应的视在距离xi;对视在距离门限si对应的最大模糊度fi进行更新,得到更新后的视在距离门限si对应的最大模糊度fxi;出距离单元序号集合距离单元序号集合中第q个元素riq为:riq=q*qi+xi,q=0,1,2...fxi;令查找表中第i行第Ri列的数值更新为1,Ri=ri0,ri1,ri2...,令j依次取1至N,如查找表中第j列的所有数值之和yj大于或等于L,则得出目标与脉冲多普勒雷达的距离A。
Description
技术领域
本发明属于脉冲多普勒解距离模糊技术领域,特别涉及一种基于一维集方法的脉冲多普勒雷达改进测距方法,适用于实际工程应用。
背景技术
一直以来,解距离模糊是脉冲多普勒雷达特有的功能。脉冲多普勒雷达是利用信号频域特性分辨和检测目标的脉冲雷达,脉冲多普勒雷达的这种信号处理方式可获得近于最佳的信噪比,及较精确的目标距离和径向速度数据,有良好的杂波抑制性能,在雷达上尤其在机载雷达上得到了越来越广泛的应用。
脉冲多普勒雷达由于在工作周期内交替发送多种频率的脉冲,使得每种脉冲的最大不模糊距离都远小于雷达探测距离,导致低、中、高重频的脉冲多普勒雷达在速度和距离上产生不同程度的模糊。解模糊技术是脉冲多普勒雷达的关键技术,采用合理而快速的方法对保证目标的检测概率和降低虚警率都有重要意义。
解决脉冲多普勒解模糊问题,主要有两种方法:
1)一维集方法。一维集方法的实质是利用穷举法解脉冲重复频率的同余方程组。即列举出对应当前目标视在距离的所有可能的真实距离,按照脉冲重复频率的不同将计算的真实距离分组,并从小到大依次排列,从每组结果各抽出一个组成所有可能的真实距离集合,分别计算各集合的均值和方差,方差值最小的集合即为最优解,其均值则为真实距离解。一维集方法可以准确地求出目标的真实距离值,但由于涉及求均值和方差,在数据量大的时候,尤其当目标的模糊度大时,计算量相当大,在工程中难以满足实时信号处理的要求。
2)查表法。首先确定脉冲重复频率及其个数、测距精度和测距范围,再在雷达威力范围内按不模糊距离单元由小到大的顺序建立表格,表格的内容为根据中国余数定理算出的所有不模糊距离单元所对应的各个脉冲重复频率的视在距离单元号。查表法简单易行,但随着目标的测距范围增大,存储空间需求快速膨胀,在查表过程中会根据距离单元的划分程度,可能会浪费大量时间在不必要的距离值上。
上述方法中,一维集方法随着测距范围的增大计算量呈指数增加,时间消耗量大;查表法随着测距范围增大存储空间需求量呈指数增大,对硬件存储能力要求高,当存储空间不足时测距精度就会相应降低。因此在处理距离解模糊问题的过程中,由于时间和存储空间的限制,现有的各种测距方法实效性差,不再适用。
发明内容
本发明的目的在于提出一种基于一维集方法的脉冲多普勒雷达改进测距方法,以进一步降低运算量和数据存储量,同时保证测距精度。
为实现上述目的,本发明包括以下步骤:
步骤1,根据脉冲多普勒雷达发射信号的每个脉冲重复频率和脉冲多普勒雷达的每个距离单元,建立查找表;所述查找表的每行对应脉冲多普勒雷达发射信号的一个脉冲重复频率,所述查找表的每列对应脉冲多普勒雷达的一个距离单元;令查找表中的每个数值取0;
步骤2,令i=1,2,...,K,K表示脉冲多普勒雷达发射信号的脉冲重复频率的个数;得出脉冲多普勒雷达发射信号的第i个脉冲重复频率Pi对应的视在距离门限si和视在距离门限si对应的最大模糊度fi;
步骤3,利用脉冲多普勒雷达向目标发射电磁波,获得对应的目标回波信号;根据目标回波信号得出脉冲多普勒雷达发射信号的第i个脉冲重复频率Pi对应的视在距离xi;
步骤4,根据视在距离门限si和视在距离门限si对应的最大模糊度fi,得到目标的视在距离单元的序号xi所对应的最大模糊度fxi:
得出距离单元序号集合距离单元序号集合中第q个元素riq为:riq=q*qi+xi,q=0,1,2...fxi;
令查找表中第i行第Ri列的数值更新为1,步骤5,令j依次取1至N,如果查找表中第j列的所有数值之和yj大于或等于L,则得出目标与脉冲多普勒雷达的距离A,A=j*c',c'表示脉冲多普勒雷达的测距精度,L为设定的整数且K>L≥K/2。
本发明的有益效果为:
1)本发明对查表法的存储方式进行优化,并用于解决了一维集方法计算量大的问题,简化了计算流程,有效提高了方法的性能。
2)本发明在测距范围大的时候更具有优势,由于查找表所需的存储空间非常小,大大缩短了查表的运算周期,也可以避免存储空间对查表精度的限制,为工程实现时,尤其是存储空间较小和工作周期较短的情况,提供了一种有效的解决方法。
附图说明
图1为本发明的一种基于一维集方法的脉冲多普勒雷达改进测距方法的流程图;
图2,为仿真实验中采用两种方法得出目标与脉冲多普勒雷达的距离时在不同测距范围内对存储空间的需求量示意图。
具体实施方式
下面结合附图对本发明作进一步说明:
参照图1,为本发明的一种基于一维集方法的脉冲多普勒雷达改进测距方法的流程图。该基于一维集方法的脉冲多普勒雷达改进测距方法包括以下步骤:
步骤1,根据脉冲多普勒雷达发射信号的每个脉冲重复频率和脉冲多普勒雷达的每个距离单元,建立查找表;所述查找表的每行对应脉冲多普勒雷达发射信号的一个脉冲重复频率,所述查找表的行数为K,K表示脉冲多普勒雷达发射信号的脉冲重复频率的个数;脉冲多普勒雷达发射信号的第i个脉冲重复频率表示Pi,i=1,2,...,K;所述查找表的每列对应脉冲多普勒雷达的一个距离单元,所述查找表的列数为N,N表示脉冲多普勒雷达的距离单元数;令查找表第i行第j列的数值yij取0,j=1,2,...,N。
本发明实施例中,脉冲多普勒雷达的距离单元数N的计算公式为:N=D/c',D表示脉冲多普勒雷达的测距范围,c'表示脉冲多普勒雷达的测距精度。
步骤2,令i=1,2,...,K,K表示脉冲多普勒雷达发射信号的脉冲重复频率的个数;得出脉冲多普勒雷达发射信号的第i个脉冲重复频率Pi对应的视在距离门限si和视在距离门限si对应的最大模糊度fi。
其具体子步骤为:
(2.1)求取脉冲多普勒雷达发射信号的第i个脉冲重复频率Pi对应的最大不模糊距离Qi:
其中,c表示光速,i=1,2,...,K。
(2.2)得出脉冲多普勒雷达发射信号的第i个脉冲重复频率Pi对应的不模糊距离单元的个数qi,qi=Qi/c',c'表示脉冲多普勒雷达的测距精度,i=1,2,...,K。
对于脉冲多普勒雷达发射信号的第i个脉冲重复频率Pi,其对应的视在距离门限si和其对应的最大模糊度fi有如下关系:
N=fi*qi+si
其中,N表示脉冲多普勒雷达的距离单元数。
因此脉冲多普勒雷达发射信号的第i个脉冲重复频率Pi对应的视在距离门限si和视在距离门限si对应的最大模糊度fi分别为:
si=N%qi,
其中,N%qi表示用N除以qi所得余数,N表示脉冲多普勒雷达的距离单元数,qi表示脉冲多普勒雷达发射信号的第i个脉冲重复频率Pi对应的不模糊距离单元的个数,表示进行向下取整。
步骤3,利用脉冲多普勒雷达向目标发射电磁波,获得对应的目标回波信号。然后,根据目标回波信号得出目标的视在距离(模糊距离),得出脉冲多普勒雷达发射信号的第i个脉冲重复频率Pi对应的视在距离单元的序号xi。(张代忠,洪一,邱炜.脉冲多普勒雷达中的解模糊算法及实现.雷达科学与技术[J].2004.10.2(5):293-297)
步骤4,根据视在距离门限si和视在距离门限si对应的最大模糊度fi,得到目标的视在距离单元的序号xi所对应的最大模糊度fxi:
其中,xi表示脉冲多普勒雷达发射信号的第i个脉冲重复频率Pi对应的视在距离单元的序号,si表示脉冲多普勒雷达发射信号的第i个脉冲重复频率Pi对应的视在距离门限,fi表示视在距离门限si对应的最大模糊度fi;i=1,2,...,K,K表示脉冲多普勒雷达发射信号的脉冲重复频率的个数。根据fxi的计算公式,当xi>si时,fxi等于fi,反之,当xi≤si时,fxi等于fi-1。
然后,得出距离单元序号集合距离单元序号集合的元素个数为fxi+1,距离单元序号集合中第q个元素riq为:riq=q*qi+xi,q=0,1,2...fxi。距离单元序号集合的每个元素表示一个可能的真实距离单元序号。
令查找表中第i行第Ri列的数值更新为1,i=1,2,...,K,K表示脉冲多普勒雷达发射信号的脉冲重复频率的个数,
步骤5,根据m/n检测准则得出目标距离值A,m/n检测准则的含义是:对于m个检测条件当有n个检测条件满足时就确定结果有效。
在本发明实施例中,对于第j个距离单元,j=1,2,...,N,当脉冲多普勒雷达发射信号的K个脉冲重复频率中的L个脉冲重复频率的计算结果有效,目标所在距离单元为第j个距离单元,L为设定的整数且K>L≥K/2。
其具体子步骤为:
(5.1)当j=1时,转至子步骤(5.2);
(5.2)得出查找表中第j列的所有数值之和yj,N表示脉冲多普勒雷达的距离单元数;yj为:
yj=y1j+y2j+y3j+...+yij+...+yKj,
其中,yij表示查找表第i行第j列的数值,i=1,2,...,K,K表示脉冲多普勒雷达发射信号的脉冲重复频率的个数。
(5.3)如果yj<L,则令j的值自增1,返回至子步骤(5.2);如果yj≥L,则得出目标所在距离单元,目标所在距离单元为第j个距离单元,然后得出目标与脉冲多普勒雷达的距离A,A=j*c',c'表示脉冲多普勒雷达的测距精度。
本发明的效果通过以下仿真实验进一步说明:
1)仿真条件
仿真过程中,采用FPGA作为硬件实现环境,脉冲多普勒雷达发射信号的脉冲重复频率的个数为3,脉冲多普勒雷达发射信号3个脉冲重复频率分别为1.5KHz、2.0KHz以及2.5KHz,脉冲多普勒雷达的测距精度为50m。
2)仿真内容
分别采用查表法和本发明得出目标与脉冲多普勒雷达的距离;计算两种方法在不同测距范围内对存储空间的需求量。参照图2,为仿真实验中采用两种方法得出目标与脉冲多普勒雷达的距离时在不同测距范围内对存储空间的需求量示意图。图2中,横轴表示距离单元的序号(从0到6000),纵轴表示表示存储空间的需求量,单位为105bit(位)。从图2中可以得出,本发明所需的存储空间需求量远小于查表法。
显然,本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围。这样,倘若本发明的这些修改和变型属于本发明权利要求及其等同技术的范围之内,则本发明也意图包含这些改动和变型在内。
Claims (3)
1.一种基于一维集方法的脉冲多普勒雷达改进测距方法,其特征在于,包括以下步骤:
步骤1,根据脉冲多普勒雷达发射信号的每个脉冲重复频率和脉冲多普勒雷达的每个距离单元,建立查找表;所述查找表的每行对应脉冲多普勒雷达发射信号的一个脉冲重复频率,所述查找表的每列对应脉冲多普勒雷达的一个距离单元;令查找表中的每个数值取0;
步骤2,令i=1,2,...,K,K表示脉冲多普勒雷达发射信号的脉冲重复频率的个数;得出脉冲多普勒雷达发射信号的第i个脉冲重复频率Pi对应的视在距离门限si和视在距离门限si对应的最大模糊度fi;
步骤3,利用脉冲多普勒雷达向目标发射电磁波,获得对应的目标回波信号;根据目标回波信号得出脉冲多普勒雷达发射信号的第i个脉冲重复频率Pi对应的视在距离xi;
步骤4,根据视在距离门限si和视在距离门限si对应的最大模糊度fi,得到目标的视在距离单元的序号xi所对应的最大模糊度fxi:
得出距离单元序号集合距离单元序号集合中第q个元素riq为:riq=q*qi+xi,q=0,1,2...fxi;qi表示脉冲多普勒雷达发射信号的第i个脉冲重复频率Pi对应的不模糊距离单元的个数;
令查找表中第i行第Ri列的数值更新为1,
步骤5,令j依次取1至N,如果查找表中第j列的所有数值之和yj大于或等于L,则得出目标与脉冲多普勒雷达的距离A,A=j*c’,c’表示脉冲多普勒雷达的测距精度,L为设定的整数且K>L≥K/2。
2.如权利要求1所述的一种基于一维集方法的脉冲多普勒雷达改进测距方法,其特征在于,所述步骤2的具体子步骤为:
(2.1)求取脉冲多普勒雷达发射信号的第i个脉冲重复频率Pi对应的最大不模糊距离Qi:
其中,c表示光速,i=1,2,...,K;
(2.2)得出脉冲多普勒雷达发射信号的第i个脉冲重复频率Pi对应的不模糊距离单元的个数qi,qi=Qi/c’,c’表示脉冲多普勒雷达的测距精度,i=1,2,...,K;
则脉冲多普勒雷达发射信号的第i个脉冲重复频率Pi对应的视在距离门限si和视在距离门限si对应的最大模糊度fi分别为:
si=N%qi,
其中,N%qi表示用N除以qi所得余数,N表示脉冲多普勒雷达的距离单元数,表示进行向下取整。
3.如权利要求1所述的一种基于一维集方法的脉冲多普勒雷达改进测距方法,其特征在于,所述步骤5的具体子步骤为:
(5.1)j=1,2,...,N,当j=1时,转至子步骤(5.2);
(5.2)得出查找表中第j列的所有数值之和yj,N表示脉冲多普勒雷达的距离单元数;yj为:
yj=y1j+y2j+y3j+...+yij+…+yKj,
其中,yij表示查找表第i行第j列的数值,i=1,2,...,K;
(5.3)如果yi<L,则令j的值自增1,返回至子步骤(5.2);如果yj≥L,则得出目标所在距离单元,目标所在距离单元为第j个距离单元,然后得出目标与脉冲多普勒雷达的距离A,A=j*c’,c’表示脉冲多普勒雷达的测距精度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410513115.0A CN104375139B (zh) | 2014-09-29 | 2014-09-29 | 一种基于一维集方法的脉冲多普勒雷达改进测距方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410513115.0A CN104375139B (zh) | 2014-09-29 | 2014-09-29 | 一种基于一维集方法的脉冲多普勒雷达改进测距方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104375139A CN104375139A (zh) | 2015-02-25 |
CN104375139B true CN104375139B (zh) | 2017-01-25 |
Family
ID=52554172
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410513115.0A Expired - Fee Related CN104375139B (zh) | 2014-09-29 | 2014-09-29 | 一种基于一维集方法的脉冲多普勒雷达改进测距方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104375139B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106291497B (zh) * | 2016-07-27 | 2018-07-06 | 陕西黄河集团有限公司 | 基于快速查表法的解速度模糊算法 |
CN107907884B (zh) * | 2017-10-02 | 2021-07-16 | 广东美的制冷设备有限公司 | 物体距离检测方法、装置、空调器及存储介质 |
CN109164441B (zh) * | 2018-09-30 | 2020-11-17 | 厦门大学 | 一种雷达测距的方法 |
CN111289953B (zh) * | 2020-01-14 | 2022-04-26 | 北京理工大学 | 一种基于模糊矩阵更新的天基雷达距离/速度解模糊方法 |
CN116755073B (zh) * | 2023-06-21 | 2024-03-26 | 上海雷骥电子科技有限公司 | 查找表解距离模糊的方法及应用 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5781149A (en) * | 1997-05-13 | 1998-07-14 | Long; Maurice W. | Doppler radar clutter spike rejector |
EP2343570A1 (en) * | 2010-01-11 | 2011-07-13 | BAE Systems PLC | Pulse radar range profile motion compensation |
CN102819010A (zh) * | 2012-08-14 | 2012-12-12 | 成都西科微波通讯有限公司 | 一种改进的二维恒虚警方法 |
-
2014
- 2014-09-29 CN CN201410513115.0A patent/CN104375139B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5781149A (en) * | 1997-05-13 | 1998-07-14 | Long; Maurice W. | Doppler radar clutter spike rejector |
EP2343570A1 (en) * | 2010-01-11 | 2011-07-13 | BAE Systems PLC | Pulse radar range profile motion compensation |
CN102819010A (zh) * | 2012-08-14 | 2012-12-12 | 成都西科微波通讯有限公司 | 一种改进的二维恒虚警方法 |
Non-Patent Citations (3)
Title |
---|
Medium PRF set selection: an approach through combinatorics;A.G.Thomas and M.C.Berg;《IEE Proceedings-Radar,Sonar and Navigation》;20020806;第141卷(第6期);第307-311页 * |
Range and velocity ambiguity resolution;G.Trunk and S.Brockett;《Radar Conference,1993.,Record of the 1993 IEEE National》;19931231;第146-149页 * |
脉冲多普勒雷达中的解模糊算法及实现;张代忠;《雷达科学与技术》;20041031;第2卷(第5期);第293-314页 * |
Also Published As
Publication number | Publication date |
---|---|
CN104375139A (zh) | 2015-02-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104375139B (zh) | 一种基于一维集方法的脉冲多普勒雷达改进测距方法 | |
CN109814073B (zh) | 一种mtd雷达解模糊测速方法 | |
Kim et al. | Human activity classification based on micro-Doppler signatures using a support vector machine | |
CN104215951B (zh) | 一种在海杂波背景下慢速小目标的检测***及其检测方法 | |
CN106468770B (zh) | K分布杂波加噪声下的近最优雷达目标检测方法 | |
CN104914415B (zh) | 基于目标距离像模板匹配的单脉冲雷达相干干扰方法 | |
Cain et al. | Convolutional neural networks for radar emitter classification | |
CN113221631B (zh) | 基于卷积神经网络的序列脉冲抗干扰目标检测方法 | |
CN104459661B (zh) | 检测快速火炮类微弱目标的方法 | |
CN104316914A (zh) | 依赖形状参数的雷达目标自适应检测方法 | |
CN104215959B (zh) | 一种多机动目标径向初速度和径向加速度的估计方法 | |
CN103885057A (zh) | 自适应变滑窗多目标跟踪方法 | |
CN106855941A (zh) | 基于雷达微多普勒信号稀疏优化的手势识别方法及*** | |
CN103954935A (zh) | 一种最小l1范数下的雷达信号分选方法 | |
CN104111449A (zh) | 一种改进的基于广义内积的空时二维自适应处理方法 | |
CN113504522B (zh) | 一种基于发射天线随机切换的时空解耦与超分辨测角方法 | |
CN105510895B (zh) | 一种随机脉冲间隔序列的信号处理方法 | |
CN106383344A (zh) | 基于融合准则的多站雷达运动目标检测方法 | |
CN113030932B (zh) | 一种扩展目标稳健自适应检测方法与*** | |
CN105093196A (zh) | 基于逆伽马纹理复合高斯模型下的相干检测方法 | |
CN105785330A (zh) | 一种认知型副瓣干扰抑制方法 | |
CN114117912A (zh) | 一种数据模型双驱动下的海杂波建模与抑制方法 | |
CN106772275A (zh) | 一种基于Alpha稳定分布的低秩矩阵恢复检测方法 | |
CN104375128B (zh) | 一种基于互相关函数的高机动目标快速积累检测方法 | |
CN104950292B (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: 20170125 Termination date: 20170929 |