[发明专利]一种基于子空间的高分辨宽测绘带SAR-GMTI杂波抑制方法

权利要求书

1.一种基于子空间的高分辨宽测绘带SAR-GMTI杂波抑制方法，其特征在于步骤如下：

步骤1：获取的多通道雷达系统的第k个通道回波信号是二维矩阵S_k(n,m)，所述的S_k(n,m)为nrn×nan维矩阵，k＝1,2，…,K，多通道系统共有K个通道；对矩阵S_k(n,m)列进行FFT处理，将结果保存在矩阵S_k(n,m)中；其中，nrn表示距离向点数，nan表示方位向点数；

步骤2：根据已知的雷达参数，构造参考信号向量S_ref(n)为nrn×1向量；其中，γ表示调频率，γ＝B/Tp，B表示发射信号带宽，Tp表示发射脉冲宽度，f_n表示为距离向频域坐标，B为发射信号带宽，Δf为距离频域间隔，

步骤3：取出步骤1得到S_k(n,m)的每一列，均点乘参考信号向量S_ref(n)的共轭，得到距离脉压后的第k通道数据矩阵S_k(f_n,x_m)；x_m表示方位向时域坐标，L表示为合成孔径长度，m＝0,1,...,nan-1；

步骤4：根据已知的雷达参数，构造第k个通道参考信号向量S_ref_k(m)为1×nan向量；其中，f_c表示雷达发射信号的载频，c为电磁波的传播速度，V为卫星速度，θ为卫星的斜视角,d_k为第k通道与参考通道沿卫星速度方向的距离，t_m为方位慢时间；

步骤5：取出步骤3得到的S_k(f_n,x_m)的每一行，均点乘参考信号向量S_ref_k(m)的共轭，得到距离走动校正之后的第k通道数据矩阵S′_k(f_n,x_m)；

步骤6：根据已知的雷达参数，构造第k个通道chirp傅里叶变换参考信号其中，R为卫星平台到场景中心斜距；

步骤7：取出步骤5得到的S′_k(f_n,x_m)点乘参考信号向量S_Chirp_k(n,m)的共轭，并对相乘结果行进行FFT处理，得到chirp傅里叶变换之后的第k通道数据矩阵S″_k(f_n,f_m)；其中，f_m表示为方位向频域坐标，PRF为方位采用频率，Δf_a为方位频域间隔，

步骤8：取出步骤7得到S″_k(f_n,f_m)，沿f_n方向进行IFFT处理，得到S_k″(n,f_m)，同时构造多通道系统回波矢量S(n,f_m)＝[S″₁(n,f_m) S″₂(n,f_m) … S″_k(n,f_m) … S″_K(n,f_m)]^T；其中[·]^T表示矩阵转置；

步骤9：取出步骤8得到的矢量对协方差矩阵进行估计，协方差矩阵

步骤10：取出步骤9得到协方差矩阵进行特征值分解，得到最新特征为λ_min，其对应的特征向量为杂波对应的零空间V_min；

步骤11：取出步骤10得到杂波对应的零空间V_min，并进行杂波抑制，杂波抑制过程及输出可以表示为