[发明专利]一种前向散射雷达阴影逆合成孔径成像方法

权利要求书

1.一种前向散射雷达动目标阴影逆合成孔径成像方法，其特征在于，包括如下步骤：

步骤1、对前向散射雷达回波信号采用分数阶傅里叶变换的成像过程与采用逆菲涅尔变换的成像过程进行等效，获得使上述两种成像过程等效的条件，即为：输入等效条件、变换参数等效条件以及变换结果中幅度和相位的等效条件，具体为：

传统的基于逆菲涅尔变换的阴影逆合成孔径成像过程如下式表示：

H·(x′)=γ2πvQ·exp[-j2(γx′v)2]∫-Ts/2+Ts/2E·(t)exp[-j(γ2t2+γx′tv)]dt---(5)]]>

其中，E·(t)=Asinφjλrc1rc2exp(j2πLλ+jΩ2zp2)·∫-∞+∞H·(x′)exp[jΩ2v2sin2φ(x′v+t)2]dx′,]]>表示根据菲涅尔衍射得到的前向散射雷达的回波信号；

表示目标复轮廓函数，H·(x′)=∫-∞+∞ϵ(x′,z′)exp(jΩ2z′2+2jΩ2zpz′)dz′;]]>x′y′z′表示以目标中心点为原点的局部坐标系；和分别表示目标中心点到发射机和接收机的距离，(x_p，y_p，z_p)为目标的中心点在全局坐标系中的坐标；L为基线长度；T_s表示相参累计时间；ε(x′，z′)为假定目标垂直穿越基线时刻的目标阴影轮廓S的指示函数，且ϵ(x′,z′)=1,(x′,z′)∈S0,(x′,z′)∉S;]]>

Q·=AsinφjλdTdRexp(j2πLλ+jΩ2zp2),]]>其中，Ω=πλ(1dT+1dR),]]>A为回波信号幅度，φ为目标与基线夹角，d_T和d_R分别表示目标在穿越基线时刻到发射机和接收机的距离；λ表示回波信号波长；

γ＝2Ω²v²sin²φ表示目标中心点的多普勒频率变化率，即调频斜率，其中，v表示目标运动速度；t表示时间；

分数阶傅里叶变换的标准定义如下式所示：

Xα(u)=Aαexp[j(πu2cotα)]×∫-∞+∞x(t)exp[j(πt2cotα-2πutcscα)]dt---(8)]]>

其中α为分数阶域中的变换角度，u为分数阶域中的变换域坐标，x(t)为待变换的目标回波信号；

则，所述输入等效条件为：

所述变换参数等效条件为：cotα=-γ2π,ucscα=γx′2πv---(10)]]>

所述幅度和相位的等效条件为：

|Xα(x′)Aα|=|H·(x′)2πvQ·γ|]]>

arg{Xα(x′)exp[-j(πu2cotα)]/Aα}=arg{H·(x′)Q·exp[jγ2(x′v)2]}---(11)]]>

其中，X_α(·)表示分数阶傅里叶变换结果；

消除(10)式中量纲归一化因子对调频斜率γ及空域坐标x′的影响，进而得到分数阶域中变换角度α和变换域坐标u的表达式：

u=x′γcsc(arccotγTs-2πfs)---(12,)]]>

α=arccotγTs-2πfs---(13,)]]>

其中，f_s表示采样率；

步骤2、成像时，利用步骤1中的式(12’)和(13’)分别得到变换角度α和变换域坐标u的值，代入到式(8)中，同时，根据式(9)的输入等效条件，用前向散射雷达的回波信号作为式(8)中的待变换的目标回波信号，作变换角度为α的分数阶傅里叶变换，得到变换结果X_α(x′)；

步骤3、目标阴影轮廓提取，具体为：

首先，根据步骤2得到的变换结果X_α(x′)以及(11)式，对复轮廓函数中幅度和相位分别进行估计：

|H·^(x′)|=|Xα(x′)γAα2πvQ·|---(14)]]>

arg{H·^(x′)}=arg{Xα(x′)Aαexp[Λ]}---(15)]]>

其中，表示复轮廓函数估计值；Λ=-j(πu2cotα+π2-2πLλ-Ω2zp2-γx′22v2);]]>

然后，根据式(14)和(15)分别估计得到的幅度和相位，对目标侧影轮廓高度差像h(x′)和中线像m(x′)进行提取：

h(x′)=|H·^(x′)|sinc(k(1dT-1dR)zph(x′)/2π)---(17)]]>

m(x′)=arg{H·^(x′)}kzp(1dT+1dR)---(18)]]>

其中，表示波数。