CN104459643A - 一种舰船尾迹雷达回波仿真方法及系统 - Google Patents

Abstract

本发明提供一种舰船尾迹雷达回波仿真方法及系统，所述方法包括建立海面几何模型和舰船尾迹几何模型；叠加所述海面几何模型和所述舰船尾迹几何模型，建立海面与舰船尾迹的一体化复合场景；根据所述一体化复合场景对其进行电磁散射计算，形成舰船尾迹雷达回波。所述系统包括：几何建模单元，用于建立海面几何模型和舰船尾迹几何模型；一体化建模单元，用于叠加所述几何建模单元中的海面几何模型和舰船尾迹几何模型，建立海面与舰船尾迹的一体化复合场景；电磁计算单元，用于根据所述一体化建模单元提供的一体化复合场景对其进行电磁散射计算，形成舰船尾迹雷达回波。本发明适用于目标识别与仿真技术领域。

Description

一种舰船尾迹雷达回波仿真方法及系统

技术领域

本发明涉及目标识别与仿真技术领域，特别是指一种舰船尾迹雷达回波仿真方法及系统。

背景技术

在星载和机载雷达对海面目标探测和电磁成像中，运动舰船的尾迹表现出强散射特性，不同舰船在不同运动状态时所具有舰船尾迹的电磁散射特征不同，该电磁散射特征是复杂海环境中目标精确预估、识别和信号分离技术的关键所在，其应用涉及军用遥感、雷达成像、制导和预警技术等众多领域。

近十几年来，欧洲遥感卫星(European Remote Sensing Satellite，ERS)-1和ERS-2获得了许多包含舰船尾迹的合成孔径雷达(Synthetic ApertureRadar，SAR)图像。多数国家都有关于仿真尾迹SAR回波和SAR图像研究，目的是通过不同的观测参数对舰船尾迹SAR图像进行仿真，来确定对不同类型舰船尾迹的最佳SAR观测条件。如Tunaley等对舰船湍流尾迹和开尔文(Kelvin)尾迹的SAR图像进行了仿真，但是没给出具体的模型和方法；Oumansour分别仿真了L和X波段雷达参数下的Kelvin尾迹SAR图像，对比分析了不同雷达频段对Kelvin尾迹的观测效果。用电磁模拟仿真的方法来研究SAR对舰船尾迹的成像，是得到对舰船尾迹最佳观测参数的一种重要手段。然而，对舰船尾迹雷达实测数据稀缺且理论模型过于简化，严重制约舰船尾迹检测和识别。

发明内容

本发明要解决的技术问题是提供一种舰船尾迹雷达回波仿真方法及系统，以解决现有技术所存在的舰船尾迹雷达实测数据稀缺且理论模型过于简化的问题。

为解决上述技术问题，本发明实施例提供一种舰船尾迹雷达回波仿真方法，包括：

建立海面几何模型和舰船尾迹几何模型；

叠加所述海面几何模型和所述舰船尾迹几何模型，建立海面与舰船尾迹的一体化复合场景；

根据所述一体化复合场景对其进行电磁散射计算，形成舰船尾迹雷达回波。

可选地，所述建立海面几何模型和舰船尾迹几何模型包括：

根据海谱模型的线性过滤法来实现海面几何模型；

根据舰船尺寸和运动参数来建立具体的舰船尾迹几何模型。

可选地，所述叠加所述海面几何模型和所述舰船尾迹几何模型，建立海面与舰船尾迹的一体化复合场景包括：

叠加所述海面几何模型和所述舰船尾迹几何模型；

利用三角面元建立海面与舰船尾迹的一体化复合场景。

可选地，所述根据所述一体化复合场景对其进行电磁散射计算，形成舰船尾迹雷达回波包括：

将所述的一体化复合场景大尺度轮廓划分为一个个小平面，同时将各个小平面按照海面的确定模拟样本的斜率特征来倾斜；

利用斜率叠加双尺度面元模型和基于海谱分布的单元场散射模型分别对所述小平面的面元和面元内的毛细波结构进行电磁散射计算，确定各面元的散射场；

