利用一维FDTD方法对电磁波传播及反射透射进行仿真2_第1页
利用一维FDTD方法对电磁波传播及反射透射进行仿真2_第2页
已阅读5页,还剩3页未读 继续免费阅读

下载本文档

版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领

文档简介

1、电磁场与电磁波实验报告实验项目:一维FDTD方法模拟电磁波传播班级:集成电路姓名:张超Forpersonaluseonlyinstudyandresearch;notforcommercialuse学号:同组姓名:林彬同组学号:指导老师:汤炜实验日期:一、实验目的要求1、了解数值方法的基本原理,熟悉时域有限差分方法(FDTD)的计算思路。2、复习Matlab语言,学习编程的基本技巧和编程思路。3、加强对电磁波理论的了解,理解反射系数,透射系数等基本概念。4、形象展示电磁波的传播及与介质板的作用过程。二、实验内容利用一维FDTD方法对电磁波传播及反射透射进行仿真三、实验仪器计算机Matlab编译

2、系统四、实验原理该实验的中心思想就是利用麦克斯韦方程组来建立模型,然后根据模型编写程序,对模型进行仿真实验,通过matlab的图形仿真来实现入射波、透射波、反射波的波形波形仿真。1、一维Maxwell方程:VxH"E+£竺_dtVxE=卩翌V H=0V E=££dzdtdHy_严£2、在将时间空间进行离散化处理,其核丿心思想是将计:心思想是将计算蟲的空间和时间进行划分空间:例如:三维空间划分为立方块,二维空间划分为正方柱,一维空间划分为平面板。划分的区域非常小,以至于可以认为场量在该区域是不变的。时间:将电磁波的与目标的作用时间划分为很多时间小

3、段,可以认为场量在该时间段内是不变的。时空的标定:空间的划分长度为As,维情况下用kAs表示每个场点的空间位置,并简记为k。例如:E(k)=E(kAs)表示k位置的电场时间的划分长度为At,利用nAt表示某个时刻,并简记为n。(书写时写在上标位置)空间序列:H(k1.5)yH(k.5)yAsH(k+.5)yH(k+1.5)yPE(k2)x时间序列OE(k-t)x0E(k+1)xHn1.5yAtHn+.5En2xEn1xEnxEn+1x表示A的x分量在(n+05)At时刻、(k+05)As位置的值3、一维FDTD方法中的离散原则:电场:空间位置位于整空间步长,时间位于整时间步长。磁场:空间位置位

4、于半空间步长,时间位于半时间步长。不得用于商业用途-空En(k+1)EnpAsxx故综合表示,电场和磁场分别可表示型如:4、麦克斯韦思维方程组与差分方程的结合:Hn+.5(k+0.5)=Hn-.5(k+0.5)空En(k+1)EnyypAsxxHn+.5(k+0.5)=Hn-.5(k+0.5)yy同理根据另一Maxwell方程得:En+1(k)=En(krhn.5(k+0.5)-Hn.5(k-0.5xx£心Lyy两个迭代方程中右边:后时刻场量左边:前时刻场量即如果能够得到前一时刻的电场和磁场,根据方程即可得到后一时刻的场量。5、空间步长和时间步长的设定:原则上说,空间步长和时间步长越

5、小越好。实际上,太小的步长会导致计算速度过慢,内存占用较多。常用的设定为:其中:为计算区域中的最小波长在以上约束条件下,迭代方程可写为yyHn+.5(k+0.5)-Hn+.5(k-0.5)En+1(k)=En(k)-Txx£r6、计算模型:模型说明:介质板厚度为9cm,相对介电常数为4。待求问题:介质板的反射系数和透射系数(05GHz)£9cm7、建模根据计算频段与媒质特性可知,媒质中的最小波长人=3cm根据前面描述的时间步长和空间步长的约束,可以得到:血二1.5加A/二2.5閃以上为该实验的原理及模型建立的方法步骤,除此之外,试验中还应注意初始值的设立,否则也得不到实验预

6、期的结果。五、实验结果实验用的MATLAB程序:clc;clear;k0=930;w=40;n_time=8100;kk0=1000;n_k=2000;cc=0.5;i0=990;订=1100;N=2A14;dt=2.5*10人(-12);f0=5*10A9;Ex=zeros(1,N);Hy=zeros(1,N);E_inc=zeros(1,N);E_ref=zeros(1,N);E_tra=zeros(1,N);forii=1:2k=1:1:n_k+1;Ex(k)=exp(log(0.001)*(k-k0).*(k-k0)/w/w);Hy(k)=exp(log(0.001)*(k-k0+0.

7、25).*(k-k0+0.25)/w/w);fori_time=1:1:n_timeeps(k)=1;ifii=2;tmp1=kk0:1:kk0+59;eps(tmp1)=4;endforkk=2:1:n_kEx(kk)=Ex(kk)-cc*(Hy(kk)-Hy(kk-1)./(eps(kk-1)+eps(kk)/2);endforkk=1:1:n_k-1Hy(kk)=Hy(kk)-cc*(Ex(kk+1)-Ex(kk);endifii=1E_inc(i_time)=Ex(i0);elseE_ref(i_time)=Ex(i0);E_tra(i_time)=Ex(i1);enddrawnow;

8、pause(0.000005);plot(Ex,'r-');holdon;plot(1000,1000,1,-0.4,1059,1059,1,-0.4);holdoff;axis(0,2000,-0.4,1);end仅供个人参考由于n_time=1000至到5000的入射波波形及其频谱图几乎没变化,而'当n-time变化到7000和8000时波形及其频谱图有了很大的3.入射波与放射波混合波特图依次为n_time=1000,2000,4000,5000,7000,80000.80.60.4-20idSOGO1001201401&01B02004反射波波形及其频谱图不得用于商业用途仅供个人参考仅供个人用于学习、研究;不得用于商业用途。Forpersonaluseonlyinstudyandresearch;notforcommercialuse.NurfurdenpersdnlichenfurStudien,Forschung,zukommerziellenZweckenverwendetwerden.Pourl'etudeetlarechercheuniquementadesfinspersonnelles

温馨提示

  • 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

评论

0/150

提交评论