CN113051804A - 一种用于ota仿真的复包络蛙跳隐式时域有限差分法 - Google Patents
一种用于ota仿真的复包络蛙跳隐式时域有限差分法 Download PDFInfo
- Publication number
- CN113051804A CN113051804A CN202110435889.6A CN202110435889A CN113051804A CN 113051804 A CN113051804 A CN 113051804A CN 202110435889 A CN202110435889 A CN 202110435889A CN 113051804 A CN113051804 A CN 113051804A
- Authority
- CN
- China
- Prior art keywords
- simulation
- ota
- finite difference
- space
- difference method
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000004088 simulation Methods 0.000 title claims abstract description 37
- 238000000034 method Methods 0.000 title claims abstract description 29
- 230000005684 electric field Effects 0.000 claims abstract description 17
- 230000005284 excitation Effects 0.000 claims description 20
- 230000005672 electromagnetic field Effects 0.000 claims description 8
- 238000012545 processing Methods 0.000 claims description 4
- 239000012530 fluid Substances 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 abstract description 9
- 230000004044 response Effects 0.000 abstract description 4
- 238000012360 testing method Methods 0.000 abstract description 4
- 238000004422 calculation algorithm Methods 0.000 description 8
- 238000010586 diagram Methods 0.000 description 4
- 230000006870 function Effects 0.000 description 3
- 239000011159 matrix material Substances 0.000 description 2
- 238000001228 spectrum Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/02—Reliability analysis or reliability optimisation; Failure analysis, e.g. worst case scenario performance, failure mode and effects analysis [FMEA]
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Monitoring And Testing Of Transmission In General (AREA)
Abstract
本发明公开了一种用于OTA仿真的复包络蛙跳隐式时域有限差分法。该方法采用复包络交替方向隐格式时域有限差分法对全空间进行电磁仿真来获取收发之间的响应,从而达到获取性能指标的目的。首先,在设置空间媒质分布、确定信号调制方式的基础上,通过复包络蛙跳隐式时域有限差分法的更新方程对电场和磁场进行不断迭代,直到达到预设的时间总步数后停止;然后根据调制方式对接收信号进行解调以获取参数指标。本发明方法具有无条件稳定和高精度的特点,适用性广泛,可用于进行OTA测试的仿真计算,也可用于其他复杂介质下的电磁计算问题。
Description
技术领域
本发明属于计算电磁学技术领域,尤其涉及一种用于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)值,包括以下步骤:
步骤一、设定求解频率、带宽并计算响应波长。
假设仿真中心频率fc为28GHz,带宽BW为1GHz。因此中心频率波长λc=C/fc为10.71mm,其中C为电磁波在真空中的波长。
步骤二、设定求解区域并进行空间和时间的离散化,确定时间步长。
假设用于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稳定性条件中对应的时间步长ΔtCFL=1.5396×10-12s。
由于本发明所提出的算法无条件稳定,因此,为了提高精度,根据图2所示误差随归一化频率变化曲线,取离散时间步长Δt=2ΔtCFL=3.0792×10-12s。
步骤三、确定发射信号和仿真总时间步数。
取发射信号激励为z方向的电场激励其波形为随机序列对应的QAM64调制的OFDM信号。信号星座图和频谱分别如图3和图4所示。根据所生成的随机信号长度以及设定的时间步长,取迭代总时间步数Nt=2×106步。
步骤四、申请内存空间并进行初始化。
电磁场空间包括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,Ny-1]∩Z,相邻的三个电场即构成一个方程,共构成Ny-1个方程,当j=0和j=Ny时,由计算空间的截断边界条件分别构成一个方程,共2个方程,因此共Ny+1个方程,又有Ny+1个未知数,方程可解。并且,方程系数全部位于主对角线和上下两次对角线,可以由追赶法进行快速求解。其他电场分量同理。
(2)电场激励处理
此步仅在指定的电场激励位置网格(187,187,25)的Ez分量上进行。
(3)磁场计算
采用追赶法求解式(10)(11)(12)对应的三对角方程.
(4)磁场激励处理
本实施例不进行磁场激励,故不做处理。
(5)保存当前时间下观测点的值。
以此法迭代Nt=2×106步,可得观测点随时间的信号情况。
步骤六、对观测信号进行解调。
对观测点的接收信号进行OFDM解调和QAM64解调,可得接收星座图如图5所示,通过计算可得均方根EVM=0.2487%。
本发明方法所对应的增长矩阵的特征值始终不大于1,因此在任何步长下均无条件稳定。如图2所示,给出了本方法误差随归一化频率(以中心频率归一化)变化的曲线,可以看到在所设定的中心频率附近具有极高的精度。
本发明给出了一种用于OTA仿真的复包络蛙跳隐式时域有限差分法,可适用于一维、二维及三维的环境中。该方法可对仿真空间进行电磁仿真来获取收发之间的响应,尤其适用于OTA测试的仿真,具有精度高和无条件稳定的特点。
可以理解,本发明是通过一些实施例进行描述的,本领域技术人员知悉的,在不脱离本发明的精神和范围的情况下,可以对这些特征和实施例进行各种改变或等效替换。另外,在本发明的教导下,可以对这些特征和实施例进行修改以适应具体的情况及材料而不会脱离本发明的精神和范围。因此,本发明不受此处所公开的具体实施例的限制,所有落入本申请的权利要求范围内的实施例都属于本发明所保护的范围内。
Claims (2)
1.一种用于OTA仿真的复包络蛙跳隐式时域有限差分法,其特征在于,包括以下步骤:
步骤1、设定求解的中心频率和带宽,根据频率与波长的关系计算中心频率的波长;
步骤2、根据仿真的环境设定相应的求解区域,确定所求空间的大小,确定空间边界的处理方式,确定空间内的散射物体位置及其电磁常数,确定激励点和接收点,然后根据Yee网格对空间进行离散化,确定空间的网格数、物体的网格位置和激励点、接收点所在的网格,最后由CFL稳定性条件求对应的时间步长,根据仿真精度需要为仿真设定时间步长倍数;
步骤3、指定发射信号的调制方式并生成发射信号,根据信号长度设定仿真的时间步数;
步骤4、根据Yee网格数,为各个电磁场分量申请内存空间并进行初始化,然后对更新迭代方程中需要用到的系数进行保存和初始化;
步骤5、对电磁场进行更新迭代,首先更新区域内的电场值,电场计算完成后加入电激励源,然后更新磁场值,加入磁激励源,最后保存接收点位置上的信号,作为该时刻的输出;通过电场、磁场不停的迭代,计算到设定的时间总步长后停止;最终可以得到接收点处整个时间的接收信号;
步骤6、根据信号的调制方式进行解调,并求解参数指标。
2.根据权利要求1所述的用于OTA仿真的复包络蛙跳隐式时域有限差分法,其特征在于,步骤4中所述的迭代系数采用索引方式进行存储。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110435889.6A CN113051804A (zh) | 2021-04-22 | 2021-04-22 | 一种用于ota仿真的复包络蛙跳隐式时域有限差分法 |
PCT/CN2022/071654 WO2022237218A1 (zh) | 2021-04-22 | 2022-01-12 | 一种用于ota仿真的复包络蛙跳隐式时域有限差分法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110435889.6A CN113051804A (zh) | 2021-04-22 | 2021-04-22 | 一种用于ota仿真的复包络蛙跳隐式时域有限差分法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN113051804A true CN113051804A (zh) | 2021-06-29 |
Family
ID=76519917
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110435889.6A Pending CN113051804A (zh) | 2021-04-22 | 2021-04-22 | 一种用于ota仿真的复包络蛙跳隐式时域有限差分法 |
Country Status (2)
Country | Link |
---|---|
CN (1) | CN113051804A (zh) |
WO (1) | WO2022237218A1 (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113987792A (zh) * | 2021-10-26 | 2022-01-28 | 山东大学 | 一种fdtd算法中实现模式源精确输入的方法 |
WO2022237218A1 (zh) * | 2021-04-22 | 2022-11-17 | 东南大学 | 一种用于ota仿真的复包络蛙跳隐式时域有限差分法 |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116861711B (zh) * | 2023-09-05 | 2023-11-17 | 北京航空航天大学 | 一种基于时域间断伽辽金方法的双各向同性媒质仿真方法 |
CN117236085B (zh) * | 2023-11-16 | 2024-02-13 | 芯瑞微(上海)电子科技有限公司 | 基于fdtd算法对圆柱体电磁场模拟的精度优化方法 |
CN117390935B (zh) * | 2023-12-11 | 2024-03-01 | 芯瑞微(上海)电子科技有限公司 | 一种计算fdtd电磁仿真收敛检测触发时刻的方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20120057274A (ko) * | 2010-11-26 | 2012-06-05 | (주)애니캐스팅 | Id-fdtd 방법을 이용한 전자기파 수치해석 방법 |
CN104408256A (zh) * | 2014-12-01 | 2015-03-11 | 天津工业大学 | 一种截断一维Debye介质Crank-Nicolson完全匹配层实现算法 |
CN107239586A (zh) * | 2016-03-29 | 2017-10-10 | 南京理工大学 | 对无条件稳定时域有限差分法有效的区域分解并行方法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6772076B2 (en) * | 2002-04-24 | 2004-08-03 | Ffc Limited | Electromagnetic field analysis method based on FDTD method, medium representation method in electromagnetic field analysis, simulation device, and storage medium |
CN107526887B (zh) * | 2017-08-22 | 2020-02-18 | 电子科技大学 | 一种基于GPU并行的LeapfrogADI-FDTD方法 |
CN113051804A (zh) * | 2021-04-22 | 2021-06-29 | 东南大学 | 一种用于ota仿真的复包络蛙跳隐式时域有限差分法 |
-
2021
- 2021-04-22 CN CN202110435889.6A patent/CN113051804A/zh active Pending
-
2022
- 2022-01-12 WO PCT/CN2022/071654 patent/WO2022237218A1/zh active Application Filing
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20120057274A (ko) * | 2010-11-26 | 2012-06-05 | (주)애니캐스팅 | Id-fdtd 방법을 이용한 전자기파 수치해석 방법 |
CN104408256A (zh) * | 2014-12-01 | 2015-03-11 | 天津工业大学 | 一种截断一维Debye介质Crank-Nicolson完全匹配层实现算法 |
CN107239586A (zh) * | 2016-03-29 | 2017-10-10 | 南京理工大学 | 对无条件稳定时域有限差分法有效的区域分解并行方法 |
Non-Patent Citations (6)
Title |
---|
付梅艳等: "基于三角形单元的时域有限积分方法研究", 《微波学报》, no. 1, 15 August 2010 (2010-08-15), pages 83 - 87 * |
付梅艳等: "基于三角形单元的时域有限积分方法研究", 微波学报, no. 1, 15 August 2010 (2010-08-15), pages 83 - 87 * |
王长华等: "基于时域有限差分法的谐振腔分析", 《现代电子技术》, no. 2, 15 January 2011 (2011-01-15), pages 103 - 105 * |
王长华等: "基于时域有限差分法的谐振腔分析", 现代电子技术, no. 2, 15 January 2011 (2011-01-15), pages 103 - 105 * |
董云扬: "同频率全双工系统信道仿真精确算法的研究", 《中国优秀硕士学位论文全文数据库(电子期刊)-信息科技辑》, 15 August 2016 (2016-08-15), pages 136 - 379 * |
董云扬: "同频率全双工系统信道仿真精确算法的研究", 中国优秀硕士学位论文全文数据库(电子期刊)-信息科技辑, 15 August 2016 (2016-08-15), pages 136 - 379 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2022237218A1 (zh) * | 2021-04-22 | 2022-11-17 | 东南大学 | 一种用于ota仿真的复包络蛙跳隐式时域有限差分法 |
CN113987792A (zh) * | 2021-10-26 | 2022-01-28 | 山东大学 | 一种fdtd算法中实现模式源精确输入的方法 |
CN113987792B (zh) * | 2021-10-26 | 2024-04-16 | 山东大学 | 一种fdtd算法中实现模式源精确输入的方法 |
Also Published As
Publication number | Publication date |
---|---|
WO2022237218A1 (zh) | 2022-11-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113051804A (zh) | 一种用于ota仿真的复包络蛙跳隐式时域有限差分法 | |
Glisson et al. | Simple and efficient numerical methods for problems of electromagnetic radiation and scattering from surfaces | |
Sun et al. | Efficient implementations of the Crank-Nicolson scheme for the finite-difference time-domain method | |
Kong et al. | High-order split-step unconditionally-stable FDTD methods and numerical analysis | |
Zunoubi et al. | CUDA Implementation of ${\rm TE}^{z} $-FDTD Solution of Maxwell's Equations in Dispersive Media | |
CN114528716B (zh) | 一种应用于多尺度电磁波问题分析的高效数值仿真方法 | |
Boag et al. | A fast physical optics (FPO) algorithm for double-bounce scattering | |
EP2287757A1 (en) | Multilevel-Multigrid simulation techniques | |
Omiya et al. | Design of cavity-backed slot antennas using the finite-difference time-domain technique | |
CN113987792A (zh) | 一种fdtd算法中实现模式源精确输入的方法 | |
JPH11295365A (ja) | 電磁界解析方法 | |
Steinle et al. | An implicit flux-corrected transport algorithm | |
Liu et al. | Application of Chaos Mind Evolutionary Algorithm in Antenna Arrays Synthesis. | |
Yilmaz et al. | A single-boundary implicit and FFT-accelerated time-domain finite element-boundary integral solver | |
Schobert et al. | Low-frequency surface integral equation solution by multilevel Green's function interpolation with fast Fourier transform acceleration | |
Hochman et al. | A numerical methodology for efficient evaluation of 2D Sommerfeld integrals in the dielectric half-space problem | |
Sun et al. | The transfinite-element time-domain method | |
Kim et al. | Adaptive integral method for higher order method of moments | |
Lei et al. | Study on MPI-based parallel modified conformal FDTD for 3-D electrically large coated targets by using effective parameters | |
CN106485071A (zh) | 一种多层分组结构快速近远场转换方法 | |
Zhang et al. | Study on the optimum virtual topology for MPI based parallel conformal FDTD algorithm on PC clusters | |
Bilbao et al. | Numerical analysis of electromagnetic fields | |
Zhang et al. | MPI-based parallelized locally conformal FDTD for modeling slot antennas and new periodic structures in microstrip | |
Boyd | Continuum Breathers: Forty Years After | |
Yu et al. | Massively parallel conformal FDTD on a BlueGene supercomputer |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination |