[RFEDA原创]FDTD迷你代码-1D UPML

;_nV*G.y#^  
clear; L>ruNw'-K  
clc; #~JR_oQE!  
c=3e8; W,D$=Bg  
m0=1e-7*4*pi; p]]*H2UD  
e0=1e-9/(36*pi); ,q".d =6  
eta0=120*pi; o4U[;.?c  
nz=100; &tjv.t  
na=10;%absorbing boundary ,zP.ch0K  
ex=zeros(1,nz+1); \* /R6svz  
hy=zeros(1,nz); ,,Ia4c  
lambda=1;%wavelength Cqra\  
f=c/lambda; 'CfM'f3uu  
omega=2*pi*f; (Qp53g  
nd=20; d1{%z\u a  
d=lambda/nd; -Uo"!o>x|  
s=1;  L7rEMq  
dt=s*d/c;  %&81xAt  
cb=dt/(d*e0); +\ZaVi  
cq=dt/(d*m0); RAs5<US:  
fh=@(t)(heaviside(t).*sin(t).*((heaviside(t)-heaviside(t-6*pi)).*(1-cos(pi*t/(6*pi)))/2+heaviside(t-6*pi))); >hFg,5 _l3  
fh=@(t)(fh(omega.*t)); ;,2;J3,pA  
m=1; l`b%imX  
sigmamax=(m+1)/(150*pi*d); QcN$TxU>  
sigma=sigmamax*((na-1:-1:0)/(na-1)).^m; u,88V@^  
sigmam=sigma*m0/e0; Cvgk67C=$  
kappa=1; I;5:jT`  
ca0=(e0*kappa/dt-(sigma(1:end-1)+sigma(2:end))/2)./(e0*kappa/dt+(sigma(1:end-1)+sigma(2:end))/2); M&h`uO/[  
cb0=1/d./(e0*kappa/dt+(sigma(1:end-1)+sigma(2:end))/2); 0kiV-yc   
cp0=(m0*kappa/dt-sigmam/2)./(m0*kappa/dt+sigmam/2); a*N<gId  
cq0=1/d./(m0*kappa/dt+sigmam/2); oScKL#Hu  
ca1=fliplr(ca0); hLo>jE  
cb1=fliplr(cb0); 0;X0<IV  
cp1=fliplr(cp0); Ir4M5OR\  
cq1=fliplr(cq0); uJ:SN;  
h=plot(ex); T!ik"YZ@i  
set(gca,'ylim',[-2,2]); P-LdzVt(^  
for n=0:500 _yu_Ev}R  
    ex(1)=fh(n*dt); hlO,mU  
    ex(2:end-na)=ex(2:end-na)-cb*diff(hy(1:end-na+1));  F<Y>  
    ex(end-na+1:end-1)=ca1.*ex(end-na+1:end-1)-cb1.*diff(hy(end-na+1:end)); m`lxQik  
    hy(1:end-na)=hy(1:end-na)-cq*diff(ex(1:end-na)); BK[ YX)  
    hy(end-na+1:end)=cp1 .. ;XDGlv%  
][[\!og  

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