FDTD分析报告1

1、电磁场与电磁波课程设计报告课设题目：基于时域有限差分法(FDTD的矩形谐振腔分析学院：信息与电气工程学院专业：电磁场与无线技术班级：1302701姓名：解宇亮学号：130270115指导教师：周洪娟哈尔滨工业大学(威海)2016年6月3日一.设计任务采用FDT躁值计算的方法来分析理想谐振腔中的场，谐振腔尺寸为25*12.5*60mm填充空气,采用直角坐标系下的场分量迭代公式，激励源采用高斯脉冲源，源的参数根据谐振腔的尺寸来确定。分析时间和空间离散度以及采样点数对分析结果的影响。二.方案设计(1)学习FDTD0论，并推导直角坐标系下Maxwell方程的差分方程；(2)理论学习并推导理想矩形谐振腔

2、中的时谐场，并分析其谐振频率分布；(3)激励源采用高斯脉冲源，导体采用PEC边界，利用FDTD编程求解谐振腔内的场分量；(4)对谐振腔内部分点处的采样数据进行频谱分析，提取其谐振频率分布，并与理论对比，并分析时间和空间离散度以及采样点数对分析结果的影响。三.设计原理3.1矩形谐振腔谐振频率分析当扰动频率恰使腔内的平均电能和平均磁能相等时便发生谐振， 等于某种模式的1/2波导波长整数倍时，该模式发生谐振，称为谐振模在直角坐标系中可写作六个标量方程这个频率称为谐振频率。腔长:y一z正 x：:Ez=1:z；x-Ey；：Ex.=-1.:HzZ.:Hx:z,：Hy；x-y：Hy；z.:t-：Hyjt；:

3、HzFt-；1mHy0mHz.:Exx.二Ex,:Hz.x汩 x:yFEyZLFt.二Ey二Ez(2)在矩形谐振腔中，TM 的场结构l,、l.,m二、,n二、,p二、Ez(x,y,z)=Emsin(x)sin(y)cos(z)ablEx(x,y,z)=,m二、zm-：、.,T()Emcos(x)sin(caay)sin(p-：z)n、Ey(x,y,z)=-r()Emsin(kcb、加二、一x)cos(y)sin(bz)Hz(x,y,z)=0(6)j；n_Hx(x,y,z)7()Emsin(kcb、，n二、,x)cos(y)cos(bz)Hy(x,y,z)=jm二m二j)EmCos(x)sin(

4、kcaa、Q：、y)cos(-z)l(8)谐振频率fmnpmnp_kmnp2(m)2(n)22a2bTE波推导出的谐振频率结果与TM波一致。3.2不同激励方向下的谐振频率表 1 1 不同激励方向下的谐振频率ExEyEz谐振频率12.25GHz6.5GHz13.42GHz3.3时域有限差分法分析图 1 1YeeYee 元胞的交替采样FDTD采用E、H分量的节点在空间上的交替排列和在时间上的交替抽样方式，从而可以在时问轴上逐次推进的求解空问电磁场的值。3.4FDTD在直角坐标系下的迭代公式直角坐标系下电场迭代公式E；1i；,j,k=CAmE：i2,j,kCBmn4iin1iin1iin1iiHz2

