CN109725297A - 一种基于回波极化度的有源转发式干扰识别方法 - Google Patents

一种基于回波极化度的有源转发式干扰识别方法 Download PDF

Info

Publication number
CN109725297A
CN109725297A CN201910048306.7A CN201910048306A CN109725297A CN 109725297 A CN109725297 A CN 109725297A CN 201910048306 A CN201910048306 A CN 201910048306A CN 109725297 A CN109725297 A CN 109725297A
Authority
CN
China
Prior art keywords
echo
signal
polarizability
polarization
active
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
Application number
CN201910048306.7A
Other languages
English (en)
Other versions
CN109725297B (zh
Inventor
刘加方
张云华
李东
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
National Space Science Center of CAS
Original Assignee
National Space Science Center of CAS
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by National Space Science Center of CAS filed Critical National Space Science Center of CAS
Priority to CN201910048306.7A priority Critical patent/CN109725297B/zh
Publication of CN109725297A publication Critical patent/CN109725297A/zh
Application granted granted Critical
Publication of CN109725297B publication Critical patent/CN109725297B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明公开了一种基于回波极化度的有源转发式干扰识别方法,所述方法包括:通过在雷达的两正交极化通道分别馈入两个去相干的噪声雷达信号来获得随机极化波并进行发射;雷达的两正交极化接收通道接收回波信号;计算回波信号的极化度;根据回波信号的极化度的差异进行真实目标和有源转发式干扰的识别。本发明的方法通过发射随机极化波,在回波域利用散射回波的极化度可实现对有源转发式干扰和真实目标的有效识别;本发明的方法过程简单,计算量小,适用于当前绝大多数有源转发式干扰。

Description

一种基于回波极化度的有源转发式干扰识别方法
技术领域
本发明涉及雷达极化抗干扰技术领域,特别涉及一种基于回波极化度的有源转发式干扰识别方法。
背景技术
随着雷达技术的快速发展,雷达对抗及电子战技术也得到不断进步,雷达所面临的电磁环境也日益恶劣,敏捷的电子侦查与先进的电子干扰技术使得雷达生存环境日益严峻。有源转发式干扰是一种重要的欺骗干扰,通过模拟目标特征并发射假目标回波信号来欺骗雷达,使雷达出现密集假目标扰乱雷达对真实目标探测,或对雷达实现角度欺骗以扰乱跟踪制导雷达对目标的跟踪。随着数字射频存储器(Digital Radio Frequency Memory,DRFM)、数字干扰合成技术(Digital Jamming Synthesize,DJS)等先进电子干扰技术的发展,雷达有源转发式干扰表现出高逼真度及智能化等特征,DRFM可以精确的截获、存储并转发敌方雷达信号,高逼真度的假目标不仅可以获得与真实目标回波相同的处理增益,而且在在时域、频域和空域的特征与真实目标特性十分相似,传统抗干扰方法已经无法识别出真假目标。
极化度(Degree of Polarization,DoP)是衡量电磁波极化程度的参数。实际的雷达目标与入射极化波相互作用的过程就是对其极化调制的过程。如果入射电磁波为完全极化波,那么目标的散射波也将是完全极化波,只是极化方式发生过了改变,例如由线极化变为圆极化,但其回波极化度不发生改变;若入射电磁波为随机极化波,那么目标为一个起偏器,其会改善电磁波的极化度,但极化度仍然小于1。因此在发射随机极化波的情况下,转发式有源干扰的极化状态由干扰机发射天线决定,通常为圆极化或线极化天线,转发的干扰回波信号为完全极化波,而目标的回波为部分极化波,极化度小于1,可以通过对回波极化度的测量鉴别有源转发式干扰。
发明内容
本发明目的在于提供一种基于回波极化度的有源转发式干扰鉴别方法。结合噪声雷达信号与极化雷达技术,通过在雷达正交极化天线发射两个去相干的噪声雷达信号来获得随机极化波,利用真实目标散射回波极化度与有源转发式干扰极化度的差别,能够实现对真实目标与有源转发式干扰的鉴别。
本发明的技术方案如下:
一种基于回波极化度的有源转发式干扰识别方法,所述方法包括:
通过在雷达的两正交极化通道分别馈入两个去相干的噪声雷达信号来获得随机极化波并进行发射;
雷达的两正交极化接收通道接收回波信号;
计算回波信号的极化度;根据回波信号的极化度的差异进行真实目标和有源转发式干扰的识别。
作为上述方法的一种改进,所述两个去相干的噪声雷达信号分别为:雷达的水平极化通道发射的噪声雷达信号sH(t)、雷达的垂直极化通道发射的噪声雷达信号sV(t);sH(t)和sV(t)是正交的,t为时间。
作为上述方法的一种改进,所述回波信号为目标回波信号sr和干扰回波信号sj的和;
所述目标回波信号sr为:
其中,目标极化散射矩阵为snh,snv分别为水平接收通道和垂直接收通道的噪声;
干扰回波信号sj为:
hjr为干扰接收极化Jones矢量,hjt为干扰转发式天线极化Jones矢量,G为有源转发式干扰信号转发增益,
作为上述方法的一种改进,所述回波信号的极化度的计算过程为:
其中,P为回波信号的极化度,J=[g0,g1,g2,g3]T为Stokes矢量;
在水平垂直极化基下有:
其中,aH为水平接收通道的信号强度,aV为垂直接收通道的信号强度;电磁波水平分量和垂直分量相位差θ=θhvh为电磁波水平分量的相位,θv为电磁波垂直分量的相位。
作为上述方法的一种改进,所述识别的过程为:设定极化度检测阈值Pe;当P≥Pe时,判断为有源干扰假目标,否则,为真实目标。
本发明的优势在于:
1、本发明的方法通过发射随机极化波,在回波域利用散射回波的极化度可实现对有源转发式干扰和真实目标的有效识别;
2、本发明的方法过程简单,计算量小,适用于当前绝大多数有源转发式干扰;
3、本发明的方法采用的随机极化波在时域和极化域都呈现良好的随机性,具有良好的低截获/低检测(LPI/LPD)特性,提高了雷达在复杂电磁环境下的生存能力。
附图说明
图1为本发明基于回波极化度的有源转发式干扰识别方法的流程图;
图2为本发明目标回波极化状态随时间变化;
图3为本发明目标回波极化状态在庞加莱极化球上分布;
图4为本发明目标回波极化度;
图5为本发明信噪比0dB情况下目标回波信号极化状态随时间变化;
图6为本发明不同信噪比下目标回波极化状态在庞加莱极化球上分布;
图7为本发明目标回波极化度随信噪比的变化;
图8为本发明转发式干扰回波极化状态在庞加莱极化球上分布;
图9为本发明干扰噪声比0dB情况下干扰信号极化状态随时间变化;
图10为本发明不同干扰噪声比下干扰信号极化状态在庞加莱极化球上分布;
图11为本发明干扰信号极化度随干扰噪声比的变化。
具体实施方式
下面结合附图和实施例对本发明的进行详细说明。
本发明提出一种基于回波极化度的有源转发式干扰识别方法,在发射随机极化波的情况下,利用散射回波的极化度可实现对有源转发式干扰和真实目标的有效识别。具体包括:
步骤1)采用两正交极化发射天线,分别设为H(Horizontal,水平)极化和V(Vertical,垂直)极化,雷达的H极化通道、V极化通道同时发射去相干波形矢量sH(t)、sV(t)为两正交极化通道发射的噪声雷达信号,采样率为1000MHz,脉冲宽度为20us,带宽为600MHz,发射随机极化波的极化度为0.005;
步骤2)雷达两正交极化接收通道分别对散射回波信号和有源转发式干扰信号进行接收;
我方采用随机极化波照射敌方目标,敌方电子对抗系统采用有源转发式干扰,有源干扰工作方式为直接放大转发。设目标在某姿态角下的散射矩阵为为考虑一般性,我们假定目标极化散射矩阵为任意一个随机的散射矩阵,真实环境中都是有各种各样噪声存在,雷达接收机接收到的回波除了目标回波外还有各类噪声,发射随机极化波经过目标散射后,雷达两正交极化通道接收到的回波信号为:
snh,snv分别为H,V接收通道噪声
干扰信号的回波可表示为:
hjr为干扰接收极化Jones矢量,hjt为干扰转发式天线极化Jones矢量,G为有源转发式干扰信号转发增益,sn为H,V接收通道噪声矢量。假设干扰机接收天线为圆极化天线,转发天线为任意线极化天线,那么
步骤3)计算接收回波信号的极化度。
由Stokes矢量J=[g0,g1,g2,g3]T的定义式,在水平垂直极化基下有:
电磁波Stokes矢量g0分量描述了电磁波的功率密度,其余三个元素所构成的子矢量表征了电磁波的极化状态,极化度定义为完全极化波强度与部分极化波强度之比,即
考虑一般性,随机给目标极化散射矩阵500个不同的散射矩阵,从图2可以看出目标回波极化随时间变化,为部分极化波。将500次不同的散射矩阵情况下计算得到目标回波Stokes矢量投影到Poincare极化球上,如图3所示,目标回波的Stokes矢量都位于Poincare极化球内,图4给出来500次不同散射矩阵情况下目标回波信号极化度,从图中可以看出,目标回波的极化度都低于0.6,说明真实目标回波是部分极化波。为进一步考虑系统噪声对目标回波极化度测量的影响,我们假设目标散射矩阵图5给出了信噪比为0dB情况下目标回波极化状态,可以看出目标回波为部分极化波,极化度小于1,图6给出了不同信噪比下目标回波极化状态在庞加莱极化球上分布,从图中可以看出,随着信噪比的降低,目标信号Stokes矢量渐渐靠近庞加莱极化球球心。图7给出了目标回波极化度随信噪比的变化曲线,随着信噪比的降低,目标回波极化度渐渐降低。
考虑一般性,当干扰机接收天线为圆极化天线,而发射天线为随机的线极化天线,随机给出500种线极化状态,将500次不同线极化天线转发状态下的干扰回波Stokes矢量投影到Poincare极化球如图8所示,干扰回波的Stokes矢量都位于Poincare极化球球面上,说明有源转发式干扰信号是完全极化波,干扰回波的极化度都为1。为进一步考虑系统噪声对有源转发式干扰信号极化度测量的影响,我们假设转发线极化天线为斜45°线极化天线,图9给出了干扰噪声比为0dB情况下有源转发式干扰信号极化状态,可以看出有源转发式干扰信号为部分极化波,极化度小于1,图10给出了不同信噪比下目标回波极化状态在庞加莱极化球上分布,从图中可以看出,随着信噪比的降低,有源干扰信号Stokes矢量渐渐从庞加莱极化球球面靠近庞加莱极化球球心。图11给出了目标回波极化度随信噪比的变化曲线,随着干扰噪声比的降低,有源干扰信号极化度渐渐降低,这也符合真实场景,随着噪声增强,两正交极化通道之间去相关性增大,回波的极化度也随着噪声的增强而降低。
步骤4)设定合适的极化度鉴别阈值。通过对目标回波和干扰回波极化度的分析可知,当无噪声时,干扰回波为完全极化波,极化度为1,而由随机极化波照射的目标回波为部分极化波,极化度介于[0,1]之间。而当有噪声时,目标信号回波和有源干扰信号极化度都随着噪声的增强而降低,我们可以通过设置合适的检测阈值,根据测量得到的回波极化度P的大小来检验是真实目标回波信号还是有源转发式干扰信号,即:
其中Pe为检测阈值。
通常真实场景下,目标回波信噪比一般为0dB-20dB左右,而有源干扰信号通常干扰功率较大,干扰噪声比一般都远大于10dB,结合图4、图7和图11分析可知,将极化度检测阈值设置为0.8可以很好的鉴别出是真实目标还是有源转发式干扰。
最后所应说明的是,以上实施例仅用以说明本发明的技术方案而非限制。尽管参照实施例对本发明进行了详细说明,本领域的普通技术人员应当理解,对本发明的技术方案进行修改或者等同替换,都不脱离本发明技术方案的精神和范围,其均应涵盖在本发明的权利要求范围当中。

