[发明专利]一种基于压缩感知的可变脉冲重复间隔SAR成像方法

权利要求书

1.一种基于压缩感知的可变脉冲重复间隔SAR成像方法，其特征在于，对于某个场景，记回波信号矩阵为D，距离压缩后回波信号矩阵为S，理想场景后向散射系数为X；所述方法包括：

S1、对D进行距离压缩，表达式为：

D_rc＝D_rf⊙H_r

其中，D_rc为距离压缩后的信号矩阵，符号⊙表示Hadamard积，H_r为距离压缩因子矩阵；

S2、初始化2D-FISTA参数及矩阵，记l为FISTA当前迭代次数，L为总迭代次数；

S3、2D-FISTA通过迭代重建无方位模糊的图像，步骤S3具体包括以下分步骤：

S31、计算第l次迭代的重建图像矩阵，其中，Ψ为Daubechies-4小波稀疏基算子，Γ为X稀疏因子矩阵；

S32、计算第l次迭代的误差矩阵，其中，为对应方位向非均匀离散傅里叶变换矩阵，F_a为对应方位向离散傅里叶变换矩阵，为一个三维张量，H_a为一致压缩矩阵，S为完成距离压缩后的二维时域信号矩阵，B为S丢失位置的索引矩阵，符号(·)^H为矩阵的共轭转置运算，符号°定义了三维张量与矩阵的运算；

S33、计算第l次迭代的梯度，其中表示梯度算子，

S34、计算场景图像矩阵的稀疏因子矩阵，Γ^(l+1)＝soft(U^(l)，λ_l/l_f)，其中soft(U^(l),λ_l/l_f)为软阈值函数，Z^(l)＝Γ^(l)+(σ_l-1-1)/σ_l·(Γ^(l)-Γ^(l-1))，σ_l表示第l次迭代中更新的σ，σ_l-1表示第l-1次迭代中更新的σ；表示第l次迭代的梯度，Γ为X稀疏因子矩阵，l_f为利普希茨常数，λ_l表示第l次迭代中更新的λ，用于控制软阈值函数的门限；λ、σ为FISTA迭代过程中使用的常数参数；

S35、更新FISTA固定所需的参数，

S36、更新迭代次数，l＝l+1；

S37、判断迭代次数是否小于等于给定的总迭代次数，即lL是否成立，若成立则继续执行S31-S37，否则结束迭代，输出

2.根据权利要求1所述的一种基于压缩感知的可变脉冲重复间隔SAR成像方法，其特征在于，步骤S3中符号的运算过程为依次取三维张量的前向切片矩阵对其要进行处理的矩阵按照行进行傅里叶变换，计算式为：

其中，表示矩阵为第p行所有列元素，表示矩阵第p行所有列元素，表示三维张量的第p个前向切片矩阵。