1.一种基于稀疏表示理论的合成孔径雷达成像方法，包括如下步骤：
(1)根据雷达向地面场景发射脉冲 获得地面场景的回波信号
式中， 为快时间，ta为慢时间，fc为载频，γ为调频斜率，v为载机速度，C为电磁波传播速度，wa(·)为方位窗函数，wr(·)为线性调频脉冲信号时间窗函数，P为地面场景在方位向的离散网格个数，Q为地面场景在距离向的离散网格个数，p为地面场景在方位向的第p个离散网格，q为地面场景在距离向的第q个离散网格，fpq为第(p,q)个网格处散射点的散射系数，R0为雷达到地面场景中心的垂直斜距，Rq为雷达到距离向第q个离散网格的垂直斜距，xp为方位向第p个离散网格的x轴坐标；
(2)对回波信号 进行二维离散采样，得到如下回波矩阵Y：
其中，M为方位向发射脉冲的数目，N为每个脉冲内，距离向采样点数； 为回波信号在第(m,n)个采样时刻的样本值，其表示式如下：
式中，λ为载波波长，ta,m为第m个脉冲， 为快时间第n个采样时刻值；
(3)构造地面场景的方位向基矩阵R：
其中
(4)构造地面场景的距离向基矩阵G：
其中，
(5)根据回波矩阵Y、方位向基矩阵R、距离向基矩阵G和地面场景的散射系数矩阵F，构建如下矩阵关系式：
Y＝RFG+N，
式中，N是与回波矩阵Y同维度的噪声矩阵，散射系数矩阵F的形式如下：
其中，p＝1,2,…,P，q＝1,2,…,Q；
(6)根据地面场景的散射系数矩阵F，获得中间变量矩阵V，水平全变差矩阵Dx，垂直全变差矩阵Dy：
V＝F
Dx＝DhF，
Dy＝FDv
式中，Dh为水平全变差算子矩阵，Dv为垂直全变差算子矩阵；
(7)利用步骤(5)的矩阵关系式和步骤(6)得到的参数，构造全变差框架下的优化函数f(F,V,Dx,Dy)：
其中 表示求函数最小值的运算符，λ为全变差罚因子，ν为场景散射系数罚因子，μ为回波罚因子，||·||F表示求矩阵的Frobenius范数，||·||1表示对矩阵的每一个元素取绝对值后再相加；
(8)用Split-Bregman方法对步骤(7)中的优化函数f(F,V,Dx,Dy)进行迭代求解，得*
到地面场景散射系数矩阵F的最终迭代结果F；
*
(9)对步骤(8)得到地面场景散射系数矩阵F的最终迭代结果F取模值，输出得到地面场景的图像。
2.根据权利要求1所述的基于稀疏表示理论的合成孔径雷达成像方法，其中所述步骤(1)雷达向地面场景发射脉冲 其表示如下：
式中，fc为载频，γ为调频斜率，t为全时间， 为快时间，ta为慢时间，这三个时间的关系为 表示线性调频脉冲信号的时间窗函数，Tr是脉冲持续的
时间。
3.根据权利要求1所述的基于稀疏表示理论的合成孔径雷达成像方法，其中所述步骤(8)用Split-Bregman方法对步骤(7)中的优化函数f(F,V,Dx,Dy)进行迭代求解，按如下步骤进行：
(8a)根据优化函数f(F,V,Dx,Dy)，得到Split-Bregman形式的优化函数L(F,V,Dx,Dy)：
式中，Bx为水平全变差矩阵辅助变量，By为垂直全变差矩阵辅助变量，W为散射系数矩阵辅助变量；
(8b)初始化迭代步数k＝0，第k次迭代后散射系数矩阵Fk为全1矩阵，第k次迭代k
后中间变量矩阵V为全1矩阵，第k次迭代后水平全变差矩阵 为全1矩阵，第k次迭代k
后垂直全变差矩阵 为全1矩阵；第k次迭代后散射系数矩阵辅助变量W为全0矩阵，第k次迭代后水平全变差矩阵辅助变量 为全0矩阵，第k次迭代后垂直全变差矩阵辅助变量 为全0矩阵；设置全变差罚因子λ>0，场景散射系数罚因子ν>0，回波罚因子μ>0，迭-4
代终止条件ε＝5×10 ；
k+1
(8c)获得第k+1次迭代后散射系数矩阵F ；
k+1
(8d)获得第k+1次迭代后中间变量矩阵V ；
k+1
(8e)根据第k+1次迭代后中间变量矩阵V ，第k次迭代后水平全变差矩阵辅助变量第k次迭代后垂直全变差矩阵辅助变量 按照如下公式获得第k+1次迭代后水平全变差矩阵 第k+1次迭代后垂直全变差矩阵
式中，p＝1,2,…,P，q＝1,2,…,Q，Dh为水平全变差算子矩阵，Dv为垂直全变差算子矩阵，sign(·)为符号函数；
(8f)根据第k+1次迭代后中间变量矩阵Vk+1，第k+1次迭代后散射系数矩阵Fk+1，第k+1次迭代后水平全变差矩阵 第k+1次迭代后垂直全变差矩阵 第k次迭代后水平全变差矩阵辅助变量 第k次迭代后垂直全变差矩阵辅助变量 第k次迭代后散射系数k k+1
矩阵辅助变量W，按照如下公式获得第k+1次迭代后散射系数矩阵辅助变量W ，第k+1次迭代后水平全变差矩阵辅助变量 第k+1次迭代后垂直全变差矩阵辅助变量k+1 k k+1 k+1
W ＝W +V -F
k+1 k
(8g)根据第k+1次迭代后散射系数矩阵F ，第k次迭代后散射系数矩阵F，按如下公式得到均方误差M：
式中，||·||F表示求矩阵的Frobenius范数；
(8h)判断均方误差M≤ε是否成立，若成立执行步骤(8i)；否则令k＝k+1，返回步骤(8c)继续迭代运行，其中ε为迭代终止条件；
(8i)令F*＝Fk+1，得到地面场景散射系数矩阵F的最终迭代结果F*。
4.根据权利要求3所述的基于稀疏表示理论的合成孔径雷达成像方法，其中所述步骤(8c)中获得第k+1次迭代后散射系数矩阵Fk+1，按照如下步骤进行：
(8c1)根据方位向基矩阵R，距离向基矩阵G，获得方位向对称矩阵Ψ，距离向对称矩阵Φ：
Ψ＝RTR
，
Φ＝GGT
其中，(·)T表示矩阵转置；
(8c2)根据特征分解方位向对称矩阵Ψ，距离向对称矩阵Φ，得到方位向正交矩阵UR，方位向对角矩阵LR，距离向正交矩阵UG，距离向对角矩阵LG；
(8c3)根据回波矩阵Y，方位向基矩阵R，距离向基矩阵G，第k次迭代后中间变量矩阵Vk，第k次迭代后散射系数矩阵辅助变量Wk，方位向正交矩阵UR，距离向正交矩阵UG，方位向对角矩阵LR，距离向对角矩阵LG，按照如下公式获得第k+1次迭代后散射系数辅助矩阵式中，p＝1,2,…,P，q＝1,2,…,Q；
(8c4)根据方位向正交矩阵UR，距离向正交矩阵UG，第k+1次迭代后散射系数辅助矩阵k+1
按照如下公式获得第k+1次迭代后散射系数矩阵F ：
5.根据权利要求3所述的基于稀疏表示理论的合成孔径雷达成像方法，其中所述步骤k+1
(8d)中获得第k+1次迭代后中间变量矩阵V ，按照如下步骤进行：
(8d1)根据水平全变差算子矩阵Dh，垂直全变差算子矩阵Dv，获得水平全变差矩阵E，垂直全变差矩阵Z：
T
其中，(·)表示矩阵转置；
(8d2)根据特征分解水平全变差矩阵E，垂直全变差矩阵Z，得到水平全变差正交矩阵Ux，水平全变差对角矩阵Lx，垂直全变差正交矩阵Uy，垂直全变差对角矩阵Ly；
k+1
(8d3)根据第k+1次迭代后散射系数矩阵F ,水平全变差算子矩阵Dh，垂直全变差算k
子矩阵Dv,第k次迭代后散射系数矩阵辅助变量W，第k次迭代后水平全变差矩阵辅助变量 第k次迭代后垂直全变差矩阵辅助变量 第k次迭代后水平全变差矩阵 第k次迭代后垂直全变差矩阵 水平全变差正交矩阵Ux，垂直全变差正交矩阵Uy，水平全变差对角矩阵Lx，垂直全变差对角矩阵Ly，按照如下公式获得第k+1次迭代后中间变量辅助矩阵式中，p＝1,2,…,P，q＝1,2,…,Q；
(8d4)根据水平全变差正交矩阵Ux，垂直全变差正交矩阵Uy，第k+1次迭代后中间变量k+1
辅助矩阵 按照如下公式获得第k+1次迭代后中间变量矩阵V ：