一种时域电磁散射分析的有限冲激偶序列表征方法技术

本发明专利技术公开了一种时域电磁散射分析的有限冲激偶序列表征方法，包括以下步骤：S1.给定已知频域平面电磁波，并进行时域表征；S2.对散射体表面进行离散化，并确定在平面电磁波照射下散射体表面单元感应电流；S3.确定散射体表面单元二次辐射远场，并进行时间离散与辐射场序列化；S4.在任意平面脉冲波形下的确定时域散射响应信号离散序列，即有限冲激偶序列。本发明专利技术利用时域物理光学得到天线安装于复杂几何平台表面形成的时域散射响应特性，该响应由表示曲面几何散射特性的有限冲激函数形式，再利用卷积理论能够得到在任意脉冲平面波激励条件下的散射响应曲线。条件下的散射响应曲线。条件下的散射响应曲线。

【技术实现步骤摘要】
一种时域电磁散射分析的有限冲激偶序列表征方法


[0001]本专利技术涉及电磁散射分析，特别是涉及一种时域电磁散射分析的有限冲激偶序列表征方法。

技术介绍

[0002]随着超宽带天线技术的日益发展，天线工作的频带越来越宽，天线端口的激励从点频时谐信号逐渐过度到时域脉冲信号，而且随着短脉冲通信和超宽带雷达系统的广泛应用，脉冲信号的脉冲宽度越来越窄，天线工作的有效频带也更宽，常规的基于频域的点频分析方法需要大量的频率采样点以得到宽频谱响应，因此点频分析并不适合于超宽带和短脉冲平面波的电磁散射分析，因此目标时域电磁散射的计算分析成为计算电磁学中广为关注的热点之一。

技术实现思路