将各面元贡献的散射场进行相关叠加形成总的散射贡献，最终形成舰船尾迹雷达回波。

可选地，所述舰船尾迹几何模型包括：开尔文尾迹、湍流尾迹、内波尾迹和潜艇尾迹。

本发明实施例所述的舰船尾迹雷达回波仿真方法，通过建立海面几何模型和舰船尾迹几何模型，并叠加所述海面几何模型和所述舰船尾迹几何模型，从而建立海面与舰船尾迹的一体化复合场景，并根据所述一体化复合场景对其进行电磁散射计算，形成舰船尾迹雷达回波。如此，通过建立海面和舰船尾迹的一体化复合场景能够为舰船尾迹的检测和识别提供模型参考，对所述一体化复合场景进行电磁散射计算形成的舰船尾迹雷达回波，能够为舰船尾迹的检测和识别提供完整的数据信息，同时该方法还能突破天气、目标和传感器等方面的局限。

另一方面，本发明实施例提供一种舰船尾迹雷达回波仿真系统，包括：

几何建模单元：用于建立海面几何模型和舰船尾迹几何模型；

一体化建模单元：用于叠加所述几何建模单元中的海面几何模型和舰船尾迹几何模型，建立海面与舰船尾迹的一体化复合场景；

电磁计算单元：用于根据所述一体化建模单元提供的一体化复合场景对其进行电磁散射计算，形成舰船尾迹雷达回波。

可选地，所述几何建模单元包括：

海面建模模块：用于根据海谱模型的线性过滤法来建立海面几何模型；

舰船建模模块：用于根据舰船尺寸和运动参数来建立具体的舰船尾迹几何模型。

可选地，所述一体化建模单元包括：

第一叠加模块：用于叠加所述海面几何模型和所述舰船尾迹几何模型；

一体化建模模块：用于利用三角面元建立海面与舰船尾迹的一体化复合场景。

可选地，所述电磁计算单元包括：

划分模块：用于将所述的一体化复合场景大尺度轮廓划分为一个个小平面，同时将各个小平面按照海面的确定模拟样本的斜率特征来倾斜；

电磁计算模块：用于利用斜率叠加双尺度面元模型和基于海谱分布的单元场散射模型分别对所述小平面的面元和面元内的毛细波结构进行电磁散射计算，确定各面元的散射场；

第二叠加模块：用于将各面元贡献的散射场进行相关叠加形成总的散射贡献，最终形成舰船尾迹雷达回波。

可选地，所述舰船尾迹几何模型包括：开尔文尾迹、湍流尾迹、内波尾迹和潜艇尾迹。

本发明实施例所述的舰船尾迹雷达回波仿真系统，通过几何建模单元建立海面几何模型和舰船尾迹几何模型，并通过一体化建模单元叠加所述几何建模单元中的海面几何模型和舰船尾迹几何模型，从而建立海面与舰船尾迹的一体化复合场景，并通过电磁计算单元根据所述一体化建模单元提供的一体化复合场景对其进行电磁散射计算，形成舰船尾迹雷达回波。如此，通过建立海面和舰船尾迹的一体化复合场景能够为舰船尾迹的检测和识别提供模型参考，对所述一体化复合场景进行电磁散射计算形成的舰船尾迹雷达回波，能够为舰船尾迹的检测和识别提供完整的数据信息，同时该系统还能突破天气、目标和传感器等方面的局限。

本发明的上述技术方案的有益效果如下：

上述方案中，如此，通过建立海面和舰船尾迹的一体化复合场景能够为舰船尾迹的检测和识别提供模型参考，对所述一体化复合场景进行电磁散射计算形成的舰船尾迹雷达回波，能够为舰船尾迹的检测和识别提供完整的数据信息，同时该方法或者系统还能突破天气、目标和传感器等方面的局限。

附图说明

图1为本发明实施例一提供的舰船尾迹雷达回波仿真方法流程图；

图2为全局直角坐标系与本地坐标系示意图；

