[RFEDA原创]FDTD迷你代码-1D 连接边界

#A@d;U%  
clear; / Ml d.  
clc; RPh8n4&("  
c=3e8; 8euh]+  
m0=1e-7*4*pi; Ed^F_Gg#  
e0=1e-9/(36*pi); Fb^Ae6/i  
eta0=120*pi; 8YLS/dN0 w  
nz=100; Wb7z&vj  
nc=10; h88IP:bo  
ex=zeros(1,nz+1); Cr;d !=  
hy=zeros(1,nz); q]}1/JZS  
lambda=1;%wavelength .zMM!l3  
f=c/lambda; A:JWUx  
omega=2*pi*f; flr&+=1?D  
nd=20; n?EL\B   
d=lambda/nd; @9}SHS  
s=1; "KhVS  
dt=s*d/c; QN'v]z  
cb=dt/(d*e0); [h=[@ji
B  
cq=dt/(d*m0); n@e[5f9?x  
mur1=(s-1)/(s+1); $(PWN6{\r^  
t0=(nz/2-nc+2)*d/c; L(ni6-  
z=(-nz/2:nz/2)*d; D,)^l@UP  
fh=@(t)(heaviside(t).*sin(t).*((heaviside(t)-heaviside(t-6*pi)).*(1-cos(pi*t/(6*pi)))/2+heaviside(t-6*pi))); =#Qm D=  
exi=@(t)(fh(omega.*t)); #2U4}#Mi  
hyi=@(t)(exi(t)/eta0); (46S^*  
h=plot(ex); h5>38Kd  
set(gca,'ylim',2*[-1,1]); w?db~"T  
for n=0:200 6lKM5,Oa  
    ex0=ex([2,end-1]); TXDb5ZCzM  
    ex(2:end-1)=ex(2:end-1)-cb*diff(hy); ="Edt+a)t  
    t=(n+0.5)*dt-t0-(z(1+nc)-0.5*d)/c; uJX(s6["=  
    ex(1+nc)=ex(1+nc)+cb*hyi(t);%connection rQ!X  
    t=(n+0.5)*dt-t0-(z(1+nz-nc)+0.5*d)/c;%connection VdfV5"  
    ex(end-nc)=ex(end-nc)-cb*hyi(t); ZR|cZH1}C  
    ex([1,end])=ex0([1,end])+mur1*(ex([2,end-1])-ex([1,end]));%ABC ~l@-g
Ayw  
    hy=hy-cq*diff(ex); A!\-e*+W=  
    t=(n+1)*dt-t0-(z(1+nc))/c; O:
>9yZhV  
    hy(nc)=hy(nc)+cq*exi(t); F#az&  
    t=(n+1)*dt-t0-(z(1+nz-nc))/c; C $r]]MSj  
     .. 0u7\*Iy  
)b_ GKA `  

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