[0003]本专利技术的目的在于克服现有技术的不足，提供一种时域电磁散射分析的有限冲激偶序列表征方法，根据时域散射响应特性，该响应由表示曲面几何散射特性的有限冲激函数形式，再利用卷积理论能够得到在任意脉冲平面波激励条件下的散射响应曲线，特别合适分析宽频带的散射场分析。
[0004]本专利技术的目的是通过以下技术方案来实现的：一种时域电磁散射分析的有限冲激偶序列表征方法，包括以下步骤：
[0005]S1.给定已知频域平面电磁波，并进行时域表征；
[0006]S2.对散射体表面进行离散化，并确定在平面电磁波照射下散射体表面单元感应电流；
[0007]S3.确定散射体表面单元二次辐射远场，并进行时间离散与辐射场序列化；
[0008]S4.在任意平面脉冲波形下的确定时域散射响应信号离散序列，即有限冲激偶序列。
[0009]本专利技术的有益效果是：本专利可以根据任意脉冲信号S的具体时间序列形式，利用卷积公式得到该任意脉冲信号S的散射场时域表达式，并根据时频分析得到S的散射场的频域响应，因此该专利方法特别合适分析宽频带的散射场分析，通过一次脉冲激励得到宽频带散射响应结果。
附图说明
[0010]图1为本专利技术的方法流程图；
[0011]图2为平面电磁波示意图；
[0012]图3为散射体表面离散化示意图；
[0013]图4为每个单元的表面感应电流示意；
[0014]图5为单元j的二次辐射示意图；
[0015]图6为实施例中网格曲面模型示意图；
[0016]图7为实施例中总散射场时间序列随时间的变化示意图；
[0017]图8为实施例中s(t)波形示意图；
[0018]图9为实施例中的波形示意图；
[0019]图10为实施例中得到的时域散射响应信号离散序列Sca示意图。
具体实施方式
[0020]下面结合附图进一步详细描述本专利技术的技术方案，但本专利技术的保护范围不局限于以下所述。
[0021]如图1所示，一种时域电磁散射分析的有限冲激偶序列表征方法，包括以下步骤：
[0022]S1.给定已知频域平面电磁波，如图2所示，并进行时域表征；
[0023]给定已知频域平面电磁波，表示为：
[0024][0025][0026]其中，ω＝2πf，为工作角频率，f为工作频率，单位为Hz，E
θ
为平面电磁波电场的垂直极化分量，为平面电磁波电场的水平极化分量,θ和为平面电磁波的入射角，为平面电磁波的入射方向单位矢量，为平面电磁波的入射方向单位矢量，和为电磁波的入射方向上的垂直和水平极化单位矢量，为复数量，η＝120π为自由空间电磁波波阻抗，E为电场，H为磁场，为场点坐标矢量；k0＝Jω/c0为自由空间波数；c0为自由空间电磁波速度；
[0027]对(1)和(2)进行傅里叶反变换，得到时域平面电磁波的表达式：
[0028][0029][0030]其中，为时域卷积运算符，δ为狄拉克函数。
[0031]S2.对散射体表面进行离散化，并确定在平面电磁波照射下散射体表面单元感应电流；
[0032]S201.对散射体表面进行离散化：
[0033]按照不高于0.2倍最短工作波长λ
m
＝c0/f
u
的网格长度将散射体离散成N个三角形网格单元，其中，f
u
为最高工作频率，有M个点坐标形成的列表[P
i
]和N个点连接列表[Elem
j
]，其中每个单元的点连接Elem
j
包含三个顶点的序号i0＝Elem
j
.p0，i1＝Elem
j
.p1和i2＝Elem
j
.p2，这三个顶点序号从列表[P
i
]中得到该单元的三个顶点坐标P[i0]、P[i1]和P[i2]，P[i]的点坐标可以表示为(P[i].x,P[i].y,P[i].z)；其中，i∈[0,1,2,
…
,M],j∈[0,1,2,
…
,N]，i为点序号，j为单元的点连接序号，如图3所示；
[0034]S202.确定在平面电磁波照射下散射体表面单元感应电流:
[0035]由物理光学近似，散射体表面被平面电磁波照亮区域上任意一点r在某一个时刻t的感应电流J
PO
(r,t)为：
[0036][0037]其中，为散射体表面在r位置上的单位法向矢量；
×
为矢量叉乘符号；
[0038]而对于信号为δ(t)的脉冲平面电磁波，其入射到位置上的瞬态磁场分量表示为：
[0039][0040][0041][0042]计算得到的d0为脉冲平面电磁波到达r场点位置的迟滞时间。
[0043]从N个单元中取出第j个单元，该三角形单元的点连接中的三个点坐标分别为和则有单元j中心点坐标为：由于三角形边长不高于0.2λ
m
，因此将三角形单元内的场当成均匀场来进行计算，用单元中心点的感应场作为整个三角形单元内的感应场，则有公式(7)和公式(8)计算第j个单元的表面感应电流：
[0044][0045]其中，为的坐标，为表面在的指向散射体外侧的单位法向矢量，为入射平面波在第j个单元中心位置的感应场的延迟时间,
·
为矢量点乘符号。每个单元的表面感应电流如图4所示
[0046]S3.确定散射体表面单元二次辐射远场，并进行时间离散与辐射场序列化；
[0047]S301.计算指定散射角方向上半径为R
s
的监测球面上二次辐射远场，为散射平面电磁波电场的垂直极化分量，为散射平面电磁波电场的水平极化分量,θ
s
和为平面电磁波的散射角，为平面电磁波的散射方向单位矢量，和为电磁波的散射方向上的垂直和水平极化单位矢量；
[0048]第j个单元的二次辐射场(频域)表示为：
[0049][0050]将公式(10)两边进行傅里叶逆变换，得到第j个单元的二次辐射场时域表示为：
[0051][0052]其中,为第j个单元在散射角方向上半径为R
s
的监测球面上二次辐射远场的延迟时间，单位为秒，为冲激偶本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种时域电磁散射分析的有限冲激偶序列表征方法，其特征在于：包括以下步骤：S1.给定已知频域平面电磁波，并进行时域表征；S2.对散射体表面进行离散化，并确定在平面电磁波照射下散射体表面单元感应电流；S3.确定散射体表面单元二次辐射远场，并进行时间离散与辐射场序列化；S4.在任意平面脉冲波形下的确定时域散射响应信号离散序列，即有限冲激偶序列。2.根据权利要求1所述的一种时域电磁散射分析的有限冲激偶序列表征方法，其特征在于：所述步骤S1包括：给定已知频域平面电磁波，表示为：给定已知频域平面电磁波，表示为：其中，ω＝2πf，为工作角频率，f为工作频率，单位为Hz，E
θ
为平面电磁波电场的垂直极化分量，为平面电磁波电场的水平极化分量,θ和为平面电磁波的入射角，为平面电磁波的入射方向单位矢量，为平面电磁波的入射方向单位矢量，和为电磁波的入射方向上的垂直和水平极化单位矢量，为复数量，η＝120π为自由空间电磁波波阻抗，E为电场，H为磁场，为场点坐标矢量；k0＝Jω/c0为自由空间波数；c0为自由空间电磁波速度；对(1)和(2)进行傅里叶反变换，得到时域平面电磁波的表达式：对(1)和(2)进行傅里叶反变换，得到时域平面电磁波的表达式：其中，为时域卷积运算符，δ为狄拉克函数。3.根据权利要求1所述的一种时域电磁散射分析的有限冲激偶序列表征方法，其特征在于：所述步骤S2包括：S201.对散射体表面进行离散化：按照不高于0.2倍最短工作波长λ
m
＝c0/f
u
的网格长度将散射体离散成N个三角形网格单元，其中，f
u
为最高工作频率，有M个点坐标形成的列表[P
i
]和N个点连接列表[Elem
j
]，其中每个单元的点连接Elem
j
包含三个顶点的序号i0＝Elem
j
.p0，i1＝Elem
j
.p1和i2＝Elem
j
.p2，这三个顶点序号从列表[P
i
]中得到该单元的三个顶点坐标P[i0]、P[i1]和P[i2]，P[i]的点坐标可以表示为(P[i].x,P[i].y,P[i].z)；其中，i∈[0,1,2,
…
,M],j∈[0,1,2,
…
,N]，i为点序号，j为单元的点连接序号；S202.确定在平面电磁波照射下散射体表面单元感应电流:由物理光学近似，散射体表面被平面电磁波照亮区域上任意一点r在某一个时刻t的感应电流J
PO
(r,t)为：其中，为散射体表面在r位置上的单位法向矢量；
×
为矢量叉乘符号；
而对于信号为δ(t)的脉冲平面电磁波，其入射到位置上的瞬态磁场分量表示为：量表示为：量表示为：计算得到的d0为脉冲平面电磁波到达r场点位置的迟滞时间。从N个单元中取出第j个单元，该三角形单元的点连接中的三个点坐标分别为和则有单元j中心点坐标为：由于三角形边长不高于0.2λ
m
，因此将三角形单元内的场当成均匀场来进行计算，用单元中心点的感应场作为整个三角形单元内的感应场，则有公式(7)和公式(8)计算第j个单元的表面感应电流：其中，为的坐标，为表面在的指向散射体外侧的单位法向矢量，为入射平面波在第j个单元中心位置的感...

【专利技术属性】
技术研发人员：李尧尧，杨顺川，
申请(专利权)人：东莘电磁科技成都有限公司，
类型：发明
国别省市：