急求个用matlab时域分析做的时域有限差分的界面

《MATLAB模拟的电磁学时域有限差分法》一书附带光盘 - 下载频道 - CSDN.NET
&&&&《MATLAB模拟的电磁学时域有限差分法》一书附带光盘
&《MATLAB模拟的电磁学时域有限差分法》一书附带光盘
《MATLAB模拟的电磁学时域有限差分法》一书附带光盘
若举报审核通过,可奖励20下载分
被举报人:
举报的资源分:
请选择类型
资源无法下载
资源无法使用
标题与实际内容不符
含有危害国家安全内容
含有反动色情等内容
含广告内容
版权问题,侵犯个人或公司的版权
*详细原因:
您可能还需要
Q.为什么我点的下载下不了,但积分却被扣了
A. 由于下载人数众多,下载服务器做了并发的限制。若发现下载不了,请稍后再试,多次下载是不会重复扣分的。
Q.我的积分不多了,如何获取积分?
A. 获得积分,详细见。
完成任务获取积分。
评价资源返积分。
论坛可用分兑换下载积分。
第一次绑定手机,将获得5个C币,C币可。
下载资源意味着您已经同意遵守以下协议
资源的所有权益归上传用户所有
未经权益所有人同意,不得将资源中的内容挪作商业或盈利用途
CSDN下载频道仅提供交流平台,并不能对任何下载资源负责
下载资源中如有侵权或不适当内容,
本站不保证本站提供的资源的准确性,安全性和完整性,同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
行业下载排行
你下载资源过于频繁,请输入验证码
如何快速获得积分?
你已经下载过该资源,再次下载不需要扣除积分
《MATLAB模拟的电磁学时域有限差分法》一书附带光盘
所需积分:10
剩余积分:
VIP会员,免积分下载
会员到期时间:日
剩余下载次数:1000
VIP服务公告:时域有限差分法的Matlab仿真_百度文库
两大类热门资源免费畅读
续费一年阅读会员,立省24元!
评价文档:
时域有限差分法的Matlab仿真
时​域​有​限​差​分​法​的​M​a​t​l​a​b​仿​真
阅读已结束,如果下载本文需要使用
想免费下载本文?
把文档贴到Blog、BBS或个人站等:
普通尺寸(450*500pix)
较大尺寸(630*500pix)
你可能喜欢时域有限差分法对平面TE波的MATLAB仿真71-第3页
上亿文档资料,等你来发现
时域有限差分法对平面TE波的MATLAB仿真71-3
Hx(i,j,k)?CP(i,j,k)Hx(i,;?EZ(i,j?1,k)?EZ(i,j,k)Ey;????;?y?z??;Hy(i,j,k)?CP(i,j,k)Hx(i,;?E(i,j,k?1)?Ex(i,j,k)EZ(;??x??;?z?x??HZ(i,j,k)?CP(i,j,k;?Ey(i?1,j,k)?Ey(i,j,k)Ex;????;?x?y??;
Hx(i,j,k)?CP(i,j,k)Hx(i,j,k)?CQ(i,j,k)?EZ(i,j?1,k)?EZ(i,j,k)Ey(i,j,k?1)?Ey(i,j,k)?(2-8d)?????y?z??Hy(i,j,k)?CP(i,j,k)Hx(i,j,k)?CQ(i,j,k)?E(i,j,k?1)?Ex(i,j,k)EZ(i?1,j,k)?EZ(i,j,k)?(2-8e)??x???z?x??HZ(i,j,k)?CP(i,j,k)HZ(i,j,k)?CQ(i,j,k)?Ey(i?1,j,k)?Ey(i,j,k)Ex(i,j?1,k)?Ex(i,j,k)?(2-8f)?????x?y??根据上述FDTD差分方程组可得出计算电磁场的时域推进计算方法,如图2-2所示。 图2-2
FDTD在时域的交叉半步逐步推进计算式(2-8a)~(2-8c)的等号左边的电场值是第n次循环的电场值,等号右边的电场值是第n-1次循环存储在内存中的电场值,磁场值是本次循环计算得到的磁场值;式(2-8d)~(2-8f) 等号左边的磁场值是第n次循环的磁场值,等号右边的磁场值和电场值都是第n-1次循环存储在内存中的场值。这样,就解决了半个时间步在程序中无法表示的问题,而且也没有破坏电磁场在时间上逐步推进的逻辑关系。2.3 时域有限差分法相关技术 2.3.1 数值稳定性问题上述FDTD方程是一种显式差分方程,在执行时,存在一个重要的问题:即算法的稳定性问题。这种不稳定性表现为在解显式方程时,随着时间步数的继续增加,计算结果也将无限制地增加。Taflove等[4]于1975年对Yee氏差分格式的稳定性进行了讨论,并导出了对时间步长的限制条件。数值解是否稳定主要取决于时间步长?t与空间步长?x、?y、?z的关系。对于非均匀媒质构成的计算空间选用如下的稳定性条件:?t?1(2-9)111v()2?()2?()2?x?y?z(2-9)式是空间和时间离散之间应当满足的关系,又称为Courant稳定性条件。若采用均匀立方体网格:?x??y??z??s?s?t?
(2-10)v3其中,v为计算空间中的电磁波的最大速度。2.3.2 数值色散FDTD方程组是对Maxwell旋度方程进行差分近似,在进行数值计算时,将会在计算网格中引起数字波模的色散,即在FDTD网格中,电磁波的相速与频率有关,电磁波的相速度随波长、传播方向及变量离散化的情况不同而改变。这种关系由非物理因素引起,且色散将导致非物理因素引起的脉冲波形畸变、人为的各向异性和虚假折射等现象。显然,色散与空间、时间的离散间隔有关,如下式所示:k?y?111112???t?2?kx?x?2?y2?kz?z?2?kz?z? (2-11) ??sin?sin?sin?sin?sin????????2??z2222c?t2?2??x2?2??y2??z??????式(2-11)是三维情况下在FDTD方法中的单色平面波数值色散关系的一般形式,它表明FDTD计算中波的传播速度与传播方向有关。式中kx、ky、kz分别是波矢量沿x、y、z方向的分量,?是角频率,c是被模拟的均匀介质中的光速。与数值色散关系相对应,在无耗介质中的单色平面波,色散解析关系是:??c?2?kx2?ky2?kz2
(2-12)由式(2-11)可知,当式(2-11)中的?t、?x、?y、?z均趋于零时,它就趋于式(2-12)。也就是说数值色散是由于用近似差分替代连续微分而引起的,而且在理论上可以减小到任意程度,只要此时时间步长和空间步长都足够小,但这将大大增加所需的计算机存储空间和计算时间,并使累积误差增加。因此,在实际计算中要根据问题的性质和计算机的软硬件条件来选择合适的时间步长和空间步长。为获得理想的色散关系,问题空间分割应按照小于正常网格的原则进行。一般选取的最大空间步长为?max??min20,?min为所研究范围内电磁波的最小波长。由上分析说明,数值色散在用FDTD法分析电磁场传播中的影响是不可能避免的,但我们可以尽可能的减小数值色散的影响。2.3.3 离散网格的确定无论是简单目标还是复杂目标,在进行FDTD离散时网格尺寸的确定,除了受计算资源的限制不可能取得很小外,还需要考虑以下几个因素:1.目标离散精确度的要求。网格应当足够小以便能精确模拟目标几何形状和电磁参数。2.FDTD方法本身的要求。主要是考虑色散误差的影响。设网格为立方体?x??y??z??,所关心频段的频率上限为fmax,对应波长为?min,则考虑FDTD的数值色散要求???minN(2-13)通常N?10。上式是根据已知所关心频率上限情况下来确定FDTD网格尺寸?的;反之,若给定?,则FDTD计算结果可用的上限频率也随之确定。3.入射波的要求。入射波的上限截止频率fc应包含所关心频率范围,即fc?fmax。2.4 吸收边界条件由时域有限差分法的基本原理可知,在利用时域有限差分法研究电磁场时,需在全部问题空间建立Yee氏网格空间,并存储每个单元网格上任一时间步的六个场分量用于下一时间步的计算。而在对于辐射、散射这类开放系统的实际研究中,不可能有无限大的存储空间。因此,必须在某处将网格空间截断,且在截断边界网格点处运用特殊的场分量计算方法,使得向边界面行进的波在边界处保持“外向行进”特性、无明显的反射现象,并且不会使内部空间的场产生畸变,从而用有限网格空间模拟电磁波在无界空间中传播的情况。具有这种功能的边界条件称之为吸收边界条件,或辐射边界条件,或网格截断条件[29~31],如图2-3所示。 图 2-3 附加截断边界使计算区域变为有限域从FDTD的基本差分方程组可以看出,在截断边界面上切向场分量的计算需要利用计算空间以外的电磁场分量,因此FDTD基本差分方程对这些截断边界面上的场分量失效。如何处理截断边界上的场分量,使之与需要考虑的无限空间有尽量小的差异,是FDTD中必须很好解决的一个重要问题。实际上,这是要求在误差可容忍的范围内,计算空间中的外向波能够顺利通过截断边界面而不引起波的明显反射,使有限计算空间的数值模拟与实际情况趋于一致,对外向波而言,就像在无限大空间中传播一样。所以,需要一种截断边界网格处的特殊计算方法,它不仅要保证边界场计算的必要精度,而且还要大大消除非物理因素引起的波反射,使得用有限的网格空间就能模拟电磁波在无限空间中的传播。但是如果处理不当,截断边界面可能造成较大反射,构成数值模拟误差的一部分,甚至可能造成算法不稳定。加于截断边界场分量符合上述要求的算法就称为吸收边界条件(Absorbing Boundary Conditions)。2.4.1 一阶和二阶近似吸收边界条件在截断边界附近通常没有激励源。考虑齐次波动方程??2?2?21?2????x2??y2??z2?c2?t2??f?0
(2-14) ??式中,f表示直角坐标系下任意电磁场分量。B.Engquist和A.Majda[15]利用偏微分算子对式(2-13)作因式分解,并分别取其Taylor级数展开式中的第一项和前两项近似,导出了适合直角坐标系下FDTD吸收边界条件的单向波动方程,这就是Engquist-Majda吸收边界条件。设三维长方体FDTD区域0&x&a,0&y&b,0&z&d,这时有六个截断边界,其一阶和二阶解 析吸收边界条件的具体形式见表2-1,其中f代表任一直角场分量。G. Mur[14]对表中的吸收边界条件引入了一种简单有效的差分数值算法,即对时间和空间的偏微分取二阶中心差分近似,将单向波方程离散化,便形成了著名的G.Mur的一阶和二阶吸收边界条件,其总体虚假反射在1%~5%。图2-5给出了反射系数(反射波与入射波)与入射角的关系,由图2-4可看出,仅当入射角较小时其反射系数较小。表2-1 三维长方体FDTD区域的一阶和二阶吸收边界条件 包含各类专业文献、高等教育、幼儿教育、小学教育、各类资格考试、应用写作文书、中学教育、时域有限差分法对平面TE波的MATLAB仿真71等内容。 
 时域有限差分法对平面TE波的MATLAB仿真_数学_自然科学_专业资料。天线与电波传播时域有限差分法对平面 TE 波的 MATLAB 仿真摘 要时域有限差分法是由有限差分法发...  时域有限差分法对平面 TE 波的 MATLAB 仿真姓名:王云璐 学号: (西北工业大学 电子信息学院 ,陕西 西安,710072) 摘要:本文分析 FDTD 算法的基本...  分) 时域有限差分方法作为一种典型的全波时域分析...利用 matlab 对一些实际的电磁场问题进行仿真研究。 ...[4] 杨臻颖.三维散射的 FDTD 模拟及平面入射波...  时域有限差分法仿真一维TE波在分裂场完全匹配层【含源码】_电子/电路_工程科技_专业资料。时域有限差分法仿真一维TE波在分裂场完全匹配层【含源码】时域...  时域有限差分法对平面TE... 35页 1下载券 时域有限差分法发展综述 7页 1下载...i , j (18) 为了编写统一的 TE 和 TM 波二维 FDTD 程序, 可将描述 TE...  二维FDTD仿真 暂无评价 3页 1下载券 基于MATLAB的二维TM波FD... 暂无评价 5...时域有限差分法对平面TE... 35页 1下载券 时域有限差分法论文 33页 免费 ...  时域有限差分法对平面TE波... 35页 1财富值 33.时域有限差分2 暂无评价 8...另外,上传者有关于导体电容,电荷分布,天线性能分析,微带线等MATLAB分析程序仿真...  时域有限差分法对平面TE... 35页 1下载券 题目16:离散系统Z域分析... 暂无...有限差分法的Matlab程序 2页 免费一维时域有限差...时域有限差分法的Matlab仿真_百度文库
两大类热门资源免费畅读
续费一年阅读会员,立省24元!
评价文档:
时域有限差分法的Matlab仿真
阅读已结束,如果下载本文需要使用
想免费下载本文?
把文档贴到Blog、BBS或个人站等:
普通尺寸(450*500pix)
较大尺寸(630*500pix)
你可能喜欢

我要回帖

更多关于 matlab时域转频域程序 的文章

 

随机推荐