- 1、本文档共5页,可阅读全部内容。
- 2、有哪些信誉好的足球投注网站(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
时域有限差分法的Matlab仿真
时域有限差分法的Matlab仿真
春波
(华北电力大学 电子与通信工程系河北 保定071003)
摘 要:介绍了时域有限差分法的基本原理,并利用Matlab仿真,对矩形波导谐振腔中的电磁场作了模拟和分析。关键词:时域有限差分法;Matlab;矩形波导;谐振腔
目前,电磁场的时域计算方法越来越引人注目。时域有限差分(Finite Difference TimeDomain,FDTD)法[1]作为一种主要的电磁场时域计算方法,最早是在1966年由K. S. Yee提出的。这种方法通过将Maxwell旋度方程转化为有限差分式而直接在时域求解,通过建立时间离散的递进序列,在相互交织的网格空间中交替计算电场和磁场。经过三十多年的发展,这种方法已经广泛应用到各种电磁问题的分析之中。作为一种工程仿真工具得到了广泛应用[2]。用于时域有限差分法,可以简化编程,使研究者的研究重心放在FDTD法本身上,而不必在编程上花费过多的时间。下面将采用FDTD法,利用Matlab仿真来分析矩形波导谐振腔的电磁场,说明了将二者结合起来的优越性。
1FDTD法基本原理时域有限差分法的主要思想是把Maxwell方程在空间、时间上离散化,用差分方程代替一阶偏微分方程,求解差分方程组,从而得出各网格单元的场值。FDTD空间网格单元上电场和磁场各分量的分布如图1所示。
电场和磁场被交叉放置,电场分量位于网格单元每条棱的中心,磁场分量位于网格单元每个面的中心,每个磁场(电场)分量都有4个电场(磁场)分量环绕。这样不仅保证了介质分界面上切向场分量的连续性条件得到自然满足,而且还允许旋度方程在空间上进行中心差分运算,同时也满足了法拉第电磁感应定律和安培环路积分定律,也可以很恰当地模拟电磁波的实际传播过程。1.1Maxwell方程的差分形式旋度方程为:将其标量化,并将问题空间沿3个轴向分成若干网格单元,用Δx,Δy和Δz分别表示每个网格单元沿3个轴向的长度,用Δt表示时间步长。网格单元顶点的坐标(x,y,z)可记为:其中:i,j,k和n为整数。同时利用二阶精度的中心有限差分式来表示函数对空间和时间的偏导数,即可得到如下FDTD基本差分式:由于方程式里出现了半个网格和半个时间步,为了便于编程,将上面的差分式改写成如下形式:其中:ε和μ分别为介质的介电常数和磁导率。1.2数值色散及稳定性条件为了减小数值色散,在选取空间网格尺寸时,应满足λmin≥10Δ,Δ=min(Δx,Δy,Δz),λmin是被研究媒质空间的最小波长值。由此可以看出:减小网格尺寸可以减小数值色散,但是会引起计算存储量的增大,因此需综合考虑,权衡处理。为了使数值计算稳定,时间步长的选择应满足:2矩形波导谐振腔模型及仿真分析矩形波导谐振腔是由两端短路的一段金属波导构成,如图2所示。本文选用的是WJB100型矩形波导,即宽度a=22.86 mm,高度b=10.16 mm,选取谐振腔长度l=50.80 mm。由微波理论[3]可知,当bAL时,TE101模的谐振波长最长,是矩形波导谐振腔的主模。TE101模在矩形腔的3个方向都不传输能量,呈驻波分布,并且电场分量只有Ey分量,在腔体中央最强;磁场有Hx和Hz两个分量,在腔壁附近最强,腔体中央为0。由于矩形波导谐振腔有着广泛的应用,并且大多工作在主模状态,这里选择微分高斯脉冲以激励起TE101模。
决定网格单元的尺寸和时间步长取Δx=Δy=Δz=Δ=1.27mm,则空间网格数为18×8×40,取Δt=Δ/(2c)2.118ps。设置激励源微分高斯脉冲的表达式为:其优点是不含零频率分量,时域波形和频谱如图3所示。为了在谐振腔中激励起TE101模,并且抑制其他高次模,选择线源脉冲,使之在腔内xz平面中心处沿y轴方向分布,并选择合适的t0和τ值。经过反复试验,取t0=95.316ps,τ=285.948ps。
(3)设置边界条件对于矩形波导谐振腔,腔体的6个面都是金属,为方便起见,这里假设都是理想导体,即腔内导体边界上的所有切向电场分量为0,所有法向磁场分量为0。矩形谐振腔的谐振波长为:则TE101模的谐振波长为λ0=41.69mm,对应的谐振频率f0=7.19GHz。图4是Matlab仿真计算出的矩形波导谐振腔的谐振频率,很明显,与理论值非常接近,由此可以确定确实激励起了TE101模,并且较好地抑制了高次模。图5是程序运行N=3 000步时,腔体中央xz平面上的各分量瞬时分布图,可以看出,符合电磁场理论。
3结语以上结合FDTD和Matlab对矩形波导谐振腔做了仿真分析,所编Matlab程序简洁明了,运行效率也较高。FDTD法在电磁场数值分析方面有很大的优越性,而Matlab具有强大的数据处理和图形处理功能,
您可能关注的文档
- 数据库中的模糊逻辑规则发现-计算机工程与应用.PDF
- 数据库理论与技术.DOC
- 数据库管理系统-忻州师范学院.DOC
- 数据库系统原理____第九章.PPT
- 数据挖掘-中国数字科技馆.PDF
- 数据挖掘与数据仓库概述.PDF
- 数据挖掘与检索.DOC
- 数据挖掘与检索-MediaandCommunicationLab.PDF
- 数据挖掘-学者网.PPT
- 数据挖掘初探.PDF
- 《GBT 43565-2023中小学合成材料面层篮球场地》必威体育精装版解读.pptx
- 《GBT 12910-2023纸和纸板 二氧化钛含量的测定》必威体育精装版解读.pptx
- 《GBT 43603.1-2023镍铂靶材合金化学分析方法 第1部分:铂含量的测定 电感耦合等离子.pptx
- 《GBT 9364.6-2023小型熔断器 第6部分:小型熔断体用熔断器支持件》必威体育精装版解读.pptx
- 《GBT 36545-2023移动式电化学储能系统技术规范》必威体育精装版解读.pptx
- 《GBT 42737-2023电化学储能电站调试规程》必威体育精装版解读.pptx
- 《GBT 43397-2023婴儿救生衣》必威体育精装版解读.pptx
- 《GBT 17626.39-2023电磁兼容 试验和测量技术 第39部分:近距离辐射场抗扰度试验》.pptx
- 《GBT 43559-2023蜂胶生产技术规范》必威体育精装版解读.pptx
- 《GBT 43623-2023大黄栽培技术规程》必威体育精装版解读.pptx
文档评论(0)