5、(i+，j+.,k)Hz2(i+复，j/)_Hy2(i+-,j,k+)-Hy2(i+/，j,k一)(必国豆LJ式中m=i1,j,kE；*(i,j+1,k)=CA(m)E；(i,j+k)+CB(m)In411n411nnHx2(i,j+；,k+白Hx2(i,j+2,k；)Hz2(iY,j+k)Hz2(iT,k)zL一(ii)式中m=i,j；,kEn1i,j,k1=CAmEni,j,k；CBm1111nF11n711n力11n差11Hy2(i+2,j,k+；)-Hy2(i-2,j,k+；)Hx2(i,j+；,k+1)-Hxi,j；,k+；)式中m=i,j,k2t直角坐标系下磁场迭代公式(12)CA

6、m=CBm=m；m1-mt.：t2t2；m=11.二mn一2；m:t;m：t(13)(14)式中式中式中n.1Hx2i,j1,k1=CPmHx2i,j-CQm11=i,J2,k2n：1Hy2iT,j,kE；(i,j+1,kFgk；_n1_n1_n)-Ez(i,j,k十1Eyi,j+,k+1)-Eyi,jLZ(15)1nF112)=CP(m)Hy2(i+j,k十1-CQgFn111Hz2i1,j；,k+1)-E：(i+；,j,k)EZ(i+1,j,k+)EZ(i,j,k+J(16)1n-k=CPmHz2i-CQmz11，jJkg。+1,j+Jk)-E；(i,j+kE；iJj1,k-E：i-j,k

7、(17)CPmt二 mm；：t2m二 mm，t2-；:：m(18)3.5FDTD数值稳定性分析在FDTN,时间增量却和空间增量Ax、Ay、Az之间不是相互独立的，它们的取值必须满足一定的关系，以避免数值结果的不稳定一一表现为随着时间步数的增加，计算结果发散。1、时间步长稳定性要求：t=-(20)71其中T为波动周期。2、时间步长与空间步长的关系在非均匀区域，v取最大值。真空中v=c(光速)。若是立方体Yee元胞，6=Ax=Ay=Az,那若电磁波所在空间的介质特性与频率有关，则电磁波的传播速度也将是频率的函数，这种现象称色散。而FDTD方程是原Maxwell方程的一种近似，所以当计算机对电磁波在

8、空间的传播进行模拟时，在非色散空问中也会出现色散现象，且电磁波的相速度随波长、传播方向及变量离散化的情况发生变化，这种非物理性的色散称为数值色散。数值色散会导致脉冲波形的破坏，出现人为的各向异性及虚假的折射现象。数值色散是由于近似差商替代连续微商引起的，这种影响可以通过减小离散化过程所取空间和时间步长而无限减小， 但计算空间的总网格数目的增加也会相应增加对计算机存储空间和计算时间的要求，所以要根据实际条件来选择合适的步长。为了减小数值色散，可选择更高的要求3.5谐振腔中激励源的选择设置激励源高斯脉冲的表达式为)(22)为了在谐振腔中激励起模， 并且抑制其他高次模， 选择线源脉冲， 使之在腔内x

9、z平面中心处沿y轴方向分布，并选择合适的和 p 值。经过反复试验，取=138.462ps,r=46.154ps图 2 2 高斯源的时域与频域分布:t1CQm:m二 mm.:t21二 mm.:tm2m(19)(21)四、实验结果分析.ey方向激励下的场结构图 3 3eyey 方向激励下的电场结构图4 4eyey 方向激励下的磁场结构上图显示了在谐振腔内部由高斯脉冲源激励，由谐振腔自身参数选择出的谐振频率，在(12.5,6.25,12.5)mm处所包含的频率分量。容易发现谐振发生在频率6.6Ghz处，与理论分析结果相差0.1GHz通过计算可知，此条彳下的频率分辨率为0.06GHz,实验结果在此频率

10、分辨率条件下与理论值相差不大。波形图中的其他波峰是由于高斯激励脉冲元不是理想脉冲所引起的，同时在9.6GHz产生谐振。.其他方向激励的电场分布图5 5e ez z方向激励下的电场结构图6 6e ex x方向激励下的电场结构由图所示，ez激励方向下的谐振频率为13.2GHz与理论分析结果相差0.26GHzex激励方向下的谐振频率为12GHz与理论分析结果相差0.26GHz不同方向激励下的谐振频率与理论值相差不大，同时也会产生高次模的谐振。在二次谐振后其余频率分量均为0.空间离散度对电场结构影响图 7 7 空间离散度对电场结构影响上图显示了空间离散度对场结构的影响，由结果对比可知，满足稳定性条件要

11、求下，空间离散度越小，谐振频率分辨率越高，高模抑制效果越好。.时间离散度对电场结构影响图 8 8 时间离散度对电场结构影响上图显示了时间离散度对场结构分布的影响，在满足时间离散度要求下，dt越大观察时间越长，频域的分辨率越高。.时间步长对电场结构影响图 9 9 时间步长对电场结构影响采样有效时间越长，即在相同的采样频率(时间步长dt的倒数)下，采样点数nmax增力口。那么频域中一个周期内频率的间隔fo(1/dt/nmax)减小，对两个频率相近序列的频谱峰值分辨能力增强，谱峰更尖锐更细，频率分辨力越好，时域波形密集。五、结论麦克斯韦方程组是支配宏观电磁现象的一组基本方程，FDTD方法是从微分形式

