利用一维FDTD方法对电磁波传播及反射透射进行仿真2

电磁场与电磁波实验报告实验项目：一维FDTD方法模拟电磁波传播班 级： 集成电路 姓名： 学 号： 1015251041 同组姓名： 林彬同组学号： 1015251017 指导老师： 汤炜 实验日期： 一、实验目的要求1、 了解数值方法的基本原理，熟悉时域有限差分方法(FDTD)的计算思路。2、 复习Matlab语言，学习编程的基本技巧和编程思路。3、 加强对电磁波理论的了解，理解反射系数，透射系数等基本概念。4、形象展示电磁波的传播及与介质板的作用过程。二、实验内容利用一维FDTD方法对电磁波传播及反射透射进行仿真三、实验仪器计算机Matlab编译系统四、实验原理该实验的中心思想就是利用麦克斯韦方程组来建立模型，然后根据模型编写程序，对模型进行仿真实验，通过matlab的图形仿真来实现入射波、透射波、反射波的波形波形仿真。1、一维Maxwell方程:VxH=oE+8竺dzdidzdl_dtdzdidzdlVxE=W半V・H=0v・E=?2、在将时间空间进行离散化处理其核心思想是将计算区域的空间和时间进行划分空间：例如：三维空间划分为立方块，二维空间划分为正方柱，一维空间划分为平面板。划分的区域非常小，以至于可以认为场量在该区域是不变的。时间：将电磁波的与目标的作用时间划分为很多时间小段，可以认为场量在该时间段内是不变的。时空的标定：空间的划分长度为As，一维情况下用kAs表示每个场点的空间位置，并简记为k。例如:E(k)=E(kAs)表示k位置的电场4+°*+0.5)时间的划分长度为At，利用nAt表示某个时刻，并简记为n。(书写时写在上标位置)表示A的x分量在(n+0.5)At时刻、(k+0・5)As位置的值4+°*+0.5)3、一维FDTD方法中的离散原则:电场：空间位置位于整空间步长，时间位于整时间步长。磁场：空间位置位于半空间步长，时间位于半时间步长。空间序列： H(k5)H(k-15) ，△， H(k+.5) H(k+1.5)TOC\o"1-5"\h\zy y y zHD • O • O • O • O>E(k-2) E(k-1) E(k) E(k+1)X时间序列： '人X xHn-1.5 也y -^Hn+.5-n • o • o • o • tEn-2 En-1 Hn-.5 En E〃+1x x y x xf A \故综合表示，电场和磁场分别可表示型如：玲(左)| IW+0.5)[En([En(k+1)-En(k)]Hn+.5(k+0.5)=Hn-.5(k+0.5)一-^Hn+.5(k+0.5)=Hn-.5(k+0.5)y y同理根据另一Maxwell方程得：\En+1(k)=En(k)-~^S[Hn+.5(k+0.5)-Hn+.5(k-0.5)]l两个迭代方程中右边：后时刻场量左边：前时刻场量即如果能够得到前一时刻的电场和磁场，根据方程即可得到后一时刻的场量。5、空间步长和时间步长的设定：原则上说，空间步长和时间步长越小越好。实际上，太小的步长会导致计算速度过慢，内

常用的设定为:存占用较多。常用的设定为:Av=—2。其中： 为计算区域中的最小波长在以上约束条件下，迭代方程可写为:7、建模根据计算频段与媒质特性可知，媒质中的最小波长E根据前面描述的时间步长和空间步长的约束，可以得到：公二1.5砌 =以上为该实验的原理及模型建立的方法步骤，除此之外，试验中还应注意初始值的设立，否则也得不到实验预期的结果。五、实验结果实验用的MATLAB程序：clc;clear;k0=930;w=40;n_time=8100;kk0=1000;n_k=2000;cc=0.5;i0=990;i1=1100;N=2A14;dt=2.5*10A(-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.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;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]);end1、空间电场分布（迭代步数为：）注迭代步数为n_timen_time=1000n_time=2000n_time=4000

n_time=5000n_time=7000n_time=80002、入射波波形及其频谱图n_time=1000n_time=7000n_time=8000由于n_time=1000到5000的入射波波形及其频谱图几乎没变化，而当n-time变