[发明专利]基于IFFT频谱相除法和梯度自适应格型滤波的直达波抑制方法

权利要求书

1.基于IFFT频谱相除法和梯度自适应格型滤波的直达波抑制方法，其特征在于：该方法的步骤如下：

步骤(1)、频谱相除并进行傅里叶逆变换，估计直达波的时间延迟和幅度；在已知信号形式的前提下，采用IFFT频谱相除法对其进行时间延迟和幅度估计；为了除去噪声的影响，在频域采用加窗函数的方法，将接收信号的傅里叶变换除以归一化后发射信号的傅里叶变换可近似得到：

$DIV\left(f\right)\·L\left(f\right)=\frac{X\left(f\right)}{X_0\left(f\right)}\·L\left(f\right)\≈A_d{e}^{-{\mathrm{jt}}_{d}f}+A_s{e}^{-{\mathrm{jt}}_{ds}f}+\underset{k=1}{\overset{M}{\Σ}}{A}_{tk}{e}^{-{\mathrm{jt}}_{dtk}f}---\left(1\right)$

其中X(f)表示接收信号的傅里叶变换；X₀(f)表示归一化后发射信号的傅里叶变换；A_d、A_s、A_tk分别表示直达波信号、其他干扰信号、目标回波信号的幅度；t_d、t_ds、t_dtk分别表示直达波信号、其他干扰信号、目标回波信号到达接收机的时间延迟；M表示目标个数；L(f)表示采用的矩形窗；然后对DIV(f)进行傅里叶逆变换，就可得到：

$F^{-1}\[DIV\left(f\right)\]=A_d\mathrm{\δ}(t-t_d)+A_s\mathrm{\δ}(t-t_{ds})+\underset{k=1}{\overset{M}{\Σ}}{A}_{tk}\mathrm{\δ}(t-t_{dtk})---\left(2\right)$

其中DIV(f)的逆傅里叶变换为一组冲激函数，分别代表接收信号中不同成分到达接收机的时延脉冲；对逆傅里叶变换的结果进行峰值检测，进而对接收信号中的直达波信号的时间延迟t_d以及幅度A_d进行估计；

步骤(2)、生成参考信号；根据直达波信号的时间延迟和幅度估计结果，生成直达波信号x_d(t)；

步骤(3)、自适应滤波处理；采用估计的直达波信号作为参考信号d(t)，对接收信号进行梯度自适应格型滤波处理，将滤波结果y(t)与参考信号d(t)相减，得到直达波抑制结果e(t)；

该基于IFFT频谱相除法和梯度自适应格型滤波的直达波抑制方法不需要采用单独的接收天线接收直达波参考信号，并能实现直达波的有效抑制，提高了系统的工作性能；

该基于IFFT频谱相除法和梯度自适应格型滤波的直达波抑制方法所采用的IFFT频谱相除方法进行时间延迟估计时具有较高的分辨率，不同于AR、MUSIC等功率谱估计方法，可以同时对直达波信号的时间延迟和幅度进行估计，便于生成参考信号；

该基于IFFT频谱相除法和梯度自适应格型滤波的直达波抑制方法所采用的梯度自适应格型滤波算法与归一化最小均方(LMS)滤波算法相比，收敛更快，具有更好的干扰杂波抑制效果。