[发明专利]一种阴影逆合成孔径雷达无线电全息信号提取方法

权利要求书

1.一种阴影逆合成孔径雷达无线电全息信号提取方法，其特征在于：

设定前向散射雷达FSR系统基线为L，工作波长为λ，载波角频率为ω₀，直达波信号幅度为U_d，直达波初始相位为初始频率为f₀；以系统基线方向为X轴，水平面上垂直于系统基线的方向为Y轴，以发射机为原点建立坐标系：目标运动速度为V，运动方向与基线夹角为目标沿基线方向的空间位置坐标为x，目标垂直于基线方向的空间位置坐标为y，目标运动速度在x方向和y方向的分量为V_x和V_y；目标从Q点穿越基线，Q点到发射机和接收机的距离分别为d_T和d_R，雷达观测时间为T_s，目标回波幅度为U_t(t)，目标回波多普勒信号频率为f_d(t)，多普勒相位为目标回波随机相位为通过发射机与接收机的位置设置使U_d＞＞U_t(t)；包括如下步骤：

步骤一，发射机的接收信号为目标回波与直达波的叠加：对所述的接收信号进行包络检波，在U_d＞＞U_t(t)的情况下，检波结果为目标回波分量与直达波分量的叠加，具体表示如下：

其中：

fd(t)=-1λ[xVx+yVyx2+y2+yVy-(L-x)Vx(L-x)2+y2]]]>

目标在x方向和y方向的速度分量和空间位置坐标分别为：

x＝d_T+V_xt，y＝V_yt

步骤二，所述的检波信号U(t)减去直达波分量：

忽略相位即为目标RHS的实部；

步骤三，对所述目标RHS的实部S_R进行分段Hilbert变换得到目标RHS的虚部S_I，具体步骤如下：

1）、S_R简化表示为S_R＝cos(t²)，进行时频分析，取零频率点对应的时刻t₀，观测时间在t₀前后对称分布，令t₀＝0；

2）、以所述的零频率点对应的时刻t₀为基准，将S_R分为t＜t₀、t＞t₀两段，对所述t＜t₀、t＞t₀两段信号分别进行Hilbert变换，变换后的信号t＜t₀部分乘以-1，两段信号组成初步的变换信号S_T，可以表示为：

ST(t)=-1π∫SR(τ)t-τdτ,t<t01π∫SR(τ)t-τdτ,t≥t0]]>

S_T在t₀周围存在幅度失真；

3）、对所述变换后的信号S_T求差分，选择距离零频率对应时刻t₀最近的两个极值点所对应的时刻t₁',t₂'，选定t₂'＞t₁'，则[t₁',t₂']之间的信号即为主瓣范围，对主瓣范围内的S_T进行曲线拟合，利用拟合曲线代替S_T原主瓣范围内曲线，重构出信号虚部S_I；

步骤四，利用信号实部S_R和信号虚部S_I构成复信号，即为目标RHS。

2.如权利要求1所述的一种阴影逆合成孔径雷达无线电全息信号提取方法，其特征在于，步骤三中的对主瓣范围内的S_T进行曲线拟合的方法为：

选择主瓣内的两个时间点：

t1=t1′+t02,]]>t2=t2′+t02]]>

在[t₁',t₂′]范围内选择一点t₃，得到三组参数对(t₁,S_T1),(t₂，S_T2),(t₃,S_T3)；利用所述的三组参数对，对主瓣范围内的目标RHS失真部分进行二次曲线拟合；

假设二次函数为y＝ax²+bx+c，将(t₁,S_T1),(t₂,S_T2),(t₃,S_T3)带入到二次函数y＝ax²+bx+c中，可以得到方程组：

at12+bt1+c=ST1at22+bt2+c=ST2at32+bt3+c=ST3]]>

求解方程组，得到：

a=(ST1-ST3)(t2-t3)-(ST2-ST3)(t1-t3)(t1-t3)(t1-t2)(t2-t3)b=(ST1-ST3)-a(t12-t32)t1-t3c=ST3-at32-bt3]]>

将上式中的参数a,b,c代入二次函数，得到拟合曲线。