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

 本人是一新手，混迹于本论坛，从各位老大的帖子里面学到很多。 t|?eNKVV9'  
刚刚开始学习FDTD，这个是我自己编写的程序，希望可以对像我这样的新手有一定启发哈。 d=eIsP'h  
  FSD~Q&9&  
程序题目要求见附件 ((AIrE>Rr  
;:1o|>mX
 
clear  %设置参数 |crm{]7X  
mu=4*pi*10^(-7);%真空中的磁导率 ?)O!(=6%'  
eps=8.8542*10^(-12);% 真空介电常数 Y6RbRcJw  
T=20*10^(-12); N1yx|g:  
t0=4*T; [=]LR9c4  
c=3*10^8; ;iO5 8S3  
f=1/T;%频率 LhC
wZ1  
longth=c/f;%波长 BFswqp:  
ds=longth/15;%网格分辨率 h|%a}])G)  
dt=ds/(2*c);%时间步长 e]*@|e
4b  
%吸收边界条件迭代方程中的系数 E RjMe'q4  
a=(c*dt-ds)/(c*dt+ds); ?+G /5,e  
aa=(c^2*eps*dt)/(2*(c*dt+ds)); ZtGtJV"H  
%角点迭代方程中的系数 '1"vwXJ"  
bb=(c*dt-sqrt(2)*ds)/(c*dt+sqrt(2)*ds); >Vph_98|  
nt=500; !:3X{)4  
nx=200; 821;;]H  
ny=nx; C;I:?4  
   Oh5aJ)"D  
%---------------------------------------------------初始化--------------------------------------------------- j@.^3:  
        Ex=zeros(nx,ny); c&'5r OY~  
        Ey=zeros(nx,ny); 61Wh %8-  
        Hz=zeros(nx,ny); #<DS-^W!  
        Hz(100,100)=exp(-16); "0aJE1)p:  
        D9JHx+Xf>  
        #CRd@k?  
        dz=zeros(nx,ny);% n-1时刻的Hz yYiu
69v  
        ex=zeros(nx,ny);% n-1/2时刻的Ex u$<>8aMei  
        ey=zeros(nx,ny);% n-1/2时刻的Ey k!py*noy  
    &3f^]n!@  
        ;T :]?5W!  
  4SJb\R)XK  
%------------------------------------------------电场和磁场迭代---------------------------------------------------% |_mN:(3  
for t=1:nt  %时间步数 o&MOcy D  
   2" u,f  
  %磁场迭代 R1~wzy  
    nLY(%):(P  
          B^@X1EE  
          *~
kHH
 
           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)); f~:wI
9  
       ]([^(&2  
     6,~]2H'zq  
           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)); ?izl#?  
           |3G;Rh9w,  
           pulse=exp(-(((t*dt-t0)/T)^2));%激励源高斯脉冲 u;9a/RI  
 %电场迭代 ~EM(*k._  
  +1nzyD_E  
        n;LjKE  
      for i=2:nx |ZG0E  
        for j=2:ny-1 LRqlK\  
            if (i==nx*1/2&j==ny*1/2) .'bhRQY  
            Hz(i,j)=dz(i,j)+(dt/(mu*ds))*(Ex(i,j)-Ex(i,j-1)-Ey(i,j)+Ey(i-1,j))+pulse; "t%Jj89a\  
            else nL+y"O  
            Hz(i,j)=dz(i,j)+(dt/(mu*ds))*(Ex(i,j)-Ex(i,j-1)-Ey(i,j)+Ey(i-1,j)); C."\ a_p  
            end d}
<-G.&_  
        end +!mEP>  
           ?N=`
}}Ky-  
      end  ;r}yeISf  
 %---------------------------------------------------设置Mur的吸收边界条件---------- R(f6uO!m  
 %------------------------------% <72q^w  
  R$bDj>8  
     for j=2:ny-1 I
Xpn(vX  
        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)); O>d [;Q  
    end g(dReC  
 %边界x=nx H'}6Mw%ra  
    for j=2:ny-1 l4ru0V8s7  
       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)); O=}d:yZb!  
    end rE%HNPO  
  %边界y=1  hv*XuT/  
    for i=2:nx-1 -I[KIeF  
        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)); 5c$\DZ(  
    end S]Ye`  
  %边界y=ny 9qvl9,*g  
    for i=2:nx-1 :!wl/X ~  
        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)); v8K`cijSS  
    end G7%f| Y  
  D2I|Z  
    ;tC$O~X  
  %---------------------------------------------------------角点处理------------------------------------------------------------% 3F32 /_`  
    Hz(1,1)=dz(2,2)+bb*(Hz(2,2)-dz(1,1));%角点x=1,y=1 .'b|pd  
   ?\$6"c<G  
    Hz(nx,1)=dz(nx-1,2)+bb*(Hz(nx-1,2)-dz(nx,1));%角点x=nx,y=1 Jwpc8MQ  
     bnZ H  
    Hz(1,ny)=dz(2,ny-1)+bb*(Hz(2,ny-1)-dz(1,ny) .. uC%mG
Za  
?5;N=\GQ  

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

不错，支持

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

这个的确要支持下

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

恩，贵在坚持啊，：）

不错，支持下

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

相互学习 共同进步啊  

赞扬~ .j}]J:{%  
qng ~,m  
学习中~