一种合成孔径雷达高精度运动补偿方法
【专利摘要】本发明提出了一种合成孔径雷达高精度运动补偿方法，它是首先用雷达平台匀速直线运动轨迹进行BP成像，获得粗聚焦的SAR图像，然后以SAR图像强度作为目标函数，然后以图像强度最优为准则迭代调整天线相位中心误差，迭代结束后获得天线相位中心误差估计，天线相位中心(Antenna?Phase?Center，APC)位置误差，最后将估计出的APC位置误差加到匀速直线运动轨迹上，获得APC绝对位置，并用BP算法进行最终的高精度成像。本发明方法与现有的基于图像强度的自聚焦BP算法相比，能获得更高的运动补偿精度，且大大降低了内存消耗。
【专利说明】一种合成孔径雷达高精度运动补偿方法
【技术领域】
[0001]本发明属于合成孔径雷达(Synthetic Aperture Radar, SAR)高分辨率成像的【技术领域】，它特别涉及到了 SAR高精度运动补偿的【技术领域】。
【背景技术】
[0002]合成孔径雷达(SAR)是一种高分辨率微波成像雷达，具有全天时和全天候工作的优点，已被广泛应用各个领域，如地形测绘、制导、环境遥感和资源勘探等。SAR应用的重要前提和信号处理的主要目标是通过成像算法获取高分辨、高精度的微波图像。但是诸多环境因素(如风场、湍流等)将导致雷达载体平台的运动轨迹偏离设计的理想轨迹，从而严重影响SAR图像的质量(包括聚焦深度、对比度等)。因此，运动补偿技术成为了 SAR成像过程中的关键技术。
[0003]后向投影(BP)算法是一种精确的SAR时域成像算法，它首先将合成孔径雷达原始数据沿距离向进行距离压缩(脉冲压缩)，然后通过选择不同慢时间观测空间中任意像素点在距离压缩后SAR数据空间中的数据，补偿方位向多普勒相位，并进行相干积累，最终获得各像素点散射系数的成像算法。由于在精确已知天线相位中心(Antenna Phase Center,APC)的前提下，BP算法可以有效补偿运动误差，因而已被广泛应用。详见“师君.双基地SAR与线阵SAR原理及成像技术研究[D].电子科技大学博士论文.2009”。但是，当平台轨迹精度较低或者未知时，BP算法的成像精度会大大降低。因此，研究针对BP算法的运动补偿技术具有重要意义。
[0004]自聚焦BP算法是一类基于空域图像质量的自聚焦算法，也可以看作是一类针对BP算法的运动补偿方法，其主要过程是根据图像质量指标优化方位向相位补偿误差向量，当图像质量指标达到最优时，SAR图像聚焦最好。目前主要的自聚焦BP算法有基于最小图像熵的自聚焦 BP 算法(详见“M.Liu, C.S.Li, X.H.Shi, A back-project1n fastautofocus algorithm based on minimum entropy for SAR imaging[C].3rd APSARConference.2011:1_4”)、结合自聚焦和快速BP的高精度成像算法(详见“L.Zhang，H.L.Li, Z.J.Qiao, M.D.Xing, Z.Bao, Integrating autofocus techniques with fastfactorized back-project1n for high-resolut1n spotlight SAR imaging[J].1EEEGeoscience and Remote Sensing Letters.2013,10 (6): 1394-1398”)和基于图像锐度的自聚焦 BP 算法(详见“J.N.Ash, An autofocus method for backproject1n imagery insynthetic aperture radar[J].1EEE Geoscience and Remote Sensing Letters.2012,9(1):104-108”)。其中基于图像锐度的自聚焦BP算法成像效果最好。但是，基于图像锐度的自聚焦BP算法对整个场景进行统一的相位补偿，降低了大部分像素点的相位补偿精度，从而限制了成像精度。
【发明内容】
[0005]为了对场景中的每个像素点进行高精度的相位补偿，本发明提出了一种合成孔径雷达高精度运动补偿方法，其特点是:首先用雷达平台匀速直线运动轨迹进行BP成像，获得粗聚焦的SAR图像，然后以SAR图像强度作为目标函数，利用最优化技术估计出天线相位中心(Antenna Phase Center,APC)位置误差,最后将估计出的APC位置误差加到勻速直线运动轨迹上，获得APC绝对位置，并用BP算法进行最终的高精度成像。本发明方法与现有的基于图像强度的自聚焦BP算法相比，能获得更高的运动补偿精度，且大大降低了内存消耗。
[0006]为了方便描述本发明的内容，首先作以下术语定义:
[0007]定义1、脉冲压缩
[0008]脉冲压缩是一种现代雷达信号处理技术，简单来说就是雷达发射宽脉冲，然后再接收端“压缩”为窄脉冲，从而改善雷达的两种性能:作用距离和距离分辨率。详见“皮亦鸣，杨建宇，付毓生，杨晓波.合成孔径雷达成像原理.第一版.电子科技大学出版社.2007.3”。
[0009]定义2、升采样
[0010]升采样是一种在离散信号域提高信号采样率的方法，有时域升采样和频域升采样两种实现方式。
[0011]定义3、快速傅里叶变换
[0012]计算离散傅里叶变换的一种快速算法，简称FFT。采用这种算法能使计算机计算离散傅里叶变换所需要的乘法次数大为减少，特别是被变换的抽样点数越多，FFT算法计算量的节省就越显著。FFT的逆变换叫做逆傅里叶变换，简称IFFT。详见“程乾生.数字信号处理.北京大学出版社，北京，2003”。
[0013]定义4、后向投影算法
[0014]后向投影算法，简称BP算法。BP算法首先利用雷达平台的轨迹信息求出雷达平台与场景像素点的距离历史，然后通过距离历史找出回波数据中对应的复数据，再进行相位补偿并相干累加，从而得到该像素点的复图像值。详见“师君.双基地SAR与线阵SAR原理及成像技术研究[D].电子科技大学博士论文.2009”。
[0015]定义5、方位向、距离向
[0016]将雷达平台运动的方向叫做方位向，将垂直于方位向的方向叫做距离向。
[0017]定义6、快时间、慢时间、慢时刻
[0018]快时间是距离向采样的时间，慢时间是方位向采样的时间，将离散的慢时间从I开始编号，每一个编号叫做一个慢时刻。
[0019]定义7、天线相位中心
[0020]天线相位中心，简称APC，是雷达天线发射信号的位置，精确的天线相位中心是BP算法能够精确成像的前提。
[0021]定义8、图像强度
[0022]图像强度是指一幅复图像中每个像素点幅度平方之和，可以用来表征SAR图像聚焦好坏，SAR图像聚焦越好，图像强度越大。
[0023]定义9、共轭梯度法
[0024]共轭梯度法是一种最优化方法，用迭代点处的负梯度向量为基础产生一组共轭方向。共轭梯度法的收敛速度快，且不用求矩阵的逆，在使用计算机求解时，所需存储量较小。详见“何坚勇.最优化方法.第一版.清华大学出版社.北京.2007.1”。
[0025]定义10、Armijo 算法
[0026]Armijo算法是一种一维搜索算法，可以保证目标函数在搜索方向上充分下降。详见“ARMIJ0.Minimizat1n of funct1ns having Lipschitz continuous first partialderivatives [J].Pacific Journal of Mathematics, vol.16，n0.1，pp.1-3.1966，，。
[0027]本发明提供了一种合成孔径雷达高精度运动补偿方法，它包括如下步骤(如附图1所示):
[0028]步骤1、相关参数的初始化
[0029]初始化的参数均为已知，且初始化的参数如下:光速为C ;雷达发射线性调频信号，载频为ω ;雷达发射脉冲的带宽为B;雷达发射脉冲的时宽为Tp;雷达脉冲重复周期为T ;雷达回波距离向采样频率为fs ;雷达回波数据矩阵为SK>a ;雷达回波数据矩阵SK>a的方位向点数和距离向点数分别为K和L(K和L均为正整数)，K也称为慢时刻数；慢时刻向量为 ts = [-K/2，1-K/2，, K/2-1] XT ;升采样数据矩阵为 SSKXP，SSkxp 的大小为 K 行、8 X L列，K为雷达回波数据矩阵的方位向点数，L为雷达回波数据的距离向点数；几何坐标系为三维笛卡尔坐标系O-XYZ ;雷达平台速度向量为V，V的大小为I行3列；雷达平台在零时刻的位置向量为Pttl, Pto的大小为I行3列；0-ΧΥ平面内的矩形场景为Θ ;将Θ离散化为像素点网格，记为ΩΜΧΝ;像素点网格ΩΜΧΝ中X方向和Y方向的像素点点数分别为M和N;像素点网格Ωμχν中X方向和Y方向的像素点间隔分别为dx和dy ;场景中心位置向量为P。，P。的大小为I行3列；共轭梯度算法最大迭代次数为Q，Q为正整数；初始APC误差向量为ΔΡ2，ΔΡ3 = Olx诚，Olx3i^ I行、3XK列的零矩阵，K为雷达回波数据矩阵SKXl的方位向点数。
[0030]步骤2、对雷达回波数据矩阵的每一行进行脉冲压缩
[0031]取出步骤I中所有雷达回波数据SKXy采用传统的脉冲压缩方法对SK>a的每一行进行脉冲压缩，得到脉冲压缩后的数据矩阵PSK>a。
[0032]步骤3、对脉冲压缩后的数据矩阵的每一行进行频域升采样
[0033]对步骤2得到的脉冲压缩后的数据矩阵PSK>a的每一行统一做如下8倍频域升采样处理:
[0034]步骤3.1、取出步骤2中脉冲压缩后的数据矩阵PSkx1j的第k行向量,记为sk, k =1，2，...，K，K为步骤I定义的慢时刻数。
[0035]步骤3.2、对向量Sk做传统的快速傅里叶变换(FFT)，得到向量fk。
[0036]步骤3.3、从向量fk的L/2+1位置开始插入7XL个零，得到向量zk，Zk= [fk⑴，fk(2), , fk (L/2), 01X7L, fk (L/2+1), , fk(L)], fk(l)为向量 fk 中的第 I 个元素，fk(2)为向量fk中的第2个元素，fk(L/2)为向量fk中的第L/2个元素,01Χι为I行、7XL列的零向量，fk(L/2+l)为向量fk中的第L/2+1个元素，fk(L)为向量fk中的第L个元素，L为步骤I定义的雷达回波数据矩阵距离向点数。
[0037]步骤3.4、对向量Zk做传统的快速傅里叶逆变换(IFFT)，得到向量ssk，将向量ssk存放到矩阵SSkxp的第k行，SSkxp为步骤I定义的升采样数据矩阵。
[0038]步骤4、计算初始搜索方向
[0039]步骤4.1、采用公式Pa% = Pto + V X ts(/c)计算第k个匀速直线APC，记为PSjc, Pttl
为步骤I定义的雷达平台在零时刻的位置向量，V为步骤I定义的雷达平台速度向量，ts(k)为向量ts的第k个元素，k = 1,2,..., K, K为步骤I定义的慢时刻数，ts为步骤I定义的慢时刻向量。
[0040]步骤4.2、采用公式Pli = Pc + Km - M/2) X dx, (η - N/2) X dy, O ]计算像素点网格ΩΜΧΝ中第mSn列像素点的位置向量，记为P。 P。为步骤I定义的场景中心位置向量，m = 1,2,..., M, M为步骤I定义的像素点网格Ωμχν中X方向的像素点点数，η = I,
2,..., N, N为步骤I定义的像素点网格Ωμχν中Y方向的像素点点数，dx为步骤I定义的像素点网格ΩΜΧΝ中X方向的像素点间隔，dy为步骤I定义的像素点网格ΩΜΧΝ中Y方向的像素点间隔。
[0041]步骤4.3、米用公式
[0042]APlk = [ΔΡπ°(3(/? -1) + I), AP°(3(/c -1)4- 2), AP°(3(/c -1) + 3)j
计算第第k个APC误差向量，记为ΔΡ&，APS力步骤I定义的初始APC误差向量，AP°(3(fe -1) + I)为 ΔΡ3中第 3(k-l)+l 个元素,APa°(3(/c -1) + 2)为 ΔΡ^ 中第 3 (k_l)+2个元素，AP°(3(k -1) + 3)为ΔΡ^中第3 (k-1) +3个元素，k = 1，2，...，K，K为步骤I定义的
pO _j_^p0 _p0
慢时刻数；采用公式= |pf+^tpri计算中间单位向量，记为.eo K0为步骤4.1
I CLlJi ClJi TiXTl ITYtTtfk,9 (X
计算出的第k个匀速直线APC P 勾步骤4.2计算出的像素点网格ΩΜΧΝ中第mSn列像素点的位置向量，m = 1,2,..., M, M为步骤I定义的像素点网格Ωμχν中X方向的像素点点数，η = 1,2,…，N, N为步骤I定义的像素点网格Ωμχν中Y方向的像素点点数。
[0043]步骤4.4、采用公式私n,fc = 2|P°fc - Ρ,°.|/?:计算第k个APC与像素点网格ΩΜΧΝ中第m行η列像素点的回波延时，记为fgn,k，Ρα%为步骤4.1计算出的第k个匀速直线APC，
为步骤4.2计算出的像素点网格ΩΜΧΝ中第m行η列像素点的位置向量,m = 1,2,...,
M，M为步骤I定义的像素点网格ΩΜΧΝ中X方向的像素点点数，η = 1,2,…，N，N为步骤I定义的像素点网格ΩΜΧΝ中Y方向的像素点点数，k = 1,2,…，K，K为步骤I定义的慢时刻数，C为步骤I定义的光速。
[0044]步骤4.5、利用步骤4.4中第k个APC与像素点网格Ω MXN中第m行η列像素点的回波延时^^^，在SSkxp第k行数据里找到对应的回波信号值，记为A0mntk, 为步骤4.4计算出的第k个APC与像素点网格Ωμχν中第m行η列像素点的回波延时，SSkxp为步骤3计算出的升采样数据矩阵，m= 1，2，…，Μ，Μ为步骤I定义的像素点网格ΩΜΧΝ中X方向的像素点点数，η = 1,2,…，N, N为步骤I定义的像素点网格Ωμχν中Y方向的像素点点数，k = 1,2,…，K，K为步骤I定义的慢时刻数，C为步骤I定义的光速。
[0045]步骤4.6、采用公式= 2|P°fe+- P二|/C计算第k个APC与像素点网格Ωμχν中第m行η列像素点矫正后的回波延时，记为τ二,fc，为步骤4.1计算出的第k个匀速直线APC，ΔΡ&为步骤4.3计算出的第k个APC误差向量，Pgn为步骤4.2计算出的像素点网格Ωμχν中第m行η列像素点的位置向量，m = 1,2,..., M, M为步骤I定义的像素点网格Ωμχν中X方向的像素点点数，η = 1,2,…，N, N为步骤I定义的像素点网格ΩΜΧΝ中Y方向的像素点点数，k= 1,2,…，K，K为步骤I定义的慢时刻数，C为步骤I定义的光速。
[0046]步骤4.7、采用公式
[0047]RL = Σ?=ι Re(A%,>k)cos(ojTfm, k) -1m{A°Uhk)sm^i°mi k)计算像素点网格
Ω MXN中第m行η列像素点后向投影值的实部，记为i?二’ fc为步骤4.5计算出的却对应的回波信号值，为步骤4.4计算出的第k个APC与像素点网格Ωμχν中第m行η列像素点的回波延时，ω为步骤I定义的雷达发射信号载频，Tfnra为步骤4.6计算出的第k
个APC与像素点网格ΩΜΧΝ中第m行η列像素点矫正后的回波延时，m = 1,2,…，M, M为步骤I定义的像素点网格ΩΜΧΝ中X方向的像素点点数，η = 1,2,…，N, N为步骤I定义的像素点网格Ωμχν中Y方向的像素点点数，k = 1,2,…，K，K为步骤I定义的慢时刻数，Re(.)表示取实部运算，Im(.)表示取虚部运算。
[0048]步骤4.8、采用公式
[0049]Ign = Σ|=ι Re(Afnn k)sh^Tlnk) —计算像素点网格
Ωμχν中第m行η列像素点后向投影值的虚部,记为为步骤4.5计算出的对应的回波信号值，却xra为步骤4.4计算出的第k个APC与像素点网格Ωμχν中第m行η列像素点的回波延时，ω为步骤I定义的雷达发射信号载频，4liMc为步骤4.6计算出的k个
APC与像素点网格ΩΜΧΝ中第m行η列像素点矫正后的回波延时,m = 1,2,...,]?,]?为步骤I定义的像素点网格ΩΜΧΝ中X方向的像素点点数，η = 1,2,…，Ν，Ν为步骤I定义的像素点网格ΩΜΧΝ中Y方向的像素点点数，k= 1,2,…，K，K为步骤I定义的慢时刻数，Re(.)表示取实部运算，Im(.)表示取虚部运算。
[0050]步骤4.9、采用公式
[0051 ] JRmrjJl — —o)Re^jsin— <y/yM(^4m7I,fc)cos(itJTmri大)计算 R^in 关于
tU勺偏导数，记为fc，为步骤4.7计算出的素点网格ΩΜΧΝ中第!!!行]!列像素点后向投影值实部，T^irhk为步骤4.6计算出的第k个APC与像素点网格ω MXN中第m行η
列像素点矫正后的回波延时，Agira为步骤4.5计算出的#对应的回波信号值，?二來为步骤4.4计算出的第k个APC与像素点网格ΩΜΧΝ中第m行η列像素点的回波延时，ω为步骤I定义的雷达发射信号载频，m = 1,2,…，M, M为步骤I定义的像素点网格ΩΜΧΝ中X方向的像素点点数，η = 1,2,…，Ν，Ν为步骤I定义的像素点网格ΩΜΧΝ中Y方向的像素点点数，k = 1,2,…，K，K为步骤I定义的慢时刻数，Re(.)表示取实部运算，Im(.)表示取虚部运算。
[0052]步骤4.10、采用公式
[0053]^mn,k —太)S01(6JTmn,k)计算/^^关于7冊认的偏导数，记为d/二,fc’ 坊m为步骤4.8计算出的像素点网格Ωμχν中第m行η列像素点后向投影值虚部，为步骤4.6计算出的第k个APC与像素点网格Ω MXN中第m行η列像素点矫正后的回波延时，为步骤4.5计算出的.^m,fc对应的回波信号值，fLijc为步骤
4.4计算出的第k个APC与像素点网格ΩΜΧΝ中第m行η列像素点的回波延时，ω为步骤I定义的雷达发射信号载频，m= 1,2,…，M，M为步骤I定义的像素点网格ΩΜΧΝ中X方向的像素点点数，η = 1,2,…，N，N为步骤I定义的像素点网格ΩΜΧΝ中Y方向的像素点点数，k = 1,2,..., K, K为步骤I定义的慢时刻数，Re (.)表示取实部运算，Im(.)表示取虚部运算。
[0054]步骤4.11、采用公式Ca = 4(Rmn ‘ dR^k + Ign.dlgn k)/c计算初始中间常量，记为《Sm为步骤4.7计算出的像素点网格ΩΜΧΝ中第mSn列像素点后向投影值的实部，dRgnik为步骤4.9计算出的R^m关于的偏导数，为步骤4.6计算出的第k个APC与像素点网格Ωμχν中第m行η列像素点矫正后的回波延时，为步骤4.8计算出的像素点网格Ωμχν中第m行η列像素点后向投影值的虚部，di二,k为步骤4.10计算出的^m关于的偏导数，4m,k为步骤4.6计算出的第k个APC与像素点网格Ω MXN中第
m行η列像素点矫正后的回波延时,m = 1,2,...,]?,]?为步骤I定义的像素点网格ΩΜΧΝ中X方向的像素点点数，η = 1,2,…，Ν，Ν为步骤I定义的像素点网格ΩΜΧΝ中Y方向的像素点点数，k= 1,2,..., K, K为步骤I定义的慢时刻数。
[0055]步骤4.12、采用公式Vf = -Σ^=?Σξ=ΛηΜ.计算第k个慢时刻对应的中间向量，记为Vf，为步骤4.11计算出的初始中间常量，e匕为步骤4.3计算出的中间单位向量，k = 1，2，…，K，K为步骤I定义的慢时刻数，M为步骤I定义的像素点网格Ωμχν中X方向的像素点点数，N为步骤I定义的像素点网格Ωμχν中Y方向的像素点点数。
[0056]步骤4.13、重复步骤4.1到步骤4.12，对所有的k = 1，2，…，K，K为步骤I定义的慢时刻数，计算出vf, V20,…，V益’ vf为第I个慢时刻对应的中间向量ν?力第2个慢时刻对应的中间向量，<为第K个慢时刻对应的中间向量。
[0057]步骤4.H、采用公式Vf0 = vf'vf7",.*.,ν|Τ计算初始梯度向量，记为.VfO，为步骤4.13计算出的第I个慢时刻对应的中间向量，v|为步骤4.13计算出的第2个慢时刻对应的中间向量，v|为步骤4.13计算出的第K个慢时刻对应的中间向量，(.)τ表示向量转置运算。
[0058]步骤4.15、采用公式dfl = -W0计算初始搜索方向，记为d°，为步骤4.14计算出的初始梯度向量。
[0059]步骤5、初始化迭代变量
[0060]步骤5.1、定义当前迭代次数为q，q = 0，1，2，…，Q，Q为步骤I定义的共轭梯度算法最大迭代次数,初始化q = O。
[0061]步骤5.2、定义第q次迭代APC误差向量力ΔΡ|，并初始化为ΔΡ^ = ΔΡ°, ΔΡ°为步骤I定义的初始APC误差向量，q为步骤5.1定义的当前迭代次数。
[0062]步骤5.3、定义第q次迭代搜索方向为(f，并初始化为(f = d°，d°为步骤4.15计算出的初始搜索方向，Q为步骤5.1定义的当前迭代次数。
[0063]步骤5.4、定义第q次迭代梯度向量为VF，并初始化为= W0 Vf0为步骤4.14计算出的初始梯度向量，q为步骤5.1定义的当前迭代次数。
[0064]步骤6、判断迭代是否结束
[0065]如果当前迭代次数q满足q > Q，Q为步骤I定义的共轭梯度算法最大迭代次数，则结束迭代，输出ΔΡ|, ΔΡΙ为步骤5.2定义的第q次迭代APC误差向量；如果当前迭代次数q满足q < Q，则继续执行步骤7。
[0066]步骤7、计算最佳搜索步长
[0067]利用Armijo算法计算第q次迭代最佳搜索步长，记为λ % λ ?的大小为I行、3XK列。
[0068]步骤8、计算第q+Ι次迭代APC误差向量
[0069]采用公式ΔΡΓ1 二 ΔΡ^ + 十算第q+Ι次迭代APC误差向量，记为ΔΡα9+1，AP^
为步骤5.2定义的第q次迭代APC误差向量，λ ?为步骤6计算出的第q次迭代最佳搜索步长，Cf为步骤5.3计算出的第q次迭代搜索方向，q为步骤5.1定义的当前迭代次数。
[0070]步骤9、计算第q+Ι次迭代搜索方向
[0071]步骤9.1、计算第q+Ι次迭代梯度向量，具体方法是:
[0072]步骤9.1.1、采用公式Pg1 = Ρω +V X ts0c)计算第q+Ι次迭代、第k个匀速直线
APC,记为，Pto为步骤I定义的雷达平台在零时刻的位置向量，V为步骤I定义的雷达平台速度向量，ts(k)为向量ts的第k个元素，k = 1,2,..., K, K为步骤I定义的慢时刻数，ts为步骤I定义的慢时刻向量，q为步骤5.1定义的当前迭代次数。
[0073]步骤9.1.2、采用公式 P=1 = Pc + [(m - M/2) X dx, (η — N/2) x dy,0 ]计算第
q+1次迭代、像素点网格Ωμχν中第m行η列像素点的位置向量,记为P=1, P。为步骤I定义的场景中心位置向量，m= 1，2，…，M，M为步骤I定义的像素点网格ΩΜΧΝ中X方向的像素点点数，η = 1,2,…，N，N为步骤I定义的像素点网格ΩΜΧΝ中Y方向的像素点点数，dx为步骤I定义的像素点网格ΩΜΧΝ中X方向的像素点间隔，dy为步骤I定义的像素点网格Ωμχν中Y方向的像素点间隔，q为步骤5.1定义的当前迭代次数。
[0074]步骤9.1.3、采用公式
[0075]
ΔΡ= = y^+1(3(k -1) + I), APf1(SC^-1)+ 2), APaq+1(3(/c _ I) + 3)]
计算第q+1次迭代、第k个APC误差向量，记为ΔΡ^1，ΔΡΓ1为步骤8计算出的第q+Ι次迭代APC误差向量，ΔΡ^+1(3(/? — I) + I)为Δρ:_+1中第3 (k_l)+1个元素，+1 (3 (k -1) + 2)为 ΔΡ$+1 中第 3 (k-1) +2 个兀素,APa (3(k — 1)+3)为ΔΡΓ1中第3(k-l)+3个元素，k= 1,2,…，K，K为步骤I定义的慢时刻数；采用公式
O =计算第q+Ι次迭代中间单位向量，记为C, Pg1为步骤
9.1.1计算出的第q+Ι次迭代、第k个匀速直线APC，P=1为步骤9.1.2计算出的第q+Ι次迭代、像素点网格ΩΜΧΝ中第m行η列像素点的位置向量,m = 1,2,...,]?, M为步骤I定义的像素点网格Ωμχν中X方向的像素点点数，η = 1，2，…，N, N为步骤I定义的像素点网格Ωμχν中Y方向的像素点点数，q为步骤5.1定义的当前迭代次数。
[0076]步骤9.1.4、采用公式= 2]?1 - P^l/C计算第q+Ι次迭代、第k个APC与像素点网格Ωμχν中第m行η列像素点的回波延时,记为为步骤9.1.1计算出的第q+Ι次迭代、第k个匀速直线APC，k = 1，2，...，K，K为步骤I定义的慢时刻数，P:1
为步骤9.1.2计算出的第q+Ι次迭代、像素点网格Ωμχν中第m行η列像素点的位置向量，m = 1,2, “.，Μ,M为步骤I定义的像素点网格ΩΜΧΝ中X方向的像素点点数，η = 1,2,...,N，N为步骤I定义的像素点网格Ωμχν中Y方向的像素点点数，C为步骤I定义的光速，q为步骤5.1定义的当前迭代次数。
[0077]步骤9.1.5、利用第q+Ι次迭代、第k个APC与像素点网格ΩΜΧΝ中第mRn列像素点的回波延时,在SSkxp第k行数据里找到对应的回波信号值，记为fi+i
THTZffC TfViXtH.为步骤9.1.4计算出的第q+Ι次迭代、第k个APC与像素点网格ΩΜΧΝ中第m行η列像素点的回波延时，SSkxp为步骤3计算出的升采样数据矩阵，m = 1，2，…，M, M为步骤I定义的像素点网格ΩΜΧΝ中X方向的像素点点数，η = I, 2,…，N, N为步骤I定义的像素点网格Ωμχν中Y方向的像素点点数，k = 1，2，...，K，K为步骤I定义的慢时刻数，q为步骤5.1定义的当前迭代次数。
[0078]步骤9.1.6、采用公式= 2jPa^1 + APjI1 — Pg^/C计算第q+Ι次迭代、第k个APC与像素点网格Ωμχν中第m行η列像素点矫正后的回波延时,记为d+Y, P=为
IiLi V^tx-
步骤9.1.1计算出的第q+Ι次迭代、第k个匀速直线APC，为步骤9.1.3计算出的第
q+Ι次迭代、第k个APC误差向量，P^1为步骤9.1.2计算出的第q+Ι次迭代、像素点网格
ΩΜΧΝ中第m行η列像素点的位置向量,m = 1，2,…，M,M为步骤I定义的像素点网格ΩΜΧΝ中X方向的像素点点数，η = 1，2，…，N, N为步骤I定义的像素点网格ΩΜΧΝ中Y方向的像素点点数，k = 1，2，...，K，K为步骤I定义的慢时刻数，C为步骤I定义的光速，q为步骤5.1定义的当前迭代次数。
[0079]步骤9.1.7、采用公式
[008。] c = zLi MCi )cos (Ci) - 二)咖(肌.二)计算像素点网格
Ω MXN中第q+Ι次迭代、第m行η列像素点后向投影值的实部,记为为步骤9.L 5
计算出的对应的回波信号值，为步骤9.1.4计算出的第q+Ι次迭代、第k个APC
与像素点网格ΩΜΧΝ中第mSn列像素点的回波延时，ω为步骤I定义的雷达发射信号载频,为步骤9.1.6计算出的第q+Ι次迭代、第k个APC与像素点网格ΩΜΧΝ中第m行η
列像素点矫正后的回波延时，m= 1，2，…，M，M为步骤I定义的像素点网格ΩΜΧΝ中X方向的像素点点数，η = 1,2,…，N，N为步骤I定义的像素点网格Ωμχν中Y方向的像素点点数，k = 1,2,..., K, K为步骤I定义的慢时刻数，Re(.)表示取实部运算，Im(.)表示取虚部运算，q为步骤5.1定义的当前迭代次数。
[0081]步骤9.1.8、采用公式
[0082]C1 = Σ?-=ι- lm{A<l；llk)cos{(m(j^k)计算像素点网格ΩΜΧΝ中第q+Ι次迭代、第m行η列像素点后向投影值的虚部，记为=，为步骤
9.1.5计算出的?ΓΛ对应的回波信号值，为步骤9.1.4计算出的第q+Ι次迭代、第k
个APC与像素点网格ΩΜΧΝ中第mSn列像素点的回波延时，ω为步骤I定义的雷达发射信号载频，τ?为步骤9.1.6计算出的第q+Ι次迭代、k个APC与像素点网格ΩΜΧΝ中第m
行η列像素点矫正后的回波延时，m = 1,2,...,]?,]?为步骤I定义的像素点网格ΩΜΧΝ中X方向的像素点点数，η = 1,2,…，N, N为步骤I定义的像素点网格Ωμχν中Y方向的像素点点数，k = 1,2,…，K，K为步骤I定义的慢时刻数，Re(.)表示取实部运算，Im(.)表示取虚部运算，q为步骤5.1定义的当前迭代次数。
[0083]步骤9.1.9、采用公式
[0084]dR:k = —ω/?￡>(Οη(ωΟ —计算 ^
关于的偏导数，记为—ΓΛ，为步骤9.1.7计算出的第q+Ι次迭代、素点网格Ωμχν中第m行η列像素点后向投影值实部，为步骤9.1.6计算出的第q+Ι次迭代、第k个APC与像素点网格ΩΜΧΝ中第m行η列像素点矫正后的回波延时为步骤9.1.5计算出的对应的回波信号值，为步骤9.1.4计算出的第q+Ι次迭代、第k个APC与像素点网格ΩΜΧΝ中第mSn列像素点的回波延时，ω为步骤I定义的雷达发射信号载频，m = 1,2, “.，Μ,M为步骤I定义的像素点网格ΩΜΧΝ中X方向的像素点点数，η = 1,2,...,N，N为步骤I定义的像素点网格ΩΜΧΝ中Y方向的像素点点数，k = 1,2,…，K，K为步骤I定义的慢时刻数，Re(.)表示取实部运算，Im(.)表示取虚部运算，q为步骤5.1定义的当前迭代次数。
[0085]步骤9.1.10、采用公式
[0086]dl=k = MRe(A^k)COS(ωτ^)—論n(0(6>d)计算 /=1 关于
C的偏导数，记为I=1为步骤9.1.8计算出的第q+Ι次迭代、像素点网格Ω_中第m行η列像素点后向投影值虚部，为步骤9.1.6计算出的第q+Ι次迭代、k个APC与像素点网格Ωμχν中第m行η列像素点矫正后的回波延时，All1k为步骤9.1.5计算出的对应的回波信号值，d为步骤9.1.4计算出的第q+Ι次迭代、第k个APC与像素点网格Ωμχν中第m行η列像素点的回波延时，ω为步骤I定义的雷达发射信号载频，m = 1，2，…，M，M为步骤I定义的像素点网格ΩΜΧΝ中X方向的像素点点数，η = 1，2，…，Ν，Ν为步骤I定义的像素点网格ΩΜΧΝ中Y方向的像素点点数，k = 1,2,…，K，K为步骤I定义的慢时刻数，Re(.)表示取实部运算，Im(.)表示取虚部运算，q为步骤5.1定义的当前迭代次数。
[0087]步骤9.1.11、采用公式^1fc = 4(i?^.dR^+I^1.d/=y/C计算第 q+1 次迭代中间常量，记为为步骤9.1.7计算出的像素点网格Ωμχν中第q+Ι次迭代、第m行η列像素点后向投影值的实部，为步骤9.1.9计算出的关于的偏导数，为步骤9.1.6计算出的第q+Ι次迭代、第k个APC与像素点网格ΩΜΧΝ中第m行η列像素点矫正后的回波延时，/If为步骤9.1.8计算出的第q+Ι次迭代、像素点网格Ωμχν中第m行η列像素点后向投影值的虚部，d!=k为步骤9.1.10计算出的/^l1关于"4?的偏导数，为步骤9.1.6计算出的第q+Ι次迭代、第k个APC与像素点网格ΩΜΧΝ中第m行
η列像素点矫正后的回波延时，m = 1,2,…，M，M为步骤I定义的像素点网格ΩΜΧΝ中X方向的像素点点数，η = 1,2,…，N，N为步骤I定义的像素点网格Ωμχν中Y方向的像素点点数，k = 1，2，…，K，K为步骤I定义的慢时刻数，q为步骤5.1定义的当前迭代次数。
[0088]步骤9.1.12、采用公式vf+1 = —计算第k个慢时刻对应的第q+Ι次迭代中间向量，记为Vf+1，为步骤9.1.11计算出的第q+Ι次迭代中间常量，<ι+Λ为步骤9.1.3计算出的第q+Ι次迭代中间单位向量，k = 1，2，…，K，K为步骤I
定义的慢时刻数，M为步骤I定义的像素点网格Ωμχν中X方向的像素点点数，N为步骤I定义的像素点网格ΩΜΧΝ中Y方向的像素点点数，q为步骤5.1定义的当前迭代次数。
[0089]步骤9.1.13、重复步骤9.1.1到步骤9.1.12，对所有的k = 1,2,…，K，K为步骤
I定义的慢时刻数，计算出<+1，Vf1,…，Vf1, Vf1为第I个慢时刻对应的第q+Ι次迭代中间向量，Vf+1为第2个慢时刻对应的第q+Ι次迭代中间向量，V|+1为第K个慢时刻对应的第q+Ι次迭代中间向量，Q为步骤5.1定义的当前迭代次数。
[0090]步骤9.1.14、采用公式VF+i = [vflT,v|+lT,.",Vf1Y计算第q+Ι次迭代梯度向量，记为VF+Svf1为步骤9.1.13计算出的第I个慢时刻对应的第q+Ι次迭代中间向量，Vf1为步骤9.1.13计算出的第2个慢时刻对应的第q+Ι次迭代中间向量，vf1为步骤
9.1.13计算出的第K个慢时刻对应的第q+Ι次迭代中间向量，(.)Τ表示向量转置运算，q为步骤5.1定义的当前迭代次数。
[0091]步骤9.2、采用公式& = IIWWI丨2/丨IWqII2计算中间参数，vri+i为步骤9.1.14计算的第q+Ι次迭代梯度向量，VP:为步骤5.4定义的第q次迭代梯度向量，q为步骤5.1定义的当前迭代次数，I I.I I表示求向量二范数运算。
[0092]步骤9.3、采用公式必+1 = -Vfi?+1 +A7CltH十算第q+Ι次迭代搜索方向，记为cT1，
▽f#1为步骤9.1.14计算出的第q+Ι次迭代梯度向量，β q为步骤9.2计算出的中间参数，Cf为步骤5.3定义的第q次迭代搜索方向，q为步骤5.1定义的当前迭代次数。
[0093]步骤10、更新迭代次数，进入下一次迭代
[0094]采用公式q — q+1更新下一次迭代次数，运算符一表示赋值操作，重复步骤6?步骤10，直到迭代结束。
[0095]本发明的创新点在于:基于图像强度最优进行合成孔径雷达天线相位中心误差估计，从而进行高精度运动误差补偿。本发明首先利用匀速直线天线相位中心进行粗聚焦后向投影成像，然后以图像强度最优为准则迭代调整天线相位中心误差，迭代结束后获得天线相位中心误差估计，最后利用匀速直线天线相位中心和估计出的天线相位中心误差进行高精度后向投影成像。
[0096]本发明的优点:与基于图像锐度最优的自聚焦BP算法相比，本发明考虑了相位误差的空变性，通过有效估计天线相位中心误差，对场景中每个像素点都能进行高精度的相位补偿；除此之外，由于不需要存储像素点在每个慢时刻的后向投影值，本发明所需内存容量比基于图像锐度最优的自聚焦BP算法大大减少，因此更适用于大场景、长孔径、高精度SAR成像。
【专利附图】
【附图说明】
[0097]图1是本发明流程图。
【具体实施方式】
[0098]本发明主要采用计算机仿真的方法进行验证，所有步骤、结论都在MATLAB_R2010b上验证正确。具体实施步骤如下:
[0099]步骤1、相关参数的初始化
[0100]初始化的参数均为已知，且初始化的参数如下:光速为C = 3 XlOVs ;雷达发射线性调频信号，载频为ω = 2X X 109rad/s ;雷达发射脉冲的带宽为B = 300X 16Hz ；雷达发射脉冲的时宽为Tp = I X 1^6s ;雷达脉冲重复周期为T = 0.002s ;雷达回波距离向采样频率为fs = 390X 16Hz ;雷达回波数据矩阵为SK>a ;雷达回波数据矩阵SK>a的方位向点数和距离向点数分别为K = 512和L = 512, K也称为慢时刻数；慢时刻向量为ts =[-512/2,1-512/2，...，512/2-1] X 0.002s ;升采样数据矩阵为 SSKXP，SSkxp 的大小为 K 行、8XL列，K为雷达回波数据矩阵的方位向点数，L为雷达回波数据的距离向点数；几何坐标系为三维笛卡尔坐标系O-XYZ ;雷达平台速度向量为V = [0，100，0]m/s，V的大小为I行3列；雷达平台在零时刻的位置向量为Pto = [0，0，4000]m，Pto的大小为I行3列;0_ΧΥ平面内的矩形场景为? ;将?离散化为像素点网格，记为ΩΜΧΝ;像素点网格ΩΜΧΝ中X方向和Y方向的像素点点数分别为M = 60和N = 60 ;像素点网格ΩΜΧΝ中X方向和Y方向的像素点间隔分别为dx = 0.5m和dy = 0.5m ;场景中心位置向量为P。= [3000，0，0]m，Pc的大小为I行3列；共轭梯度算法最大迭代次数为Q = 50，Q为正整数；初始APC误差向量为ΔΡ|，ΔΡ| = O1x3k, Oix3k为I行、3XK列的零矩阵，K为雷达回波数据矩阵Skxl的方位向点数。
[0101]步骤2、对雷达回波数据矩阵的每一行进行脉冲压缩
[0102]取出所有雷达回波数据SK>a，利用脉冲压缩方法对SK>a的每一行进行脉冲压缩，得到脉冲压缩后的数据矩阵psK>a。
[0103]步骤3、对脉冲压缩后的数据矩阵的每一行进行频域升采样
[0104]对步骤2得到的脉冲压缩后的数据矩阵PSK>a的每一行统一做如下8倍频域升采样处理:
[0105]步骤3.1、取出步骤2中脉冲压缩后的数据矩阵PSkx1j的第k行向量,记为sk, k =1，2，...，K，K为步骤I定义的慢时刻数。
[0106]步骤3.2、对向量Sk作快速傅里叶变换(FFT)，得到向量fk。
[0107]步骤3.3、从向量fk的L/2+1位置开始插入7XL个零，得到向量zk，Zk= [fk⑴，fk(2), , fk (L/2), 01X7L, fk (L/2+1), , fk(L)], fk(l)为向量 fk 中的第 I 个元素，fk(2)为向量fk中的第2个元素，fk(L/2)为向量fk中的第L/2个元素,01Χι为I行、7XL列的零向量，fk(L/2+1)为向量fk中的第L/2+1个元素，fk(L)为向量fk中的第L个元素，L为步骤I定义的雷达回波数据矩阵距离向点数。
[0108]步骤3.4、对向量Zk作快速傅里叶逆变换(IFFT)，得到向量ssk，将向量ssk存放到矩阵SSkxp的第k行，SSkxp为步骤I定义的升采样数据矩阵。
[0109]步骤4、计算初始搜索方向
[0110]步骤4.1、采用公式巧^ = Pm + V X ts(fc)计算第k个匀速直线APC，记为Pttl为步骤I定义的雷达平台在零时刻的位置向量，V为步骤I定义的雷达平台速度向量，ts(k)为向量ts的第k个元素，k = 1,2,..., K, K为步骤I定义的慢时刻数，ts为步骤I定义的慢时刻向量。
[0111]步骤4.2、采用公式P二 = Pc + [O — M/2) X dx, (η — N/2) x dy, O ]计算像素点网格ΩΜΧΝ中第mSn列像素点的位置向量，记为Pfnn，P。为步骤I定义的场景中心位置向量，m = 1,2,..., M, M为步骤I定义的像素点网格Ωμχν中X方向的像素点点数，η = I,
2,..., N, N为步骤I定义的像素点网格Ωμχν中Y方向的像素点点数，dx为步骤I定义的像素点网格ΩΜΧΝ中X方向的像素点间隔，dy为步骤I定义的像素点网格ΩΜΧΝ中Y方向的像素点间隔。
[0112]步骤4.3、米用公式
[0113]ΔΡα% = [ΔΡ?(3(/? -1) + I), ΔΡ0°(3(/? -1) + 2), ΔΡ°(3(^ -1) + 3)]
计算第第k个APC误差向量，记为ΔΡα%，ΔΡ公为步骤I定义的初始APC误差向量，AP°(3(fc -1) + I)为 ΔΡ5 中第 3(k-l)+l 个元素，AP°(3(fc -1) + 2)为 ΔΡ| 中第 3 (k_l)+2
个元素，6?2(3伏-1) + 3)为4卩^中第3(1^-1)+3个元素，1^=1，2，...，1(，1(为步骤1定义的
po +ApO —po
慢时刻数;采用公式=计算中间单位向量，记为e0 SO为步骤4.1
I'fc衫1Vfc Fmnlemn,k，ra,k
计算出的第k个匀速直线APC, Ptl为步骤4.2计算出的像素点网格Ωμχν中第m行η列像素点的位置向量，m = 1,2,..., M, M为步骤I定义的像素点网格Ωμχν中X方向的像素点点数，η = 1,2,…，N, N为步骤I定义的像素点网格Ωμχν中Y方向的像素点点数。
[0114]步骤4.4、采用公式fLu = 2lP°fe - P二|/C计算第k个APC与像素点网格Ωμχν中第m行η列像素点的回波延时，记为巧?,为步骤4.1计算出的第k个匀速直线APC，
P?m为步骤4.2计算出的像素点网格ΩΜΧΝ中第m行η列像素点的位置向量,m = 1,2,...,
M，M为步骤I定义的像素点网格ΩΜΧΝ中X方向的像素点点数，η = 1,2,…，N，N为步骤I定义的像素点网格ΩΜΧΝ中Y方向的像素点点数，k = 1,2,…，K，K为步骤I定义的慢时刻数，C为步骤I定义的光速。
[0115]步骤4.5、利用在SSkxp第k行数据里找到对应的回波信号值，记为为步骤4.4计算出的第k个APC与像素点网格Ωμχν中第m行η列像素点的回波延时，SSkxp为步骤3计算出的升采样数据矩阵，m = 1，2，…，M，M为步骤I定义的像素点网格Ωμχν中X方向的像素点点数，η = 1,2,…，N, N为步骤I定义的像素点网格ΩΜΧΝ中Y方向的像素点点数，k= 1,2,…，K，K为步骤I定义的慢时刻数，C为步骤I定义的光速。
[0116]步骤4.6、采用公式= 2|P°fe+AP°fc - P;|/C计算第k个APC与像素点网格ΩΜΧΝ中第mSn列像素点矫正后的回波延时，记为<?太，P^k为步骤4.1计算出的第k个匀速直线ΑΡ(:，ΔΡ^为步骤4.3计算出的第k个APC误差向量P0nri力步骤4.2计算出的像素点网格Ωμχν中第m行η列像素点的位置向量，m = 1,2,..., M, M为步骤I定义的像素点网格Ωμχν中X方向的像素点点数，η = 1,2,…，N, N为步骤I定义的像素点网格ΩΜΧΝ中Y方向的像素点点数，k= 1,2,…，K，K为步骤I定义的慢时刻数，C为步骤I定义的光速。
[0117]步骤4.7、采用公式
[0118]R^ut = ZLi "?η大)cos(wtU -1m(A^nn k)si7i(MT°m k)计算像素点网格
Ω MXN中第m行η列像素点后向投影值的实部，记为为步骤4.5计算出的[^^对应的回波信号值，为步骤4.4计算出的第k个APC与像素点网格Ωμχν中第m行η列像素点的回波延时，ω为步骤I定义的雷达发射信号载频，为步骤4.6计算出的第k
个APC与像素点网格ΩΜΧΝ中第m行η列像素点矫正后的回波延时，m = 1,2,…，M, M为步骤I定义的像素点网格ΩΜΧΝ中X方向的像素点点数，η = 1,2,…，N, N为步骤I定义的像素点网格Ωμχν中Y方向的像素点点数，k = 1,2,…，K，K为步骤I定义的慢时刻数，Re(.)表示取实部运算，Im(.)表示取虚部运算。
[0119]步骤4.8、采用公式
[0120]I^in = ILi Re(A%inJi)Sin(OJTfmik) -计算像素点网格
Ω MXN中第m行η列像素点后向投影值的虚部，记为I’ 为步骤4.5计算出的〒^^对应的回波信号值，AiLa为步骤4.4计算出的第k个APC与像素点网格ΩΜΧΝ中第mSn列像素点的回波延时，ω为步骤I定义的雷达发射信号载频，T^irhk为步骤4.6计算出的k个
APC与像素点网格ΩΜΧΝ中第m行η列像素点矫正后的回波延时,m = 1,2,...,]?,]?为步骤I定义的像素点网格ΩΜΧΝ中X方向的像素点点数，η = 1,2,…，Ν，Ν为步骤I定义的像素点网格ΩΜΧΝ中Y方向的像素点点数，k= 1,2,…，K，K为步骤I定义的慢时刻数，Re(.)表示取实部运算，Im(.)表示取虚部运算。
[0121]步骤4.9、米用公式
[0122]dRmn Jc — - (O R emn s in〔co ? ~ ω/?τι(./1.7ηη 火)cos(cot 冊j,女)计算关于
^mn,k的偏导数，记为dhc，为步骤4.7计算出的素点网格Ω MXN中第m行η列像素点后向投影值实部Wtnii为步骤4.6计算出的第k个APC与像素点网格ΩΜΧΝ中第mSn列像素点矫正后的回波延时，AcLra为步骤4.5计算出的#对应的回波信号值，T0nnk为步骤4.4计算出的第k个APC与像素点网格ΩΜΧΝ中第m行η列像素点的回波延时，ω为步骤I定义的雷达发射信号载频，m= 1,2,…，M，M为步骤I定义的像素点网格ΩΜΧΝ中X方向的像素点点数，η = 1,2,…，N，N为步骤I定义的像素点网格Ωμχν中Y方向的像素点点数，k = 1,2,..., K, K为步骤I定义的慢时刻数，Re(.)表示取实部运算，Im(.)表示取虚部运算。
[0123]步骤4.10、采用公式
[0124],紀)siji(ciiTmn,釔)计算/^m 关于的偏导数，记为d/^fc，为步骤4.8计算出的像素点网格ΩΜΧΝ中第m行η列像素点后向投影值虚部，为步骤4.6计算出的第k个APC与像素点网格Ωμχν中第m行η列像素点矫正后的回波延时为步骤4.5计算出的Pn ifc对应的回波信号值，ftin,k为步骤4.4
计算出的第k个APC与像素点网格ΩΜΧΝ中第m行η列像素点的回波延时，ω为步骤I定义的雷达发射信号载频，m= 1，2，…，M，M为步骤I定义的像素点网格ΩΜΧΝ中X方向的像素点点数，η = 1,2,…，N，N为步骤I定义的像素点网格ΩΜΧΝ中Y方向的像素点点数，k = 1,2,..., K, K为步骤I定义的慢时刻数，Re (.)表示取实部运算，Im(.)表示取虚部运算。
[0125]步骤4.11、采用公式Ca = ^iRmn.dRQmn_k + fgn.U/C计算初始中间常量，记为RSm为步骤4.7计算出的像素点网格ΩΜΧΝ中第mSn列像素点后向投影值的实部，d/Cg为步骤4.9计算出的/?二关于勺偏导数彳_lk为步骤4.6计算出的第k个APC与像素点网格Ωμχν中第m行η列像素点矫正后的回波延时，Pnn力步骤4.8计算出的像素点网格ΩΜΧΝ中第m行η列像素点后向投影值的虚部，为步骤4.10计算出的從Tl关于Pnfc的偏导数，'dn.fc为步骤4.6计算出的第k个APC与像素点网格Ωμχν中第
m行η列像素点矫正后的回波延时,m = 1,2,...,]?,]?为步骤I定义的像素点网格ΩΜΧΝ中X方向的像素点点数，η = 1,2,…，Ν，Ν为步骤I定义的像素点网格ΩΜΧΝ中Y方向的像素点点数，k= 1,2,..., K, K为步骤I定义的慢时刻数。
[0126]步骤4.12、采用公式# = -Σ^=ιΣ^ι^0?α.egm,fc计算第k个慢时亥Ij对应的中间向量，记为vf，为步骤4.11计算出的初始中间常量，eQmnik为步骤4.3计算出的中间单位向量，k= 1,2,…，K，K为步骤I定义的慢时刻数，M为步骤I定义的像素点网格Ωμχν中X方向的像素点点数，N为步骤I定义的像素点网格ΩΜΧΝ中Y方向的像素点点数。
[0127]步骤4.13、重复步骤4.1到步骤4.12，对所有的k = 1，2，…，K，K为步骤I定义的慢时刻数，计算出vf, V20,…，v|, Vl°为第I个慢时刻对应的中间向量，v|为第2个慢时刻对应的中间向量，V益为第K个慢时刻对应的中间向量。
[0128]步骤4.14、采用公式Vf0 二 [vf,vf,…,vff计算初始梯度向量，记为可为步骤4.13计算出的第I个慢时刻对应的中间向量vf力步骤4.13计算出的第2个慢时刻对应的中间向量，站为步骤4.13计算出的第K个慢时刻对应的中间向量，(.)τ表示向量转置运算。
[0129]步骤4.15、采用公式dQ = -VfH十算初始搜索方向，记为d°，VfD为步骤4.14计算出的初始梯度向量。
[0130]步骤5、初始化迭代变量
[0131]步骤5.1、定义当前迭代次数为q，q = 0，1，2，…，Q，Q为步骤I定义的共轭梯度算法最大迭代次数,初始化q = O。
[0132]步骤5.2、定义第q次迭代APC误差向量为ΔΡ^并初始化为ΔΡ| = ΔΡ°, ΔΡ°为步骤I定义的初始APC误差向量，q为步骤5.1定义的当前迭代次数。
[0133]步骤5.3、定义第q次迭代搜索方向为(f，并初始化为(f = d°，d°为步骤4.15计算出的初始搜索方向，Q为步骤5.1定义的当前迭代次数。
[0134]步骤5.4、定义第q次迭代梯度向量为VfS并初始化为Vff? = Wq, Vf0为步骤
4.14计算出的初始梯度向量，q为步骤5.1定义的当前迭代次数。
[0135]步骤6、判断迭代是否结束
[0136]如果当前迭代次数q满足q > Q，Q为步骤I定义的共轭梯度算法最大迭代次数，贝IJ结束迭代，输出ΔΡ|，ΔΡ|为步骤5.2定义的第q次迭代APC误差向量；如果当前迭代次数q满足q < Q，则继续执行步骤7。
[0137]步骤7、计算最佳搜索步长
[0138]利用Armijo算法计算第q次迭代最佳搜索步长，记为λ % λ ?的大小为I行、3XK列。
[0139]步骤8、计算第q+Ι次迭代APC误差向量
[0140]采用公式ΔΡ^+1 = ΔΡ^ +^dtH十算第q+1次迭代APC误差向量，记为ΔΡ!+1, ΔΡ^
为步骤5.2定义的第q次迭代APC误差向量，λ ?为步骤6计算出的第q次迭代最佳搜索步长，Cf为步骤5.3计算出的第q次迭代搜索方向，q为步骤5.1定义的当前迭代次数。
[0141]步骤9、计算第q+Ι次迭代搜索方向
[0142]步骤9.1、计算第q+Ι次迭代梯度向量
[0143]步骤9.1.1、采用公式^1 = Pt() + V X ts(fc)计算第q+Ι次迭代、第k个匀速直线
APC,记为Pt。为步骤I定义的雷达平台在零时刻的位置向量，V为步骤I定义的雷达平台速度向量，ts(k)为向量ts的第k个元素，k = 1,2,..., K, K为步骤I定义的慢时刻数，ts为步骤I定义的慢时刻向量，q为步骤5.1定义的当前迭代次数。
[0144]步骤9.1.2、采用公式P=1 = Pc + [(?π — M/2) X dx,(n — N/2) X d}/,0 ]计算第q+1次迭代、像素点网格Ωμχν中第m行η列像素点的位置向量,记为P=1, P。为步骤I定义的场景中心位置向量，m = 1,2,…，M，M为步骤I定义的像素点网格Ωμχν中X方向的像素点点数，η = 1,2,..., N, N为步骤I定义的像素点网格ΩΜΧΝ中Y方向的像素点点数，dx为步骤I定义的像素点网格ΩΜΧΝ中X方向的像素点间隔，dy为步骤I定义的像素点网格ΩΜΧΝ中Y方向的像素点间隔，q为步骤5.1定义的当前迭代次数。
[0145]步骤9.1.3、采用公式
[0146]ΔΡ^；1 = [AP^+1(3(k -1) + I), ΔΡ?+1(3(Λ -1) + 2), AP^+1(3(k -1) + 3)]
计算第q+1次迭代、第k个APC误差向量，记为ΔΡ^\ ΔΡΓ1为步骤8计算出的第q+Ι次迭代APC误差向量，把+i(3(k — I) + I)为ΔΡΓ1中第3(k_l)+l个元素，
AP^+1(3(k — I) + 2)为 AP^+1 中第 3 (k-1) +2 个元素，ΔΡ兰+1(3〔k — I) + 3)为 ΔΡ^+1
中第3(k-l)+3个元素，k = 1，2，...，K，K为步骤I定义的慢时刻数；采用公式
^fnnlk =丨计算第q+Ι次迭代中间单位向量，记为Z+1 pq+ι为步骤9.1.1
\ a,k a,k mn |mn,k， a,k
计算出的第q+1次迭代、第k个匀速直线APC，P=1为步骤9.1.2计算出的第q+Ι次迭代、
像素点网格ΩΜΧΝ中第m行η列像素点的位置向量,m = 1,2,…，M,M为步骤I定义的像素点网格Ωμχν中X方向的像素点点数，η = 1,2,…，N, N为步骤I定义的像素点网格ΩΜΧΝ中Y方向的像素点点数，q为步骤5.1定义的当前迭代次数。
[0147]步骤9.1.4、采用公式= 2|Ρ^ - P^^l/C计算第q+1次迭代、第k个APC与像素点网格Ωμχν中第m行η列像素点的回波延时，记为PU1为步骤9.1.1计算出的第q+Ι次迭代、第k个匀速直线APC，k= 1,2,...,K,K为步骤I定义的慢时刻数，P^1
为步骤9.1.2计算出的第q+Ι次迭代、像素点网格Ωμχν中第m行η列像素点的位置向量，m = 1,2, “.，Μ,M为步骤I定义的像素点网格ΩΜΧΝ中X方向的像素点点数，η = 1,2,...,N，N为步骤I定义的像素点网格Ωμχν中Y方向的像素点点数，C为步骤I定义的光速，q为步骤5.1定义的当前迭代次数。
[0148]步骤9.1.5、利用在SSkxp第k行数据里找到对应的回波信号值，记为
A^nnik > 为步骤9.1.4计算出的第q+Ι次迭代、第k个APC与像素点网格ΩΜΧΝ中第m
行η列像素点的回波延时，SSkxp为步骤3计算出的升采样数据矩阵，m = 1,2,…，M, M为步骤I定义的像素点网格ΩΜΧΝ中X方向的像素点点数，η = 1,2,…，N, N为步骤I定义的像素点网格Ωμχν中Y方向的像素点点数，k = 1,2,...，K，K为步骤I定义的慢时刻数，q为步骤5.1定义的当前迭代次数。
[0149]步骤9.1.6、采用公式= 2\Ρ^ + APJJ1 — Pgtl/C计算第q+1次迭代、第k个APC与像素点网格Ωμχν中第m行η列像素点矫正后的回波延时，记为朽^1为步骤9.1.1计算出的第q+Ι次迭代、第k个匀速直线APC，Mt=1为步骤9.1.3计算出的第q+Ι次迭代、第k个APC误差向量，为步骤9.1.2计算出的第q+Ι次迭代、像素点网格ΩΜΧΝ
中第m行η列像素点的位置向量，m = 1,2,...,]?,]?为步骤I定义的像素点网格ΩΜΧΝ中X方向的像素点点数，η = 1,2,…，N, N为步骤I定义的像素点网格Ωμχν中Y方向的像素点点数，k = 1，2，...，K，K为步骤I定义的慢时刻数，C为步骤I定义的光速，q为步骤5.1定义的当前迭代次数。
[0150]步骤9.1.7、采用公式
[0151]R +1 = Σ?=ι ReUl^k)cos(MT^k) -5?η(ωτ.^)计算像素点网格
Ωμχν中第q+Ι次迭代、第m行η列像素点后向投影值的实部，记为β=1，为步骤9.1.5
计算出的对应的回波信号值，d为步骤9.1.4计算出的第q+Ι次迭代、第k个APC
与像素点网格ΩΜΧΝ中第mSn列像素点的回波延时，ω为步骤I定义的雷达发射信号载频为步骤9.!.6计算出的第q+Ι次迭代、第k个APC与像素点网格ΩΜΧΝ中第m行η
列像素点矫正后的回波延时，m = 1，2，…，M，M为步骤I定义的像素点网格Ωμχν中X方向的像素点点数，η = 1,2,…，N，N为步骤I定义的像素点网格Ωμχν中Y方向的像素点点数，k = 1,2,..., K, K为步骤I定义的慢时刻数，Re(.)表示取实部运算，Im(.)表示取虚部运算，q为步骤5.1定义的当前迭代次数。
[0152]步骤9.1.8、采用公式
[0153]C,1 = ILi 如(Oin(d) — /m(C)cos(w4:y 计算像素点网格
Ωμχν中第q+Ι次迭代、第m行η列像素点后向投影值的虚部,记为J=1，为步骤9.1.5
计算出的^^夂对应的回波信号值，为步骤9.1.4计算出的第q+Ι次迭代、第k个APC
与像素点网格ΩΜΧΝ中第mSn列像素点的回波延时，ω为步骤I定义的雷达发射信号载频，为步骤9.1.6计算出的第q+Ι次迭代、k个APC与像素点网格Ω MXN中第m行η列像素点矫正后的回波延时，m = 1,2,…，M，M为步骤I定义的像素点网格ΩΜΧΝ中X方向的像素点点数，η = 1,2,…，N，N为步骤I定义的像素点网格ΩΜΧΝ中Y方向的像素点点数，k = 1,2,..., K, K为步骤I定义的慢时刻数，Re (.)表示取实部运算，Im(.)表示取虚部运算，q为步骤5.1定义的当前迭代次数。
[0154]步骤9.1.9、采用公式
[0155]dR^k = -6>i?e(0(wd) - w/m(OsOO 计算关于的偏导数，记为/Ctl1为步骤9.1.7计算出的第q+1次迭代、素点网格Ωμχν中第m行η列像素点后向投影值实部，为步骤9.1.6计算出的第q+Ι次迭代、第k个APC与像素点网格Ωμχν中第m行η列像素点矫正后的回波延时，为步骤9.1.5计算出对应的回波信号值，为步骤9.1.4计算出的第q+Ι次迭代、第k个APC与像素点网格ΩΜΧΝ中第mSn列像素点的回波延时，ω为步骤I定义的雷达发射信号载频，m= 1,2,..., M, M为步骤I定义的像素点网格ΩΜΧΝ中X方向的像素点点数，η = 1,2,...,N，N为步骤I定义的像素点网格ΩΜΧΝ中Y方向的像素点点数，k = 1,2,…，K，K为步骤I定义的慢时刻数，Re(.)表示取实部运算，Im(.)表示取虚部运算，q为步骤5.1定义的当前迭代次数。
[0156]步骤9.1.10、采用公式
[0157]= MR^lm^C0K(0T^计算 4?1 关于
的偏导数，记为力步骤9.1.8计算出的第q+1次迭代、像素点网格Ωμχν中第m行η列像素点后向投影值虚部，为步骤9.1.6计算出的第q+Ι次迭代、第k个APC与像素点网格Ωμχν中第m行η列像素点矫正后的回波延时，A^k为步骤9.1.5计算出的对应的回波信号值，为步骤9.1.4计算出的第q+Ι次迭代、第k个APC与像素点网格ΩΜΧΝ中第mSn列像素点的回波延时，ω为步骤I定义的雷达发射信号载频，m= 1,2,..., M, M为步骤I定义的像素点网格ΩΜΧΝ中X方向的像素点点数，η = 1,2,...,N，N为步骤I定义的像素点网格ΩΜΧΝ中Y方向的像素点点数，k = 1,2,…，K，K为步骤I定义的慢时刻数，Re(.)表示取实部运算，Im(.)表示取虚部运算，q为步骤5.1定义的当前迭代次数。
[0158]步骤9.1.11、采用公式= A{Rltn.dRq:k + C1.d/=；y/C 计算第 q+l 次迭代中间常量，记为为步骤9.1.7计算出的像素点网格Ωμχν中第q+Ι次迭代、第m行η列像素点后向投影值的实部，为步骤9.1.9计算出的关于的偏导数，为步骤9.1.6计算出的第q+Ι次迭代、第k个APC与像素点网格ΩΜΧΝ中第m行η列像素点矫正后的回波延时，/^1为步骤9.1.8计算出的第q+Ι次迭代、像素点网格Ωμχν中第m行η列像素点后向投影值的虚部，dim为步骤9.L 10计算出的关于τΙι"?的偏导数，τ?1ι"?为步骤9.1.6计算出的第q+Ι次迭代、第k个APC与像素点网格ΩΜΧΝ中第m行
η列像素点矫正后的回波延时，m = 1,2,…，M，M为步骤I定义的像素点网格ΩΜΧΝ中X方向的像素点点数，η = 1,2,…，N，N为步骤I定义的像素点网格Ωμχν中Y方向的像素点点数，k = 1，2，…，K，K为步骤I定义的慢时刻数，q为步骤5.1定义的当前迭代次数。
[0159]步骤9.1.12、采用公式Vf+1 = -Ση=?Ση=Λη%.计算第k个慢时刻对应的第q+Ι次迭代中间向量，记为4+1，为步骤9.1.11计算出的第q+Ι次迭代中间常量，Ci为步骤9.1.3计算出的第q+Ι次迭代中间单位向量，k = 1,2,…，K，K为步骤I定义的慢时刻数，M为步骤I定义的像素点网格ΩΜΧΝ中X方向的像素点点数，N为步骤I定义的像素点网格ΩΜΧΝ中Y方向的像素点点数，q为步骤5.1定义的当前迭代次数。
[0160]步骤9.1.13、重复步骤9.1.1到步骤9.1.12，对所有的k = 1,2,…，K，K为步骤
I定义的慢时刻数，计算出vf1，v|+1，…，V§+1’ <+1为第I个慢时刻对应的第q+1次迭代中间向量，Vf1为第2个慢时刻对应的第q+Ι次迭代中间向量Vf1为第K个慢时刻对应的第q+Ι次迭代中间向量，Q为步骤5.1定义的当前迭代次数。
[0161]步骤9.1.14、采用公式VF+i = [4+1'4+1'.",4+1丫计算第q+Ι次迭代梯度向量，记为Vf#1，#+1为步骤9.1.13计算出的第I个慢时刻对应的第q+Ι次迭代中间向量，Vf1为步骤9.1.13计算出的第2个慢时刻对应的第q+Ι次迭代中间向量vf1为步骤
9.1.13计算出的第K个慢时刻对应的第q+Ι次迭代中间向量，(.)Τ表示向量转置运算，q为步骤5.1定义的当前迭代次数。
[0162]步骤9.2、采用公式尖=IIWi?+1丨|2/||Vff计算中间参数，vfc?+1为步骤9.1.14计算的第q+Ι次迭代梯度向量，Vfi?为步骤5.4定义的第q次迭代梯度向量，q为步骤5.1定义的当前迭代次数，I I.I I表示求向量二范数运算。
[0163]步骤9.3、采用公式必+1 = —VF+1 +禹CR十算第q+Ι次迭代搜索方向，记为cT1，
Vfi?+1为步骤9.1.14计算出的第q+Ι次迭代梯度向量，β q为步骤9.2计算出的中间参数，Cf为步骤5.3定义的第q次迭代搜索方向，q为步骤5.1定义的当前迭代次数。
[0164]步骤10、更新迭代次数，进入下一次迭代
[0165]采用公式q — q+1更新下一次迭代次数，运算符一表示赋值操作，重复步骤6?步骤10，直到迭代结束。
[0166]通过本发明的具体实施可以看出，相比于基于图像锐度最优的自聚焦BP算法，本发明通过最优化方法估计天线相位中心误差，实现了对场景中每个像素点都进行高精度运动补偿，从而大大提高了成像精度；由于本发明中没有存储像素点在每个慢时刻的后向投影值，所以内存要求大大降低。
【权利要求】
1.一种合成孔径雷达高精度运动补偿方法，其特征是它包括如下步骤: 步骤1、相关参数的初始化 初始化的参数均为已知，且初始化的参数如下:光速为C ;雷达发射线性调频信号，载频为ω ;雷达发射脉冲的带宽为B ;雷达发射脉冲的时宽为Tp ;雷达脉冲重复周期为T ;雷达回波距离向采样频率为fs ;雷达回波数据矩阵为SK>a ;雷达回波数据矩阵SK>a的方位向点数和距离向点数分别为K和L(K和L均为正整数)，K也称为慢时刻数；慢时刻向量为ts=?-Κ/2,1~Κ/2, , Κ/2-1]ΧΤ ;升采样数据矩阵为 SSKXP，SSkxp 的大小为 K 行、8 X L 列，K为雷达回波数据矩阵的方位向点数，L为雷达回波数据的距离向点数；几何坐标系为三维笛卡尔坐标系O-XYZ ;雷达平台速度向量为V，V的大小为I行3列；雷达平台在零时刻的位置向量为Pt(l，Pt(l的大小为I行3列；0ΧΥ平面内的矩形场景为Θ ;将Θ离散化为像素点网格，记为ΩΜΧΝ;像素点网格ΩΜΧΝ中X方向和Y方向的像素点点数分别为M和N;像素点网格Ωμχν中X方向和Y方向的像素点间隔分别为dx和dy ;场景中心位置向量为P。，P。的大小为I行3列；共轭梯度算法最大迭代次数为Q，Q为正整数；初始APC误差向量为ΔΡ|，^P| = O1XBiT- Oixsk为I行、3ΧΚ列的零矩阵，K为雷达回波数据矩阵SKM的方位向点数； 步骤2、对雷达回波数据矩阵的每一行进行脉冲压缩 取出步骤I中所有雷达回波数据SKXy采用传统的脉冲压缩方法对SK>a的每一行进行脉冲压缩，得到脉冲压缩后的数据矩阵PSK>a ； 步骤3、对脉冲压缩后的数据矩阵的每一行进行频域升采样 对步骤2得到的脉冲压缩后的数据矩阵PSK>a的每一行统一做如下8倍频域升采样处理: 步骤3.1、取出步骤2中脉冲压缩后的数据矩阵PSloa的第k行向量，记为sk，k = 1，.2，...，K，K为步骤I定义的慢时刻数； 步骤3.2、对向量Sk做传统的快速傅里叶变换(FFT)，得到向量fk ； 步骤3.3、从向量fk的L/2+1位置开始插入7XL个零，得到向量zk，zk = [fk(l)，fk(2),..., fk(L/2) ,01X7L, fk(L/2+l),..., fk (L)], fk(l)为向量 fk 中的第 I 个元素，fk(2)为向量fk中的第2个元素，fk(L/2)为向量fk中的第L/2个元素,01Χι为I行、7XL列的零向量，fk(L/2+l)为向量fk中的第L/2+1个元素，fk(L)为向量fk中的第L个元素，L为步骤I定义的雷达回波数据矩阵距离向点数； 步骤3.4、对向量Zk做传统的快速傅里叶逆变换(IFFT)，得到向量ssk，将向量SSk存放到矩阵SSkxp的第k行，SSkxp为步骤I定义的升采样数据矩阵； 步骤4、计算初始搜索方向 步骤4.1、采用公式= Pr0+VX 计算第k个匀速直线APC，记为Pto为步骤I定义的雷达平台在零时刻的位置向量，V为步骤I定义的雷达平台速度向量，ts(k)为向量ts的第k个元素，k = 1,2,..., K, K为步骤I定义的慢时刻数，ts为步骤I定义的慢时刻向量； 步骤4.2、采用公式PiL = Pc + [(；m — M/2) K dx, (η — N/2) X #,0 ]计算像素点网格Ωμχν中第m行η列像素点的位置向量，记为P。为步骤I定义的场景中心位置向量，m= 1,2,..., M, M为步骤I定义的像素点网格ΩΜΧΝ中X方向的像素点点数，η = 1,2,...,N, N为步骤I定义的像素点网格Ωμχν中Y方向的像素点点数，dx为步骤I定义的像素点网格ΩΜΧΝ中X方向的像素点间隔，dy为步骤I定义的像素点网格ΩΜΧΝ中Y方向的像素点间隔； 步骤4.3、采用公式 ?Ρ|Λ = [ΔΡ|(3(Ι?-?)-?), AP|(3(fc-l) + 2)i APfle(3(fc — i) + 3〕]计算第第 k个APC误差向量，记为ΔΡ乙，ΔΡΠ0为步骤I定义的初始APC误差向量，APc°(3(fe -1)十I〕为APc0中第3(k_l)+l个元素，ΔΡα0(3^-1)十2)为ΔΡα°中第3(k_l)+2个元素，ΔΡβ°(3(& _ 1) + 3)为ΔΡ:中第3 (k-1) +3个元素，k = 1，2，...，K，K为步骤I定义的慢时
刻数；采用公式=计算中间单位向量，记为#“，为步骤4.1计算出的第k个匀速直线APC, 为步骤4.2计算出的像素点网格ΩΜΧΝ中第!11行11列像素点的位置向量，m = 1,2,..., M, M为步骤I定义的像素点网格Ωμχν中X方向的像素点点数，η = 1,2,…，N, N为步骤I定义的像素点网格Ωμχν中Y方向的像素点点数； 步骤4.4、采用公式= 2|Ρ^ — Ρ5|β|/€计算第k个APC与像素点网格ΩΜΧΝ中第m行η列像素点的回波延时，记为为步骤4.1计算出的第k个匀速直线APC，PL为步骤4.2计算出的像素点网格ΩΜΧΝ中第m行η列像素点的位置向量,m = 1,2, *..,Μ,Μ为步骤I定义的像素点网格Ωμχν中X方向的像素点点数,η = 1,2,…,N,N为步骤I定义的像素点网格Ωμχν中Y方向的像素点点数，k = 1,2,…，K，K为步骤I定义的慢时刻数，C为步骤I定义的光速； 步骤4.5、利用步骤4.4中第k个APC与像素点网格Ωμχν中第m行η列像素点的回波延时’在SSkxp第k行数据里找到+f二k对应的回波信号值，记为为步骤.4.4计算出的第k个APC与像素点网格ΩΜΧΝ中第!11行11列像素点的回波延时，SSkxp为步骤3计算出的升采样数据矩阵，m= 1,2,…，Μ，Μ为步骤I定义的像素点网格ΩΜΧΝ中X方向的像素点点数，η = 1,2,…，N，N为步骤I定义的像素点网格Ωμχν中Y方向的像素点点数，k = 1,2,…，K，K为步骤I定义的慢时刻数，C为步骤I定义的光速； 步骤4.6、采用公式= 2lP|,fe + &Plk — PiM|/C计算第k个APC与像素点网格ΩΜΧΝ中第!!!行]!列像素点矫正后的回波延时,记为τ0 pt>为步骤4.1计算出的第k个
，* a,k匀速直线APC，AP|,k为步骤4.3计算出的第k个APC误差向量，Pl?为步骤4.2计算出的像素点网格ΩΜΧΝ中第m行η列像素点的位置向量,m = 1,2,…，M,M为步骤I定义的像素点网格ΩΜΧΝ中X方向的像素点点数，η = 1，2，…，Ν，Ν为步骤I定义的像素点网格Ωμχν中Y方向的像素点点数，k= 1,2,…，K，K为步骤I定义的慢时刻数，C为步骤I定义的光速；步骤 4.7、采用公式Ulw = Σ|=ιJfe)cos(o>[U -1m(i4U0i(a>rU计算像素点网格Ωμχν中第m行η列像素点后向投影值的实部,记为，Λ『Μ.ι|￡为步骤4.5计算出的对应的回波信号值，为步骤4.4计算出的第k个APC与像素点网格Ωμχν中第m行η列像素点的回波延时，ω为步骤I定义的雷达发射信号载频，为步骤4.6计算出的第k个APC与像素点网格ΩΜΧΝ中第!11行11列像素点矫正后的回波延时，m= 1，.2，…，M，M为步骤I定义的像素点网格ΩΜΧΝ中X方向的像素点点数，η = 1，2，…，Ν，Ν为步骤I定义的像素点网格ΩΜΧΝ中Y方向的像素点点数，k = 1,2,…，K，K为步骤I定义的慢时刻数，Re(.)表示取实部运算，Im(.)表示取虚部运算； 步骤 4.8、采用公式,/^ =计算像素点网格Ωμχν中第m行η列像素点后向投影值的虚部,记为Cn，aLiJc为步骤4.5计算出的1^?,Ir对应的回波信号值，为步骤4.4计算出的第k个APC与像素点网格ΩΜΧΝ中第m行η列像素点的回波延时，ω为步骤I定义的雷达发射信号载频，为步骤4.6计算出的k个APC与像素点网格ΩΜΧΝ中第!11行11列像素点矫正后的回波延时，m= 12，…，M，M为步骤I定义的像素点网格ΩΜΧΝ中X方向的像素点点数，η = 1,2,…，N，N为步骤I定义的像素点网格ΩΜΧΝ中Y方向的像素点点数，k = 1,2,…，K，K为步骤I定义的慢时刻数，Re(.)表示取实部运算，Im(.)表示取虚部运算； 步骤 4.9、采用公式- aim)cos(￡OT,|.Hj￡)计算关于T,U勺偏导数，记为, J^j3为步骤4.7计算出的素点网格ΩΜΧΝ中第m行η列像素点后向投影值实部，为步骤4.6计算出的第k个APC与像素点网格Ω MXN中第m行η列像素点矫正后的回波延时，AU步骤4.5计算出的对应的回波信号值，为步骤4.4计算出的第k个APC与像素点网格ΩΜΧΝ中第mSn列像素点的回波延时，ω为步骤I定义的雷达发射信号载频，m= 1，2，…，M，M为步骤I定义的像素点网格Ωμχν中X方向的像素点点数，η = 1,2,…，N, N为步骤I定义的像素点网格ΩΜΧΝ中Y方向的像素点点数，k= 1,2,..., K, K为步骤I定义的慢时刻数,Re(.)表示取实部运算，Im(.)表示取虚部运算； 步骤 4.10、采用公式rU - —?GCnA)#.计算关于Wmjc的偏导数，记为&1^， iL.为步骤4.8计算出的像素点网格ΩΜΧΝ中第m行η列像素点后向投影值虚部，+An.k为步骤4.6计算出的第k个APC与像素点网格Ωμχν中第m行η列像素点矫正后的回波延时，为步骤4.5计算出的对应的回波信号值，fLi,k为步骤4.4计算出的第k个APC与像素点网格Ωμχν中第m行η列像素点的回波延时，ω为步骤I定义的雷达发射信号载频，m = 1,2,…，M，M为步骤I定义的像素点网格Ωμχν中X方向的像素点点数，η = 1,2,…，N, N为步骤I定义的像素点网格ΩΜΧΝ中Y方向的像素点点数，k = 1,2,..., K, K为步骤I定义的慢时刻数，Re(.)表示取实部运算，Im(.)表示取虚部运算； 步骤4.11、采用公式逸计算初始中间常量，记为KL+Jc, 为步骤4.7计算出的像素点网格ω ΜΧΝ中第m行η列像素点后向投影值的实部，为步骤4.9计算出的关于'tU勺偏导数，T7LJr为步骤4.6计算出的第k个APC与像素点网格Ωμχν中第m行η列像素点矫正后的回波延时，为步骤4.8计算出的像素点网格Ωμχν中第m行η列像素点后向投影值的虚部，dlhk为步骤4.10计算出的C关于的偏导数，为步骤4.6计算出的第k个APC与像素点网格ΩΜΧΝ中第m行η列像素点矫正后的回波延时,m = 1,2,...,]?,]?为步骤I定义的像素点网格ΩΜΧΝ中X方向的像素点点数，η = 12，…，N, N为步骤I定义的像素点网格ΩΜΧΝ中Y方向的像素点点数，k = 1,2,..., K, K为步骤I定义的慢时刻数； 步骤4.12、采用公式= —Σ==1计算第k个慢时刻对应的中间向量，记为1<，为步骤4.11计算出的初始中间常量，为步骤4.3计算出的中间单位向量，k = 1，2，…，K，K为步骤I定义的慢时刻数，M为步骤I定义的像素点网格Ωμχν中X方向的像素点点数，N为步骤I定义的像素点网格ΩΜΧΝ中Y方向的像素点点数； 步骤4.13、重复步骤4.1到步骤4.12，对所有的k=l，2，…，K，K为步骤I定义的慢时刻数，计算出vf，v|,…，i4, vf为第I个慢时刻对应的中间向量，Vf为第2个慢时刻对应的中间向量，为第K个慢时刻对应的中间向量； 步骤4.14、采用公式Ffe=…,vff计算初始梯度向量，记为wo，为步骤.4.13计算出的第I个慢时刻对应的中间向量，为步骤4.13计算出的第2个慢时刻对应的中间向量，4为步骤4.13计算出的第K个慢时刻对应的中间向量，(.)τ表示向量转置运算； 步骤4.15、采用公式# = 一VfH十算初始搜索方向，记为d0, Vfe为步骤4.14计算出的初始梯度向量； 步骤5、初始化迭代变量 步骤5.1、定义当前迭代次数为q，q = 0,1,2,…，Q，Q为步骤I定义的共轭梯度算法最大迭代次数，初始化q = O ; 步骤5.2、定义第q次迭代APC误差向量为ΔΡ=，并初始化为APf = ΔΡο0, ΔΡ|为步骤I定义的初始APC误差向量，q为步骤5.1定义的当前迭代次数； 步骤5.3、定义第q次迭代搜索方向为d%并初始化为(f = d°，d°为步骤4.15计算出的初始搜索方向，q为步骤5.1定义的当前迭代次数； 步骤5.4、定义第q次迭代梯度向量为？F,并初始化为= Vf0, ?f°.为步骤4.14计算出的初始梯度向量，q为步骤5.1定义的当前迭代次数； 步骤6、判断迭代是否结束 如果当前迭代次数q满足q > Q，Q为步骤I定义的共轭梯度算法最大迭代次数，则结束迭代，输出—I，APa5为步骤5.2定义的第q次迭代APC误差向量；如果当前迭代次数q满足q < Q，则继续执行步骤7 ； 步骤7、计算最佳搜索步长 利用Armijo算法计算第q次迭代最佳搜索步长，记为λ % λ ?的大小为I行、3XK列； 步骤8、计算第q+Ι次迭代APC误差向量 采用公式ΔΡ^Γ1 = APJ+FdH十算第q+Ι次迭代APC误差向量，记为ΔΡ^1，ΔΡ:为步骤5.2定义的第q次迭代APC误差向量，λ ?为步骤6计算出的第q次迭代最佳搜索步长，Cf为步骤5.3计算出的第q次迭代搜索方向，q为步骤5.1定义的当前迭代次数； 步骤9、计算第q+Ι次迭代搜索方向 步骤9.1、计算第q+Ι次迭代梯度向量，具体方法是:步骤9.1.1、采用公式P^1 =Pfe+VXts (fe〕计算第q+Ι次迭代、第k个匀速直线APC，记为PqI' Pto为步骤I定义的雷达平台在零时刻的位置向量，V为步骤I定义的雷达平台速度向量，ts(k)为向量ts的第k个元素，k = 1,2,..., K, K为步骤I定义的慢时刻数，ts为步骤I定义的慢时刻向量，q为步骤5.1定义的当前迭代次数； 步骤9.1.2、采用公式P=1 =Pc + [Cm — M/2) X dx, {η-N/2) X dy, O }计算第 q+Ι 次迭代、像素点网格Ωμχν中第m行η列像素点的位置向量,记为P=1，P。为步骤I定义的场景中心位置向量，m = 1,2,…，M，M为步骤I定义的像素点网格Ωμχν中X方向的像素点点数，η = 12,..., N, N为步骤I定义的像素点网格Ω_中Y方向的像素点点数，dx为步骤I定义的像素点网格ΩΜΧΝ中X方向的像素点间隔，dy为步骤I定义的像素点网格Ωμχν中Y方向的像素点间隔，q为步骤5.1定义的当前迭代次数； 步骤9.1.3、采用公式 ΔΡ"；1= [APr1CSCfe-1) +I), AP^Hsik-1) + 2% ΔΡ3+1(3(λ — I) + 3)]计算第q+Ι次迭代、第k个APC误差向量，记为ΛΡ】;1, ΔΡΓ1为步骤8计算出的第q+Ι次迭代APC 误差向量，-1) + 1)为啤外中第 3 (k-1) +1 个元素，ΔΡ|+1(3(Λ-1) + 2}为APf1 中第 3 (k-1)+2 个元素，APm I) + 3〕为APf+1 中第 3 (k-1)+3 个元素，k =.1，2，...，K，K为步骤I定义的慢时刻数;采用公式.==丨丨计算第q+Ι次迭代中间单位向量，记为P=1为步骤9.1.1计算出的第q+Ι次迭代、第k个匀速直线APC，PJ1=1为步骤9.1.2计算出的第q+Ι次迭代、像素点网格ΩΜΧΝ中第m行η列像素点的位置向量,m= 1,2,…，M,M为步骤I定义的像素点网格ΩΜΧΝ中X方向的像素点点数，η = 1,2,..., N, N为步骤I定义的像素点网格Ωμχν中Y方向的像素点点数，q为步骤5.1定义的当前迭代次数； 步骤9.1.4、采用公式= 21 P! J1 — P=|/C计算第q+1次迭代、第k个APC与像素点网格ΩΜΧΝ中第m行η列像素点的回波延时,记为为步骤9.1.1计算出的第q+Ι次迭代、第k个匀速直线APC，k = 1，2，...，K，K为步骤I定义的慢时刻数，P=1为步骤9.1.2计算出的第q+Ι次迭代、像素点网格Ωμχν中第m行η列像素点的位置向量，m =.12，…，M，M为步骤I定义的像素点网格ΩΜΧΝ中X方向的像素点点数，η= 1，2，…，N，N为步骤I定义的像素点网格ΩΜΧΝ中Y方向的像素点点数，C为步骤I定义的光速，q为步骤5.1定义的当前迭代次数； 步骤9.1.5、利用第q+Ι次迭代、第k个APC与像素点网格ΩΜΧΝ中第m行η列像素点的回波延时?.?，在SSkxp第k行数据里找到对应的回波信号值，记为为步骤9.1.4计算出的第q+Ι次迭代、第k个APC与像素点网格Ω MXN中第m行η列像素点的回波延时，SSkxp为步骤3计算出的升采样数据矩阵，m = 1,2,…，M，M为步骤I定义的像素点网格ΩΜΧΝ中X方向的像素点点数，η = 12，…，Ν，Ν为步骤I定义的像素点网格ΩΜΧΝ中Y方向的像素点点数，k = 1,2,...,K,K为步骤I定义的慢时刻数，q为步骤5.1定义的当前迭代次数； 步骤9.1.6、采用公式.^ = 2]?1 + ΔΡ^1 — Ρ,ΙΠ/C计算第q+i次迭代、第k个APC与像素点网格Ωμχν中第m行η列像素点矫正后的回波延时,记为T^1k，P=1为步骤.9.1.1计算出的第q+Ι次迭代、第k个匀速直线APC，ΔΡ=1为步骤9.1.3计算出的第q+Ι次迭代、第k个APC误差向量，P=1为步骤9.1.2计算出的第q+Ι次迭代、像素点网格ΩΜΧΝ中第m行η列像素点的位置向量= 1,2,...,]?,]?为步骤I定义的像素点网格ΩΜΧΝ中X方向的像素点点数，η = 1,2,…，N, N为步骤I定义的像素点网格Ωμχν中Y方向的像素点点数，k = 1,2,…，K，K为步骤I定义的慢时刻数，C为步骤I定义的光速，q为步骤5.1定义的当前迭代次数； 步骤 9.1.7、采用公式 --=1 = Σξ=1^(Α^)εο3(ωτ^Ι) -1m(Al^k)sin(iaT^k)计算像素点网格Ωμχν中第q+Ι次迭代、第m行η列像素点后向投影值的实部，记为为步骤9.1- 5计算出的C1fe对应的回波信号值，为步骤9.1.4计算出的第q+Ι次迭代、第k个APC与像素点网格ΩΜΧΝ中第m行η列像素点的回波延时，ω为步骤I定义的雷达发射信号载频，为步骤9.1.6计算出的第q+Ι次迭代、第k个APC与像素点网格Ωμχν中第m行η列像素点矫正后的回波延时，m = 1,2,..., M, M为步骤I定义的像素点网格Ωμχν中X方向的像素点点数，η = 1，2，…，N, N为步骤I定义的像素点网格Ωμχν中Y方向的像素点点数，k = 1,2,..., K, K为步骤I定义的慢时刻数，Re(.)表示取实部运算，Im(.)表示取虚部运算，q为步骤5.1定义的当前迭代次数；步骤 9.1.8、采用公式,/,I=1 = if_=i MOinGd)—imCd1Xd)'计算像素点网格Ωμχν中第q+Ι次迭代、第m行η列像素点后向投影值的虚部，记为为步骤9.1.5计算出的?=对应的回波信号值，为步骤9.1.4计算出的第q+Ι次迭代、第k个APC与像素点网格ΩΜΧΝ中第m行η列像素点的回波延时，ω为步骤I定义的雷达发射信号载频，为步骤9.1.6计算出的第q+Ι次迭代、k个APC与像素点网格ΩΜΧΝ中第m行η列像素点矫正后的回波延时，m = 1,2,...,]?, M为步骤I定义的像素点网格Ωμχν中X方向的像素点点数，η = 1，2，…，N, N为步骤I定义的像素点网格Ωμχν中Y方向的像素点点数，k = 12,..., K, K为步骤I定义的慢时刻数，Re(.)表示取实部运算，Im(.)表示取虚部运算，q为步骤5.1定义的当前迭代次数； 步骤 9.1.9、采用公式=~ ω^?ΛΜ)εθ4ωΤ?2?)计算Λ=关于C的偏导数，记为战^，Λ=1为步骤9.1.7计算出的第q+Ι次迭代、素点网格ΩΜΧΝ中第!!!行]!列像素点后向投影值实部，d.为步骤9.1.6计算出的第q+Ι次迭代、第k个APC与像素点网格ΩΜΧΝ中第mSn列像素点矫正后的回波延时，为步骤.9.1.5计算出的对应的回波信号值，为步骤9.1.4计算出的第q+Ι次迭代、第k个APC与像素点网格ΩΜΧΝ中第mSn列像素点的回波延时，ω为步骤I定义的雷达发射信号载频，m= 1,2,…，M，M为步骤I定义的像素点网格ΩΜΧΝ中X方向的像素点点数，η=1,2,…，N，N为步骤I定义的像素点网格ΩΜΧΝ中Y方向的像素点点数，k = 1，2，…，K，K为步骤I定义的慢时刻数，Re (.)表示取实部运算，Im(.)表示取虚部运算，q为步骤.5.1定义的当前迭代次数； 步骤 9.1.10、采用公式dC = Me(Os(C) - ω/.(0.η(?^=)计算I=1关于'Ci的偏导数，记为为步骤9.1.8计算出的第q+Ι次迭代、像素点网格Ωμχν中第m行η列像素点后向投影值虚部，为步骤9.1.6计算出的第q+1次迭代、k个APC与像素点网格ΩΜΧΝ中第mSn列像素点矫正后的回波延时，为步骤.9.1.5计算出的*?二对应的回波信号值,C1fc为步骤9.1.4计算出的第q+Ι次迭代、第k个APC与像素点网格ΩΜΧΝ中第mSn列像素点的回波延时，ω为步骤I定义的雷达发射信号载频，m= 1,2,…，M，M为步骤I定义的像素点网格ΩΜΧΝ中X方向的像素点点数，η=1,2,…，N，N为步骤I定义的像素点网格ΩΜΧΝ中Y方向的像素点点数，k = 1，2，…，K，K为步骤I定义的慢时刻数，Re (.)表示取实部运算，Im(.)表示取虚部运算，q为步骤.5.1定义的当前迭代次数； 步骤9.1.H、采用公式.C + &.计算第q+1次迭代中间常量，记为】m 为步骤9.1.7计算出的像素点网格Ωμχν中第q+Ι次迭代、第m行η列像素点后向投影值的实部，O步骤9.1.9计算出的I?=关于的偏导数，rH1Jc为步骤9.1.6计算出的第q+Ι次迭代、第k个APC与像素点网格ΩΜΧΝ中第m行η列像素点矫正后的回波延时，4Γ为步骤9.1.8计算出的第q+Ι次迭代、像素点网格Ωμχν中第m行η列像素点后向投影值的虚部，为步骤9.1.10计算出的巧1关于的偏导数，为步骤9.L 6计算出的第q+Ι次迭代、第k个APC与像素点网格Ωμχν中第m行η列像素点矫正后的回波延时，m = 1，2，…，M，M为步骤I定义的像素点网格Ωμχν中X方向的像素点点数，η = 1,2,…，N，N为步骤I定义的像素点网格Ωμχν中Y方向的像素点点数，k = 1，2，…，K，K为步骤I定义的慢时刻数，q为步骤5.1定义的当前迭代次数；步骤9.1.12、采用公式< + i = ~Σ^=1Σξ=1Υ^.计算第k个慢时刻对应的第q+Ι次迭代中间向量，记为<+1，为步骤9.1.11计算出的第q+Ι次迭代中间常量， 为步骤9.1.3计算出的第q+Ι次迭代中间单位向量，k= 1，2，…，K，K为步骤I定义的慢时刻数，M为步骤I定义的像素点网格ΩΜΧΝ中X方向的像素点点数，N为步骤I定义的像素点网格ΩΜΧΝ中Y方向的像素点点数，q为步骤5.1定义的当前迭代次数； 步骤9.1.13、重复步骤9.1.1到步骤9.1.12，对所有的k = 1，2，…，K，K为步骤I定义的慢时刻数，计算出vf1，…，vf1，vf1为第I个慢时刻对应的第q+1次迭代中间向量，1为第2个慢时刻对应的第q+Ι次迭代中间向量，为第K个慢时刻对应的第q+Ι次迭代中间向量，Q为步骤5.1定义的当前迭代次数； 步骤9.1.14、采用公式？=+.“,<.+lTf计算第q+1次迭代梯度向量，记为+ 1为步骤9.1.13计算出的第I个慢时刻对应的第q+Ι次迭代中间向量，vf1为步骤9.1.13计算出的第2个慢时刻对应的第q+Ι次迭代中间向量，v|+1为步骤9.1.13计算出的第K个慢时刻对应的第q+Ι次迭代中间向量，(.)τ表示向量转置运算，q为步骤.5.1定义的当前迭代次数；步骤9.2、采用公式& =计算中间参数，为步骤9.1.14计算的第q+Ι次迭代梯度向量，为步骤5.4定义的第q次迭代梯度向量，q为步骤5.1定义的当前迭代次数，I I.I I表示求向量二范数运算； 步骤9.3、采用公式《1^1 = -7:^1十馬#计算第q+1次迭代搜索方向，记为cT1，▽F+1为步骤9.1.14计算出的第q+Ι次迭代梯度向量，β q为步骤9.2计算出的中间参数，Cf为步骤5.3定义的第q次迭代搜索方向，q为步骤5.1定义的当前迭代次数； 步骤10、更新迭代次数，进入下一次迭代 采用公式q — q+1更新下一次迭代次数，运算符一表示赋值操作，重复步骤6?步骤.10，直到迭代结束。
【文档编号】G01S7/40GK104181514SQ201410407342
【公开日】2014年12月3日 申请日期:2014年8月18日 优先权日:2014年8月18日
【发明者】张晓玲, 胡克彬, 熊海进, 赵韩星, 何蜀丰 申请人:电子科技大学