12、的麦克斯韦旋度方程出发进行差分离散从而得到一组时域推进公式。使用matlab对已知尺寸填充空气理想导体边界的矩形谐振腔进行三维的maxwell旋度方程分析，使用高斯脉冲激励源激励，验证了物体的电尺寸决定物体的电磁谐振频率（固有频率）和腔内存在的主模。另外，在离散空间中，当时间和空间步长确定后，随着采样点数的增加，幅频响应途中频率分辨率变高。六、参考文献.电磁场与电磁波，Matlab相关书籍等.葛德彪，闫玉波，电磁波时域有限差分方法，西安电子科技大学出版.基于FDTD法的电磁仿真的优化算法谭文泉（南京工业大学信息科学与工程学院，江苏南京210009）七、附录functiondx=DX1(inpu

13、t1,input2,input3)%UNTITLEDSummaryofthisfunctiongoeshere%Detailedexplanationgoesherec=3e8;muz=4.0*pi*1.0e-7;epsz=1.0/(cA2*muz);dx=input1;X=25e-3;Y=12.5e-3;Z=60e-3;ie=X/dx;je=Y/dx;ke=Z/dx;ib=ie+1;jb=je+1;kb=ke+1;%源及观察的位置is=ie/2+1;js=je/2+1;ks=ke/2+1;is=ie/2+1;js=je/2+1;ks=ke/2+1;却寸间步长dt=dx/(2*c);f0=13

14、.42e9;%总步数nmax=1000；%选择f=1;nmax=nmax(f);利水冲的系数rtau=0.3/f0;%决定脉宽的时间长度tau=rtau/dt;%3dt化ndelay=tau*3;%寸延%介质的参数%选择腔内任意一点进行验证i=10;j=5;k=10;%介质的参数epsr=1.0;sig=0.0;mur=1.0;sigm=0.0;%迭代系数ca=(1.0-(dt*sig)/(2.0*epsz*epsr)/(1.0+(dt*sig)/(2.0*epsz*epsr);cb=(dt/epsz/epsr/dx)/(1.0+(dt*sig)/(2.0*epsz*epsr);da=(1.0

15、-(dt*sigm)/(2.0*muz*mur)/(1.0+(dt*sigm)/(2.0*muz*mur);db=(dt/muz/mur/dx)/(1.0+(dt*sigm)/(2.0*muz*mur);%初始化ex=zeros(ie,jb,kb);ey=zeros(ib,je,kb);ez=zeros(ib,jb,ke);hx=zeros(ib,je,ke);hy=zeros(ie,jb,ke);hz=zeros(ie,je,kb);forn=1:nmaxex(1:ie,2:je,2:ke)=ca*ex(1:ie,2:je,2:ke)+cb*(hz(1:ie,2:je,2:ke)-hz(1:

16、ie,1:je-1,2:ke)+hy(1:ie,2:je,1:ke-1)-hy(1:ie,2:je,2:ke);ey(2:ie,1:je,2:ke)=ca*ey(2:ie,1:je,2:ke)+cb*(hx(2:ie,1:je,2:ke)-hx(2:ie,1:je,1:ke-1)+hz(1:ie-1,1:je,2:ke)-hz(2:ie,1:je,2:ke);ez(2:ie,2:je,1:ke)=ca*ez(2:ie,2:je,1:ke)+cb*(hx(2:ie,1:je-1,1:ke)-hx(2:ie,2:je,1:ke)+hy(2:ie,2:je,1:ke)-hy(1:ie-1,2:je,

17、1:ke);%高斯线源脉冲，电流源沿着y方向激励ey(is,1:je,ks)=ey(is,1:je,ks)+8*exp(-(n-ndelay)A2/tauA2);hx(2:ie,1:je,1:ke)=da*hx(2:ie,1:je,1:ke)+db*(ey(2:ie,1:je,2:kb)-ey(2:ie,1:je,1:ke)+ez(2:ie,1:je,1:ke)-ez(2:ie,2:jb,1:ke);hy(1:ie,2:je,1:ke)=da*hy(1:ie,2:je,1:ke)+db*(ex(1:ie,2:je,1:ke)-ex(1:ie,2:je,2:kb)+ez(2:ib,2:je,1:ke)-ez(1:ie,2:je,1:ke);hz(1:ie,1:je,2:ke)=da*hz(1:ie,1:je,2:ke)+db*(ex(1:ie,2:jb,2:ke)-ex(1:ie,1:je,2:ke)+ey(1:ie,1:je,2:ke)-ey(2:ib,1:je,2:ke);%等某点处的值取出来o(1,n)=ex(i,j,k);o(2,n)=ey(i,j,