一种获取瞬态电磁脉冲阵列天线辐射场的方法与流程

1.本发明属于时域阵列天线技术领域，具体涉及一种获取瞬态电磁脉冲阵列天线辐射场的方法。


背景技术：

2.瞬态电磁脉冲技术是高功率微波领域的前沿课题。在瞬态电磁脉冲辐射技术研究中，由于单一辐射单元的功率容量受器件的物理性能限制且时域方向图波束宽度不能灵活调节，因此多个辐射单元组成阵列进行空间合成是一个重要发展方向。为设计和优化瞬态电磁脉冲阵列天线的规模和布局，使其达到特定的场强分布，需要计算和分析瞬态电磁脉冲阵列天线的辐射场及其均匀性。
3.上述问题的求解计算通常依赖于数值计算方法，但该方法存在建模复杂、计算量过大、计算周期长、硬件资源要求高等问题，且对于大尺寸大规模阵列天线远距离辐射场的计算，数值方法难以实现，所以瞬态电磁脉冲阵列天线辐射场的快速计算方法具有研究价值。
4.因此，需要一种获取瞬态电磁脉冲阵列天线辐射场的方法，在不依赖数值计算的情况下，能够快速且较为准确地计算任意观察点的辐射场。


技术实现要素：

5.本发明的目的是提供一种获取瞬态电磁脉冲阵列天线辐射场的方法，解决现有技术中依靠数值方法获取阵列天线辐射场困难，且大尺寸大规模阵列天线远距离辐射场的获取数值方法难以实现的技术问题。
6.为达到上述目的，解决上述技术问题，本发明的技术方案如下：
7.一种获取瞬态电磁脉冲阵列天线辐射场的方法，包括以下步骤：
8.步骤1：设定平面阵列天线的阵元间距d
x
、dy，其中：d
x
、dy分别是x、y方向的阵元间距；搭建a
×
b元子阵，a≥3、b≥3，并采集阵列中心某一阵元在远场距离r处的时域方向图信息；
9.远场是指场强距离乘积近似不变的区域；
10.时域方向图信息为球坐标系下单阵元在远场距离r处不同方向上的电场主要包括h面不同方向的辐射场和e面不同方向的辐射场其中m、n为正整数，具体数值由所需方向图的角度分辨率确定；将各方向辐射场的时间起点对齐以消除延时，得到h面、e面辐射场分别为
11.步骤2：计算阵列天线各阵元s(i,j)辐射场在观察点p(x,y,z)处空间合成的相对延时t
ij
；
12.i，j为阵元在阵列中的行数和列数；
13.p点距离阵列口面中心为ro，且ro应不小于步骤1中的r，确保待计算观察点位于阵元的远场；
14.步骤3：计算各阵元s(i,j)在观察点p(x,y,z)处的近似辐射场
15.步骤3.1：计算观察点p(x,y,z)相对阵元s(i,j)的h面偏角和e面偏角；由于阵元s(i,j)与观察点p(x,y,z)的距离远大于阵元s(i,j)的尺寸，其波程可近似为阵元s(i,j)末端中心(x
ij
,y
ij
,z
ij
)与观察点p(x,y,z)的距离；已知观察点p坐标为(x,y,z)，令则与阵元s(i,j)主轴的h面偏角α
ij
和e面偏角β
ij
[0016][0017][0018]
步骤3.2：假设各阵元的时域方向图相同，根据观察点p(x,y,z)与阵元s(i,j)的距离和方位角数据，结合单阵元在主轴r处h面和e面不同偏角方向上的辐射场，计算出各阵元在观察点p(x,y,z)的近似辐射场；
[0019]
由于各阵元的不一定等于r，根据时域远场的场强距离乘积近似不变的规律，采用公式(3)、(4)分别对h面偏角、e面偏角辐射场进行幅度修正，得到阵元s在|r
ij
|距离上的等效时域方向图信息
[0020][0021][0022]
根据h面偏角α
ij
近似计算出处的辐射场；当且与之间角度较小时，将与之间的波形近似为线性渐变，则
[0023][0024]
同理，根据e面偏角β
ij
近似计算出处的辐射场
[0025][0026]
在h面偏角α
ij
近似计算得到的处辐射场基础上，再考虑e面偏角β
ij
对辐射场修正，得到阵元s(i,j)在观察点p处的近似辐射场为
[0027][0028]
其中，为阵元s(i,j)在其主轴距离处的辐射场；
[0029]
步骤4：将所有阵元s(i,j)的近似辐射场延时t
ij
后叠加合成，计算出观察点p(x,y,z)点的辐射场；
[0030]
根据步骤3得到所有阵元在观察点的近似辐射场结合步骤2得到的各阵元在观察点p(x,y,z)处空间合成的相对延时t
ij
，对各阵元在观察点的辐射场进行延时矢量合成，由公式(8)
[0031][0032]
计算得到瞬态电磁脉冲阵列天线在观察点p(x,y,z)的辐射场。
[0033]
进一步的，计算不同阵元s(i,j)的辐射场到达观察点p因波程差导致的延时，即为各阵元s(i,j)辐射场在观察点p(x,y,z)处空间合成的相对延时t
ij
。
[0034]
进一步的，步骤2～4中的阵列天线能够采用与步骤1阵列间距不同、规模更大的阵列天线，但若采用与步骤1阵列间距相同的阵列天线，计算准确性将更高。
[0035]
本发明的有效收益如下：
[0036]
1、本发明提出的一种瞬态电磁脉冲阵列天线辐射场计算方法，基于辐射场空间矢量合成的原理，采用单阵元辐射场参数数据与理论近似相结合的混合方法，解决了数值计算方法对于大尺寸大规模阵列天线辐射场计算的速度慢、甚至无法计算的问题，且计算结果比较准确。
[0037]
2、本发明基于小型阵列天线单元时域方向图数据，实现更大规模阵列天线的合成场计算，为分析瞬态电磁脉冲阵列天线合成场场强及均匀性、设计和优化阵列天线的规模布局提供有效支撑。
[0038]
3、本发明提出的一种瞬态电磁脉冲阵列天线辐射场计算方法，通过小型阵列的单个阵元时域方向图计算大规模阵列天线的合成场，避免了数值计算方法对大尺寸大规模阵列天线计算的复杂建模过程和庞大的数据量，同时极大地拓展了计算的空间区域。
附图说明
[0039]
图1为本发明实现的流程示意图；
[0040]
图2为本发明单个阵元时域方向图数据采集示意图；
[0041]
图3为本发明实施例示意图。
具体实施方式
[0042]
下面结合附图和具体实施例对本发明进行详细的阐述和说明。
[0043]
图1为本发明提供的计算方法流程图，结合图2、图3，一种瞬态电磁脉冲阵列天线辐射场计算方法，包括以下步骤：
[0044]
步骤1：设定平面阵列天线的阵元间距d
x
、dy(d
x
、dy分别是x、y方向的阵元间距，天
线口面位于xoy平面，天线主轴方向为z轴方向)，搭建a
×
b元子阵(a、b≥3)并采集阵列中心某一阵元在远场距离r处的时域方向图信息，为步骤3提供数据支撑；
[0045]
图2所示为3
×
3平面天线阵列，通过实测或者数值计算等方法采集阵元(2,2)的时域方向图信息，即采集不同方向远场距离r处的电场主要包括h面不同方向的辐射场和e面不同方向的辐射场将主波形的时间起点对齐以消除延时，得到h面、e面各辐射场分别为波形的时间起点对齐以消除延时，得到h面、e面各辐射场分别为
[0046]
步骤2：计算m
×
n元平面阵列天线任一阵元s(i,j)在观察点p(x,y,z)处空间合成的相对时间差t
ij
；
[0047]
以5
×
6阵列(或m为奇数、n为偶数的阵列)为例，如图3所示，其中阵元为平面天线(同图2)，e面尺寸为w
antenna
，厚度忽略不计。计算不同阵元s(i,j)的辐射场到达p点因波程差产生的时间差t
sij
，即可得到p点处空间合成的相对时间差t
ij
＝t
sij
；由于阵元s与p点的距离远大于阵元的尺寸，其波程可近似为阵元末端中心(x
ij
,y
ij
,z
ij
)与p的距离，其中
[0048][0049]
已知p点坐标为(x,y,z)，则阵元的波程可由勾股定理算出。
[0050]
步骤3：计算阵元s(i,j)在观察点p处的近似辐射场
[0051]
步骤3.1：计算观察点p相对阵元s(i,j)的h面偏角α
ij
和e面偏角β
ij
。如图3所示，令在以为半径、(x
ij
,y
ij
,0)为球心的球面上，将分别在h面和e面进行投影，投影延长线分别交于球面上的ph和pe点，则α
ij
＝∠phsa，β
ij
＝∠pesa。计算得到α
ij
和β
ij
分别为
[0052][0053][0054]
步骤3.2：假设各阵元的时域方向图相同，根据观察点p与阵元s(i,j)的距离和方位角数据，结合单元在主轴r处h面和e面不同偏角方向上的辐射场，计算出各阵元在观察点的近似辐射场。
[0055]
由于各阵元的不一定为r，根据时域远场的场强距离乘积近似不变的规律，采用公式(4)、(5)分别对h面偏角、e面偏角辐射场进行幅度修正，可得到阵元s在距离上的等效时域方向图信息
[0056][0057][0058]
根据h面偏角α
ij
近似计算出处的辐射场。当且与之间角度较小时，将与之间的波形近似为线性渐变，则
[0059][0060]
同理，当时，根据e面偏角β
ij
近似计算出处的辐射场波形
[0061][0062]
在h面偏角α
ij
近似计算得到的处辐射场基础上，再考虑e面偏角β
ij
对波形修正，得到阵元s(i,j)在观察点p处的近似辐射场为
[0063][0064]
其中，为阵元s(i,j)在其主轴距离处的辐射场。
[0065]
步骤4：将阵元s(i,j)的辐射场延时t
ij
后叠加合成，计算出阵列在p点的辐射场；
[0066]
根据步骤3得到所有阵元在p点的近似辐射场并通过步骤2得到的各阵元辐射电磁波到达p点的合成相对时间差t
ij
，对各阵元在p点的辐射场进行延时平移和矢量合成，由公式(9)计算得到瞬态电磁脉冲阵列天线在观察点p的辐射场
[0067][0068]
以上是本发明的较佳实施例，并非是对本发明作其他形式的限制，任何熟悉本专业的技术人员可能利用上述揭示的技术内容加以变更或改型为等同变化的等效实施例。凡依本发明技术方案所作的任何简单修改、等同变化与改型，仍属于本发明的保护范围。