CN103364772A - 基于实数域广义多重信号分类算法的目标低仰角估计方法 - Google Patents
基于实数域广义多重信号分类算法的目标低仰角估计方法 Download PDFInfo
- Publication number
- CN103364772A CN103364772A CN2013102965300A CN201310296530A CN103364772A CN 103364772 A CN103364772 A CN 103364772A CN 2013102965300 A CN2013102965300 A CN 2013102965300A CN 201310296530 A CN201310296530 A CN 201310296530A CN 103364772 A CN103364772 A CN 103364772A
- Authority
- CN
- China
- Prior art keywords
- theta
- real number
- number field
- matrix
- angle
- 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
Images
Landscapes
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种基于实数域广义多重信号分类算法的目标低仰角估计方法。其实现步骤是:1)对雷达接收回波进行采样;2)利用采样数据计算采样协方差矩阵;3)对采样协方差矩阵进行空间平滑和酉变换,得到实数域协方差矩阵;4)对实数域协方差矩阵进行特征值分解,获得噪声投影矩阵;5)构造实数域导向矢量流形;6)利用噪声投影矩阵和实数域导向矢量流形构造空间谱,对空间谱二维角度搜索,获得初次角度估计值;7)利用初次角度估计值估计多径衰减系数,构造二次空间谱,通过二维角度搜索得到角度估计值;8)比较估计值中两个角度大小,将最大角度作为目标仰角值。本发明降低了运算复杂度,提高了雷达在低信噪比下的角度估计性能。
Description
技术领域
本发明属于雷达信号处理技术领域,特别涉及目标的低仰角估计方法,可用于雷达在多径环境下对目标的探测与跟踪。
背景技术
在雷达对低空目标的探测和跟踪中,目标仰角值是一个非常重要的测量参数,只有准确估计出目标仰角才能实现对目标的稳定跟踪。雷达在测量和跟踪低空目标时,除了接收目标返回的直射信号外,还会接收到地面或海面产生的反射信号,从而形成多径效应,导致测角误差,甚至跟踪目标丢失。由于直射信号与反射信号是相干信号,且对于远距低空目标,两个信号的波达方向夹角往往小于阵列的半功率波束宽度,这就要求目标低仰角估计方法既能分辨相干信号又具有较高的角度分辨能力,因此具有较高分辨力的超分辨技术得到了广泛关注。但是由于许多超分辨测角技术运算量大,难以进行实时处理,如最大似然估计、多重信号分类方法等,因此,如何降低运算量也是一个亟待解决的重要问题。
目前,已经有很多学者利用阵列超分辨技术解决目标低仰角估计问题。张文俊等人在文章“广义MUSIC算法在米波雷达测高中的应用及改进”中提出了能够直接处理相干信号的广义多重信号分类测角方法。这种广义多重信号分类测角方法在保证天线的有效孔径不变的前提下对相干信号具有高分辨力,可用于低空目标俯仰角估计,其基本思想在于根据数据空间中的信号子空间与噪声子空间的正交性,利用多径衰减系数形成合成导向矢量,构造二维空间谱函数获得低仰角估计值。该方法虽然对低空目标仰角的估计精度高,但它需要对空间谱函数进行复数域二维角度搜索,计算量大。尤其对于低空高速目标,如反舰导弹、作战飞机等,该方法难以实现实时信号处理,不利于工程应用。
发明内容
本发明的目的在于针对上述已有技术的不足,提出了一种基于实数域广义多重信号分类算法的目标低仰角估计方法,在保证低空测角精度的同时减小运算复杂度,提高雷达在低信杂比下的角度估计性能,满足雷达对远距低空目标测角精度的要求。
实现本发明目的的技术思路是:利用酉变换将数据相关矩阵转换到实数域,并估计多径衰减系数,结合多径衰减系数和广义多重信号分类算法实现对低空目标的俯仰角估计,其实现步骤包括如下:
(1)在相控阵雷达的俯仰方向上设计含有N个阵元的等距离线阵,接收目标回波信号;
(2)对线阵接收到的回波信号进行采样,获得第n个阵元第l次快拍在目标距离单元处的采样数据xn(l),其中,n=1,2,…N,l=1,2,…L,L为采样数据的总快拍数;
(4)对采样协方差矩阵进行双向空间平滑,获得平滑后的协方差矩阵
(5)构造酉变换矩阵U:
当阵元个数N为偶数时:
式中,j表示虚部,IN/2为(N/2)×(N/2)维的单位矩阵,JN/2为(N/2)×(N/2)维的置换矩阵;
当阵元个数N为奇数时:
式中,I(N-1)/2为((N-1)/2)×((N-1)/2)维的单位矩阵,J(N-1)/2为((N-1)/2)×((N-1)/2)维的置换矩阵;
式中,[·]H表示矩阵共轭转置;
(7)对实数域协方差矩阵RU进行特征值分解,得到N个特征值及与特征值对应的特征矢量,对N个特征值从小到大排列,取前(N-2)个小特征值对应的特征矢量v1,v2…vN-2,构成实数域噪声子空间矩阵:
EnU=[v1 v2 … vN-2];
(8)利用实数域噪声子空间矩阵EnU,获得实数域噪声投影矩阵PU:
(9)根据已知的雷达角度搜索范围[α,β],构造(θ1,θ2)的实数域导向矢量流形AU(θ1,θ2),其中,α为最小搜索角度,β为最大搜索角度,(θ1,θ2)为两个不同的搜索角度θ1,θ2构成的向量,且θ1,θ2∈[α,β];
本发明与现有技术相比具有以下优点:
1)减小了角度估计的运算复杂度。现有的广义多重信号分类测角方法中的特征值分解、空间谱搜索等都是在复数域中进行的,运算量大。本发明将酉变换与广义多重信号分类算法相结合,将特征值分解和空间谱搜索过程转换到实数域进行,有效地减小了运算复杂度。
2)提高了低信杂比下低仰角目标波达方向估计精度。本发明充分利用多径衰减系数进行二次空间谱搜索,减小了低信噪比下角度估计值的误差,改善了雷达在低信杂比下对低仰角目标波达方向估计精度。
附图说明
图1是本发明的实现流程图;
图2是在多径环境下,本发明与广义多重信号分类测角方法的运算量对比图;
图3是在多径环境下,利用本发明与广义多重信号分类测角方法分别进行测角,得到的角度均方根误差变化曲线。
具体实施方式
参照图1,本发明的实现步骤如下:
步骤1,接收目标回波信号。
在相控阵雷达的俯仰方向上设计等距离线阵,该线阵由N个阵元组成,这些个阵元在俯仰方向上等间隔排列,且均为无方向性阵元,用于接收所有方向上的目标回波信号。
步骤2,获取目标所在距离单元的采样数据。
对线阵接收到的回波信号进行数字采样,通过目标检测技术确定目标所在距离单元,获得第n个阵元第l次快拍在目标所在距离单元处的采样数据xn(l),其中,n=1,2,…N,l=1,2,…L,L为采样数据的总快拍数;
(3a)将N个阵元的采样数据x1(l),x2(l),…,xN(l),排列成N×1维采样数据矢量X(l):
X(l)=[x1(l),x2(l),…,xN(l)]T,
其中,[·]T表示矩阵转置;
式中,[·]H表示矩阵的共轭转置。
式中,[·]*表示矩阵的共轭,置换矩阵JN为:
步骤5,根据阵元个数N的奇偶性,构造酉变换矩阵U。
酉变换矩阵U的构造与阵元个数N的奇偶性及单位矩阵和置换矩阵有关,通过阵元个数N确定单位矩阵和置换矩阵的维数,以保证酉变换矩阵U的正交性,为此本实例构造的酉变换矩阵U如下:
当阵元个数N为偶数时,
式中,j表示虚部,IN/2为(N/2)×(N/2)维的单位矩阵,JN/2为(N/2)×(N/2)维的置换矩阵;
当阵元个数N为奇数时:
式中,I(N-1)/2为((N-1)/2)×((N-1)/2)维的单位矩阵,J(N-1)/2为((N-1)/2)×((N-1)/2)维的置换矩阵。
步骤6,计算实数域协方差矩阵。
式中,[·]H表示矩阵的共轭转置。
步骤7,产生噪声子空间矩阵。
对实数域协方差矩阵RU进行特征值分解,得到N个特征值及与特征值对应的特征矢量,将这N个特征值从小到大排列,取前(N-2)个小特征值对应的特征矢量v1,v2…vN-2,构成实数域噪声子空间矩阵:
EnU=[v1 v2 … vN-2]。
步骤8,利用噪声子空间矩阵EnU,获得实数域噪声投影矩阵PU:
步骤9,根据已知的雷达角度搜索范围[α,β],利用酉变换矩阵U构造角度向量(θ1,θ2)的实数域导向矢量流形AU(θ1,θ2)。
(9a)根据已知的雷达载频f和阵元间距d,分别计算两个不同的搜索角度θ1和θ2方向上的信号导向矢量a(θ1)和a(θ2):
式中,(θ1,θ2)为两个不同的搜索角度θ1,θ2构成的角度向量,且θ1,θ2∈[α,β],α为最小搜索角度,β为最大搜索角度,j表示虚部,c表示光速;
(9b)利用上述两个信号导向矢量a(θ1)和a(θ2),构造角度向量(θ1,θ2)的N×2维导向矢量流形A(θ1,θ2):
A(θ1,θ2)=[a(θ1)a(θ2)];
(9c)对导向矢量流形A(θ1,θ2)进行酉变换,获得实数域导向矢量流形AU(θ1,θ2):
AU(θ1,θ2)=UHA(θ1,θ2),
其中,U为酉变换矩阵。
(10a)利用实数域导向矢量流形AU(θ1,θ2)和实数域噪声投影矩阵PU,得到空间谱函数F0(θ1,θ2):
式中,|·|表示行列式;
步骤11,估计多径衰减系数矩阵。
步骤12,计算二次空间谱的搜索区间。
式中,j表示虚部,f表示雷达载频,d表示雷达阵元间距,c表示光速;
其中,U为酉变换矩阵。
步骤14,对二次空间谱进行二维角度搜索。
步骤15,获得目标仰角估计值。
本发明的效果通过以下仿真实验进一步说明:
1.仿真条件:
假设一个俯仰方向上的均匀线阵,阵元数为34,阵元间距为半波长,雷达载频为14GHz,雷达天线高度为10m,接收到的信号均为窄带信号,噪声为与信号不相关的高斯白噪声,快拍数为30,多径镜面反射系数为-0.8。目标高度为100m,目标与雷达距离为6km,角度搜索间隔为0.05°,角度的均方根误差值是通过100次蒙特卡洛实验得到,均方根误差计算公式为:其中,Monte为蒙特卡洛实验次数,θreal为目标真实俯仰角度,表示第i次蒙特卡洛实验估计角度值。
2.仿真内容:
仿真1,改变雷达阵元个数,分别利用本发明和广义多重信号分类测角方法进行3000次实验,获得两者在不同阵元个数下核心步骤所消耗的时间,作出关于阵元个数的消耗时间曲线,如图2所示。
仿真2,分别利用本发明和广义多重信号分类测角方法进行100次角度估计,获得两者在不同单元信噪比下的角度均方根误差曲线,如图3所示。
3.仿真分析:
从图2中可以看出,本发明的运算量要小于广义多重信号分类测角方法,可降低到原来运算量的1/3左右,阵元个数越多,本发明在运算量上的优势越明显,有利于工程应用;
从图3中可以看出,本发明与广义多重信号分类测角方法相比,在高信噪比下,两者测角精度差不多,而在低信噪比下,本发明的测角精度更高,当单元信噪比为0dB时,与广义多重信号分类测角方法相比,精度可提高1-2mrad。
综上,本发明的角度估计性能明显优于广义多重信号分类测角方法,降低了运算复杂度,提高了雷达在低信噪比下的低仰角估计性能。
Claims (6)
1.一种基于实数域广义多重信号分类算法的目标低仰角估计方法,包括如下步骤:
(1)在相控阵雷达的俯仰方向上设计含有N个阵元的等距离线阵,接收目标回波信号;
(2)对线阵接收到的回波信号进行采样,获得第n个阵元第l次快拍在目标距离单元处的采样数据xn(l),其中,n=1,2,…N,l=1,2,…L,L为采样数据的总快拍数;
(5)构造酉变换矩阵U:
当阵元个数N为偶数时:
式中,j表示虚部,IN/2为(N/2)×(N/2)维的单位矩阵,JN/2为(N/2)×(N/2)维的置换矩阵;
当阵元个数N为奇数时:
式中,I(N-1)/2为((N-1)/2)×((N-1)/2)维的单位矩阵,J(N-1)/2为((N-1)/2)×((N-1)/2)维的置换矩阵;
式中,[·]H表示矩阵共轭转置;
(7)对实数域协方差矩阵RU进行特征值分解,得到N个特征值及与特征值对应的特征矢量,对N个特征值从小到大排列,取前(N-2)个小特征值对应的特征矢量v1,v2…vN-2,构成实数域噪声子空间矩阵:
EnU=[v1 v2 … vN-2];
(8)利用实数域噪声子空间矩阵EnU,获得实数域噪声投影矩阵PU:
(9)根据已知的雷达角度搜索范围[α,β],构造(θ1,θ2)的实数域导向矢量流形AU(θ1,θ2),其中,α为最小搜索角度,β为最大搜索角度,(θ1,θ2)为两个不同的搜索角度θ1,θ2构成的向量,且θ1,θ2∈[α,β];
4.根据权利要求1所述的基于实数域广义多重信号分类算法的目标低仰角估计方法,其中,步骤(9)所述的构造(θ1,θ2)的实数域导向矢量流形AU(θ1,θ2),按如下步骤进行:
(9a)分别计算搜索角度θ1和搜索角度θ2方向上的信号导向矢量a(θ1)和a(θ2):
式中,j表示虚部,f表示雷达载频,d表示雷达阵元间距,c表示光速;
(9b)利用上述两个信号导向矢量a(θ1)和a(θ2),构造搜索角度(θ1,θ2)的N×2维导向矢量流形A(θ1,θ2):
A(θ1,θ2)=[a(θ1) a(θ2)];
(9c)利用导向矢量流形A(θ1,θ2),获得实数域导向矢量流形AU(θ1,θ2):
AU(θ1,θ2)=UHA(θ1,θ2)
其中,U为酉变换矩阵。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310296530.0A CN103364772B (zh) | 2013-07-14 | 2013-07-14 | 基于实数域广义多重信号分类算法的目标低仰角估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310296530.0A CN103364772B (zh) | 2013-07-14 | 2013-07-14 | 基于实数域广义多重信号分类算法的目标低仰角估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103364772A true CN103364772A (zh) | 2013-10-23 |
CN103364772B CN103364772B (zh) | 2015-04-22 |
Family
ID=49366548
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310296530.0A Expired - Fee Related CN103364772B (zh) | 2013-07-14 | 2013-07-14 | 基于实数域广义多重信号分类算法的目标低仰角估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103364772B (zh) |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103913731A (zh) * | 2014-02-25 | 2014-07-09 | 西安电子科技大学 | 基于先验信息的目标仰角估计方法 |
CN105913044A (zh) * | 2016-05-04 | 2016-08-31 | 大连理工大学 | 一种基于Sigmoid协方差矩阵的多重信号分类方法 |
CN106125056A (zh) * | 2016-06-13 | 2016-11-16 | 西安电子科技大学 | 基于修正因子的最小方差谱估计方法 |
CN106371079A (zh) * | 2016-08-19 | 2017-02-01 | 西安电子科技大学 | 基于谱锐化的多重信号分类谱估计方法 |
CN108459312A (zh) * | 2018-03-26 | 2018-08-28 | 西安电子科技大学 | 基于复合多径因子估计的加权多频最大似然仰角估计方法 |
CN108549059A (zh) * | 2018-03-26 | 2018-09-18 | 西安电子科技大学 | 一种复杂地形条件下的低空目标仰角估计方法 |
CN108957389A (zh) * | 2018-05-25 | 2018-12-07 | 哈尔滨工程大学 | 一种实数域多通道信号目标方位估计方法 |
CN110231617A (zh) * | 2019-05-30 | 2019-09-13 | 深圳市华讯方舟微电子科技有限公司 | 目标障碍物方位探测方法、装置、车载雷达及存储介质 |
CN110673086A (zh) * | 2019-10-31 | 2020-01-10 | 上海无线电设备研究所 | 一种基于数字阵列雷达的二维角度超分辨方法 |
CN110940973A (zh) * | 2019-12-10 | 2020-03-31 | 成都纳雷科技有限公司 | 一种用于雷达目标检测的角度测量方法及装置 |
CN113341372A (zh) * | 2021-05-31 | 2021-09-03 | 电子科技大学 | 一种基于去相干的酉Root-MUSIC算法的DOA估计方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP3681986B2 (ja) * | 2001-02-15 | 2005-08-10 | 三菱電機株式会社 | 測角装置及び測角方法 |
CN101349742A (zh) * | 2008-08-29 | 2009-01-21 | 西安电子科技大学 | 用于宽带相控阵阵元间距优化和频域多目标测频测向方法 |
CN101526995A (zh) * | 2009-01-19 | 2009-09-09 | 西安电子科技大学 | 基于对角子类判决分析的合成孔径雷达目标识别方法 |
WO2010092421A1 (en) * | 2009-02-16 | 2010-08-19 | The B.T.I. Company | Wireless communication systems and methods with source localization and self-calibration |
-
2013
- 2013-07-14 CN CN201310296530.0A patent/CN103364772B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP3681986B2 (ja) * | 2001-02-15 | 2005-08-10 | 三菱電機株式会社 | 測角装置及び測角方法 |
CN101349742A (zh) * | 2008-08-29 | 2009-01-21 | 西安电子科技大学 | 用于宽带相控阵阵元间距优化和频域多目标测频测向方法 |
CN101526995A (zh) * | 2009-01-19 | 2009-09-09 | 西安电子科技大学 | 基于对角子类判决分析的合成孔径雷达目标识别方法 |
WO2010092421A1 (en) * | 2009-02-16 | 2010-08-19 | The B.T.I. Company | Wireless communication systems and methods with source localization and self-calibration |
Cited By (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103913731B (zh) * | 2014-02-25 | 2016-04-06 | 西安电子科技大学 | 基于先验信息的目标仰角估计方法 |
CN103913731A (zh) * | 2014-02-25 | 2014-07-09 | 西安电子科技大学 | 基于先验信息的目标仰角估计方法 |
CN105913044A (zh) * | 2016-05-04 | 2016-08-31 | 大连理工大学 | 一种基于Sigmoid协方差矩阵的多重信号分类方法 |
CN105913044B (zh) * | 2016-05-04 | 2019-05-24 | 大连理工大学 | 一种基于Sigmoid协方差矩阵的多重信号分类方法 |
CN106125056A (zh) * | 2016-06-13 | 2016-11-16 | 西安电子科技大学 | 基于修正因子的最小方差谱估计方法 |
CN106125056B (zh) * | 2016-06-13 | 2018-07-06 | 西安电子科技大学 | 基于修正因子的最小方差谱估计方法 |
CN106371079A (zh) * | 2016-08-19 | 2017-02-01 | 西安电子科技大学 | 基于谱锐化的多重信号分类谱估计方法 |
CN108459312B (zh) * | 2018-03-26 | 2022-05-03 | 西安电子科技大学 | 基于复合多径因子估计的加权多频最大似然仰角估计方法 |
CN108459312A (zh) * | 2018-03-26 | 2018-08-28 | 西安电子科技大学 | 基于复合多径因子估计的加权多频最大似然仰角估计方法 |
CN108549059A (zh) * | 2018-03-26 | 2018-09-18 | 西安电子科技大学 | 一种复杂地形条件下的低空目标仰角估计方法 |
CN108549059B (zh) * | 2018-03-26 | 2022-05-03 | 西安电子科技大学 | 一种复杂地形条件下的低空目标仰角估计方法 |
CN108957389A (zh) * | 2018-05-25 | 2018-12-07 | 哈尔滨工程大学 | 一种实数域多通道信号目标方位估计方法 |
CN110231617B (zh) * | 2019-05-30 | 2021-06-18 | 深圳市华讯方舟微电子科技有限公司 | 目标障碍物方位探测方法、装置、车载雷达及存储介质 |
CN110231617A (zh) * | 2019-05-30 | 2019-09-13 | 深圳市华讯方舟微电子科技有限公司 | 目标障碍物方位探测方法、装置、车载雷达及存储介质 |
CN110673086A (zh) * | 2019-10-31 | 2020-01-10 | 上海无线电设备研究所 | 一种基于数字阵列雷达的二维角度超分辨方法 |
CN110940973A (zh) * | 2019-12-10 | 2020-03-31 | 成都纳雷科技有限公司 | 一种用于雷达目标检测的角度测量方法及装置 |
CN113341372A (zh) * | 2021-05-31 | 2021-09-03 | 电子科技大学 | 一种基于去相干的酉Root-MUSIC算法的DOA估计方法 |
Also Published As
Publication number | Publication date |
---|---|
CN103364772B (zh) | 2015-04-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103364772B (zh) | 基于实数域广义多重信号分类算法的目标低仰角估计方法 | |
CN103353596B (zh) | 基于压缩感知的波束空间域米波雷达测高方法 | |
CN106772225B (zh) | 基于压缩感知的波束域doa估计 | |
CN106526530B (zh) | 基于传播算子的2-l型阵列二维doa估计算法 | |
CN103353595B (zh) | 基于阵列内插压缩感知的米波雷达测高方法 | |
CN103901417B (zh) | L型阵列mimo雷达低复杂度空间目标二维角度估计方法 | |
CN103901395B (zh) | 一种冲击噪声环境下相干信号波达方向动态跟踪方法 | |
CN111123192B (zh) | 一种基于圆形阵列和虚拟扩展的二维doa定位方法 | |
CN103760546B (zh) | 一种雷达用低空目标波达方向估计方法 | |
CN103744076B (zh) | 基于非凸优化的mimo雷达动目标检测方法 | |
CN108535698B (zh) | 基于波束空间的米波雷达低仰角估计方法 | |
CN106021637A (zh) | 互质阵列中基于迭代稀疏重构的doa估计方法 | |
CN105259550A (zh) | 基于压缩感知的多输入多输出雷达二维角度估计方法 | |
CN105182325B (zh) | 基于秩1约束的米波mimo雷达低仰角目标测高方法 | |
CN105403874A (zh) | 非均匀阵列欠定波达方向估计方法 | |
CN111650556B (zh) | 一种宽带辐射源参数估计方法 | |
CN103364762B (zh) | 任意阵列流形的单基地mimo雷达波达方向估计方法 | |
Zhu et al. | Altitude measurement based on terrain matching in VHF array radar | |
CN108398669A (zh) | 一种基于无需预延迟处理的空时宽带自适应单脉冲测角方法 | |
CN105699950A (zh) | 基于自适应迭代前后向平滑共轭梯度的雷达杂波抑制方法 | |
CN109696657A (zh) | 一种基于矢量水听器的相干声源定位方法 | |
CN107450046B (zh) | 低仰角多径环境下的波达角估计方法 | |
CN109557504A (zh) | 一种近场窄带信号源的定位方法 | |
CN108318855A (zh) | 基于均匀圆阵的近场和远场混合信号源定位方法 | |
CN106970348B (zh) | 电磁矢量传感器阵列解相干二维music参数估计方法 |
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: 20150422 Termination date: 20210714 |
|
CF01 | Termination of patent right due to non-payment of annual fee |