一种用于OTA仿真的复包络蛙跳隐式时域有限差分法

技术领域

本发明属于计算电磁学技术领域，尤其涉及一种用于OTA仿真的复包络蛙跳隐式时域有限差分法。

背景技术

随着计算电磁学的发展，时域有限差分法(FDTD)被广泛的应用于电磁问题的求解。然而在一些计算问题中，传统的FDTD算法需要计算极长的时间步数，而传统的FDTD方法受到稳定性条件的限制，无法加大时间步长，从而使得整体计算效率低下。比如对收发机的OTA测试，为了获得较好的结果，要发送不少于一帧的数据，如果利用传统FDTD算法，则需要运行百万千万量级的时间步长。为了突破时间步长的限制，人们开始研究无条件稳定的FDTD算法，交替方向隐格式(ADI)FDTD算法和蛙跳式(leapfrog)ADI-FDTD算法也随之出现。但由于时间步长变大后，造成了精度的下降。因此需要提高长时间仿真下带限系统的仿真精度。

发明内容

本发明目的在于提供一种用于OTA仿真的复包络蛙跳隐式时域有限差分法,以解决传统时域有限差分法时间步长过小，计算效率低，交替方向隐格式算法和蛙跳式算法精度低的技术问题。

为解决上述技术问题，本发明的具体技术方案如下：

一种用于OTA仿真的复包络蛙跳隐式时域有限差分法，包括以下步骤：

步骤1、设定求解的中心频率和带宽，并根据频率与波长的关系计算相应的波长；

步骤2、根据仿真的环境设定相应的求解区域，确定所求空间的大小，确定空间边界的处理方式，确定空间内的散射物体位置及其电磁常数，确定激励点和接收点，然后根据Yee网格对空间进行离散化，确定空间的网格数、物体的网格位置和激励点、接收点所在的网格，最后由CFL稳定性条件求对应的时间步长，根据仿真精度需要为仿真设定时间步长倍数；

步骤3、指定发射信号的调制方式并生成发射信号，根据信号长度设定仿真的时间步数；

步骤4、根据Yee网格数，为各个电磁场分量申请内存空间并进行初始化，然后对更新迭代方程中需要用到的系数进行保存和初始化；

步骤5、对电磁场进行更新迭代，首先更新区域内的电场值，电场计算完成后加入电激励源，然后更新磁场值，加入磁激励源，最后保存接收点位置上的信号，作为该时刻的输出；通过电场、磁场不停的迭代，计算到设定的时间总步长后停止；最终可以得到接收点处整个时间的接收信号；

步骤6、根据信号的调制方式进行解调，并求解参数指标。

进一步的，步骤4中迭代系数采用索引方式进行存储，节省了内存空间，提高了内存的利用效率。

本发明的一种用于OTA仿真的复包络蛙跳隐式时域有限差分法具有以下优点：

本发明对仿真空间进行电磁仿真来获取收发之间的响应，可适用于OTA测试的仿真，具有精度高和无条件稳定的特点。

附图说明

图1为本发明的仿真计算流程图；

图2为本发明的方法误差随归一化频率变化的图像示意图；

图3为发射端信号的星座图；

图4为发射端信号的频谱图；

图5为接收端信号的星座图。

具体实施方式

为了更好地了解本发明的目的、结构及功能，下面结合附图，对本发明一种用于OTA仿真的复包络蛙跳隐式时域有限差分法做进一步详细的描述。

如图1所示，是本发明的仿真计算流程图。

本发明用于OTA仿真的复包络蛙跳隐式时域有限差分法，包括以下步骤：

步骤1、设定求解的中心频率和带宽，并根据频率与波长的关系计算相应的波长；

步骤2、根据仿真的环境设定相应的求解区域，确定所求空间的大小，确定空间边界的处理方式，确定空间内的散射物体位置及其电磁常数，确定激励点和接收点，然后根据Yee网格对空间进行离散化，确定空间的网格数、物体的网格位置和激励点、接收点所在的网格，最后由CFL稳定性条件求对应的时间步长，根据仿真精度需要为仿真设定时间步长倍数；

步骤3、指定发射信号的调制方式并生成发射信号，根据信号长度设定仿真的时间步数；