Claims (5)

1.一种基于回波极化度的有源转发式干扰识别方法,所述方法包括:
通过在雷达的两正交极化通道分别馈入两个去相干的噪声雷达信号来获得随机极化波并进行发射;
雷达的两正交极化接收通道接收回波信号;
计算回波信号的极化度;根据回波信号的极化度的差异进行真实目标和有源转发式干扰的识别。
2.根据权利要求1所述的基于回波极化度的有源转发式干扰识别方法,其特征在于,所述两个去相干的噪声雷达信号分别为:雷达的水平极化通道发射的噪声雷达信号sH(t)、雷达的垂直极化通道发射的噪声雷达信号sV(t);sH(t)和sV(t)相互正交,t为时间。
3.根据权利要求2所述的基于回波极化度的有源转发式干扰识别方法,其特征在于,所述回波信号为目标回波信号sr和干扰回波信号sj的和;
所述目标回波信号sr为:
其中,目标极化散射矩阵为snh,snv分别为水平接收通道和垂直接收通道的噪声;
干扰回波信号sj为:
hjr为干扰接收极化Jones矢量,hjt为干扰转发式天线极化Jones矢量,G为有源转发式干扰信号转发增益,
4.根据权利要求3所述的基于回波极化度的有源转发式干扰识别方法,其特征在于,所述回波信号的极化度的计算过程为:
其中,P为回波信号的极化度,J=[g0,g1,g2,g3]T为Stokes矢量;
在水平垂直极化基下有:
其中,aH为水平接收通道的信号强度,aV为垂直接收通道的信号强度;电磁波水平分量和垂直分量相位差θ=θhv,其中,θh为电磁波水平分量的相位,θv为电磁波垂直分量的相位。
5.根据权利要求4所述的基于回波极化度的有源转发式干扰识别方法,其特征在于,所述识别的过程为:设定极化度检测阈值Pe;当P≥Pe时,判断为有源干扰假目标,否则,为真实目标。
CN201910048306.7A 2019-01-18 2019-01-18 一种基于回波极化度的有源转发式干扰识别方法 Expired - Fee Related CN109725297B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910048306.7A CN109725297B (zh) 2019-01-18 2019-01-18 一种基于回波极化度的有源转发式干扰识别方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910048306.7A CN109725297B (zh) 2019-01-18 2019-01-18 一种基于回波极化度的有源转发式干扰识别方法

