[RFEDA原创]新手学习FDTD，自己编写的程序

 本人是一新手，混迹于本论坛，从各位老大的帖子里面学到很多。 8!35 K  
刚刚开始学习FDTD，这个是我自己编写的程序，希望可以对像我这样的新手有一定启发哈。 +mBS&FK  
  1.@{5f3T  
程序题目要求见附件 L *\[;.mk  
'V%w{ZiiV  
clear  %设置参数 GmNCw5F  
mu=4*pi*10^(-7);%真空中的磁导率 #m{(aa9;  
eps=8.8542*10^(-12);% 真空介电常数 .*6Nq
X$  
T=20*10^(-12); @BLB.=  
t0=4*T; @nwVl8  
c=3*10^8; )0E_Y@  
f=1/T;%频率 T\Xf0|y  
longth=c/f;%波长 -cUbIbW  
ds=longth/15;%网格分辨率 ~6t<`&f  
dt=ds/(2*c);%时间步长 tYTl-c  
%吸收边界条件迭代方程中的系数 Qa/1*Mb  
a=(c*dt-ds)/(c*dt+ds); MPnMLUB$\  
aa=(c^2*eps*dt)/(2*(c*dt+ds)); KsIHJr7-  
%角点迭代方程中的系数 0.+Eo.AX4M  
bb=(c*dt-sqrt(2)*ds)/(c*dt+sqrt(2)*ds); 8Ckd.HKpQ  
nt=500; +a,#
BSt  
nx=200; dpE^BWv3  
ny=nx; {.oz^~zs]g  
   ^! h3#4  
%---------------------------------------------------初始化--------------------------------------------------- q
Pz
_PRje  
        Ex=zeros(nx,ny); W2-l_{  
        Ey=zeros(nx,ny); R%r25_8  
        Hz=zeros(nx,ny); 00IW9B-  
        Hz(100,100)=exp(-16); 4P?`<K'  
        %cDGs^lgA  
        pvl];w  
        dz=zeros(nx,ny);% n-1时刻的Hz .n_Z0&i/w  
        ex=zeros(nx,ny);% n-1/2时刻的Ex 6@lZVM)E  
        ey=zeros(nx,ny);% n-1/2时刻的Ey -}4CY\d6'  
    v(0ujfSR0  
        LsI@_,XW<  
  cg^~P-i@*  
%------------------------------------------------电场和磁场迭代---------------------------------------------------% " ]k}V2l  
for t=1:nt  %时间步数 2?LPr  
   XBX`L"0  
  %磁场迭代 =/j!S|P  
    whe%o  
          QsKnaRT  
          bHm/ZZx  
           Ex(1:nx,1:ny-1)=ex(1:nx,1:ny-1)+(dt/(eps*ds))*(Hz(1:nx,2:ny)-Hz(1:nx,1:ny-1)); Chl^LEN:  
       !oi {8X@  
     o_=4Ex "  
           Ey(1:nx-1,1:ny)=ey(1:nx-1,1:ny)-(dt/(eps*ds))*(Hz(2:nx,1:ny)-Hz(1:nx-1,1:ny)); -wH0g^Ed  
           >q?{'#i /  
           pulse=exp(-(((t*dt-t0)/T)^2));%激励源高斯脉冲 bbS,pid1  
 %电场迭代 sa<\nH$_X  
  ;=@O.iF;H  
        7)}_'p  
      for i=2:nx AW\uE[kg  
        for j=2:ny-1 _l] 0V g`  
            if (i==nx*1/2&j==ny*1/2) xX<T5Ls  
            Hz(i,j)=dz(i,j)+(dt/(mu*ds))*(Ex(i,j)-Ex(i,j-1)-Ey(i,j)+Ey(i-1,j))+pulse; zQO 1%g  
            else ;uc3_J]  
            Hz(i,j)=dz(i,j)+(dt/(mu*ds))*(Ex(i,j)-Ex(i,j-1)-Ey(i,j)+Ey(i-1,j)); 28M^F~0  
            end Muq~p~m}  
        end )%p46(]  
           kz}Bc
F  
      end  )$1j"mV  
 %---------------------------------------------------设置Mur的吸收边界条件---------- j0oto
6z~b  
 %------------------------------% 8[,R4@  
  UxB3/!<5g3  
     for j=2:ny-1 lmUCrs37  
        Hz(1,j)=dz(2,j)+a*(Hz(2,j)-dz(1,j))+aa*(ex(1,j)-ex(1,j-1)+ex(2,j)-ex(2,j-1)); ,_V/W'  
    end %OHWGac"i  
 %边界x=nx I+W,%)vb  
    for j=2:ny-1 vq9O|E3  
       Hz(nx,j)=dz(nx-1,j)+a*(Hz(nx-1,j)-dz(nx,j))+aa*(ex(nx-1,j)-ex(nx-1,j-1)+ex(nx,j)-ex(nx,j-1)); ;2`t0#J$]  
    end x ]}'H  
  %边界y=1  zN5};e}^v  
    for i=2:nx-1 {4SaSv^/  
        Hz(i,1)=dz(i,2)+a*(Hz(i,2)-dz(i,1))-aa*(ey(i,1)+ey(i,2)-ey(i-1,2)-ey(i-1,1)); 0iC5,  
    end dm_Pz\*  
  %边界y=ny r/':^Ex  
    for i=2:nx-1 _IEbRVpb  
        Hz(i,ny)=dz(i,ny-1)+a*(Hz(i,ny-1)-dz(i,ny))-aa*(ey(i,ny-1)-ey(i-1,ny-1)+ey(i,ny)-ey(i-1,ny)); ?Qd`Vlp7  
    end 9 4bDJy1  
  p;7wH\c  
    6NyUGGRq  
  %---------------------------------------------------------角点处理------------------------------------------------------------% E0HE@pqr  
    Hz(1,1)=dz(2,2)+bb*(Hz(2,2)-dz(1,1));%角点x=1,y=1 "+{2!  
   `=V1w4J  
    Hz(nx,1)=dz(nx-1,2)+bb*(Hz(nx-1,2)-dz(nx,1));%角点x=nx,y=1 [vuikJP>1k  
     /EQ^-4yr  
    Hz(1,ny)=dz(2,ny-1)+bb*(Hz(2,ny-1)-dz(1,ny) .. ro?.w  
/W f.Gt9[  

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

不错，支持

呵呵，多发帖子，多交流，进步会快很多的。

这个的确要支持下

同感……交流才能共同进步。

恩，贵在坚持啊，：）

不错，支持下

程序只是手段，结果才是目的。加油。

相互学习 共同进步啊  

赞扬~ R7O<>kt  
q_98=fyE6  
学习中~