CN113376601B - 基于clean算法的捷变频雷达旁瓣抑制方法 - Google Patents
基于clean算法的捷变频雷达旁瓣抑制方法 Download PDFInfo
- Publication number
- CN113376601B CN113376601B CN202110506114.3A CN202110506114A CN113376601B CN 113376601 B CN113376601 B CN 113376601B CN 202110506114 A CN202110506114 A CN 202110506114A CN 113376601 B CN113376601 B CN 113376601B
- Authority
- CN
- China
- Prior art keywords
- echo signal
- frequency
- pulse
- radar
- 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.)
- Active
Links
- 238000004422 calculation algorithm Methods 0.000 title claims abstract description 30
- 238000000034 method Methods 0.000 title claims abstract description 27
- 230000001629 suppression Effects 0.000 title claims abstract description 16
- 239000011159 matrix material Substances 0.000 claims abstract description 32
- 238000007906 compression Methods 0.000 claims abstract description 23
- 230000006835 compression Effects 0.000 claims abstract description 21
- 238000006243 chemical reaction Methods 0.000 claims abstract description 11
- 230000001427 coherent effect Effects 0.000 claims abstract description 10
- 238000009825 accumulation Methods 0.000 claims abstract description 6
- 230000002401 inhibitory effect Effects 0.000 claims abstract 2
- 230000000875 corresponding effect Effects 0.000 claims description 10
- 238000005070 sampling Methods 0.000 claims description 4
- 230000005540 biological transmission Effects 0.000 claims description 3
- 230000002596 correlated effect Effects 0.000 claims description 2
- 230000021615 conjugation Effects 0.000 claims 1
- 238000001514 detection method Methods 0.000 abstract description 7
- 238000004364 calculation method Methods 0.000 abstract description 3
- 238000004088 simulation Methods 0.000 description 7
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000001228 spectrum Methods 0.000 description 2
- 230000008859 change Effects 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
Images
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/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明属于雷达信号处理技术领域,具体公开了一种基于CLEAN算法的捷变频雷达旁瓣抑制方法,包括:建立捷变频雷达的发射信号模型和回波信号模型;对每个脉冲的回波信号依次进行下变频和脉冲压缩处理,得到脉冲压缩后的回波信号;构建捷变频雷达回波信号基矩阵,采用相关算法对脉冲压缩后的回波信号进行脉间相参积累;采用CLEAN算法抑制所述相关处理结果中的目标旁瓣,得到旁瓣抑制后的结果。本发明通过CLEAN算法抑制相关处理结果中目标的旁瓣,从而提高对弱小目标的检测概率,计算复杂度低,运算量小。
Description
技术领域
本发明涉及雷达信号处理技术领域,尤其涉及一种基于CLEAN算法的捷变频雷达旁瓣抑制方法,可应用于捷变频雷达相参积累后的目标旁瓣抑制。
背景技术
捷变频雷达是指在相邻脉冲间发射信号载频在一定带宽范围内快速跳变,各个脉冲的载频可按一定规律变化,也可随机跳变。捷变频雷达具有一系列优势:(1)强电子对抗能力。雷达载频捷变可以主动规避干扰覆盖频段,能有效对抗窄带瞄准式干扰、跨脉冲重复周期干扰及部分前拖干扰。(2)良好的目标探测能力。脉间频率捷变降低了目标长期处于雷达散射截面积衰落区的可能性,有利于目标检测概率的提高。(3)优异的电磁兼容性能。捷变频雷达在单频点驻留时间短,可以有效减少对邻近设备的干扰。此外,也可以通过频谱管理,实现设备之间的频谱共享。因此,在日益复杂的电磁环境中,捷变频雷达具有广阔的应用前景。
对于捷变频雷达距离-速度二维联合相参处理,可以采用匹配滤波方法,即相关处理,该方法的优点是计算复杂度低,计算结果比较稳定,但是积累后目标的旁瓣较高。在检测过程中,随机起伏的旁瓣不仅会掩盖弱小目标,也会导致虚警。针对旁瓣的抑制,可以优化设计载频序列或者采用基于稀疏重构理论的信号处理算法。但是这两类方法存在如下问题:(1)采用预先设计载频序列的方法来降低旁瓣水平的效果是有限的。(2)在电子对抗中,固定的载频跳变序列容易被干扰机识别,抗干扰性能较差。(3)虽然基于稀疏重构理论的信号处理算法能有效抑制旁瓣,但是算法的计算复杂度较高,稀疏重构结果不稳定,且受观测场景稀疏性以及信噪比的限制。
发明内容
针对现有技术存在的问题,本发明的目的在于提供一种基于CLEAN算法的捷变频雷达旁瓣抑制方法,通过CLEAN算法抑制相关处理结果中目标的旁瓣,从而提高对弱小目标的检测概率,计算复杂度低,运算量小。
为了达到上述目的,本发明采用以下技术方案予以实现。
基于CLEAN算法的捷变频雷达旁瓣抑制方法,包括以下步骤:
步骤1,建立捷变频雷达的发射信号模型和回波信号模型;
步骤2,对每个脉冲的回波信号依次进行下变频和脉冲压缩处理,得到脉冲压缩后的回波信号;
步骤3,构建捷变频雷达回波信号基矩阵,采用相关算法对脉冲压缩后的回波信号进行脉间相参积累,得到对应的相关处理结果;
步骤4,采用CLEAN算法抑制所述相关处理结果中的目标旁瓣,得到旁瓣抑制后的结果。
与现有技术相比,本发明的有益效果为:
(1)本发明选用捷变频雷达,相对于传统脉冲多普勒雷达,具有强电子对抗能力,良好的目标探测能力以及优异的电磁兼容性能等诸多优点;
(2)与预先设计载频序列的方法相比,本发明没有对载频捷变序列提出要求,不会导致捷变频雷达的抗干扰性能下降;
(3)相对于基于稀疏重构理论的信号处理算法,本发明不涉及复杂的矩阵求逆运算,计算复杂度低,运行速度快,且不受观测场景稀疏性的限制。
附图说明
下面结合附图和具体实施例对本发明做进一步详细说明。
图1是本发明基于CLEAN算法的捷变频雷达旁瓣抑制方法的实现流程图;
图2是本发明实施例仿真实验中的相关处理结果图;
图3是本发明实施例仿真实验中旁瓣抑制后的结果图。
具体实施方式
下面将结合实施例对本发明的实施方案进行详细描述,但是本领域的技术人员将会理解,下列实施例仅用于说明本发明,而不应视为限制本发明的范围。
参考图1,本发明提供的一种基于CLEAN算法的捷变频雷达旁瓣抑制方法,包括以下步骤:
步骤1,建立捷变频雷达的发射信号模型和回波信号模型;
其中,为快时间,tn为慢时间,为发射信号的复包络,为矩形窗函数,Tp为脉冲宽度,γ为调频率,fn=f0+dn·Δf为第n个脉冲发射信号的载频,f0为起始频率,dn∈{0,1,…,M-1}为频率调制码字,M为跳频点数,Δf为跳频间隔。
步骤2,对每个脉冲的回波信号依次进行下变频和脉冲压缩处理,得到脉冲压缩后的回波信号;
步骤3,构建捷变频雷达回波信号基矩阵,采用相关算法对脉冲压缩后的回波信号进行脉间相参积累,得到对应的相关处理结果;
3a)构建捷变频雷达距离维基矩阵R,表示为:
构建捷变频雷达速度维基矩阵D,表示为:
Λ=[diag(r1)D diag(r2)D … diag(rM)D]
yq=[spc(q,t1) spc(q,t2) … spc(q,tN)]T
其中,yq表示第q个距离单元的回波信号。
3c)将第q个距离单元的回波信号yq与捷变频雷达回波信号基矩阵Λ做相关处理,得到相关处理结果为:
sMF=ΛH·yq
其中,sMF表示相关处理结果,上标H表示矩阵的共轭转置。
步骤4,采用CLEAN算法抑制所述相关处理结果中的目标旁瓣,得到旁瓣抑制后的结果。
4a)寻找相关处理结果sMF的最大值及其对应的位置,分别记为smax和i;
其中,sMF中的每个元素从前往后按顺序进行位置标记,分别标记为1,2,…,i,…MN;
4b)将捷变频雷达回波信号基矩阵Λ中第i列与基矩阵Λ做相关,并把得到的相关结果记为sest,同时将向量sest中第i个元素置零;
sest=ΛH·Λi;
其中,Λi表示基矩阵Λ的第i列;
4c)从相关处理结果sMF中减去sest,并更新sMF为:
SMF=SMF-Sest
计算sMF的总能量:
4d)循环执行步骤4a)-步骤4c),直到相邻两次循环对应的总能量Etotal之差小于预设阈值,停止循环;当前对应的sMF为旁瓣抑制后的结果,记为ssupp。
仿真实验
本发明的效果可以通过以下具体实例进一步说明:
1、仿真条件:
捷变频雷达发射波形采用线性调频信号,信号带宽为12MHz,脉宽为4us,一个相参处理间隔内的脉冲数为64个,雷达跳频的起始频率为12GHz,调频间隔为9MHz,调频点数为96,脉冲重复频率为25KHz。假设雷达观测场景中存在两个运动点目标,与雷达的径向距离分别为5756.5m和5752.6m,与雷达的径向速度分别为58.6m/s和96.2m/s,并且两个目标的后向散射系数分别为1和0.4。
2、仿真内容:
在上述仿真参数下,采用本发明基于CLEAN算法的捷变频雷达旁瓣抑制方法,对捷变频雷达相关处理结果中目标的旁瓣进行抑制。仿真结果如图2和图3所示,图2为相关处理结果,从图2中可以看出,相关处理只能检测到后向散射系数为1的目标,而后向散射系数为0.4的目标被另外一个目标的旁瓣所掩盖,无法准确检测。图3为本发明的旁瓣抑制后的结果,从图3中可以看出,后向散射系数为0.4的目标也能被检测。可见,本发明能够有效抑制捷变频雷达相关处理结果中的目标旁瓣。
以上仿真实验验证了本发明的正确性、有效性。
虽然,本说明书中已经用一般性说明及具体实施方案对本发明作了详尽的描述,但在本发明基础上,可以对之作一些修改或改进,这对本领域技术人员而言是显而易见的。因此,在不偏离本发明精神的基础上所做的这些修改或改进,均属于本发明要求保护的范围。
Claims (5)
1.基于CLEAN算法的捷变频雷达旁瓣抑制方法,其特征在于,包括以下步骤:
步骤1,建立捷变频雷达的发射信号模型和回波信号模型;
步骤2,对每个脉冲的回波信号依次进行下变频和脉冲压缩处理,得到脉冲压缩后的回波信号;
步骤3,构建捷变频雷达回波信号基矩阵,采用相关算法对脉冲压缩后的回波信号进行脉间相参积累,得到对应的相关处理结果;
步骤4,采用CLEAN算法抑制所述相关处理结果中的目标旁瓣,得到旁瓣抑制后的结果;
所述采用CLEAN算法抑制所述相关处理结果中的目标旁瓣,具体步骤为:
4a)寻找相关处理结果sMF的最大值及其对应的位置,分别记为smax和i;
其中,sMF中的每个元素从前往后按顺序进行位置标记,分别标记为1,2,L,i,L MN;
4b)将捷变频雷达回波信号基矩阵Λ中第i列与基矩阵Λ做相关,并把得到的相关结果记为sest,同时将向量sest中第i个元素置零;
sest=ΛH·Λi;
其中,Λi表示基矩阵Λ的第i列;
4c)从相关处理结果sMF中减去sest,并更新sMF为:
sMF=sMF-sest
计算sMF的总能量:
4d)循环执行步骤4a)-4c),直到相邻两次循环对应的总能量Etotal之差小于预设阈值,停止循环;则当前对应的sMF为旁瓣抑制后的结果,记为ssupp。
2.根据权利要求1所述的基于CLEAN算法的捷变频雷达旁瓣抑制方法,其特征在于,步骤1包含以下子步骤:
其中,j表示虚数单位,为快时间,tn为慢时间,为发射信号的复包络,为矩形窗函数,Tp为脉冲宽度,γ为调频率,fn=f0+dn·Δf为第n个脉冲发射信号的载频,f0为起始频率,dn∈{0,1,L,M-1}为频率调制码字,M为跳频点数,Δf为跳频间隔;
4.根据权利要求1所述的基于CLEAN算法的捷变频雷达旁瓣抑制方法,其特征在于,所述构建捷变频雷达回波信号基矩阵,具体为:
首先,构建捷变频雷达距离维基矩阵R,表示为:
其次,构建捷变频雷达速度维基矩阵D,表示为:
最后,通过捷变频雷达距离维基矩阵R和速度维基矩阵D,构建回波信号基矩阵Λ:
Λ=[diag(r1)D diag(r2)D L diag(rM)D]
5.根据权利要求2所述的基于CLEAN算法的捷变频雷达旁瓣抑制方法,其特征在于,所述采用相关算法对脉冲压缩后的回波信号进行脉间相参积累,具体过程为:
则第q个距离单元的回波信号表示为:
yq=[spc(q,t1) spc(q,t2) L spc(q,tN)]T
其中,上标T表示转置操作,yq表示第q个距离单元的回波信号;
然后,将第q个距离单元的回波信号yq与捷变频雷达回波信号基矩阵Λ做相关处理,得到相关处理结果为:
sMF=ΛH·yq
其中,sMF表示相关处理结果,上标H表示矩阵的共轭转置。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110506114.3A CN113376601B (zh) | 2021-05-10 | 2021-05-10 | 基于clean算法的捷变频雷达旁瓣抑制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110506114.3A CN113376601B (zh) | 2021-05-10 | 2021-05-10 | 基于clean算法的捷变频雷达旁瓣抑制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113376601A CN113376601A (zh) | 2021-09-10 |
CN113376601B true CN113376601B (zh) | 2022-11-01 |
Family
ID=77572452
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110506114.3A Active CN113376601B (zh) | 2021-05-10 | 2021-05-10 | 基于clean算法的捷变频雷达旁瓣抑制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113376601B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114578351B (zh) * | 2022-02-28 | 2023-09-29 | 北京城建设计发展集团股份有限公司 | 一种雷达、雷达旁瓣抑制方法及*** |
CN114325633B (zh) * | 2022-03-16 | 2022-06-17 | 西安电子科技大学 | 一种双波段雷达的相参处理方法及处理装置 |
CN116027280B (zh) * | 2023-03-30 | 2023-06-09 | 西安电子科技大学 | 一种低峰值旁瓣频率编码雷达波形设计方法 |
CN116366093B (zh) * | 2023-06-01 | 2023-08-25 | 南京邮电大学 | 分块捷变跳频方法 |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5726657A (en) * | 1996-03-22 | 1998-03-10 | Lockheed Martin Corporation | Phase coherent radar system using fast frequency agile waveform synthesis |
US7042386B2 (en) * | 2001-12-11 | 2006-05-09 | Essex Corporation | Sub-aperture sidelobe and alias mitigation techniques |
CN102033227B (zh) * | 2010-11-30 | 2013-01-16 | 哈尔滨工程大学 | 以gps导航卫星为外辐射源的无源雷达微弱目标检测方法 |
CN104459667B (zh) * | 2014-12-01 | 2017-02-01 | 西安电子工程研究所 | 一种基于clean的稀疏阵列波达方向doa估计方法 |
CN111781574B (zh) * | 2020-07-13 | 2022-04-29 | 西安电子科技大学 | 基于子空间正交投影的捷变频雷达目标参数估计方法 |
CN112014818B (zh) * | 2020-08-25 | 2022-09-16 | 北京理工大学 | 一种基于批处理clean的杂波抑制方法 |
CN112346019B (zh) * | 2020-09-11 | 2023-05-16 | 中国人民解放军国防科技大学 | 噪声雷达脉冲波形和低旁瓣脉冲压缩相参积累处理方法 |
-
2021
- 2021-05-10 CN CN202110506114.3A patent/CN113376601B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN113376601A (zh) | 2021-09-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113376601B (zh) | 基于clean算法的捷变频雷达旁瓣抑制方法 | |
Colone et al. | Cancellation of clutter and multipath in passive radar using a sequential approach | |
CN111693983B (zh) | 互补波形构建方法及模块、认知雷达***及波形发射方法 | |
CN113721216B (zh) | 一种捷变相参雷达的目标检测波形优化与处理方法 | |
CN107192991B (zh) | 一种射频辐射式仿真宽带脉冲雷达成像等效模拟方法 | |
CN112881982B (zh) | 一种捷变频雷达抑制密集转发干扰方法 | |
CN113640761B (zh) | 一种频率捷变雷达抗脉间假目标的信号处理方法 | |
CN106199549B (zh) | 一种采用谱减法提升lfmcw雷达信噪比的方法 | |
CN112436905B (zh) | 一种通信雷达联合*** | |
CN113884992A (zh) | 一种频率捷变雷达的自适应抗干扰方法 | |
CN116930907A (zh) | 一种ddma-mimo雷达的目标相参积累检测方法 | |
CN116953683A (zh) | 基于脉内波形认知优化设计的间歇采样转发干扰抑制方法 | |
CN110109075B (zh) | 基于白化滤波的捷变频雷达抗干扰方法 | |
CN116819458A (zh) | 一种多站频率捷变雷达抗干扰及协同增效方法 | |
CN112014807B (zh) | 一种频率捷变雷达的自适应杂波抑制方法 | |
CN109061626B (zh) | 一种步进频相参处理检测低信杂比动目标的方法 | |
Li et al. | Random-frequency-Coded waveform optimization and signal coherent accumulation against compound deception jamming | |
Elgamel et al. | Empirical mode decomposition-based monopulse processor for enhanced radar tracking in the presence of high-power interference | |
De Martín et al. | Sidelobe mitigation in noise radar using sparse signal processing | |
CN113640754B (zh) | 一种基于低秩稀疏模型的抗主瓣间歇采样转发干扰方法 | |
CN113376607B (zh) | 机载分布式雷达小样本空时自适应处理方法 | |
US9035824B2 (en) | System and method of radar location | |
CN111239698A (zh) | 基于互补离散混沌相位编码波形设计的拖曳式干扰抑制方法 | |
CN112881984B (zh) | 一种雷达信号抗干扰处理方法、装置及存储介质 | |
CN116449313B (zh) | 基于脉内多载频信号雷达抗主瓣压制噪声干扰方法与装置 |
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 |