图3为本发明实施例提供的海面和舰船尾迹的一体化复合场景示意图；

图4为本发明实施例提供的HH极化下的面元后向散射系数分布图；

图5为本发明实施例提供的HH极化下的场景后向散射系数分布图；

图6为本发明实施例提供的VV极化下的场景后向散射系数分布图；

图7为本发明实施例提供的VV极化下的SAR图像；

图8为本发明实施例提供的HH极化下的SAR图像；

图9为本发明实施例二提供的舰船尾迹雷达回波仿真系统的结构示意图；

图10为图9中几何建模单元单元101的结构示意图；

图11为图9中一体化建模单元102的结构示意图；

图12为图9中电磁计算单元103的结构示意图。

具体实施方式

为使本发明要解决的技术问题、技术方案和优点更加清楚，下面将结合附图及具体实施例进行详细描述。

本发明针对现有的存在的舰船尾迹雷达实测数据稀缺且理论模型过于简化的问题，提供一种舰船尾迹雷达回波仿真方法及系统。

实施例一

参看图1所示，本发明实施例提供的舰船尾迹雷达回波仿真方法的具体实施方式，该方法包括以下步骤：

S101：建立海面几何模型和舰船尾迹几何模型。

S102：叠加所述海面几何模型和所述舰船尾迹几何模型，建立海面与舰船尾迹的一体化复合场景。

S103：根据所述一体化复合场景对其进行电磁散射计算，形成舰船尾迹雷达回波。

以下对图1中各步骤的具体实施方法进行详细说明。

S101中，所述建立海面几何模型和舰船尾迹几何模型的步骤包括：

根据海谱模型的线性过滤法来实现海面几何模型，同时根据舰船尺寸和运动参数来建立具体的舰船尾迹几何模型，其中，所述舰船尾迹几何模型包括：开尔文尾迹、湍流尾迹、内波尾迹和潜艇尾迹。

S102中，叠加S101中所述海面几何模型和所述舰船尾迹几何模型，建立海面与舰船尾迹的一体化复合场景的步骤包括：

叠加所述海面几何模型和所述舰船尾迹几何模型，并利用三角面元建立海面与舰船尾迹的一体化复合场景。

如图3所示为海面和舰船尾迹的一体化复合场景示意图，该基于三角面元的一体化复合场景能够逼真模拟含有舰船尾迹的海面，通过设置海面几何模型和舰船尾迹几何模型的参数可以模拟不同海况和不同目标的舰船尾迹。

如图4所示，展示的是与图3所示的一体化复合场景对应的HH极化下的面元后向散射系数分布图，可见二者具有很高的相关性。由此可见，上述实施例通过建立海面和舰船尾迹的一体化复合场景能为舰船尾迹的检测和识别提供模型参考。

其中，所述建立海面与舰船尾迹的一体化复合场景中舰船尾迹以Kelvin尾迹为例，Kelvin尾迹是由船体对水面的扰动直接产生的一种表面波尾迹，通常由扩散波(Divergent Waves)和横断波(Transverse Waves)组成。Kelvin尾迹波属于长重力波，尾迹表面起伏最大的区域位于尾迹角为16°至19.5°之间，在这个区域，扩散波和横断波发生相互干涉形成尖头波(Cusp Waves)。

假定舰船以速度U_s沿着x方向航行，产生的自由表面波沿着与x轴夹角为θ的不同角度传播，相应的Kelvin尾迹波高可以用这些平面波的叠加表示为式(1)：

$h_k(x,y)=\mathrm{Re}{\∫}_{-\mathrm{\π}/2}^{\mathrm{\π}/2}A\left(\mathrm{\θ}\right)\exp \left[k_0{\sec }^2\mathrm{\θ}(i\cos \mathrm{\θx}+i\sin \mathrm{\θy})\right]\mathrm{d\θ}---\left(1\right)$

式(1)中，k₀sec²θ是沿着与x轴夹角θ传播的波分量的波数，A(θ)是复数，表示船的特征参数。对于有限水深，对于无限水深，即(H→∞)， $k_0=\frac{g}{U_s^2}.$

