[发明专利]一种基于形态分量分析的外辐射源雷达风场杂波抑制方法

摘要

本发明公开了一种基于形态分量分析的外辐射源雷达风场杂波抑制方法，利用雷达接收回波中目标信号与风场杂波具有不同的形态特性，且都可在相应的变换域进行稀疏表示，从而可利用形态分量分析的方法对信号进行分离；本发明的方法使风场杂波得到较好抑制，同时具有稳定性好、计算简便等优点，提高了外辐射源雷达目标检测性能。

主权项

一种基于形态分量分析的外辐射源雷达风场杂波抑制方法，其特征在于，包括下述步骤：步骤1，对参考通道、监测通道分别进行参考信号重构和多径杂波抑制，然后计算参考通道与监测通道二维互相关函数；步骤2，是否存在沿多普勒轴与零多普勒对称的、呈周期出现的多次谐波副峰；如果存在，则需进行风场杂波抑制；如果不存在，则表明无风场杂波；步骤3，当存在沿多普勒轴与零多普勒对称的呈周期出现的多次谐波副峰时，通过求解以下约束优化问题得到稀疏系数和{c^1,c^2}=argminc1,c2λ1||c1||1+λ2||c2||1]]>s.t.y＝A1c1+A2c2其中，c1和c2是稀疏系数；和分别是c1和c2的最佳估计；λ1、λ2为参数，且λ1＝1‑λ2；y为雷达接收到的复合信号；A1、A2分别为傅里叶变换及短时傅里叶变换；然后根据估计的得到信号分量从而实现信号分离；所述的步骤3中，最佳稀疏系数和通过以下子步骤进行求解：步骤3.1，初始化变量：迭代次数初始值k＝0，惩罚参数μ＞0，待求解稀疏系数初始值c1,0、c2,0，迭代过程中间变量初始值d1,0、d2,0；步骤3.2，更新变量：vi,k+1＝soft(ci,k+di,k,λi/μ)‑di,k，i＝1,2gk+1＝y‑A1v1,k+1‑A2v2,k+1di,k+1=12AiHgk+1,i=1,2]]>ci,k+1＝di,k+1+vi,k+1，i＝1,2其中，soft(.)为软阈值函数，soft(x,T)＝max(1‑T/|x|,0)，x,T分别软阈值函数中的变量，|x|为变量x的绝对值；vi,k+1、gk+1、di,k+1、ci,k+1分别表示第i个信号分量在进行第k+1次迭代时的中间变量；A1、A2分别为傅里叶变换及短时傅里叶变换；H为转置复共轭；步骤3.3，判断是否满足终止条件，满足则停止；否则，k＝k+1。