[发明专利]一种球面波成像数学模型及近场效应的补偿方法

权利要求书

1.一种球面波成像数学模型及近场效应的补偿方法，其特征在于：包括如下步骤：

步骤1，对现有转台目标二维成像数学模型作改进如下：令

其中，g(x，y)为目标二维散射密度函数，目标与x-y坐标系绕O点顺时钟旋转，θ为u轴相对于x轴的旋转角度；x-y和u-v坐标系的原点均为O，R₀为雷达到目标旋转中心O的距离，R_θ为目标在相对于u轴转过角度θ时，目标上任一点到达雷达天线的距离；φ为R_θ与R₀之间的夹角；为目标x-y各点相对于坐标系的极坐标；弧线S代表到达雷达天线等距离的散射点的连线；

将R_θ做二阶泰勒级数做近似可得：

步骤2，由于实际情况一般满足小角度成像的条件，θ非常小，于是步骤1得到的R_θ可以化简为：

Rθ=R0+x+yθ+y2-2xyθ2(R0+x)]]>

步骤3，在步骤2的基础上，对项进一步做近似：

Rθ2≈(R0+x+y22R0+2x)2]]>

步骤4，在步骤3的基础上，对于频率步进体制测量雷达，满足k＝k₀+nΔk(Δk远小于k₀)，(其中k＝2f/c，f为电磁波的频率，c为光速，n＝0，1，…N-1)；θ＝mΔθ，m＝0，1，…M-1；设P为目标区域独立的散射点个数(P＜N，M)；对步骤1所述的现有转台目标二维成像数学模型进行离散化，并进一步忽略二阶小量，得：

z(nΔk,mΔθ)=∑p=0P-1g(xp,yp)ΔxΔy]]>

/(R0+xp+yp22R0+2xp)2]]>

·exp[-j2π(k₀+nΔk)R₀]

·exp[-j2π(k0+nΔk)·(xp+yp22R0+2xp)]]]>

·exp[-j2πk0mΔθ(yp-xpypR0+xp)]]]>

其中，x_p，y_p为第P个目标散射点的横纵坐标；上式即为改进后的近场目标二维成像的数学模型；

步骤5，进一步变形步骤4得到的数学模型；令：

z′(n,m)=z(nΔk,mΔθ)exp[j2π(k0+nΔk)R0]xp′=xp+yp22R0+2xpyp′=yp-xpypR0+xp]]>

由于ΔxΔy为常数量，在数学变换中将被合并为一系数，因此在模型中可以省略，则步骤4所述的改进后的近场目标二维成像的数学模型变为：

z′(n,m)=∑p=0P-1g(xp,yp)/(R0+xp′)2]]>

·exp(-j2πk₀x′_p)

·exp(-j2πnΔkx′_p)

·exp(-j2πk₀mΔθy′_p)

步骤6，在步骤5的基础上，进行子空间矩阵分裂；

由于z′(n，m)与z′(n+1，m)及z′(n，m+1)分别相差一个固定相位，则将z′(n，m)按一定方式排列，并形成三个子空间矩阵X，Y，Z；

根据步骤5所得z′(n，m)的表示式，可以将X，Y，Z表示成如下形式：

X=AS+NxY=AΦS+NyZ=AΘS+Nz]]>

X，Y，Z均为(N-1)×(M-1)维矩阵，A为(N-1)×(M-1)×P维矩阵，N_i(i＝x，y，z)为测量噪声矩阵；而且有：

A(n,m,p)=exp[-j2π(k0+nΔk)xp′]·exp(-j2πk0mΔθyp′)S=[g(x0,yp)/(R0+xp′)2,g(x1,y1)/(R0+x1′)2,···,g(xP,yP)/(R0+xP′)2]Φ=diag(exp(-j2πΔkx0′),exp(-j2πΔkx1′),···,exp(-j2πΔkxP′)Θ=diag(exp(-j2πk0Δθy0′),(exp(-j2πk0Δθy1′),···,(exp(-j2πk0ΔθyP′)]]]>

步骤7，在步骤6的基础上进行奇异值分解，得到目标散射点的个数；

步骤8，通过步骤7中奇异值矩阵U₁，V₁计算的特征值，得到步骤6中Φ和Θ的对角元素的值；

取步骤7得到的奇异值矩阵U₁，V₁中与P个散射点对应的奇异值的列矢量构成从而可以得到三个新矩阵：

Ex=UΛ1HXVΛ2Ey=UΛ1HYVΛ2Ez=UΛ1HZVΛ2]]>

并定义Eθ=Ex-1Ez]]>

计算E_θ的特征值就可以得到Φ和Θ的对角元素的值；

步骤9，在步骤8的基础上，确定变换矩阵Q使得对角化，并将这种变换应用于E_θ，即：

Q^HE_θQ＝T_θ

这里T_θ是一个近似上三角矩阵，和T_θ的主对角元素分别等于Φ和Θ的主对角元素；由于采用了相同的变换矩阵Q，和T_θ的主对角元素是一一对应的，故形成了参数的自动配对；

实际上T_θ已经非常接近上三角矩阵，再做进一步的变换反而可能造成参数失配情况，这里直接取T_θ的对角元素为特征值；

步骤10，在步骤9的基础上，得到目标散射点的位置估计坐标x_p和y_p；

步骤11，完成步骤10后，计算得到目标散射点的散射强度。