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

说明书

技术领域

本发明属于航天和微波遥感的交叉技术领域，特别涉及一种星载SAR（Synthetic Aperture Radar，合成孔径雷达）数据中射频干扰抑制方法。

背景技术

星载SAR系统可以完成宽测绘带、高分辨率对地观测，是进行军事测绘和战场监视的有效手段。为了满足日益增长的军事和民用需求，特别是更高分辨率及隐蔽目标侦察的迫切需求，星载SAR的工作频段往更高和更低两个方向逐渐扩展，面临的电磁环境将会越来越复杂，特别是低波段星载SAR系统，工作频谱范围内存在大量的广播电视和通信系统，它们处于连续工作状态，而且发射功率较高，甚至有些系统的工作带宽较宽，比如DVB-T(Digital Video Broadcasting--Terrestrial,数字电视地面广播)系统，带宽有8MHz，这些干扰信号进入到星载SAR系统接收机，将严重影响星载SAR成像质量，导致图像解译和判读困难。

Feng Zhou提出的基于特征子空间分解的SAR射频干扰抑制方法对窄带射频干扰抑制效果比较好。但对于宽带射频干扰，由于带宽宽、功率高，在进行矩阵特征值分解时会产生数目较多的大特征值，现有的特征子空间分解法采用的主特征值数目判定方法会导致主特征值数目判定的误差。基于小波分解在信号奇异点检测方面的优势，本发明提出一种小波分解和特征子空间分解联合的星载SAR射频干扰抑制方法。本方法对窄带射频干扰和宽带射频干扰均有良好的抑制效果，且性能稳健。

发明内容

本发明的目的是：提出一种小波分解和特征子空间分解联合的星载SAR宽带射频干扰抑制方法，该方法充分利用了小波分解在信号奇异点检测方面的优势，对窄带射频干扰和宽带射频干扰均有良好的抑制效果，且性能稳健。

本发明技术方案的思路是：首先构造原始数据的协方差相关矩阵，对其进行特征分解，利用小波分解的方法进行主特征值判定，主特征值对应的即为干扰分量，利用主特征值对应的特征向量构造干扰子空间，然后将原始信号在特征干扰子空间上的映射分量减去，即可将干扰信号抑制掉。

本发明技术方案是：

假设x表示某一方位时刻星载SAR接收的数据，x＝[x₁,…,x_m,…,x_M]^T，其中M表示距离向采样点数，上标T表示转置操作。

第一步、构建数据矩阵

将数据x分为K个L维的子向量，L根据实际需要确定，通常32≤L≤128，其中K＝M-L+1。定义第k个子向量x_k为：

x_k＝[x_k,…x_k+L-1]^T,k＝1,2,…,K

然后，构建L×K的数据矩阵X：

X＝[x₁,…,x_k]

第二步、估计协方差相关矩阵

利用数据矩阵X估计其协方差相关矩阵

R^=1MΣk=1KxkxkT]]>

第三步、特征分解

对进行特征值分解，将特征值按照从大到小顺序排列，得到特征值序列λ₁,λ₂,…,λ_L，其相应的特征向量为u₁,u₂,…,u_L。

第四步、估计主特征值数目，构造干扰子空间

对特征值序列λ₁,λ₂,…,λ_L进行多尺度一维离散小波分解，得到的第i层的小波系数c_i为