低秩矩阵分解的多维域联合SAR宽带干扰抑制方法技术

本发明专利技术提供了一种低秩矩阵分解的多维域联合SAR宽带干扰抑制方法，当前回波数据多个脉冲内存在宽带干扰信号，分别将脉冲回波信号的短时傅里叶变换矩阵向量化，进行RPCA分解，得到低秩矩阵和稀疏矩阵两个矩阵，将分解后的稀疏矩阵每一行重排成短时傅里叶矩阵形式，并对重排后的短时傅里叶矩阵进行短时傅里叶逆变换，将原始回波信号与重构出的干扰信号相减，即可实现宽带干扰抑制，利用现有的成像算法对干扰抑制后的数据进行成像，即可得到高分辨的图像。本发明专利技术避免了时频滤波引起的有用信号损失的问题，与传统的基于能量特性差异的方法相比，在抑制干扰的同时可有效地保留有用信号信息，能较大程度提升宽带干扰抑制后的图像质量。

Wideband interference suppression of multi-dimensional joint SAR Based on low rank matrix decomposition

【技术实现步骤摘要】
低秩矩阵分解的多维域联合SAR宽带干扰抑制方法
本专利技术涉及信号领域，尤其是一种时-频-脉冲多维域联合的干扰抑制方法，适用于在时频平面内具有时频聚集特性的宽带干扰抑制。
技术介绍
由于合成孔径雷达(SyntheticApertureRadar,SAR)具有全天时、全天候工作和广域监视等特点，在战场侦察和资源勘探等方面的运用受到越来越广泛的重视。近几年，干扰技术的快速发展对SAR精确遥感技术造成巨大的威胁。根据干扰信号的带宽可以把干扰分为窄带干扰(Narrow-BandInterference,NBI)和宽带干扰(Wide-BandInterference,WBI)。目前窄带干扰抑制方法较为成熟，例如频域限波法、子空间投影法、独立分量分析法、经验模态分解法等，都能够获得比较好的干扰抑制效果。窄带干扰由于其带宽小于有用信号带宽的1％，并且雷达回波建模相对简单，因此通过时域或频域的分析手段能够实现有效的干扰抑制。而对于宽带干扰，其干扰信号带宽一般大于有用信号带宽的10％，并且宽带干扰信号在时域和频域内均与有用信号高度重合，因此仅通过一维的时域或频域分析工具难以实现宽带干扰的有效抑制。现有的关于宽带干扰抑制方法基本是针对某种特定类型的干扰形式。现有的宽带干扰抑制方法大致可以分为两大类：参数化方法和非参数化方法。参数化方法是假设WBI为幅度恒定的二次调频信号、三次调频甚至高阶调频信号的叠加，通过Wigner分布、分数阶傅里叶变换、高阶模糊函数等时频分析工具估计干扰信号的频率、二次调频率、高阶调频率等参数，然后在时域重构并对消干扰信号。在实际的干扰信号模型与假设的模型匹配且模型参数精确估计条件下，理论上参数化干扰抑制方法是最优的。但是在实际情况中，干扰信号的数学建模比较复杂，而且模型参数估计误差会导致WBI重构不准确。另一种方法是非参数化方法，将原始回波的信号变换映射到时频二维域上，利用干扰和信号在时频域内的聚集特性差异，实现干扰与有用信号分离，最终实现干扰抑制的目的。非参数化方法不需要干扰模型的先验信息，也避免了高维的参数搜索过程，因此现有的大多数算法均采用非参数方法。短时傅里叶变换(ShortTimeFrequencyTransform，STFT)是最常用的一种时频分析工具，由于该时频分布式一种线性时频分布，且具有逆变换形式，现有的大多数的宽带干扰抑制方法均是基于STFT工具设计的。但这类方法在区分干扰和有用信号是基于信号能量在时频域内的能量特性差异，将能量强的信号认为是干扰信号，并采用限波或提取的方式，实现干扰抑制。但这类方法在抑制干扰的同时，会损失与干扰重合的有用信号，造成图像质量的下降。
技术实现思路
为了克服现有技术的不足，本专利技术提供一种低秩矩阵分解的多维域联合SAR宽带干扰抑制方法。为了克服现有基于短时傅里叶变换导致的有用信号能量损失的问题，本专利技术将每个脉冲回波数据变换到短时傅里叶变换域，对于有用信号，在相邻的若干个脉冲回波内，其回波信号是不变的或是缓变的，因此可以认为有用信号的时频分布在多个脉冲内满足低秩特性；而对于宽带干扰信号，其干扰信号发射脉冲重复周期一般情况下不会与有用信号发射周期严格保持一致，表现为在有用信号的脉冲重复周期内，起始频率、调频率、调制指数、初始相位等参数会发生比较大的变化，因此在相邻的脉冲重复周期内，干扰信号不具备低秩特性，如图1所示。本专利技术主要利用干扰信号和有用信号的在多个脉冲间的低秩特性差异，结合鲁棒主分量分析方法(RobustPrincipalComponentAnalysis,RPCA)，实现有用信号和宽带干扰信号的有效分离，在抑制干扰的同时尽可能地保留有用回波信号。本专利技术解决其技术问题所采用的技术方案采用以下步骤：步骤一：假设当前回波数据多个脉冲内存在宽带干扰信号，回波信号表示为X(n,q)，其中n代表快时间采样点，q表示脉冲数；将Q个脉冲的回波信号分别进行短时傅里叶变换，表示为：STFTX,q＝STFTI,q+STFTS,q+STFTN,q,q＝1,2,..,Q(1)其中，STFTX,q、STFTI,q、STFTS,q、STFTN,q分别表示为第q个脉冲的回波信号、宽带干扰信号、有用信号、噪声信号的短时傅里叶变换矩阵，矩阵STFTX,q、STFTI,q、STFTS,q、STFTN,q的维度均为N×M；步骤二：分别将Q个脉冲回波信号的短时傅里叶变换矩阵向量化，即将N×M维度的短时傅里叶变换矩阵拉成维度1×NM的短时傅里叶变换向量，如图2所示。然后，将Q个向量构成一个新的矩阵STFTX：其中，STFTX(q)表示第q个脉冲拉成维度1×NM的短时傅里叶变换向量，则新矩阵STFTX由以低秩的有用信号矩阵STFTS、稀疏的宽带干扰信号矩阵STFTI以及噪声矩阵构成STFTN：步骤三：对STFTX进行RPCA分解，由于信号满足低秩矩阵、稀疏矩阵和噪声矩阵之和的形式，因此经过RPCA分解之后，得到低秩矩阵和稀疏矩阵两个矩阵；步骤四：将分解后的稀疏矩阵每一行重排成短时傅里叶矩阵形式，即将1×NM维度的短时傅里叶向量重排成维度为N×M的短时傅里叶变换矩阵，如图2所示，并对重排后的短时傅里叶矩阵进行短时傅里叶逆变换，即可得到Q个脉冲内重构后的干扰信号步骤五：将原始回波信号与重构出的干扰信号相减，即可实现宽带干扰抑制最后，利用现有的成像算法对干扰抑制后的数据进行成像，即可得到高分辨的图像。本专利技术步骤三中采用GoDecomposition(GoDec)算法进行低秩矩阵分解：其中，||·||F表示Frobenius范数，rank(·)表示求秩操作，card(·)表示稀疏矩阵的基数，r表示秩，ε表示稀疏度。本专利技术的有益效果在于所利用的是宽带干扰信号和有用回波信号在多个脉冲间的低秩特性差异，在多个脉冲内同时实现干扰信号与有用信号的分离。由于本专利技术避免了时频滤波引起的有用信号损失的问题，与传统的基于能量特性差异的方法相比，在抑制干扰的同时可有效地保留有用信号信息，能较大程度提升宽带干扰抑制后的图像质量。附图说明图1是本专利技术宽带干扰信号和有用信号的短时傅里叶变换示意图。图2是本专利技术矩阵化与向量化相互转换关系示意图。图3是本专利技术基于低秩矩阵分解的时-频-脉冲多维域联合SAR宽带干扰抑制方法的流程图。图4是频率限波、时频限波以及时-频-脉冲多维域联合宽带干扰抑制结果。具体实施方式下面结合附图和实施例对本专利技术进一步说明。步骤一：假设当前回波数据多个脉冲内存在宽带干扰信号，回波信号表示为X(n,q)，其中n代表快时间采样点，q表示脉冲数；将Q个脉冲的回波信号分别进行短时傅里叶变换，表示为：STFTX,q＝STFTI,q+STFTS,q+STFTN,q,q＝1,2,..,Q(1)其中，STFTX,q、STFTI,q、STFTS,q、STFTN,q分别表示为第q个脉冲的回波信号、宽带干扰信号、有用信号、噪声信号的短时傅本文档来自技高网...

