[发明专利]一种对非协同目标高精度的ISAR运动补偿方法

权利要求书

1.一种对非协同目标高精度的ISAR运动补偿方法，其特征在于实现步骤如下：

（1）由步进频率逆合成孔径雷达获得大小为M*N的非协同目标回波数据E[m,n]，其中M表示传感器接收的总脉冲串组数，N表示每一组脉冲串中步进频率脉冲的个数；

（2）对目标回波数据沿距离向进行脉冲压缩，得到M组长度为N的距离像RP_m(n)，m＝0,1,…,M-1;n＝0,1,…,N-1，其中m表示脉冲串数，n表示脉冲数；

（3）定义第一组对准距离像式中RP₀(n)表示第一个组未对准距离像；

（4）由第m组已对准距离像和参考距离像RP_refm(n)，按如下公式计算出下一组参考距离像RP_refm+1(n)：

RPrefm+1(n)=mm+1RPrefm(n)+1m+1|RPm(n)‾|---(1)]]>

（5）将RP_m(n)右移τ_m+1，计算RP_refm+1(n)与RP_m+1(n-τ_m+1)之间的互相关函数：

EC(τm+1)=Σn=0N-1|RPrefm+1(n)|·|RPm+1(n-τm+1)|---(2)]]>

其中τ_m+1为[0,1,…,N-1]中的整数；

（6）计算EC(τ_m+1)取得最大值时，τ_m+1在[0,1,…,N-1]中的值，记为τ_m+1，0；

（7）根据Nelder-Mead法则，将τ_m+1，0作为迭代技术的初始猜测值，得到最佳距离徙动值

τ^m+1=argmaxτm+1EC(τm+1)---(3)]]>

不一定为整数；

（8）计算如下式所示：

RPm+1(n)‾=RPm+1(n-τ^m+1)---(4)]]>

当不是整数时，由傅里叶变换的平移性质，按如下方法实现上述过程：

RPm+1(n-τ^m+1)=FFT{ej(2π/N)τ^m+1qIFFT{RPm+1(n)}}---(5)]]>

q表示向量[0,1,…,N-1]^T，[·]^T表示向量转置；

（9）令m＝m+1，如果m＜M-1，重复以上（4）-（8）进行下一个距离像的对准；

（10）所有距离像对准完毕后，进行相位校正：

式中c为光速，f_n为第n个脉冲的频率，为矫正后相位。

2.根据权利要求1所述的一种对非协同目标高精度的ISAR运动补偿方法，其特征在于：所述步骤（7）的具体实现过程是：

（a）定义最优点：B＝τ_m+1，0，比较（2）中定义的互相关函数在τ_m+1，0+1和τ_m+1，0-1两点的取值，EC(τ_m+1，0±1)较大的定义为次优点G，较小的定义为差点W，令f(x)＝EC(x)表示在x点RP_refm+1(n)与RP_m+1(n-x)的互相关函数，x为任意点；

（b）取中心点：反射点：R＝2O-W，R代表反射点，O代表中心点；

（c）若f(R)＞f(G)，转到（i）；若f(R)≤f(G)，转到（ii）；

（i）反射或者扩展

如果f(B)＞f(R)，以R代替W；如果f(B)≤f(R)，计算扩展点E＝2R-O和f(E)，如果f(E)＞f(B)以E代替W；如果f(E)≤f(B)以R代替W；

（ii）压缩或者收缩

如果f(R)＞f(W)，以R代替W；f(R)≤f(W)则W不变，计算压缩点或者以及f(C)，如果f(C)＞f(W)，以C代替W；如果f(C)≤f(W)，计算收缩点S和f(S)，以S代替W,O代替G；

（d）通过以上步骤（i）-（ii）进行迭代，得到