Publications (2)

Publication Number Publication Date
CN109725297A true CN109725297A (zh) 2019-05-07
CN109725297B CN109725297B (zh) 2021-07-30

Family

ID=66298396

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910048306.7A Expired - Fee Related CN109725297B (zh) 2019-01-18 2019-01-18 一种基于回波极化度的有源转发式干扰识别方法

Country Status (1)

Country Link
CN (1) CN109725297B (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113589248A (zh) * 2021-07-14 2021-11-02 中国科学院国家空间科学中心 一种随机极化雷达系统
CN114637001A (zh) * 2022-05-17 2022-06-17 艾索信息股份有限公司 雷达的目标检测方法、装置、设备及存储介质
CN116930905A (zh) * 2023-09-07 2023-10-24 中国人民解放军海军工程大学 单通道变极化雷达的全极化测量方法及系统
CN117452342A (zh) * 2023-12-26 2024-01-26 湖南力研光电科技有限公司 一种基于极化特征的箔条干扰检测方法

Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0701144A1 (en) * 1994-09-06 1996-03-13 Rockwell International Corporation Radar terrain bounce jamming detector
US20070152867A1 (en) * 2006-01-03 2007-07-05 Mitch Randall Polarization and frequency diverse radar system for complete polarimetric characterization of scatterers with increased scanning speed
UA89945U (ru) * 2013-10-14 2014-05-12 Запорізький Національний Технічний Університет Способ помехозащиты когерентно-импульсной рлс от комбинированных помех
CN103941239A (zh) * 2014-05-08 2014-07-23 西安电子科技大学 网络化雷达极化对消抑制主瓣压制式干扰方法
CN104215938A (zh) * 2014-09-30 2014-12-17 中国人民解放军63892部队 一种对地形反弹干扰的极化检测识别方法和装置
GB2517661A (en) * 1995-10-24 2015-03-04 Thomson Csf An anti-jamming antenna
CN104765020A (zh) * 2015-04-12 2015-07-08 西安电子科技大学 有源假目标干扰的极化鉴别方法
CN106125053A (zh) * 2016-06-17 2016-11-16 中国电子科技集团公司第十研究所 脉冲多普勒雷达极化抗干扰方法
CN106249212A (zh) * 2016-08-23 2016-12-21 中国人民解放军国防科学技术大学 主瓣压制干扰背景下有源假目标的极化鉴别方法
CN106338715A (zh) * 2016-08-23 2017-01-18 中国人民解放军国防科学技术大学 一种全极化有源假目标的鉴别方法
CN106918804A (zh) * 2016-11-27 2017-07-04 西南电子技术研究所(中国电子科技集团公司第十研究所) 雷达干扰识别信号处理系统
CN108896969A (zh) * 2018-07-20 2018-11-27 北京遥感设备研究所 一种基于极化对消的抗有源诱饵干扰系统及其方法

Patent Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0701144A1 (en) * 1994-09-06 1996-03-13 Rockwell International Corporation Radar terrain bounce jamming detector
GB2517661A (en) * 1995-10-24 2015-03-04 Thomson Csf An anti-jamming antenna
US20070152867A1 (en) * 2006-01-03 2007-07-05 Mitch Randall Polarization and frequency diverse radar system for complete polarimetric characterization of scatterers with increased scanning speed
UA89945U (ru) * 2013-10-14 2014-05-12 Запорізький Національний Технічний Університет Способ помехозащиты когерентно-импульсной рлс от комбинированных помех
CN103941239A (zh) * 2014-05-08 2014-07-23 西安电子科技大学 网络化雷达极化对消抑制主瓣压制式干扰方法
CN104215938A (zh) * 2014-09-30 2014-12-17 中国人民解放军63892部队 一种对地形反弹干扰的极化检测识别方法和装置
CN104765020A (zh) * 2015-04-12 2015-07-08 西安电子科技大学 有源假目标干扰的极化鉴别方法
CN106125053A (zh) * 2016-06-17 2016-11-16 中国电子科技集团公司第十研究所 脉冲多普勒雷达极化抗干扰方法
CN106249212A (zh) * 2016-08-23 2016-12-21 中国人民解放军国防科学技术大学 主瓣压制干扰背景下有源假目标的极化鉴别方法
CN106338715A (zh) * 2016-08-23 2017-01-18 中国人民解放军国防科学技术大学 一种全极化有源假目标的鉴别方法
CN106918804A (zh) * 2016-11-27 2017-07-04 西南电子技术研究所(中国电子科技集团公司第十研究所) 雷达干扰识别信号处理系统
CN108896969A (zh) * 2018-07-20 2018-11-27 北京遥感设备研究所 一种基于极化对消的抗有源诱饵干扰系统及其方法

Non-Patent Citations (8)

* Cited by examiner, † Cited by third party
Title
HUANYAO DAI; XUESONG WANG; YONGZHEN LI; YONG LIU; SHUNPING XIAO: "《Main-Lobe Jamming Suppression Method of using Spatial Polarization Characteristics of Antennas》", 《 IEEE TRANSACTIONS ON AEROSPACE AND ELECTRONIC SYSTEMS》 *
JAE-WON RIM ; IL-SUEK KOH ; SUENG-HO CHOI: "《A3-Jamming performance analysis for repeater-type active decoy against ground tracking radar considering dynamics of platform and decoy》", 《 2017 18TH INTERNATIONAL RADAR SYMPOSIUM (IRS)》 *
宋广; 徐振海: "《部分极化干扰条件下极化阵列性能分析》", 《火控雷达技术》 *
张明胜: "《极化信息与变极化技术的研究》", 《中国优秀硕士学位论文全文数据库 信息科技辑》 *
徐建平: "《基于极化的雷达抗干扰技术研究》", 《中国优秀硕士学位论文全文数据库 信息科技辑》 *
施龙飞: "《雷达极化抗干扰技术研究》", 《中国博士学位论文全文数据库 信息科技辑》 *
李永祯 等: "《有源假目标的高分辨极化鉴别研究》", 《兵工学报》 *
王睿智: "《末制导雷达极化抗干扰技术研究》", 《中国优秀硕士学位论文全文数据库 信息科技辑》 *

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113589248A (zh) * 2021-07-14 2021-11-02 中国科学院国家空间科学中心 一种随机极化雷达系统
CN114637001A (zh) * 2022-05-17 2022-06-17 艾索信息股份有限公司 雷达的目标检测方法、装置、设备及存储介质
CN116930905A (zh) * 2023-09-07 2023-10-24 中国人民解放军海军工程大学 单通道变极化雷达的全极化测量方法及系统
CN116930905B (zh) * 2023-09-07 2023-12-19 中国人民解放军海军工程大学 单通道变极化雷达的全极化测量方法及系统
CN117452342A (zh) * 2023-12-26 2024-01-26 湖南力研光电科技有限公司 一种基于极化特征的箔条干扰检测方法
CN117452342B (zh) * 2023-12-26 2024-03-19 湖南力研光电科技有限公司 一种基于极化特征的箔条干扰检测方法

Also Published As

Publication number Publication date
CN109725297B (zh) 2021-07-30

Similar Documents

Publication Publication Date Title
CN109725297A (zh) 一种基于回波极化度的有源转发式干扰识别方法
US11714199B2 (en) System and method for detecting spoofing of GNSS signals
Daneshmand et al. A GNSS structural interference mitigation technique using antenna array processing
CN113534067B (zh) 一种多干扰源下多站雷达抗欺骗式干扰方法及系统
Ma et al. Angle estimation of extended targets in main-lobe interference with polarization filtering
Lukacs et al. “RF-DNA” fingerprinting for antenna classification
Xu et al. Research on active jamming recognition in complex electromagnetic environment
Earp et al. Ultra-wideband ground-penetrating radar for the detection of buried metallic mines
Sangston et al. New results on coherent radar target detection in heavy-tailed compound-Gaussian clutter
CN113484838B (zh) 一种多基地雷达有源假目标鉴别方法及系统
Pearce et al. Sharpening sidescan sonar images for shallow-water target and habitat classification with a vertically stacked array
Lukacs et al. Device identification using active noise interrogation and RF-DNA" fingerprinting" for non-destructive amplifier acceptance testing
Son et al. Passive polarimetrie multistatic radar for ground moving target
Tao et al. Characterization and mitigation of radio frequency interference in PolSAR data
Xiang et al. Statistical method with dual-polarized MIMO array for target discrimination
CN112859015B (zh) 一种结构新颖的全极化雷达诱饵装置及方法
Heimann et al. Non-asymptotic performance bounds of eigenvalue based detection of signals in non-Gaussian noise
Zhao et al. Backdoor Attacks Against Deep Learning-Based Massive MIMO Localization
Shi et al. Entropy-based multipath detection model for MIMO radar
CN115267696B (zh) 一种抑制雷达多径干扰的方法
Zhou et al. Adaptive two-step Bayesian generalized likelihood ratio test algorithm for low-altitude detection
Abdullah et al. Target prediction in forward scattering radar
Arisesa et al. The effect of external interference on the performance of navigation radar
Tao et al. Effects of Interference on Synthetic Aperture Radar Measurements: An Illustrative Example
Zhang et al. Distance estimation and pairing for low angle group targets based on radar time reversal

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for 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: 20210730

Termination date: 20220118

CF01 Termination of patent right due to non-payment of annual fee