[发明专利]一种非匀速太赫兹雷达平台SAR成像方法

权利要求书

1.一种非匀速太赫兹雷达平台SAR成像方法，包括以下步骤：

步骤1：雷达发射太赫兹频段0.1THz～10THz线性调频信号

s(τ,t)=rect(τTp)exp(j2πfcτ+jπγτ2)]]>

其中，τ为快时间，t为慢时间，T_p为脉冲宽度，f_c为中心频率，选取在太赫兹频段，γ为调频斜率，表示虚数单位；

点目标坐标为(x_p,y_p,0)，在慢时间时刻t，雷达的位置为(x_a(t),0,H)，其中方位向坐标v为理想状态下的雷达速度，为雷达速度不均匀导致的位置偏离；则t时刻，点目标与雷达的瞬时距离R(t)为

R(t)=(vt+x~(t)-xp)2+yp2+H2]]>

雷达接收到的回波信号为

以雷达运动轨迹与场景中心的垂直距离作为参考距离R_ref，做去调频处理，输出信号为

sIF(τ,t)=rect(τ-2R(t)/cTp)exp{j(-4πγc(τ-2Rrefc)ΔR-4πfccΔR+4πγc2(ΔR)2)}]]>

其中ΔR＝R(t)-R_ref，式中相位项是去调频处理产生的剩余视频相位，在距离压缩前将其去除；

对去除剩余视频相位后的信号做距离向傅立叶变换，实现距离向的压缩，得到在距离频域上各点目标的sinc状脉冲信号；

步骤2：通过回波信号计算雷达在方位向各个采样时刻的位置

对步骤1中距离向傅立叶变换后的数据进行截取，获得特显点目标对应的数据矩阵，对所得数据矩阵在一个方位向上搜索，得到该方位向上的最大值以及最大值对应的位置，对每个方位向重复以上搜索得到所有方位向上的最大值和其对应的位置；

由前述所有方位向上的最大值对应的位置信息得到点目标与雷达的瞬时距离取其最小值为点目标与雷达的最近距离由公式

x^a(t)=(R^(t))2-Rmin2]]>

计算出每个方位向采样时刻t₀，t₁，…，t_N对应的方位向位置，进而得到对应的采样序列号k

k=[x^a(t)-x^a(t0)]·N-1x^a(tN)-x^a(t0)]]>

分别表示慢时间t取值为t₀、t_N时的方位向位置；

步骤3：基于步骤2得到的采样序列号，对步骤1中得到的脉冲信号数据x(ε)做方位向，即慢时间域的非均匀傅立叶变换，得到方位向频域均匀分布的信号，ε＝-N/2,…,N/2-1表示方位向序列号；

步骤3-1：构造矩阵

F=Nω-N/2-ωN/21-ω...ω-qN/2-ωqN/21-ωqωN/2-ω-N/21-ωN...ω-(q-1)N/2-ω(q-1)N/21-ωq-1............ωqN/2-ω-qN/21-ω-qω(q-1)N/2-ω-(q-1)N/21-ω-(q-1)...N]]>

其中，N为方位向数据长度，ω＝e^j2π/mN，m为信号过采样率，q为插值点数；

构造向量其中{mk}＝mk-[mk]，[mk]是mk的整数部分，n＝0,…,q；

s_ε为定标因子，取余弦定标因子，即

由方程组Fz_n(k)＝a_n(k)，求得在不同采样位置k处对应的q+1个插值系数z_n(k)；

步骤3-2：对方位向数据x(ε)，由公式计算傅立叶变换系数再利用快速傅立叶变换方法计算

通过定标因子进行数据变换，得到方位向频域数据

X(ϵ)=Φϵ·sϵ-1]]>

通过以上变换将回波信号的方位向变换到频域，即得到回波信号在方位向的频域形式；

步骤4：信号做STOLT插值处理并进行二维傅立叶反变换得到SAR图像；

将步骤3得到的数据在对应的波数域分布为(K_x,K_R)域，其中，K_x表示方位向波数，K_R表示点目标与雷达的视线波数，二者并非垂直关系，stolt变换公式为

Ky=KR2-Kx2]]>

通过距离向sinc插值即可得到变换结果S(K_x,K_y)；

步骤5：将步骤4中的结果S(K_x,K_y)做二维逆傅立叶变换，即得到最终的成像结果。