一种前向散射雷达地面运动目标信号建模及成像方法

摘要

本发明涉及一种运动目标前向散射信号建模及成像的信号处理方法，属于雷达信号处理技术领域。本发明所提出的方法考虑地面反射、大衍射角以及目标运动方向较为随机等实际情况，建立了精确的地面运动目标前向散射信号模型，修正了SISAR成像算法，并对目标像进行了更为精确的估计。修正后的成像算法在大衍射角、目标斜穿基线和多径干涉情况下依然可以获得精确的目标像，这些目标像与目标真实轮廓形状非常接近，可有效地对不同形状的目标进行分类识别。

主权项

1.一种前向散射雷达地面运动目标信号建模及成像方法，其特征在于包括如下具体步骤：步骤1、建立地面前向散射雷达系统的几何结构模型；xoy平面是地平面，发射机置于坐标系(x，y，z)的原点O，接收机在y轴B(0，L，0)处，发射机和接收机之间的距离为基线长度L，假设目标以速度v平行于xoy平面运动，与基线夹角为α_h，β_h分别表示从发射机和接收机方向观测的目标方位角，发射机到目标穿越基线位置的距离记为d_T，接收机到目标穿越基线位置的距离记为d_R，坐标系(x′，y′，z′)与坐标系(x，y，z)各轴相互平行，原点P(x_p，y_p，z_p)为目标的中心，T，R，T′，R′分别表示发射机、接收机、镜像发射机和镜像接收机天线相位中心；发射机和镜像发射机天线相位中心的高度是h_T；接收机和镜像接收机的天线高度分别为h_R，M₁和M₂是地面反射点；θ是反射波擦地角，假设P′为目标上任意一点，r₁和r₂分别表示发射机和接收机天线相位中心到该点的斜距，r′₁和r′₂分别表示镜像发射机和镜像接收机天线相位中心到该点的斜距；步骤2、对接收信号建模，在地面前向散射雷达系统中，接收到的信号认为是四条路径上传播信号的矢量和，这四条路径分别为：发射机到目标到接收机(以下称为路径1)；发射机到地面反射点M1到目标再到接收机(以下称为路径2)；发射机到目标到地面反射点M2到接收机(以下称为路径3)；发射机到M1到目标到M2到接收机(以下称为路径4)，四条路径的长度分别为：r₁+r₂，r′₁+r₂，r₁+r′₂，r′₁+r′₂，接收信号表示如公式(1)所示：$\stackrel{·}{E}={\stackrel{·}{E}}_1+\Gamma {\stackrel{·}{E}}_2+\Gamma {\stackrel{·}{E}}_3+{\Gamma }^2{\stackrel{·}{E}}_4$$=\frac{A({\cos \alpha }_h+{\cos \beta }_h)}{2\mathrm{j\lambda }{r}_{c1}{r}_{c2}}{\int \int }_s{e}^{\mathrm{jk}(r_1+r_2)}+{\mathrm{\Gamma e}}^{\mathrm{jk}(r_1^{\prime }+r_2)}+{\mathrm{\Gamma e}}^{\mathrm{jk}(r_1+r_2^{\prime })}+{\Gamma }^2{e}^{\mathrm{jk}(r_1^{\prime }+r_2^{\prime })}\mathrm{dS}---\left(1\right)$其中：$r_1\approx r_{c1}+\frac{x^{\prime 2}+{2x}^{\prime }{x}_P}{{2r}_{c1}}+\frac{z^{\prime 2}+{2z}^{\prime }{z}_P}{{2r}_{c1}}+\frac{h_T^2-{2h}_T(z_p+z^{\prime })}{{2r}_{c1}}---\left(2\right)$$r_2\approx r_{c2}+\frac{x^{\prime 2}+{2x}^{\prime }{x}_P}{{2r}_{c2}}+\frac{z^{\prime 2}+{2z}^{\prime }{z}_P}{{2r}_{c2}}+\frac{h_R^2-{2h}_R(z_p+z^{\prime })}{{2r}_{c2}}---\left(3\right)$$r_1^{\prime }\approx r_{c1}+\frac{x^{\prime 2}+{2x}^{\prime }{x}_P}{{2r}_{c1}}+\frac{z^{\prime 2}+{2z}^{\prime }{z}_P}{{2r}_{c1}}+\frac{h_T^2+{2h}_T(z_p+z^{\prime })}{{2r}_{x1}}---\left(4\right)$$r_2^{\prime }\approx r_{c2}+\frac{x^{\prime 2}+{2x}^{\prime }{x}_P}{{2r}_{c2}}+\frac{z^{\prime 2}+{2z}^{\prime }{z}_P}{{2r}_{c2}}+\frac{h_R^2-{2h}_R(z_p+z^{\prime })}{{2r}_{c2}}---\left(5\right)$其中，和分别表示反射系数等于1时，四条路径上的接收信号；Γ为复地面反射系数，与地面特性、反射波擦地角、波长和天线的极化方式有关；A是发射信号幅度，λ是连续波信号的波长，k＝2π/λ是波数，S是阴影孔径，r_c1，r_c2是目标中心点分别到发射机和接收机之间的距离，当目标以速度v平行于xoy平面运动，与基线夹角为时，目标的运动方程可以表示为公式(6)和公式(7)：步骤3、对接收信号进行相位补偿和傅里叶变换得到目标复轮廓函数，如公式(8)所示：$\stackrel{·}{H}\left(x^{\prime }\right)=\frac{\gamma }{2\mathrm{\pi v}\stackrel{·}{Q}}{e}^{-j\frac{\gamma x^{\prime 2}}{{2v}^2}}{\int }_{-T_s/2}^{T_s/2}\stackrel{·}{E}\left(t\right)\exp (-\mathrm{j\gamma }{x}^{\prime }t/v)\mathrm{dt}---\left(8\right)$其中，T_s为相参积累时间，运动补偿参数和γ可分别表示如公式(9)和(10)：步骤4、多径干涉下目标复轮廓函数可表示为：$\stackrel{·}{H}\left(x^{\prime }\right)={\stackrel{·}{H}}_1\left(x^{\prime }\right)+\Gamma {\stackrel{·}{H}}_2\left(x^{\prime }\right)+\Gamma {\stackrel{·}{H}}_3\left(x^{\prime }\right)+\Gamma {\stackrel{·}{H}}_4\left(x^{\prime }\right)---\left(11\right)$其中，分别表示每条传播路径上，目标的目标复轮廓函数，当忽略目标高度二次项引入相位时，可分别近似表示为：${\stackrel{·}{H}}_1\left(x^{\prime }\right)\approx {\stackrel{·}{Q}}_1\exp \left[\mathrm{jk}{\alpha }_{1}m\left(x^{\prime }\right)\right]·h\left(x^{\prime }\right)\sin c[k{\alpha }_{1}h\left(x^{\prime }\right)/2\pi ]---\left(12\right)$${\stackrel{·}{H}}_2\left(x^{\prime }\right)\approx {\stackrel{·}{Q}}_1\exp \left\{\mathrm{jk}\frac{{2h}_T[z_P+m\left(x^{\prime }\right)]}{r_{c1}}\right\}\exp \left[\mathrm{jk}{\alpha }_{1}m\left(x^{\prime }\right)\right]h\left(x^{\prime }\right)\sin c(k{\alpha }_{2}h\left(x^{\prime }\right)/2\pi )---\left(13\right)$${\stackrel{·}{H}}_3\left(x^{\prime }\right)\approx {\stackrel{·}{Q}}_1\exp \left\{\mathrm{jk}\frac{{2h}_T[z_P+m\left(x^{\prime }\right)]}{r_{c2}}\right\}\exp \left[\mathrm{jk}{\alpha }_{1}m\left(x^{\prime }\right)\right]h\left(x^{\prime }\right)\sin c(k{\alpha }_{3}h\left(x^{\prime }\right)/2\pi )---\left(14\right)$${\stackrel{·}{H}}_4\left(x^{\prime }\right)\approx {\stackrel{·}{Q}}_1\exp \left\{\mathrm{jk}\right[\frac{{2h}_T[z_P+m\left(x^{\prime }\right)]}{r_{c1}}+\frac{{2h}_R[z_P+m\left(x^{\prime }\right)]}{r_{c2}}\left\}\exp \right[\mathrm{jk}{\alpha }_{1}m\left(x^{\prime }\right)]h\left(x^{\prime }\right)\sin c(k{\alpha }_{4}h\left(x^{\prime }\right)/2\pi )---\left(15\right)$其中：${\stackrel{·}{Q}}_1=\exp \left[\mathrm{jk}(\frac{h_T^2-{2z}_P{h}_T}{{2r}_{c1}}+\frac{h_R^2-{2z}_P{h}_R}{{2r}_{c2}})\right]---\left(16\right)$是路径1上由发射天线和接收天线高度的一次项和二次项引入的前向散射信号相位，α₁，α₂，α₃，α₄分别如公式(17)至(20)所示，可近似认为是四条路径上目标中心点俯仰向观测角α_v，β_v之和，${\alpha }_1=\frac{z_P-h_T}{r_{c1}}+\frac{z_P-h_R}{r_{c2}}---\left(17\right)$${\alpha }_2=\frac{z_P+h_T}{r_{c1}}+\frac{z_P-h_R}{r_{c2}}---\left(18\right)$${\alpha }_3=\frac{z_P-h_T}{r_{c1}}+\frac{z_P+h_R}{r_{c2}}---\left(19\right)$${\alpha }_4=\frac{z_P+h_T}{r_{c1}}+\frac{z_P+h_R}{r_{c2}}---\left(20\right)$步骤5、在多径条件下估计目标轮廓高度差和中线，当地面为理想地面，即反射系数Γ≈-1时：$\stackrel{·}{H}\left(x^{\prime }\right)={\stackrel{·}{H}}_1\left(x^{\prime }\right){\stackrel{·}{-H}}_2\left(x^{\prime }\right)-{\stackrel{·}{H}}_3\left(x^{\prime }\right)+{\stackrel{·}{H}}_4\left(x^{\prime }\right)$$=h\left(x^{\prime }\right)\exp \left[\mathrm{jk}\left(\frac{h_T(h_T-{2z}_P)+2m\left(x^{\prime }\right)(z_P-h_T)}{r_{c1}}\right)\right]\{1-2\exp [\mathrm{jk}\frac{{2h}_T(z_P+m\left(x^{\prime }\right))}{r_{c1}}]+\exp [\mathrm{jk}\frac{{4h}_T(z_P+m\left(x^{\prime }\right))}{r_{c1}}\left]\right\}---\left(21\right)$由目标复轮廓函数的相位信息提取目标的中线高度：$m\left(x^{\prime }\right)=[\frac{r_{c1}\psi \left(x^{\prime }\right)}{k}-h_T^2]/{2z}_p---\left(22\right)$由目标中线高度及目标复轮廓函数幅度可以获得目标轮廓的高度差信息，$h\left(x^{\prime }\right)\frac{H\left(x^{\prime }\right)}{4{\sin }^2\left(\frac{\psi \left(x^{\prime }\right)r_{c1}{h}_T-k{h}_T^3+2k{z}_p^2{h}_T}{{2z}_p{r}_{c1}}\right)}---\left(23\right)$当地面为真实地面时，假设地面反射系数(其中：ρ为反射系数的幅度，为反射系数的相位)为简化复杂性，假设发射机和接收机的天线具有相同的高度，目标垂直穿越基线中点，且经路径4到接收机相位中心的前向散射信号，由于经两次地面反射幅度衰减很大，忽略不计，仅考虑经一次反射到达到接收机相位中心的多径信号，此时，可表示如下：$\stackrel{·}{H}\left(x^{\prime }\right)={\stackrel{·}{H}}_1\left(x^{\prime }\right)+\Gamma {\stackrel{·}{H}}_2\left(x^{\prime }\right)+\Gamma {\stackrel{·}{H}}_3\left(x^{\prime }\right)$的幅度可表示为：其中：由于在真实地面的情况下，目标复轮廓函数幅度受多径信号的影响相对较小，目标复轮廓函数幅度与目标高度差的相关性很高，因此，初步粗略的认为然后再由中线高度的估计值及地面反射率根据公式(25)修正目标轮廓高度差的估计值，通过迭代的方法最终可以获得精确的目标轮廓高度差和中线高度估计值。