1.一种基于低秩分解与矩阵填充的SAR抗射频干扰方法，具体步骤如下：S1、对接收到的回波进行预处理，通过短时傅里叶变换至时频域，获得时频谱；对接收到的回波进行脉冲压缩处理，得到脉压后的回波矩阵S；然后对矩阵S做短时傅里叶变换得到时频谱S tf ，表达式如下：S tf ＝STFT{S}其中，STFT表示短时傅里叶变换；S tf 表示时频谱；S2、基于步骤S1，在时频域使用一维陷波滤波器去除受干扰的时频谱，将陷波后的时频谱变换回距离频域获得成分缺失的频谱；首先将时频谱S tf 输入自适应一维陷波滤波器进行滤波，表达式如下：其中， 表示陷波后的时频谱，ANF 1 表示自适应一维陷波滤波器；然后将陷波后的时频谱进行逆短时傅里叶变换，表达式如下：其中， 表示陷波后的回波，ISTFT表示逆短时傅里叶变换；再对陷波后的回波做距离向傅里叶变换，表达式如下：其中， 表示陷波后的频谱，FFT ran 表示距离向傅里叶变换；S3、基于步骤S2，将频谱每个方位向脉冲重塑为Toeplitz矩阵；将频谱矩阵 沿方位向划分为N个脉冲向量 脉冲向量 表达式如下：其中，M表示脉冲向量的长度；然后定义Toeplitz变换算子Γ: 表达式如下：其中， 表示复数域，n 1 和n 2 分别表示变换后矩阵的行数和列数；h是一个超参数，通过调整其大小改变变换后Toeplitz矩阵的形状；最后对N个脉冲向量分别应用Toeplitz变换算子Γ，得到低秩矩阵 表达式如下：S4、基于步骤S3，利用基于低秩分解填充矩阵的方法恢复全部Toeplitz矩阵；S41、初始化参数，利用奇异值分解初始化两个矩阵因子，对步骤S3得到的所有低秩矩阵均进行矩阵补全恢复；先初始化参数，包括：正则化参数α，缩放因子β，惩罚因子ρ 0 ，维度参数d，迭代次数k＝0和最大迭代次数K；其中，矩阵的上标表示迭代次数，惩罚因子ρ 0 表示其初值；再初始化矩阵，即利用奇异值分解初始化两个矩阵因子，对矩阵 进行奇异值分解 且A 0 ＝U(:,1:d)S(1:d,1:d) 1/2 ,Z 0 ＝zero(d,n 2 ),Y 0 ＝0, 其中，上标H表示共轭转置；矩阵A和Z表示对 的分解，Y表示拉格朗日乘子矩阵，zero表示生成零矩阵， 表示每次迭代恢复的变量 的初值；S42、基于步骤S41，对待优化矩阵和参数迭代更新；(1)更新Z k+1 ；通过对增广拉格朗日函数求偏导等于零，可得表达式如下：其中，I表示单位矩阵；(2)更新A k+1 ；通过对增广拉格朗日函数求偏导等于零，可得表达式如下：(3)更新 通过对增广拉格朗日函数求偏导等于零，可得表达式如下：(4)更新拉格朗日乘子Y k+1 ；(5)更新惩罚因子ρ k+1 ；ρ k+1 ＝βρ kS43、判断是否满足迭代终止条件，当k＝K时，输出 作为 的恢复结果，进入步骤S44，否则返回步骤S42，令k＝k+1，继续迭代；S44、遍历方位向，对每个构造的Toeplitz矩阵 均进行步骤S41～S43的优化恢复，得到N个恢复出的频谱矩阵 其中，上标n表示每个脉冲的恢复结果；S5、基于步骤S4，将恢复出的全部Toeplitz矩阵重塑为脉冲向量，拼接得到完整频谱；对每个恢复出的频谱矩阵 按第一行和第一列重构为一个脉冲向量，然后按脉冲顺序拼接重塑为恢复的频谱矩阵S F ，并利用距离向逆快速傅里叶变换将频谱矩阵变换为回波矩阵ζ est ，表达式如下：ζ est ＝IFFT ran {S F }其中，IFFT ran 表示距离向逆快速傅里叶变换；S6、基于步骤S5，对回波矩阵ζ est 采用后向投影法对目标成像；首先进行场景网格划分，将成像区域分割为多个网格，最终生成的SAR图像中每个像素点对应一个网格单元；然后在每个方位时刻，计算雷达与每个网格点之间的距离，并确定双程时延，从当前方位时刻获取的数据中，依据双程时延找到对应的距离数据，该数据即为该网格点在此方位时刻接收到的信号；最后进行相干叠加，对每个方位时刻的数据进行相位校正后，再进行叠加，最终生成SAR图像。
