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

j[.nk  
clear; gKQV99  
clc; .6xMLo,R  
c=3e8; aU^6FI  
m0=1e-7*4*pi; <fM}Kk  
e0=1e-9/(36*pi); _TGs .t  
eta0=120*pi; ^'a#FbMtt  
nz=100; =*"8N-FU  
ex=zeros(1,nz+1); O1S7t)ag  
hy=zeros(1,nz); Od)y4nr3~  
lambda=1;%wavelength 5doi4b>]!  
f=c/lambda; &*%x]fQ@  
omega=2*pi*f; 8'%m
!  
nd=20; Jv+w{"&  
d=lambda/nd; cR_pC 9z  
s=1; @XXPJq;J  
dt=s*d/c; IWq#W(yM  
cb=dt/(d*e0); se`Eez}  
cq=dt/(d*m0); n\3#69VY  
mur1=(s-1)/(s+1); J=k=cFUX  
z=(-nz/2:nz/2)*d; h<TZJCt  
fh=@(t)(heaviside(t).*sin(t).*((heaviside(t)-heaviside(t-6*pi)).*(1-cos(pi*t/(6*pi)))/2+heaviside(t-6*pi)));  @OkoT:  
fh=@(t)(fh(omega.*t)); DA.k8M  
h=plot(ex); W\NC3]  
set(gca,'ylim',[-2,2]); `;J`O02  
for n=0:200 ".N+nM~  
    Ex0=ex(end-1); &J
catI  
    ex(2:end-1)=ex(2:end-1)-cb*diff(hy); v"W*@7<`S  
    t=n*dt; b^Re947{g  
    ex(1)=fh(t); 
zQG{j\  
    ex(end)=Ex0+mur1*(ex(end-1)-ex(end));%ABC /s} "0/Y\  
    hy=hy-cq*diff(ex); V~ORb1  
    set(h,'y .. <|jh3Hlp  
vg1JN"S[  

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

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