实际情况下，舰船尾迹由许多波长成分的海浪组成，其相应的最大波长为λ_m＝2π/k₀，这就为在无限水深和有限水深的海面上叠加相应的尾迹提供了理论基础。根据船理论，舰船周围的水流由船的中心平面分布的源产生，水流源强度σ(x，y，z)与船体的局部斜度成正比，A(θ)表示舰船参数特征的自由波谱，由式(2)计算：

$A\left(\mathrm{\θ}\right)=\frac{4k_0}{U_s}\left({\sec }^3\mathrm{\θ}\right)H(k_0{\sec }^2\mathrm{\θ},\mathrm{\θ})---\left(2\right)$

$H(K,\mathrm{\θ})=\underset{S_H}{\∫}\∫\mathrm{\σ}(x,y,z)\·\exp \left(K(\mathrm{ix}\cos \mathrm{\θ}+\mathrm{iy}\sin \mathrm{\θ}+z)\right)\mathrm{dxdz}---\left(3\right)$

式(3)中，H(K，θ)是Kochin函数，S_H是舰船表面，水流源强度σ(x，y，z)由式(4)计算：

$\mathrm{\σ}(x,y,z)=-\frac{{2U}_s}{4\mathrm{\π}}\frac{\∂}{\∂x}f(x,z)---\left(4\right)$

式(4)中，f是舰船的船体方程，若船体形状为简单的抛物型威格利(Wigley)船形式，则有

