[转载]二维FDTD算法的matlab程序

这是二维的程序，上下两个平行板之间有两个长度为的d/3的挡板，位于L/2处。 F 09DV<j  
function FDTDtwodimension(L,d,T) w6dFb6~R  
%version1.0 caizhongheng 终端短路 raVA?|'g~  
%FDTDtwodimension(6,0.18,0.5e-9) %owsBO+  
t0=3*T; e pCLM_yA  
c=3e8; I V%VU  
u=4*pi*1e-7; Wjh/M&,  
e=8.8541878e-12; ajRSMcKb7i  
dz=T*c/10; Z}A%=Z\/3  
Nz=round(L/dz); mV73 \P6K  
dt=dz/4/c; ./j,
Z$|  
dx=dz; {K{EOB_u  
Nx=round(d/dx); ZlYPoOq  
Ex=zeros(Nx,Nz+1); @4(k(  
Hy=zeros(Nx,Nz); ;|p$\26S)%  
Ez=zeros(Nx+1,Nz); Hd7,ZHj3^  
Nt=4*Nz; Ch()P.n?  
[OTJVpC  
Sw`RBN[ yo  
for n=0:Nt _sE#)@p  
%       初始脉冲加边界条件 I\?9+3 XnQ  
    t=n*dt; /WV7gO&L1  
    F=exp(-(t-t0).^2./T^2); <WXzh5D2  
    Ex(1:Nx,1)=F; (C] SH\  
    if t>2*t0 Ll4bdz,  
        Ex(1:Nx,1)=0; 8l?piig#  
    end _{vkX<s  
%   终端短路 +QM@VQ  
    Ex(1:Nx,Nz+1)=0; d8SE,A&  
   p47S^gW  
%   蛙跳格式 SjjI
r ^  
    for k=1:Nz iGDLZE+
?  
        for i=1:Nx 'evj,zFhW  
            Hy(i,k)=Hy(i,k)+dt./u.*((Ez(i+1,k)-Ez(i,k))./dx+(Ex(i,k)-Ex(i,k+1))./dz); }ZSQ>8a  
        end dUgrKDNyA  
    end MC((M,3L  
    for k=2:Nz-1 G'zF)0oD  
        for i=1:Nx R8L_J6Kpa  
            Ex(i,k)=Ex(i,k)+dt./e.*(Hy(i,k-1)-(Hy(i,k)))./dz; /:6Wzj  
        end uV-'~8  
    end Kz<@x`0   
%   挡板所处的边界条件
.O*bILU  
    Ex(1:round(Nx/3),round(Nz/2))=0; CSbI85F  
    Ex(round(2*Nx/3):Nx,round(Nz/2))=0; .I VlEG0  
   9jp:k><\(c  
    for k=1:Nz @\oz4^  
        for i=2:Nx WD;Y~|  
            Ez(i,k)=Ez(i,k)+dt./e.*(Hy(i,k)-Hy(i-1,k))./dx; 0c.s -  
        end 0U/K7sZ  
    end ]Fvm 7V  
%   上下平行板的边界条件 b96%")  
    Ez(1,1:Nz)=0; Bx"7%[  
    Ez(Nx+1,1:Nz)=0; VN0mDh?E  
   to8X=80-3  
    if mod(n,10)==0 YI-O{U  
        subplot(221); i`/+,<  
        mesh(Ex); Iq+>qX  
        subplot(222); .,~(%#Wl$  
  &n .. 
UY2X  
fl Jp4-nx  

未注册仅能浏览部分内容，查看全部内容及附件请先 登录 或 注册