【技术保护点】
1.一种低秩矩阵分解的多维域联合SAR宽带干扰抑制方法，其特征在于包括下述步骤：/n步骤一：假设当前回波数据多个脉冲内存在宽带干扰信号，回波信号表示为X(n,q)，其中n代表快时间采样点，q表示脉冲数；将Q个脉冲的回波信号分别进行短时傅里叶变换，表示为：/nSTFT

【技术特征摘要】
1.一种低秩矩阵分解的多维域联合SAR宽带干扰抑制方法，其特征在于包括下述步骤：
步骤一：假设当前回波数据多个脉冲内存在宽带干扰信号，回波信号表示为X(n,q)，其中n代表快时间采样点，q表示脉冲数；将Q个脉冲的回波信号分别进行短时傅里叶变换，表示为：
STFTX,q＝STFTI,q+STFTS,q+STFTN,q,q＝1,2,..,Q(1)
其中，STFTX,q、STFTI,q、STFTS,q、STFTN,q分别表示为第q个脉冲的回波信号、宽带干扰信号、有用信号、噪声信号的短时傅里叶变换矩阵，矩阵STFTX,q、STFTI,q、STFTS,q、STFTN,q的维度均为N×M；
步骤二：分别将Q个脉冲回波信号的短时傅里叶变换矩阵向量化，即将N×M维度的短时傅里叶变换矩阵拉成维度1×NM的短时傅里叶变换向量，然后，将Q个向量构成一个新的矩阵STFTX：



其中，STFTX(q)表示第q个脉冲拉成维度1×NM的短时傅里叶变换向量，则新矩阵STFTX由以低秩的有用信号矩阵STFTS、稀疏的...

【专利技术属性】
技术研发人员：粟嘉，李浩江，陶明亮，崔国楠，范一飞，谢坚，王伶，张兆林，韩闯，宫延云，
申请(专利权)人：西北工业大学，
类型：发明
国别省市：陕西;61