$f(x,z)=\left\{\begin{array}{cc}b(1-x^2/l^2)& (-d<z<0,-l<x<l)\\ 0& (z<-d)\end{array}\right.---\left(5\right)$

式(5)中，x为偏移位置，z为吃水深度，b为半船宽，l为半船长，d为侧壁吃水深度。

S103中，根据S102提供的用三角面元描述的一体化复合场景对其进行电磁散射计算，形成舰船尾迹雷达回波。

所述根据S102提供的一体化复合场景对其进行电磁散射计算，形成舰船尾迹雷达回波的步骤包括：

将S102提供的一体化复合场景大尺度轮廓划分为一个个小平面，同时将各个小平面按照海面的确定模拟样本的斜率特征来倾斜，并利用斜率叠加双尺度面元模型和基于海谱分布的单元场散射模型分别对各个平面的面元和面元内的毛细波结构进行电磁散射计算，确定各面元的散射场，因为每个面元均可按照海面的确定模拟样本的斜率特征来倾斜，同时若忽略各个面元之间的相互作用以及多次散射，则将各面元贡献的散射场进行相关叠加即可得总的散射贡献，最终形成舰船尾迹雷达回波，其中所述毛细波结构为小尺度部分。

首先介绍一下面元的后向散射系数。

按照布拉格(Bragg)谐振假说，沿雷达视线正负方向传播的Bragg波成分均会对雷达接收器产生贡献，任意倾斜微粗糙面元的后向散射系数为式(6)：

${\mathrm{\&sigma;}}_{\mathrm{PQ}}^{\mathrm{TSPM}}({\hat{k}}_i,{\hat{k}}_s)={\mathrm{\&pi;k}}^4{|\mathrm{\&epsiv;}-1|}^2{\left|{\tilde{F}}_{\mathrm{PQ}}\right|}^2{S}_{\mathrm{capi}}\left(q_l\right)---\left(6\right)$

式(6)中，下标P，Q＝H，V分别表示入射波和散射波在全局坐标系下的极化方式。为电磁波入射散射单位矢量，表示全局坐标系下的极化因子

$\left[\begin{array}{cc}{\tilde{F}}_{\mathrm{VV}}& {\tilde{F}}_{\mathrm{VH}}\\ {\tilde{F}}_{\mathrm{HV}}& {\tilde{F}}_{\mathrm{HH}}\end{array}\right]=\left[\begin{array}{cc}{\tilde{V}}_s\&CenterDot;{\hat{v}}_s& {\hat{H}}_s\&CenterDot;{\hat{v}}_s\\ {\hat{V}}_s\&CenterDot;{\hat{h}}_s& {\hat{H}}_s\&CenterDot;{\hat{h}}_s\end{array}\right]\left[\begin{array}{cc}{F}_{\mathrm{vv}}& F_{\mathrm{vh}}\\ F_{\mathrm{hv}}& F_{\mathrm{hh}}\end{array}\right]\left[\begin{array}{cc}{\hat{V}}_i\&CenterDot;{\hat{v}}_i& {\hat{H}}_i\&CenterDot;{\hat{v}}_i\\ {\hat{V}}_i\&CenterDot;{\hat{h}}_i& {\hat{H}}_i\&CenterDot;{\hat{h}}_i\end{array}\right]---\left(7\right)$

式(7)中，和分别表示如图2所示的全局直角坐标系和本地坐标系对应的单位极化矢量，

$\begin{array}{c}{F}_{\mathrm{vv}}=\frac{1}{\mathrm{\&epsiv;}}[1+R_v\left({\mathrm{\&theta;}}_i^l\right)][1+R_v\left({\mathrm{\&theta;}}_s^l\right)]\sin {\mathrm{\&theta;}}_i^l\sin {\mathrm{\&theta;}}_s^l\\ -[1-R_v\left({\mathrm{\&theta;}}_i^l\right)][1-R_v\left({\mathrm{\&theta;}}_s^l\right)]\cos {\mathrm{\&theta;}}_i^l\cos {\mathrm{\&theta;}}_s^l{\cos \mathrm{\&phi;}}_s^l\end{array}---\left(8\right)$

$\begin{array}{c}{F}_{\mathrm{vh}}=[1-R_v\left({\mathrm{\&theta;}}_i^l\right)][1+R_h\left({\mathrm{\&theta;}}_s^l\right)]\cos {\mathrm{\&theta;}}_i^l{\sin \mathrm{\&phi;}}_s^l---\left(9\right)\end{array}$

$\begin{array}{c}{F}_{\mathrm{hv}}=[1+R_h\left({\mathrm{\&theta;}}_i^l\right)][1-R_v\left({\mathrm{\&theta;}}_s^l\right)]\cos {\mathrm{\&theta;}}_i^l{\sin \mathrm{\&phi;}}_s^l---\left(10\right)\end{array}$

$\begin{array}{c}{F}_{\mathrm{hh}}=[1+R_h\left({\mathrm{\&theta;}}_i^l\right)][1+R_h\left({\mathrm{\&theta;}}_s^l\right)]{\sin \mathrm{\&phi;}}_s^l---\left(11\right)\end{array}$

式(8)-(11)中，为本地双站角，R_v和R_h分别为垂直(VV)和水平(HH)极化下的菲涅尔反射系数。

式(6)中，S_capi(q_l)表示为式(12)，表示对雷达有贡献的Bragg毛细波(小尺度)部分：

$S_{\mathrm{capi}}\left(q_l\right)=\frac{1}{2}[S_E^{\mathrm{capi}}\left({-q}_l\right)+S_E^{\mathrm{capi}}\left(q_l\right)]---\left(12\right)$

考虑到要用毛细波谱来描述面元的小尺度结构，引入截断波数k_cut划分大小尺度的波谱为式(13)：

$S_E^{\mathrm{capi}}\left(k\right)=\left\{\begin{array}{cc}{S}_E\left(k\right),& \left|k\right|\&GreaterEqual;k_{\mathrm{cut}}\\ 0,& \left|k\right|<k_{\mathrm{out}}\end{array}\right.---\left(13\right)$

为了使得提高面元的后向散射系数的精确度，通过引入基尔霍夫近似模型修正近垂直入射区倾斜调制的微扰系数，得到修正后的面元的后向散射系数为式(14)：

${\mathrm{\&sigma;}}_{\mathrm{PQ}}^{\mathrm{facet}}\left({\hat{k}}_i{,\hat{k}}_s\right)={\mathrm{\&sigma;}}_{\mathrm{PQ}}^{\mathrm{KAM}}({\hat{k}}_i,{\hat{k}}_s)+{\mathrm{\&sigma;}}_{\mathrm{PQ}}^{\mathrm{TSPM}}({\hat{k}}_i,{\hat{k}}_s)---\left(14\right)$

式(14)中，微扰系数由式(15)计算：

${\mathrm{\&sigma;}}_{\mathrm{PQ}}^{\mathrm{KAM}}({\hat{k}}_i,{\hat{k}}_s)=\frac{{\mathrm{\&pi;k}}^2{\left|q\right|}^2}{q_z^4}{\left|{\tilde{F}}_{\mathrm{PQ}}^{\mathrm{KAM}}\right|}^{2}P(z_x^{\tan },z_y^{\tan })---\left(15\right)$

式(15)中，k为电磁波波数，和为镜像点切平面的斜率，KAM的极化因子可表示成

${\tilde{F}}_{\mathrm{vv}}^{\mathrm{KAM}}=M_0[R_v\left({\mathrm{\&theta;}}_i^l\right)({\hat{V}}_s\&CenterDot;{\hat{k}}_i)({\hat{V}}_i\&CenterDot;{\hat{k}}_s)+R_h\left({\mathrm{\&theta;}}_i^l\right)({\hat{H}}_s\&CenterDot;{\hat{k}}_i)({\hat{H}}_i\&CenterDot;{\hat{k}}_s)]---\left(16\right)$

${\tilde{F}}_{\mathrm{vh}}^{\mathrm{KAM}}=M_0[R_v\left({\mathrm{\&theta;}}_i^l\right)({\hat{V}}_s\&CenterDot;{\hat{k}}_i)({\hat{H}}_i\&CenterDot;{\hat{k}}_s)+R_h\left({\mathrm{\&theta;}}_i^l\right)({\hat{H}}_s\&CenterDot;{\hat{k}}_i)({\hat{V}}_i\&CenterDot;{\hat{k}}_s)]---\left(17\right)$

${\tilde{F}}_{\mathrm{hv}}^{\mathrm{KAM}}=M_0[R_v\left({\mathrm{\&theta;}}_i^l\right)({\hat{H}}_s\&CenterDot;{\hat{k}}_i)({\hat{V}}_i\&CenterDot;{\hat{k}}_s)-R_h\left({\mathrm{\&theta;}}_i^l\right)({\hat{V}}_s\&CenterDot;{\hat{k}}_i)({\hat{H}}_i\&CenterDot;{\hat{k}}_s)]---\left(18\right)$

${\tilde{F}}_{\mathrm{vv}}^{\mathrm{KAM}}=M_0[R_v\left({\mathrm{\&theta;}}_i^l\right)({\hat{H}}_s\&CenterDot;{\hat{k}}_i)({\hat{V}}_i\&CenterDot;{\hat{k}}_s)+R_h\left({\mathrm{\&theta;}}_i^l\right)({\hat{V}}_s\&CenterDot;{\hat{k}}_i)({\hat{V}}_i\&CenterDot;{\hat{k}}_s)]---\left(19\right)$

式(16)-(19)中， $M_0=\left|q\right|\left|q_z\right|/\{{({\hat{H}}_s\&CenterDot;{\hat{k}}_i)}^2+{({\hat{V}}_s\&CenterDot;{\hat{k}}_i)}^2\left]{\mathrm{kq}}_z\right\}.$

最后，通过修正后的面元的后向散射系数得到面元的散射场，再将各面元贡献的散射场进行相关叠加形成总的散射贡献，最终形成舰船尾迹雷达回波。

如图5和图6所示，分别给出了HH极化和VV极化下的不同海况的Kelvin尾迹的一体化复合场景的后向散射系数分布图，由图像结果可以明确反映出不同海况的电磁散射系数不同。图7和8分别展示了HH极化和VV极化下相同海况、相同舰船尾迹的SAR成像结果，该成像结果可以正确反映HH和VV的极化差异，即VV极化下的幅度量级较HH极化要高，由此可见，上述实施例通过电磁散射计算形成的雷达回波仿真能为舰船尾迹的检测和识别提供完整的数据信息。

本发明实施例所述的舰船尾迹雷达回波仿真方法，通过建立海面几何模型和舰船尾迹几何模型，并叠加所述海面几何模型和所述舰船尾迹几何模型，从而建立海面与舰船尾迹的一体化复合场景，并根据所述一体化复合场景对其进行电磁散射计算，形成舰船尾迹雷达回波。如此，通过建立海面和舰船尾迹的一体化复合场景能够为舰船尾迹的检测和识别提供模型参考，对所述一体化复合场景进行电磁散射计算形成的舰船尾迹雷达回波，能够为舰船尾迹的检测和识别提供完整的数据信息，同时该方法还能突破天气、目标和传感器等方面的局限。

实施例二

如图9所示为本发明实施例二提供的舰船尾迹雷达回波仿真系统的具体实施方式，该系统包括：

几何建模单元101：用于建立海面几何模型和舰船尾迹几何模型；

一体化建模单元102：用于叠加所述几何建模单元101中的海面几何模型和舰船尾迹几何模型，建立海面与舰船尾迹的一体化复合场景；

电磁计算单元103：用于根据所述一体化建模单元102提供的一体化复合场景对其进行电磁散射计算，形成舰船尾迹雷达回波。

参看图10所示，在前述舰船尾迹雷达回波仿真系统的具体实施方式中，可选地，所述几何建模单元101包括：

海面建模模块1011：用于根据海谱模型的线性过滤法来建立海面几何模型；

舰船建模模块1012：用于根据舰船尺寸和运动参数来建立具体的舰船尾迹几何模型，其中，所述舰船尾迹几何模型包括：开尔文尾迹、湍流尾迹、内波尾迹和潜艇尾迹。

参看图11所示，在前述舰船尾迹雷达回波仿真系统的具体实施方式中，可选地，所述一体化建模单元102包括：

第一叠加模块1021：用于叠加所述海面几何模型和所述舰船尾迹几何模型；

一体化建模模块1022：用于利用三角面元建立海面与舰船尾迹的一体化复合场景。

参看图12所示，在前述舰船尾迹雷达回波仿真系统的具体实施方式中，可选地，所述电磁计算单元103包括：

划分模块1031：用于将所述的一体化复合场景大尺度轮廓划分为一个个小平面，同时将各个小平面按照海面的确定模拟样本的斜率特征来倾斜；

电磁计算模块1032：用于利用斜率叠加双尺度面元模型和基于海谱分布的单元场散射模型分别对所述小平面的面元和面元内的毛细波结构进行电磁散射计算，确定各面元的散射场；

第二叠加模块1033：用于将各面元贡献的散射场进行相关叠加形成总的散射贡献，最终形成舰船尾迹雷达回波。

本发明实施例所述的舰船尾迹雷达回波仿真系统，通过几何建模单元101建立海面几何模型和舰船尾迹几何模型，并通过一体化建模单元102叠加所述几何建模单元101中的海面几何模型和舰船尾迹几何模型，从而建立海面与舰船尾迹的一体化复合场景，并通过电磁计算单元103根据所述一体化建模单元102提供的一体化复合场景对其进行电磁散射计算，形成舰船尾迹雷达回波。如此，通过建立海面和舰船尾迹的一体化复合场景能够为舰船尾迹的检测和识别提供模型参考，对所述一体化复合场景进行电磁散射计算形成的舰船尾迹雷达回波，能够为舰船尾迹的检测和识别提供完整的数据信息，同时该系统还能突破天气、目标和传感器等方面的局限。

综上所述，本发明实施例通过建立海面和舰船尾迹的一体化复合场景能够为舰船尾迹的检测和识别提供模型参考，对所述一体化复合场景进行电磁散射计算形成的舰船尾迹雷达回波，能够为舰船尾迹的检测和识别提供完整的数据信息，同时该方法或者系统还能突破天气、目标和传感器等方面的局限。

以上所述是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明所述原理的前提下，还可以作出若干改进和润饰，这些改进和润饰也应视为本发明的保护范围。

Claims (10)

1.一种舰船尾迹雷达回波仿真方法，其特征在于，包括：

建立海面几何模型和舰船尾迹几何模型；

叠加所述海面几何模型和所述舰船尾迹几何模型，建立海面与舰船尾迹的一体化复合场景；

根据所述一体化复合场景对其进行电磁散射计算，形成舰船尾迹雷达回波。

2.根据权利要求1所述的舰船尾迹雷达回波仿真方法，其特征在于，所述建立海面几何模型和舰船尾迹几何模型包括：

根据海谱模型的线性过滤法来实现海面几何模型；

根据舰船尺寸和运动参数来建立具体的舰船尾迹几何模型。

3.根据权利要求1所述的舰船尾迹雷达回波仿真方法，其特征在于，所述叠加所述海面几何模型和所述舰船尾迹几何模型，建立海面与舰船尾迹的一体化复合场景包括：

叠加所述海面几何模型和所述舰船尾迹几何模型；

利用三角面元建立海面与舰船尾迹的一体化复合场景。

4.根据权利要求1所述的舰船尾迹雷达回波仿真方法，其特征在于，所述根据所述一体化复合场景对其进行电磁散射计算，形成舰船尾迹雷达回波包括：

将所述的一体化复合场景大尺度轮廓划分为一个个小平面，同时将各个小平面按照海面的确定模拟样本的斜率特征来倾斜；

利用斜率叠加双尺度面元模型和基于海谱分布的单元场散射模型分别对所述小平面的面元和面元内的毛细波结构进行电磁散射计算，确定各面元的散射场；

将各面元贡献的散射场进行相关叠加形成总的散射贡献，最终形成舰船尾迹雷达回波。

5.根据权利要求1所述的舰船尾迹雷达回波仿真方法，其特征在于，所述舰船尾迹几何模型包括：开尔文尾迹、湍流尾迹、内波尾迹和潜艇尾迹。

6.一种舰船尾迹雷达回波仿真系统，其特征在于，包括：

几何建模单元：用于建立海面几何模型和舰船尾迹几何模型；

一体化建模单元：用于叠加所述几何建模单元中的海面几何模型和舰船尾迹几何模型，建立海面与舰船尾迹的一体化复合场景；

电磁计算单元：用于根据所述一体化建模单元提供的一体化复合场景对其进行电磁散射计算，形成舰船尾迹雷达回波。

7.根据权利要求6所述的舰船尾迹雷达回波仿真系统，其特征在于，所述几何建模单元包括：

海面建模模块：用于根据海谱模型的线性过滤法来建立海面几何模型；

舰船建模模块：用于根据舰船尺寸和运动参数来建立具体的舰船尾迹几何模型。

8.根据权利要求6所述的舰船尾迹雷达回波仿真系统，其特征在于，所述一体化建模单元包括：

第一叠加模块：用于叠加所述海面几何模型和所述舰船尾迹几何模型；

一体化建模模块：用于利用三角面元建立海面与舰船尾迹的一体化复合场景。

9.根据权利要求6所述的舰船尾迹雷达回波仿真系统，其特征在于，所述电磁计算单元包括：

划分模块：用于将所述的一体化复合场景大尺度轮廓划分为一个个小平面，同时将各个小平面按照海面的确定模拟样本的斜率特征来倾斜；

电磁计算模块：用于利用斜率叠加双尺度面元模型和基于海谱分布的单元场散射模型分别对所述小平面的面元和面元内的毛细波结构进行电磁散射计算，确定各面元的散射场；

第二叠加模块：用于将各面元贡献的散射场进行相关叠加形成总的散射贡献，最终形成舰船尾迹雷达回波。

10.根据权利要求6所述的舰船尾迹雷达回波仿真系统，其特征在于，所述舰船尾迹几何模型包括：开尔文尾迹、湍流尾迹、内波尾迹和潜艇尾迹。