一种电磁散射中目标瞬态感应特征图像生成方法与流程

1.本发明涉及电磁散射，特别是涉及一种电磁散射中目标瞬态感应特征图像生成方法。


背景技术：

2.随着短脉冲通信和超宽带雷达系统的广泛应用，采用时域方法来分析目标的瞬态电磁散射已经成为计算电磁学中的热点之一。了能够快速准确分析电大尺寸目标的瞬态散射，近年来发展出与各种频域高频算法相对应的时域方法。其中，时域物理光学是物理光学法的时域版本，其主要思路是通过傅里叶变换将频域物理光学公式转化为时域公式，然后利用把目标表面的时域等效电流密度代入到时域散射电场表达式中，得到散射远场时域分布。然而众多学者从瞬态电磁散射的结果即目标的最终雷达散射截面的快速计算层面，而鲜有对瞬态电磁散射过程中平面波与散射体的时域作用瞬态过程进行分析。


技术实现要素：

3.本发明的目的在于克服现有技术的不足，提供一种电磁散射中目标瞬态感应特征图像生成方法,采用数字采样的方式将目标表面的时域响应信号离散成时间与空间的二维序列，通过数字图像表征了目标在瞬态脉冲平面波激励条件下的感应场特征图像,能够实现目标在瞬态脉冲平面波下感应场特征的可视化展示。
4.本发明的目的是通过以下技术方案来实现的：s1.给定规则曲面受脉冲平面电磁波辐射的场景信息；
5.s2.将规则曲面在xoy平面进行投影，得到投影的矩形区域r
xoy
,并在矩形区域r
xoy
内进行网格生成，得到平面二维离散矩形网格；
6.s3.遍历平面二维离散矩形网格中的每一个网格，得到每一个网格的中心点坐标和在曲面上的映射点，计算平面波波前从o点到达曲面的最小和最大光程，并进行区间分割得到时间序列和幅度序列；
7.s4.对于二维离散矩形网格中的每一个网格,确定网格中心点位置上的z
ij
、j
t
(r
ij
,t)的xyz分量，分别用z
ij
、表示，并形成列表
8.s5.生成电磁散射中目标瞬态感应特征图像。
9.本发明的有益效果是：
10.本专利给出了目标在瞬态脉冲平面波激励条件下的感应场表达式，该表达式相比时谐频域表达式的最大不同是采用了时域卷积的形式，而后采用数字采样的方式将目标表面的时域响应信号离散成时间与空间的二维序列，通过数字图像表征了目标在瞬态脉冲平面波激励条件下的感应场特征图像，能够实现目标在瞬态脉冲平面波下感应场特征的可视化展示。
附图说明
11.图1为本发明的方法流程图；
12.图2为本发明的辐射场景示意图；
13.图3为规则曲面投影示意图；
14.图4为网格剖分示意图；
15.图5为平面波波前从o点到达曲面的最小和最大光程示意图；
16.图6为实施例中得到的iz示意图；
17.图7为实施例中得到的i
x
示意图；
18.图8为实施例中得到的i
ty
示意图；
19.图9为实施例中得到的i
tz
示意图。
具体实施方式
20.下面结合附图进一步详细描述本发明的技术方案，但本发明的保护范围不局限于以下所述。
21.如图1所示，一种电磁散射中目标瞬态感应特征图像生成方法，其特征在于：包括以下步骤：
22.s1.给定规则曲面受脉冲平面电磁波辐射的场景信息：如图2所示，由脉冲平面电磁波照射到曲面表面，产生电流感应，设所述的场景信息包括：
23.(1)脉冲平面电磁波的入射方向为水平和垂直极化场强分别为和脉冲信号的形式s(t)，将脉冲离散为时间序列，其采样时间间隔为δt，脉冲的采样点数为n
t
，则脉冲的总时长为n
t
δt；；
24.(2)规则曲面的表征公式为z＝f(x,y)，参数范围为：x1≤x≤x2和y1≤y≤y2；
25.(3)网格剖分密度为d，在x轴和y轴上的网格剖分长度δx＝δy＝λ/d，其中λ＝c0/f0为工作频率对应的波长,c0为自由空间波速。
26.s2.将规则曲面在xoy平面进行投影，得到投影的矩形区域r
xoy
,如图3所示；然后在矩形区域r
xoy
内进行网格生成，得到平面二维离散矩形网格，如图4所示；
27.在矩形区域r
xoy
内进行网格生成，得到平面二维离散矩形网格的过程包括：
28.根据设置的网格密度和得到的x轴网格剖分长度δx和y轴网格剖分长度δy，计算x轴剖分数nx和y轴剖分数ny：
[0029][0030][0031]
s3.遍历平面二维离散矩形网格中的每一个网格，得到每一个网格的中心点坐标和在曲面上的映射点，计算平面波波前从o点到达曲面的最小和最大光程，并进行区间分割得到时间序列和幅度序列；
[0032]
所述步骤s3包括：
[0033]
s301：遍历x轴上的(nx-1)个网格和y轴上的(ny-1)个网格，得到x轴上的第i个，y轴上的第j个网格g
ij
,g
ij
网格中心点二维坐标r
ij
(xi,yj)为：
[0034]
xi＝(i 0.5)δx,yj＝(j 0.5)δy
[0035]
s302：网格中心点二维坐标(xi,yi)代入规则曲面表征公式得到该网格中心点在曲面上的映射点r
ij
(xi,yj,z
ij
),其在xoy平面上的投影点与网格中心点重合，其中
[0036]zij
＝f(xi,yj)。
[0037]
s303.循环步骤s301和步骤s302直到遍历完所有网格为止，并得到和分别为平面波波前从o点到达曲面的最小和最大光程；o点表示为世界坐标坐标系原点；
[0038]
s304：将区间均匀分割成ns 1段：
[0039][0040]
以使得电磁波经过每段的时间长度为δt，每段的长度δl＝(t
e-tb)/ns；
[0041]
s305：将脉冲信号s(t)在[0,n
t
δt]范围内均匀分割成n
t
1段，得到信号的时间序列以及对应的幅度序列
[0042]
s4.对于二维离散矩形网格中的每一个网格,确定网格中心点位置上的z
ij
、j
t
(r
ij
,t)的xyz分量，分别用z
ij
、表示，并形成列表
[0043]
所述步骤s4包括：
[0044]
s401.遍历所有曲面上的网格中心映射点r
ij
，计算r
ij
的瞬态感应电流j
ij
(t)；
[0045]
由时域物理光学近似，得到规则曲面被脉冲平面电磁波照亮区域上上任意一点r
′
在某一个时刻t的表面感应电流j
po
(r
′
,t)为：
[0046][0047]
而对于信号为s(t)的脉冲平面电磁波，其入射到位置上的瞬态磁场分量表示为：
[0048][0049]
其中,为平面电磁波波前抵达曲面位置r
′
上所经历的延迟时间，上所经历的延迟时间，上所经历的延迟时间，上所经历的延迟时间，为自由空间波数,η＝120π为自由空间波阻抗，c0为自由空间波速；δ为单位冲击函数；为卷积符号；
[0050]
计算第ij号网格中心r
ij
的瞬态感应电流j
ij
(t)在δt时间间隔下的幅度序列的过程包括有如下步骤：
[0051]
步骤a1:计算r
ij
与o之间的光程差对应的时间延迟
[0052]
步骤a2:计算δt
ij
在区间对应的序号k
ij
＝δt
ij
/δt；
[0053]
步骤a3:计算icur的范围区间icur∈[mincur,maxcur],由用户在[mincur,maxcur]范围内给定值，记为icur0；
[0054]
步骤a3:构造长度为ns 1，值全为0的序列并使其第k
ij
个值设置为1，
[0055]
步骤a4:计算其中为序列卷积符号，得到的序列长度为nc；
[0056]
步骤a5:完成j
ij
(t)在δt时间间隔下的幅度序列的计算；
[0057]
s402:计算t＝icur0*δt时刻某个网格中心点位置上的z
ij
、j
t
(r
ij
,t)的xyz分量，分别用z
ij
、表示；
[0058]
s403：重复步骤s401～步骤s402直到完成nx
×
ny个网格的遍历为止，得到列表
[0059]
以上步骤完成了规则曲面到xoy平面的平面投影矩形网格模型的映射，得到了矩形网格模型上每个矩形网格对应的规则曲面所受瞬态平面波激励导致的t＝icur*δt时刻瞬态表面感应场ju、,jv、j
tx
、j
ty
和j
tz
。
[0060]
s5.生成电磁散射中目标瞬态感应特征图像：
[0061]
分为为列表z、j
tx
、j
ty
和j
tz
创建四个分辨率都为nx
×
ny的数字图像数据imagz,imagjtx,imagjty,imagjtz，并将列表z、j
tx
、j
ty
和j
tz
的矩阵数据按照行列的二重循环的顺序分别逐一赋值给imagz,imagjtx,imagjty,imagjtz的各个像素点，这四个数字图像采用分辨率为nx
×
ny，每个像素代表的实际物理尺寸分别是δx和δy；设置数字图像的颜色深度和显示的通道数为1；
[0062]
分别将这四个数字图像数据分别保存至四个bmp文件，取名分别为iz.bmp、itx.bmp、ity.bmp和itz.bmp。
[0063]
在本技术的实施例中，按照本技术的方法，进行高斯曲面目标表面感应瞬态电流特征图像表征进行处理，涉及到的参数如下
[0064]
f0＝4ghz，f
bw
＝0.5ghzθ
inc
＝35,高斯曲面的表征公式为a＝3，σ＝5，参数范围为：-12.5≤x≤12.5和-12.5≤y≤12.5；脉冲采用高斯脉冲，表示为：
[0065][0066]
其中，τ＝4/f
bw
,t0＝0.8τ
[0067]
得到的iz、i
tx
、i
ty
、i
tz
如图6～9所示；从图8和图9中可以清楚且直观地观察到目标曲面在瞬态平面波映射下的散射中心即为图像中清晰显示的波动中心位置点，结合图6的
高度特征图像可以看到，波动中心位置点即与图6中所示的曲面高度的峰值点基本重合，由此可见，曲面的峰值点是构成目标散射中心的一类非常重要的因素。图7可以看到没有散射中心的存在，说明散射中心具有矢量特性，不是在各个方向都存在，而是与入射和散射方向有关，例如在本例的x方向上就不存在散射中心。
[0068]
以上所述是本发明的优选实施方式，应当理解本发明并非局限于本文所披露的形式，不应该看作是对其他实施例的排除，而可用于其他组合、修改和环境，并能够在本文所述构想范围内，通过上述教导或相关领域的技术或知识进行改动。而本领域人员所进行的改动和变化不脱离本发明的精神和范围，则都应在本发明所附权利要求的保护范围内。