[发明专利]一种SAR射频干扰抑制方法

摘要

本发明提供一种SAR射频干扰抑制方法。技术方案是：以SAR数据的一条距离线回波信号为输入，根据回波信号构建相空间矩阵，利用子空间投影方法得到参考信号的相空间矩阵；采用对角平均化方法，从参考信号相空间矩阵中恢复出一维的参考信号；根据参考信号，确定干扰信号的判别门限；根据确定的干扰判别门限，确定回波中存在的干扰频点，并用参考信号的幅度均值作为替换，完成该条距离线的射频干扰抑制。本发明可以有效提升干扰抑制后的SAR图像质量。

主权项

1.一种SAR射频干扰抑制方法，SAR是指合成孔径雷达，假设x(n),1≤n≤Nr，表示某一慢时刻SAR接收的回波信号，其中Nr表示距离向采样点数，其特征在于，包括下述步骤：第一步：提取参考信号：将雷达接收的回波信号x(n)按照时延嵌套构造成相空间矩阵S：上式中M＝Nr+1-L，L的大小取决于射频干扰的个数；令协方差矩阵R＝SSH，对协方差矩阵R进行特征值分解：R＝UΛUH其中Λ＝diag[λ1,λ2,…,λL]是一个L×L的对角矩阵，对应的特征矢量矩阵为U＝[u1,u2,…,uL]，并且λ1≥λ2≥…≥λL，λi为特征值，i＝1,2，…L；上标H表示求矩阵的共轭转置；将前k个特征值所对应的特征矢量提取出，组成参考子空间Uk＝[u1,u2,…,uk]，k的取值根据实际情况确定；令参考信号相空间矩阵采用对角平均化方法从参考信号相空间矩阵Sk中恢复出1维参考信号序列xk(n)；第二步，对回波信号进行滤波将回波信号x(n)和参考信号xk(n)分别进行快速傅里叶变换得到X和Xk：X＝[X(1),X(2),…,X(Nr)]Xk＝[Xk(1),Xk(2),…,Xk(Nr)]利用下式计算判别门限γ： γ = | X ‾ k + 3 δ | ]]> X ‾ k = 1 N r Σ i = 1 N r X k ( i ) ]]> δ 2 = 1 N r Σ i = 1 N r ( X k ( i ) - X ‾ k ) 2 ]]>利用下式对回波信号进行滤波，得到射频干扰抑制后的回波信号 X ^ ( n ) = | X ‾ k