带挡板的平行板传输线二维FDTD仿真其图如下：.DOC

带挡板的平行板传输线二维FDTD仿真，其图如下： 建立如图所示坐标，当时，高斯脉冲激励在t=0处近似为零，脉冲表达式为：。 用matlab编程求观测面的电压： clear all L=4;d=0.20; %基本参数设置 T=0.4e-9; c=3e8; u=4*pi*1e-7; %相对磁导率 e=1/u/c^2; %相对介电常数 %空间步长及z,x方向网格数 f=1/2/T; %频率 wl=c/f; %波长 dz=wl/20; dx=dz; nz=round(L/dz); %沿Z方向的网格数 nx=round(d/dz); %沿X方向的网格数 x1=d/3; x2=d*2/3; z1=L/2; nx1=round(x1/dx); %沿X方向上挡板的网格数 nx2=round(x2/dx); %沿X方向下挡板的网格数 nz1=round(z1/dz); %沿Z方向的网格数 dt=dz/2/c; %时间步长 nt=4*nz; %时间网格数 nz11=round(L/3/dz); %取L/3处观测面 nz12=round(2*L/3/dz); %取2L/3处观测面 ex=zeros(nx,nz+1); ez=zeros(nx+1,nz); hy=zeros(nx,nz); ex0=zeros(nx,nz+1); %生成空间零矩阵 for n=1:nt V1=0;V2=0; t=n*dt; t1=3*T; F=exp(-(t-t1)^2/T^2); %在输入口加高斯脉冲 ex(1:nx,1)=F; ez(1,1:nz)=0;ez(nx+1,1:nz)=0;%上下金属挡板初值 ex(1:nx1,nz1)=0;ex(nx2:nx,nz1)=0; %中间挡板情况 if (t6*T) ex(1:nx,1)=ex0(1:nx,2)+(c*dt-dz)/(c*dt+dz)*(ex(1:nx,2)-ex0(1:nx,1)); %一阶Mur吸收边界 ex(1:nx,nz+1)=ex0(1:nx,nz)+(c*dt-dz)/(c*dt+dz)*(ex(1:nx,nz)-ex0(1:nx,nz+1)); end %蛙跳 for i=1:nx for k=1:nz hy(i,k)=hy(i,k)+dt/u*((ez(i+1,k)-ez(i,k))/dx+(ex(i,k)-ex(i,k+1))/dz); %差分求磁场 end end for i=1:nx for k=1:nz-1 ex(i,k+1)=ex(i,k+1)+dt/e/dz*(hy(i,k)-hy(i,k+1));%差分求电场 end end for i=1:nx-1 for k=1:nz ez(i+1,k)=ez(i+1,k)+dt/e/dx*(hy(i+1,k)-hy(i,k));%差分求电场 end end ex0=ex; %画两观测面处的电压 for i=1:nx V1=V1+ex(i,nz11)*dx; V2=V2+ex(i,nz12)*dx; end subplot(2,1,1); plot(t,V1,+); hold on subplot(2,1,2); plot(t,V2,+); pause(0.01); hold on end 程序运行图形如下： ABC 观测面222 观测面111 Z X ABC L== d