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

S=|@L<O  
clear; V A^l+Z,d  
clc; 8MwK.H[U  
c=3e8; 6QQfQ,  
m0=1e-7*4*pi; clO,}Ph>  
e0=1e-9/(36*pi); MjL)IgT  
eta0=120*pi; mh/n.*E7  
nz=100; |bv,2uWz  
ex=zeros(1,nz+1); i'/m4 !>h  
hy=zeros(1,nz); #;KsJb)N.  
lambda=1;%wavelength W&Y"K)`  
f=c/lambda; uzr\
oj+>  
omega=2*pi*f; V&{MQWy  
nd=20; ^C)n$L>C0  
d=lambda/nd; g{k1&|  
s=1; kma?v B  
dt=s*d/c; p*QKK@C  
cb=dt/(d*e0); l,|%7-  
cq=dt/(d*m0); ~l[ra  
mur1=(s-1)/(s+1); %r{3wH#D@  
z=(-nz/2:nz/2)*d; iP@6hG`:  
fh=@(t)(heaviside(t).*sin(t).*((heaviside(t)-heaviside(t-6*pi)).*(1-cos(pi*t/(6*pi)))/2+heaviside(t-6*pi))); 8T[ 6J{|C  
fh=@(t)(fh(omega.*t)); Y-!YhWsS  
h=plot(ex);
/7c~nBU  
set(gca,'ylim',[-2,2]); BHDML.r }M  
for n=0:200 nSx8E7 |V  
    Ex0=ex(end-1); p/6zEZ*  
    ex(2:end-1)=ex(2:end-1)-cb*diff(hy); l:|D,q  
    t=n*dt; x[_=#8~.1x  
    ex(1)=fh(t); 4vnUN  
    ex(end)=Ex0+mur1*(ex(end-1)-ex(end));%ABC ,~PYt*X4  
    hy=hy-cq*diff(ex); {lG@h
N'
 
    set(h,'y .. )Q;978:  
5f'DoT  

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

提示:鼓励发贴交流加RF币！希望你将FDTD迷你代码汇总成一个帖子。再次分享之后，再给你加技术分！