CN105652249A - 一种干扰环境下的目标检测方法 - Google Patents
一种干扰环境下的目标检测方法 Download PDFInfo
- Publication number
- CN105652249A CN105652249A CN201610007990.0A CN201610007990A CN105652249A CN 105652249 A CN105652249 A CN 105652249A CN 201610007990 A CN201610007990 A CN 201610007990A CN 105652249 A CN105652249 A CN 105652249A
- Authority
- CN
- China
- Prior art keywords
- value
- target
- detection
- pulse
- confidence
- 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
- 238000001514 detection method Methods 0.000 title claims abstract description 88
- 238000000034 method Methods 0.000 claims description 20
- 238000007906 compression Methods 0.000 claims description 13
- 230000006835 compression Effects 0.000 claims description 9
- 230000008878 coupling Effects 0.000 claims description 6
- 238000010168 coupling process Methods 0.000 claims description 6
- 238000005859 coupling reaction Methods 0.000 claims description 6
- 230000001932 seasonal effect Effects 0.000 claims description 5
- 230000006870 function Effects 0.000 abstract description 27
- 230000010365 information processing Effects 0.000 abstract description 4
- 238000005259 measurement Methods 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 4
- 238000012913 prioritisation Methods 0.000 description 4
- 238000005457 optimization Methods 0.000 description 3
- 238000012360 testing method Methods 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 230000008485 antagonism Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000004927 fusion Effects 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
- 230000035485 pulse pressure Effects 0.000 description 1
- 230000009131 signaling function Effects 0.000 description 1
- 239000000758 substrate Substances 0.000 description 1
- 238000012549 training Methods 0.000 description 1
- 238000013519 translation 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
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/36—Means for anti-jamming, e.g. ECCM, i.e. electronic counter-counter measures
Abstract
本发明公开了一种干扰环境下的目标检测方法,在干扰环境下,将目标信息建模成雷达信息处理过程中的一个关于目标状态的时间序列(x1,x2,…,xk),值函数K在状态xk的数值用于反映目标和干扰在幅度、多普勒以及包括运动特性在内的其他细微特征上的差异。通过对所有可能的目标状态时间序列进行搜索,找到一组使得目标值函数K达到最大的状态序列。再通过使用DPA对多个脉冲的值函数K进行积累,最后才给出目标的检测结果,有效改善干扰环境下的目标检测性能。
Description
技术领域
本发明属于雷达目标检测技术领域,尤其涉及一种干扰环境下的目标检测方法。
背景技术
随着新的电子干扰技术的迅速发展,先进的电子干扰措施与装备的不断涌现,导致现代雷达面临的工作电磁环境日趋复杂。特别是近年来迅猛发展的基于数字射频存储器(DRFM)技术的有源欺骗干扰,以其样式多、使用灵活、对抗性强等特点,对军用雷达构成了最常见的威胁。
基于数字射频存储器的干扰机将接收到的雷达信号存储,然后根据需要对信号进行延时、相位调制等处理后再转发到雷达,产生类似于目标回波的信号作用于雷达检测和跟踪系统中,进行速度欺骗干扰和距离欺骗干扰等。转发式干扰信号与雷达信号的相关性很高,可通过雷达匹配滤波器获得高增益。另外,转发式干扰机为达到干扰雷达的目的,通常会产生功率很大的干扰信号,在所产生的干扰信号功率一定的情况下,雷达接收到的干扰信号的功率与距离的二次方成反比,而真实目标的雷达回波功率与距离的四次方成反比,所以雷达接收到的转发式干扰的功率通常远高于真实目标的功率,这就使得转发式干扰信号比真实目标回波信号更容易超过检测门限而被检测出来。所以,基于数字射频存储器转发的欺骗式干扰很容易进入雷达接收系统,造成大量虚警,影响对目标的正常检测和跟踪。
检测前跟踪(TrackBeforeDetect,TBD)是在低信噪比的情况下对目标进行检测和跟踪的一种技术。与一般的检测方法的不同之处在于,检测前跟踪在单帧内并不宣布检测结果,而是将单帧信息数字化并存储起来,对多帧数据处理后同时宣布检测结果与目标的航迹。
时间序列的动态规划方法(DPA)是穷尽搜索法的一种等效实现算法,但是它的计算效率远远高于穷尽搜索。它的主要思想是将n维优化问题转化为n个1维优化问题,对优化问题的分级处理大大降低了它的计算量。
将时间序列的动态规划方法DPA应用到雷达在干扰环境下的目标检测问题时,它把目标信息看成雷达信息处理过程中的一个关于目标状态的时间序列(x1,x2,...,xk),状态xk描述了在雷达信息处理过程中的第k个状态时关于目标的相关信息。对干扰环境下的目标检测就是N个脉冲通过对所有可能的目标状态的时间序列(x1,x2,...,xk)进行搜索,找到一组使得值函数K达到最大的状态序列(x1,x2,...,xk)。值函数K能够反映目标和干扰在幅度、多普勒以及包括运动特性在内的其他细微特征上的差异,具有预期目标特征状态可以获得较高的值函数K,相反,具有干扰特征的状态获得较低的值函数K。再通过使用DPA对多个脉冲的值函数K进行积累,可以有效改善干扰环境下的目标检测性能。
发明内容
本发明提供了一种干扰环境下的目标检测方法。在干扰环境下,本发明将目标信息建模成雷达信息处理过程中的一个关于目标状态的时间序列(x1,x2,...,xk),值函数K在状态xk的数值用于反映目标和干扰在幅度、多普勒以及包括运动特性在内的其他细微特征上的差异。通过对所有可能的目标状态时间序列进行搜索,找到一组使得目标值函数K达到最大的状态序列。再通过使用DPA对多个脉冲的值函数K进行积累,最后才给出目标的检测结果,有效改善干扰环境下的目标检测性能。
本发明为解决上述技术问题采用以下技术方案:
一种干扰环境下的目标检测方法,包括以下步骤:
步骤1),将N个脉冲数据送入雷达目标检测器;
步骤2),对送入的N个脉冲数据分别进行目标状态时间序列信息检测,形成N个脉冲值函数,具体为:
步骤2.1),对于每一个送入的脉冲数据,将其幅度与预设的幅度门限阈值进行比较,如果其幅度大于等于预设的幅度门限阈值,赋予其预设的第一置信值x1,然后将第一置信值x1以及被赋予第一置信值x1的脉冲数据送入后续的目标状态时间序列检测链路;
步骤2.2),对步骤2.1)送入的每一个脉冲数据进行脉冲匹配压缩,将其脉冲匹配压缩后的输出信噪比与预设的第一信噪比门限阈值进行比较,如果其输出信噪比值大于等于预设的第一信噪比门限阈值,再将其目标峰的宽度和强度分别和未受干扰时目标脉冲匹配压缩后目标峰的的宽度和强度进行比较,如果其目标峰的宽度和强度分别大于等于未受干扰时目标脉冲匹配压缩后目标峰的的宽度和强度,赋予其预设的第二置信值x2,然后将第二置信值x2以及脉冲匹配压缩后被赋予第二置信值x2的脉冲数据送入后续的目标状态时间序列检测链路;
步骤2.3),对步骤2.2)送入的每一个脉冲数据进行动目标显示或动目标检测,将其动目标显示或动目标检测后的输出信噪比与预设的第二信噪比门限阈值进行比较,如果其输出信噪比大于等于预设的第二信噪比门限阈值,赋予其预设的第三置信值x3,然后将第三置信值x3以及动目标显示或动目标检测后被赋予第三置信值x3的脉冲数据送入后续的目标状态时间序列检测链路;
步骤2.4),对步骤2.3)送入的每一个脉冲数据进行恒虚警检测,将其恒虚警检测输出信噪比与预设的第三信噪比门限阈值进行比较,如果其输出信噪比大于等于预设的第三信噪比门限阈值,赋予其预设的第四置信值x4,然后第四置信值x4以及恒虚警检测后被赋予第四置信值x4的脉冲数据送入后续的目标状态时间序列检测链路;
步骤2.5),计算目标状态的时间序列(x1,x2,...,xk)在xi为不同状态值下的值函数K,形成条时间检测链路,其中,k=4;i=1,2,3,4;nk为目标状态的时间序列(x1,x2,...,xk)在状态xk处的状态个数;值函数K用于反映目标和干扰在幅度、多普勒以及运动特性的差异;
步骤2.6),选取使得值函数K达到最优的检测链路;
步骤3),对N个脉冲值函数K进行联合检测,并给出点迹;
步骤4),对步骤3)送入的点迹进行点迹航迹关联处理,给出航迹。
作为本发明一种干扰环境下的目标检测方法进一步的优化方案,所述步骤2.6)中根据目标状态时间序列的最大概率准则选取使得值函数K达到最优的检测链路。
作为本发明一种干扰环境下的目标检测方法进一步的优化方案,所述步骤3)中对N个脉冲值函数K进行M/N联合检测。
作为本发明一种干扰环境下的目标检测方法进一步的优化方案,所述步骤3)中根据最大概率准则对N个脉冲值函数K进行M/N联合检测。
作为本发明一种干扰环境下的目标检测方法进一步的优化方案,进行M/N联合检测时,个数M大于等于
本发明的有益效果是:
本发明采用多脉冲目标状态的时间序列联合检测方法来解决雷达干扰环境下的目标检测问题,它避免了目标检测简单门限判决造成的目标点迹丢失问题,又能够利用目标和干扰在不同状态下的细微特征上的差异,将目标和干扰区分开来,同时,又利用多个脉冲的联合检测,提高目标的正确发现概率。
附图说明
图1为本发明的流程图;
图2为常规的雷达信号处理目标检测流程图;
图3是本发明提出的目标检测方法对实测转发假目标干扰数据的处理航迹;
图4为是经典目标检测方法对实测转发假目标干扰数据的处理航迹。
具体实施方式
下面结合附图对本发明的技术方案做进一步的详细说明:
如图1所示,本发明公开了一种干扰环境下的目标检测方法,包括以下步骤:
步骤1),将N个脉冲数据送入雷达目标检测器;
步骤2),对送入的N个脉冲数据分别进行目标状态时间序列信息检测,形成N个脉冲值函数,具体为:
步骤2.1),预处理,对于每一个送入的脉冲数据,将其幅度与预设的幅度门限阈值进行比较,如果其幅度大于等于预设的幅度门限阈值,赋予其预设的第一置信值x1,然后将第一置信值x1以及被赋予第一置信值x1的脉冲数据送入后续的目标状态时间序列检测链路;
步骤2.2),对步骤2.1)送入的每一个脉冲数据进行脉冲匹配压缩,将其脉冲匹配压缩后的输出信噪比与预设的第一信噪比门限阈值进行比较,如果其输出信噪比值大于等于预设的第一信噪比门限阈值,再将其目标峰的宽度和强度分别和未受干扰时目标脉冲匹配压缩后目标峰的的宽度和强度进行比较,如果其目标峰的宽度和强度分别大于等于未受干扰时目标脉冲匹配压缩后目标峰的的宽度和强度,赋予其预设的第二置信值x2,然后将第二置信值x2以及脉冲匹配压缩后被赋予第二置信值x2的脉冲数据送入后续的目标状态时间序列检测链路;
步骤2.3),对步骤2.2)送入的每一个脉冲数据进行动目标显示(MTI)或动目标检测(MTD),将其动目标显示或动目标检测后的输出信噪比与预设的第二信噪比门限阈值进行比较,如果其输出信噪比大于等于预设的第二信噪比门限阈值,赋予其预设的第三置信值x3,然后将第三置信值x3以及动目标显示或动目标检测后被赋予第三置信值x3的脉冲数据送入后续的目标状态时间序列检测链路;
步骤2.4),对步骤2.3)送入的每一个脉冲数据进行恒虚警检测(CFAR),将其恒虚警检测输出信噪比与预设的第三信噪比门限阈值进行比较,如果其输出信噪比大于等于预设的第三信噪比门限阈值,赋予其预设的第四置信值x4,然后第四置信值x4以及恒虚警检测后被赋予第四置信值x4的脉冲数据送入后续的目标状态时间序列检测链路;
步骤2.5),计算目标状态的时间序列(x1,x2,...,xk)在xi为不同状态值下的值函数K,形成条时间检测链路,其中,k=4;i=1,2,3,4;nk为目标状态的时间序列(x1,x2,...,xk)在状态xk处的状态个数;值函数K用于反映目标和干扰在幅度、多普勒以及运动特性的差异;
步骤2.6),选取使得值函数K达到最优的检测链路;
步骤3),对N个脉冲值函数K进行联合检测,并给出点迹;
步骤4),对步骤3)送入的点迹进行点迹航迹关联处理,给出航迹。
所述步骤2.6)中根据目标状态时间序列的最大概率准则选取使得值函数K达到最优的检测链路。
所述步骤3)中根据最大概率准则对N个脉冲值函数K进行M/N联合检测,进行M/N联合检测时,个数M大于等于
对实测转发假目标干扰数据进行处理,在预处理时,将干扰噪声基底超过雷达训练时噪声基地门限30dB、50dB和70dB的置信值分别设为60%、75%和90%。对脉冲压缩处理后的超过门限的点的脉压峰和间距设置一定得门限,并给出相应的置信值。在进行多脉冲联合检测时,将本雷达的8个脉冲按照M/N准则联合检测,只要有超过4个脉冲的值函数K超过门限,则认为该距离点上存在一个目标,并将该点的坐标和信息送至航迹融合处理。
图2为常规的雷达信号处理目标检测流程图。
图3为本发明提出的目标检测方法对实测转发假目标干扰数据的处理航迹,图4为是经典目标检测方法对实测转发假目标干扰数据的处理航迹。比较两图可知,本发明提出的方法处理结果虚假航迹点明显减少,并且能够明显改善经典方法处理时出现的航迹丢失情况。
本技术领域技术人员可以理解的是,除非另外定义,这里使用的所有术语(包括技术术语和科学术语)具有与本发明所属领域中的普通技术人员的一般理解相同的意义。还应该理解的是,诸如通用字典中定义的那些术语应该被理解为具有与现有技术的上下文中的意义一致的意义,并且除非像这里一样定义,不会用理想化或过于正式的含义来解释。
以上所述的具体实施方式,对本发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施方式而已,并不用于限制本发明,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (5)
1.一种干扰环境下的目标检测方法,其特征在于,包括以下步骤:
步骤1),将N个脉冲数据送入雷达目标检测器;
步骤2),对送入的N个脉冲数据分别进行目标状态时间序列信息检测,形成N个脉冲值函数,具体为:
步骤2.1),对于每一个送入的脉冲数据,将其幅度与预设的幅度门限阈值进行比较,如果其幅度大于等于预设的幅度门限阈值,赋予其预设的第一置信值x1,然后将第一置信值x1以及被赋予第一置信值x1的脉冲数据送入后续的目标状态时间序列检测链路;
步骤2.2),对步骤2.1)送入的每一个脉冲数据进行脉冲匹配压缩,将其脉冲匹配压缩后的输出信噪比与预设的第一信噪比门限阈值进行比较,如果其输出信噪比值大于等于预设的第一信噪比门限阈值,再将其目标峰的宽度和强度分别和未受干扰时目标脉冲匹配压缩后目标峰的的宽度和强度进行比较,如果其目标峰的宽度和强度分别大于等于未受干扰时目标脉冲匹配压缩后目标峰的的宽度和强度,赋予其预设的第二置信值x2,然后将第二置信值x2以及脉冲匹配压缩后被赋予第二置信值x2的脉冲数据送入后续的目标状态时间序列检测链路;
步骤2.3),对步骤2.2)送入的每一个脉冲数据进行动目标显示或动目标检测,将其动目标显示或动目标检测后的输出信噪比与预设的第二信噪比门限阈值进行比较,如果其输出信噪比大于等于预设的第二信噪比门限阈值,赋予其预设的第三置信值x3,然后将第三置信值x3以及动目标显示或动目标检测后被赋予第三置信值x3的脉冲数据送入后续的目标状态时间序列检测链路;
步骤2.4),对步骤2.3)送入的每一个脉冲数据进行恒虚警检测,将其恒虚警检测输出信噪比与预设的第三信噪比门限阈值进行比较,如果其输出信噪比大于等于预设的第三信噪比门限阈值,赋予其预设的第四置信值x4,然后第四置信值x4以及恒虚警检测后被赋予第四置信值x4的脉冲数据送入后续的目标状态时间序列检测链路;
步骤2.5),计算目标状态的时间序列(x1,x2,...,xk)在xi为不同状态值下的值函数K,形成条时间检测链路,其中,k=4;i=1,2,3,4;nk为目标状态的时间序列(x1,x2,...,xk)在状态xk处的状态个数;值函数K用于反映目标和干扰在幅度、多普勒以及运动特性的差异;
步骤2.6),选取使得值函数K达到最优的检测链路;
步骤3),对N个脉冲值函数K进行联合检测,并给出点迹;
步骤4),对步骤3)送入的点迹进行点迹航迹关联处理,给出航迹。
2.根据权利要求1所述的一种干扰环境下的目标检测方法,其特征在于,所述步骤2.6)中根据目标状态时间序列的最大概率准则选取使得值函数K达到最优的检测链路。
3.根据权利要求1所述的一种干扰环境下的目标检测方法,其特征在于,所述步骤3)中对N个脉冲值函数K进行M/N联合检测。
4.根据权利要求3所述的一种干扰环境下的目标检测方法,其特征在于,所述步骤3)中根据最大概率准则对N个脉冲值函数K进行M/N联合检测。
5.根据权利要求3所述的一种干扰环境下的目标检测方法,其特征在于,进行M/N联合检测时,个数M大于等于
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610007990.0A CN105652249B (zh) | 2016-01-06 | 2016-01-06 | 一种干扰环境下的目标检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610007990.0A CN105652249B (zh) | 2016-01-06 | 2016-01-06 | 一种干扰环境下的目标检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105652249A true CN105652249A (zh) | 2016-06-08 |
CN105652249B CN105652249B (zh) | 2018-03-13 |
Family
ID=56491736
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610007990.0A Expired - Fee Related CN105652249B (zh) | 2016-01-06 | 2016-01-06 | 一种干扰环境下的目标检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105652249B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107290726A (zh) * | 2017-07-03 | 2017-10-24 | 中国人民解放军海军航空工程学院 | 基于信号—数据联合处理的有源压制干扰抑制方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB1144778A (en) * | 1966-04-07 | 1969-03-12 | Hughes Aircraft Co | Automatic target detection radar system |
JPS6157872A (ja) * | 1984-08-29 | 1986-03-24 | Nec Corp | 信号処理装置 |
US6271786B1 (en) * | 1999-06-07 | 2001-08-07 | Raytheon Company | Random noise radar target detection device |
CN104007424A (zh) * | 2014-05-28 | 2014-08-27 | 电子科技大学 | 基于时频分析的机动目标检测方法 |
-
2016
- 2016-01-06 CN CN201610007990.0A patent/CN105652249B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB1144778A (en) * | 1966-04-07 | 1969-03-12 | Hughes Aircraft Co | Automatic target detection radar system |
JPS6157872A (ja) * | 1984-08-29 | 1986-03-24 | Nec Corp | 信号処理装置 |
US6271786B1 (en) * | 1999-06-07 | 2001-08-07 | Raytheon Company | Random noise radar target detection device |
CN104007424A (zh) * | 2014-05-28 | 2014-08-27 | 电子科技大学 | 基于时频分析的机动目标检测方法 |
Non-Patent Citations (1)
Title |
---|
陈华明等: "基于动态规划和置信度检验的小目标检测", 《系统工程与电子技术》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107290726A (zh) * | 2017-07-03 | 2017-10-24 | 中国人民解放军海军航空工程学院 | 基于信号—数据联合处理的有源压制干扰抑制方法 |
CN107290726B (zh) * | 2017-07-03 | 2020-07-28 | 中国人民解放军海军航空大学 | 基于信号—数据联合处理的有源压制干扰抑制方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105652249B (zh) | 2018-03-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Aprile et al. | Track-before-detect for sea clutter rejection: Tests with real data | |
CN104914415B (zh) | 基于目标距离像模板匹配的单脉冲雷达相干干扰方法 | |
CN104267379A (zh) | 一种基于波形设计的主被动雷达协同抗干扰方法 | |
CN103558595B (zh) | 雷达通用恒虚警率检测器及其数据排序方法 | |
CN104931939A (zh) | 一种基于假目标特征的弹载转发式干扰抑制方法 | |
US8976059B2 (en) | Identification and removal of a false detection in a radar system | |
CN109444820A (zh) | 杂波和干扰共存时多通道雷达先干扰抑制后目标检测方法 | |
CN105891799B (zh) | 适用于机械扫描雷达的有源干扰侦察方法 | |
CN105388464A (zh) | 一种对抗数字射频存储转发干扰的多雷达干扰抑制方法 | |
CN108614245A (zh) | 一种基于射频隐身的边跟踪边干扰方法 | |
CN104931946A (zh) | 外辐射源雷达多普勒-时延两维双门限目标检测方法 | |
CN102866392A (zh) | 一种波形捷变认知雷达系统及实现方法 | |
CN204989459U (zh) | 一种针对单脉冲跟踪制导雷达的角度欺骗干扰装置 | |
CN108445461A (zh) | 一种多径条件下雷达目标检测方法 | |
CN102879774A (zh) | 短航迹合成方法和装置 | |
Yan et al. | Intelligent and fast two-dimensional CFAR procedure | |
Grossi et al. | A two-step multi-frame detection procedure for radar systems | |
CN105652249A (zh) | 一种干扰环境下的目标检测方法 | |
CN108535704A (zh) | 一种基于自适应二维聚类的信号预分选方法 | |
CN114325599B (zh) | 一种针对不同环境的自动门限检测方法 | |
CN105137398A (zh) | 基于遗传算法的雷达抗转发式干扰的脉压滤波器优化方法 | |
Aprile et al. | An application of track-before-detect to sea-clutter rejection: experimental results based on real data | |
Aboutanios et al. | Evaluation of the single and two data set STAP detection algorithms using measured data | |
Chen et al. | Non-coherent MIMO radar in a non-Gaussian noise-plus-clutter environment | |
US3718926A (en) | Instantaneous sidelobe interrogation decode suppressor |
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: 20180313 |
|
CF01 | Termination of patent right due to non-payment of annual fee |