[发明专利]高速运动目标逆合成孔径雷达成像方法

摘要

本发明公开了一种高速运动目标逆合成孔径雷达成像方法。针对高速运动目标，首先利用最小二乘法和包络对齐估计目标运动参数，然后采用估计的运动参数对解线频调回波信号进行相干化处理，通过加权特征向量自聚焦算法对运动参数估计带来的相位误差予以消除，再利用楔形变换校正距离徙动，最后得到成像结果。现有的高速运动目标成像算法采用特显点法进行自聚焦，而本发明采用的加权特征向量自聚焦算法具有更好的聚焦效果和抗噪声性能，从而采用本发明能获得更好的成像效果。

主权项

一种高速运动目标逆合成孔径雷达成像方法，其特征在于步骤如下：第一步，雷达发射线性调频信号，该信号形式为其中Tr是脉冲时宽，fc是信号载频，γ是信号调频率，若已知信号带宽Br，则有关系γ＝Br/Tr，则目标上第i个散射点的雷达回波信号形式为：Sr(t~,tm)=Airect[t~-τTr]exp[j2πfc(t-τ)]exp[jπγ(t~-τ)2]---(1)]]>其中τ代表回波延时量，为瞬时距离，tm表示第m个慢时间，m＝1,2…M，表示快时间，Ri(tm)，VT(tm)分别表示目标上第i个散射点与雷达之间距离和径向速度；第二步，构造解线频调的参考信号为：Sref(t~,tm)=rect[t~-2Rs(tm)/cTref]exp[j2πfc(t-2Rs(tm)/c)]exp[jπγ(t~-2Rs(tm)/c)2]---(2)]]>其中Tref为参考信号的脉宽，Rs(tm)表示由窄带雷达测得的目标在tm时刻的距离；第三步，利用式(2)对回波信号进行解线频调处理，可得解线频调回波信号为：Sdechirp(t~,tm)=Sr(t~,tm)Sref*(t~,tm)=Airect[t~-2Ri(t~,tm)/cTr]exp{j[Φ1(tm)+Φ2(t~,tm)+Φ3(t~,tm)]}---(3)]]>其中Φ1(tm)＝Φ1a(tm)+Φ1b(tm)+Φ1c(tm)   (4)Φ1a(tm)=-4πcfcRΔ---(5)]]>Φ1b(tm)=-4πc[fc-VT(tm)cγ2Rs(tm)c]VT(tm)2Rs(tm)c---(6)]]>Φ1c(tm)=4πcγ2RΔcVT(tm)2Rs(tm)c---(7)]]>Φ2(t~,tm)=Φ2a(t~,tm)+Φ2b(t~,tm)+Φ2c(t~,tm)---(8)]]>Φ2a(t~,tm)=-4πcRΔγ(t~-2Rs(tm)c)---(9)]]>Φ2b(t~,tm)=-4πc{[fcγ+(1-2VT(tm)c)2Rs(tm)c]VT(tm)}γ(t~-2Rs(tm)c)---(10)]]>Φ2c(t~,tm)=2VT(tm)c4πcRΔγ(t~-2Rs(tm)c)---(11)]]>Φ3(t~,tm)=Φ3a(t~,tm)+Φ3b(t~,tm)---(12)]]>Φ3a(t~,tm)=4πc2γRΔ2---(13)]]>Φ3b(t~,tm)=4πcγ[VT(tm)c-1]VT(tm)[t~-2Rs(tm)c]2---(14)]]> RΔ＝Ri(tm)‑Rs(tm)   (15)第四步，精确估计目标运动参数第五步，利用精确估计的运动参数和对解线频调回波信号进行相干化处理，构造时域补偿函数H1(t~,tm)=exp{-j[Φ1b(tm)+Φ2b(t~,tm)+Φ3b(t~,tm)]},---(16)]]>H2(t~,tm)=exp{-j[4πc(fc+γ(t~-2Rs(tm)c))RΔso+4πc2γRΔso2]},---(17)]]>其中，RΔso＝Rs(tm)‑Ro(tm)，Ro(tm)，VT(tm)采用其精确估计值得到时域补偿后信号形式：S1(t~,tm)=Sdechirp(t~,tm)H1(t~,tm)H2(t~,tm)=Airect[t~-2Ri(tm)/cTr]exp[-j4πcγ(t~-2Rs(tm)c)RΔ′]exp(-j4πcfcRΔ′)exp(j4πcγRΔ′2)exp[j4πc2γ(-2RΔ′RΔso)]exp(jΦ1c(tm))---(18)]]>其中R′Δ＝Ri(tm)‑Ro(tm)，作傅里叶变换，则回波的距离频域表示为：S2(fr,tm)=AiTrsinc[Tr(fr+2γcRΔ′)]exp(-j4πcfcRΔ′)exp(j4πc2γRΔ′2)exp(j4πcfrRΔso)exp(jΦ1c(fr,tm))---(19)]]>其中Φ1c(fr,tm)是Φ1c(tm)在距离频域的表示，对式(19)所示信号进行相位补偿时只需要补偿fr＝‑2γR′Δ/c处的相位即可；结合(7)式，将式(19)中的Φ1c(fr,tm)改写成如下表达式：Φ1c(fr,tm)=4πc[2γRΔ′c-2γRΔsoc]VT(tm)2Rs(tm)c=-4πcfrVT(tm)2Rs(tm)c-4πcγ2RΔso(tm)cVT(tm)2Rs(tm)c---(20)]]>在距离向频域构造补偿函数其表达式如下式：H3(fr,tm)=exp(-jπfr2γ)exp(-j4πcRΔsofr)exp[-jΦ1c(fr,tm)]---(21)]]>则补偿后作距离向频域逆傅里叶变换(IFFT)可得：S3(t~0,tm)=IFFT[S2(fr,tm)·H3(fr,tm)]=Airect[t~0Tr]exp[-j4πc(fc+γt~0)RΔ′]---(22)]]>上式中此时运动目标回波已经被补偿为以目标上点O的运动轨迹Ro(tm)为参考的转台模型回波，即完成回波信号的相干化处理，并且同时补偿掉高速运动带来额外相位的影响；第六步，假设第四步精确估计的目标运动参数与真实值Ro(tm)之间存在误差RΔsoε(tm)，即则在通过上述时域和频域相干处理后，距离向作傅立叶变换有：S3′(f~0,tm)=AiTrsinc{Tr[f~0+2γcRΔ′]}exp(-j4πcfcRΔ′)exp(-j4πcfcRΔsoϵ(tm))---(23)]]>最后一个相位项就是由RΔsoε(tm)引入的误差相位，利用加权特征向量自聚焦法进行消除；第七步，对(23)式自聚焦后结果中的R′Δ进行泰勒展开，可得RΔ′=ri(0)+[vilos(0)tm+12v·ilos(0)tm2+16v··ilos(0)tm3+...]---(24)]]>上式中ri(0)和分别表示目标在等效成转台之后，目标上第i个散射点在初始时刻相对于参考点O的径向距离和径向速度，是一阶导数，此时，相干化处理之后信号变为：S4(f~0,tm)=Airect[f~0/γTr]exp[-j4πc(fc+f~0)ri(0)]exp[-j4πc(fc+f~0)vilos(0)tm]exp(-jΦi(f~0,tm))---(25)]]>上式中利用楔形变换(Keystone变换)进行距离徙动校正，定义时间变换关系作完距离徙动校正之后信号形式为：S5(f~0,τm)Airect[f~0/γTr]exp[-j4πc(fc+f~0)ri(0)]exp[-j4πcfcvilos(0)τm]exp(-jΦi′(f~0,τm))---(26)]]>其中高次相位第八步，通过加权特征向量自聚焦法消除(26)式中的高次相位，则经过自聚焦处理之后信号形式为：S6(f~0,τm)=Airect[f~0/γTr]exp[-j4πc(fc+f~0)ri(0)]exp[-j4πcfcvilos(0)τm]---(27)]]>第九步，进行二维成像；第六步加权特征向量自聚焦法的具体步骤如下：6.1、对(23)式所示各距离单元数据计算归一化幅度方差，第j个距离单元的归一化幅度方差定义为：σuj2=(uj-u‾j)2‾/uj2‾---(36)]]>其中，符号上一横表示对该距离单元内的各元素取平均值，表示第j个距离单元内各方位向数据的均值，表示第j个距离单元内各方位向数据的均方值；选取归一化幅度方差在0.2以内的距离单元记为特显点距离单元，设一共有L个特显点单元，第i个特显点距离单元的方位向信号记作xi，i＝1,2,...,L；6.2、将上述L个特显点单元的数据通过傅里叶变换变到图像域，得到各特显点单元的横向像；6.3、将特显点单元横向像的峰值圆周移位至图像中心；6.4、对6.3得到的各特显点单元数据进行加窗，窗宽的选取标准为取最大峰值下降10dB对应的宽度；6.5、对6.4得到的各特显点单元数据做逆傅立叶变换，变换到数据域，则此时忽略常数相位之后的第i个特显点单元信号xi表示为：xi=AiTrexp(-j4πcfcRΔsoϵ(tm)),m=1,2,...,M---(37)]]>其中M是方位向信号采样数，以第i(i＝1,2…L)个特显点单元的第一个慢时间采样数据为参考，该单元内第m个慢时间采样数据与第一个慢时间采样数据的相位差为ψm，则根据(37)，该单元M个方位向慢时间回波数据可表示为：gi1=ai+ηi1gi2=aiejψ2+ηi2...giM=aiejψM+ηiM---(38)]]>ai是该距离单元强散射点的信号复幅度，(38)式表示为向量形式： xi＝aiv+ηi   (39)其中ηi＝[ηi1,ηi2…ηiM]T是该距离单元的噪声序列，则表示的是不随距离单元变化的误差相位矢量；6.6、对L个特显点单元的信号进行加权，第i个特显点单元的权值为ϵi=SNR(i)Σi=1LxiHxiΣi=1L(SNR(i)xiHxi)---(40)]]>SNR(i)为第i个特显点单元的信噪比，利用(36)式所示的第i个特显点单元的归一化幅度方差来代替信噪比SNR(i)；6.7、加权之后二维信号数据域表示为：X＝[ε1x1,ε2x2,…εLxL]＝va+η  (41)a＝[ε1a1,ε2a2,…εLaL]是加权之后L个特显点单元幅度矢量，η＝[ε1η1,ε2η2,…εLηL]是加权之后L个特显点单元的噪声信号矢量，求解加权之后信号X的协方差矩阵C^=1LΣi=1Lϵi2xixiH---(42)]]>6.8、对协方差矩阵进行特征值分解，求取最大特征值对应的特征向量Φ1，则误差相位矢量v＝Φ1，利用该误差相位矢量对(23)式所示信号进行相位补偿，即可实现自聚焦操作。