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

这是二维的程序，上下两个平行板之间有两个长度为的d/3的挡板，位于L/2处。 o Hrx$>W]  
function FDTDtwodimension(L,d,T) k_1oj[O  
%version1.0 caizhongheng 终端短路 (Jk:Qz5  
%FDTDtwodimension(6,0.18,0.5e-9) yJw4!A 1!  
t0=3*T; op
|x~Thf  
c=3e8; !}HT
&N8[r  
u=4*pi*1e-7; LdyE*u_  
e=8.8541878e-12; h3;RVtS  
dz=T*c/10; ,tuZ_"?M  
Nz=round(L/dz); G[!<mh4h|  
dt=dz/4/c; a0Q\]S  
dx=dz; k)JwCt.%  
Nx=round(d/dx); dZ\T@9+j+  
Ex=zeros(Nx,Nz+1); \>4x7mF!  
Hy=zeros(Nx,Nz); rE\.[mFI  
Ez=zeros(Nx+1,Nz); ~<[]l~`  
Nt=4*Nz; hj[sxC>z5  
6dYUMqQ  
@m"P_1`*  
for n=0:Nt |R(rb-v  
%       初始脉冲加边界条件 sUsIu,1Q  
    t=n*dt; !Ly1!;<  
    F=exp(-(t-t0).^2./T^2); XPcx"zv\  
    Ex(1:Nx,1)=F;
.vXe}%  
    if t>2*t0 a4N8zDS  
        Ex(1:Nx,1)=0; RTTEAh:.  
    end I^S{V^Ty  
%   终端短路 ^Xsl
j  
    Ex(1:Nx,Nz+1)=0; 6}PoBhgSg-  
   H$]FUv8  
%   蛙跳格式 YQ5d!a.  
    for k=1:Nz uB35CRd  
        for i=1:Nx }pv<<7}|  
            Hy(i,k)=Hy(i,k)+dt./u.*((Ez(i+1,k)-Ez(i,k))./dx+(Ex(i,k)-Ex(i,k+1))./dz); 3<HPZWc  
        end ">S.~'ds  
    end H/O
v8|  
    for k=2:Nz-1 ) 2Hl\"F  
        for i=1:Nx CB?,[#r5f  
            Ex(i,k)=Ex(i,k)+dt./e.*(Hy(i,k-1)-(Hy(i,k)))./dz; V$ac}A,!  
        end tNCKL.yU  
    end SL>0_
 
%   挡板所处的边界条件 mKBPIQ+ZS  
    Ex(1:round(Nx/3),round(Nz/2))=0; u1(8a%ZC  
    Ex(round(2*Nx/3):Nx,round(Nz/2))=0; P3V}cGZ  
   !,PoH  
    for k=1:Nz pw1&WP&?3  
        for i=2:Nx _1P8rc"Dx  
            Ez(i,k)=Ez(i,k)+dt./e.*(Hy(i,k)-Hy(i-1,k))./dx; Di #Em[  
        end
Z6oA>D  
    end 5[$jrG\!  
%   上下平行板的边界条件 WK6|e[iP  
    Ez(1,1:Nz)=0; q=cH ^`<.  
    Ez(Nx+1,1:Nz)=0; )rv<"  
   -44''w?z  
    if mod(n,10)==0 o^~KAB7  
        subplot(221); AN-;*n<'  
        mesh(Ex); %ck/ Z  
        subplot(222); 6w^Fee`>]  
  &n .. +jX.::UPm  
l%$co07cX  

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