[发明专利]一种外辐射源雷达系统多径抑制方法、系统及存储介质

权利要求书

1.一种外辐射源雷达系统多径抑制方法，其特征在于，所述方法包括如下步骤：

采用OMP算法进行杂波多径时延估计，确定离散时延个数及时延位置；

根据所确定的离散时延个数及时延位置重构杂波多径信号；

将重构的杂波多径信号与监测通道接收信号相比较，得到误差信号；

根据误差信号调节自适应对消滤波器权系数，使误差信号的均方差达到最小值；

杂波多径时延估计的方法包括：

初始化残差r₀＝q(k)，q(k)表示监测通道接收信号，迭代次数m＝1，索引集Λ₀为空集，选中的列集合Q(k)为空集，重复步骤a～e，最大迭代次数K：

a、寻找匹配向量r_m-1表示第m-1次迭代后的残差，β_j表示Ψ(k)的第j列，＜r_m-1,β_j＞表示求r_m-1和β_j的内积运算，|·|表示求绝对值运算，argmax{·}表示求|＜r_m-1,β_j＞|最大值时取的j值，N_Ψ表示Ψ(k)的列数；

b、扩展索引集，Λ_m＝Λ_m-1∪{λ_m}，λ_m表示第m次迭代搜索到的列序号，Λ_m-1表示第m-1次迭代搜索到的列序号集合，∪表示求并集运算；

c、扩展选中的列集合Q_m-1(k)表示按列序号集合Λ_m-1从Ψ(k)中搜出来的列向量集合，表示按列序号λ_m从Ψ(k)中搜出来的列向量；

d、利用Q_m(k)中的向量获得新的杂波多径传递函数：h(k)表示杂波多径传递函数；argmin{·}表示求||q(k)-Q_m(k)h(k)||最小值时取的h值；

e、更新残差r_m表示第m次迭代后的残差，q(k)表示监测通道接收信号，Q_m(k)表示按列序号集合Λ_m从Ψ(k)中搜出来的列向量集合，表示杂波多径传递函数的估计；

其中：k表示离散时间，Ψ(k)表示参考通道接收的直达波数据s(k)组成的hankel矩阵；

当时，记录k和的信息；

当时，将设置为0，此即估计的稀疏多径参数；

其中：Γ为预设置门限；

中非零抽头系数的个数等于稀疏度K，等于杂波多径条数，非零抽头的位置就是杂波多径离散时延位置；

采用下述公式重构杂波多径信号：

p_i(k)＝s(k-l_i)

式中：p_i(k)表示重构的第i个杂波多径信号；k表示离散时间；l_i表示第i个时延位置；s(k-l_i)表示第i个经过时延后的直达波信号，i＝1,2,…，K，K表示稀疏度；

调节自适应对消滤波器权系数的方法包括：

令为第i个杂波多径自适应对消滤波器权值，则自适应对消滤波器权向量表达式为

式中：[·]^T表示转置运算；k表示离散时间；

自适应对消滤波器输入递归向量P(k)为：

P(k)＝[p₁(k),p₂(k),…p_i(k),…p_K(k)]^T

式中：p_i(k)表示重构的第i个杂波多径信号；i≤K，K表示稀疏度；

自适应对消滤波器输出为：

式中：[·]^H表示共轭转置运算，q(k)表示监测通道接收信号，表示自适应对消滤波器权向量，P(k)表示自适应对消滤波器输入递归向量；

采用下式对自适应对消滤波器权向量进行更新：

式中：α、β为小正数，[·]^*表示共轭运算；

自适应对消滤波器输出e(k)与直达波信号s(k)进行距离多普勒互相关处理，公式如下所示：

其中，τ是时间延迟，f是多普勒频移，|Φ(τ,f)|表示距离多普勒互相关处理结果；N为数据长度；s(k-τ)是延时τ的直达波数据。