步骤4、根据Yee网格数，为各个电磁场分量申请内存空间并进行初始化，然后对更新迭代方程中需要用到的系数进行保存和初始化，所述的迭代系数采用索引方式进行存储，节省了内存空间，提高了内存的利用效率；

步骤5、对电磁场进行更新迭代，首先更新区域内的电场值，电场计算完成后加入电激励源，然后更新磁场值，加入磁激励源，最后保存接收点位置上的信号，作为该时刻的输出；通过电场、磁场不停的迭代，计算到设定的时间总步长后停止；最终可以得到接收点处整个时间的接收信号；

步骤6、根据信号的调制方式进行解调，并求解参数指标。

在实施例中对一毫米波暗室的OTA测量环境进行分析，计算接收端信号的误差向量幅度(Error Vector Magnitude，EVM)值，包括以下步骤：

步骤一、设定求解频率、带宽并计算响应波长。

假设仿真中心频率f

步骤二、设定求解区域并进行空间和时间的离散化，确定时间步长。

假设用于28GHz测量的小型暗箱长0.3m，宽0.3m，高0.1m，六个面均布置有渐变的吸波介质。激励点位于点(0.15m，0.15m，0.02m)处，观测点分别位于(0.15m，0.15m，0.04m)处。取x，y，z方向空间离散步长为0.8mm，对求解区域按Yee网格进行离散。空间网格数为375×375×125，激励点位于网格(187，187，25)，观测点分别位于网格(187，187，50)，根据式(1)可得CFL稳定性条件中对应的时间步长Δt

由于本发明所提出的算法无条件稳定，因此，为了提高精度，根据图2所示误差随归一化频率变化曲线，取离散时间步长Δt＝2Δt

步骤三、确定发射信号和仿真总时间步数。

取发射信号激励为z方向的电场激励

步骤四、申请内存空间并进行初始化。

电磁场空间包括Ex，Ey，Ez，Hx，Hy和Hz六个分量，每个分量需要用到375×375×125的复数double类型内存空间，并初始化为0+0i。然后对更新迭代方程中需要用到的系数以索引方式进行保存和初始化。

步骤五、推导更新迭代公式，并按照公式对电磁场进行迭代。

将电场和磁场写成如式(2)的复包络形式

其中

将式(2)带入FDTD的更新方程并结合ADI方法和leapfrog方法，可得迭代更新方程为

其中

算子

将式(3)和(4)写成分量形式为

于是，在每一个时间步中依次进行如下操作：

(一)电场计算

采用追赶法求解式(7)(8)(9)对应的三对角方程。

以x方向的电场为例，等式(7)左侧为未知项，右侧全部为已知项。对于j∈[1,N

(2)电场激励处理

此步仅在指定的电场激励位置网格(187，187，25)的Ez分量上进行。

(3)磁场计算

采用追赶法求解式(10)(11)(12)对应的三对角方程.

(4)磁场激励处理

本实施例不进行磁场激励，故不做处理。

(5)保存当前时间下观测点的值。

以此法迭代N

步骤六、对观测信号进行解调。

对观测点的接收信号进行OFDM解调和QAM64解调，可得接收星座图如图5所示，通过计算可得均方根EVM＝0.2487％。

本发明方法所对应的增长矩阵的特征值始终不大于1，因此在任何步长下均无条件稳定。如图2所示，给出了本方法误差随归一化频率(以中心频率归一化)变化的曲线，可以看到在所设定的中心频率附近具有极高的精度。

本发明给出了一种用于OTA仿真的复包络蛙跳隐式时域有限差分法，可适用于一维、二维及三维的环境中。该方法可对仿真空间进行电磁仿真来获取收发之间的响应，尤其适用于OTA测试的仿真，具有精度高和无条件稳定的特点。

可以理解，本发明是通过一些实施例进行描述的，本领域技术人员知悉的，在不脱离本发明的精神和范围的情况下，可以对这些特征和实施例进行各种改变或等效替换。另外，在本发明的教导下，可以对这些特征和实施例进行修改以适应具体的情况及材料而不会脱离本发明的精神和范围。因此，本发明不受此处所公开的具体实施例的限制，所有落入本申请的权利要求范围内的实施例都属于本发明所保护的范围内。