[RFEDA原创]FDTD迷你代码-1D 色散介质（Debye模型，RC方法）

i`Lt=)@&  
clear; s)~Q@ze2  
clc; ={#r/x  
c=3e8; ]E[Mv} =  
e0=1e-9/(36*pi); ,h8)5Mj/J  
m0=4e-7*pi; -_b}b)2iYN  
einf=1; |Vd)7/LN  
nz=300; 0fi+tc30  
lambda=1;%wavelength -qid.  
f=c/lambda;%frequency sI9~TZ :  
omega=2*pi*f; NBUM* Z  
d=lambda/20; K'1rS[^>R  
s=1; ,((5|MbM/  
mur1=(s-1)/(s+1); >qcir~ &  
dt=s*d/c; ;9B:E"K?@1  
ex=zeros(nz+1,1); <aL$d7  
hy=zeros(nz,1); k\HRG@ /G  
psi=zeros(nz-1,1); Q'$aFl'NR  
nu=0.01/dt;
A3Ltk 2<  
deltaEpsilon=1; 5-=mtvA:  
expNuDt=exp(-nu*dt); &qr
H  
chi0=deltaEpsilon*(1-expNuDt); "b-6kM  
deltachi0=deltaEpsilon*(1-expNuDt)^2; XK(<N<Z@|e  
h=plot(ex); a)y8MGx?  
ylim(2*[-1,1]); (doFYF~w  
for n=1:600 Dc #iM0  
    ex0=ex(end-1); cJi5\<b  
    ex(1)=sin(omega*n*dt); 9FJU'$FN  
    ex(2:end-1)=(einf*ex(2:end-1)+psi-dt/d/e0*diff(hy,1,1))/(einf+chi0); !{Y$5)Xh`]  
    psi=ex(2:end-1)*deltachi0+expNuDt*psi; OEq8gpqY  
    ex(end)=ex0+mur1*(ex(end-1)-ex(end)); j&Z:|WniK  
    hy=hy-dt/d/m0*diff(ex,1,1); /+02BP  
  & ..
LR-op?W  
==UH)o`?8  

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