1157阅读
0回复

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

UID ：90389

注册：2012-03-16
登录：2014-09-11
发帖：17
等级：仿真新人

0楼发表于: 2014-06-11 08:53:35

clear;
clc;
c=3e8;
e0=1e-9/(36*pi);
m0=4e-7*pi;
einf=1;
nz=300;
lambda=1;%wavelength
f=c/lambda;%frequency
omega=2*pi*f;
d=lambda/20;
s=1;
mur1=(s-1)/(s+1);
dt=s*d/c;
ex=zeros(nz+1,1);
hy=zeros(nz,1);
psi=zeros(nz-1,1);
nu=0.01/dt;
deltaEpsilon=1;
expNuDt=exp(-nu*dt);
chi0=deltaEpsilon*(1-expNuDt);
deltachi0=deltaEpsilon*(1-expNuDt)^2;
h=plot(ex);
ylim(2*[-1,1]);
for n=1:600
    ex0=ex(end-1);
    ex(1)=sin(omega*n*dt);
    ex(2:end-1)=(einf*ex(2:end-1)+psi-dt/d/e0*diff(hy,1,1))/(einf+chi0);
    psi=ex(2:end-1)*deltachi0+expNuDt*psi;
    ex(end)=ex0+mur1*(ex(end-1)-ex(end));
    hy=hy-dt/d/m0*diff(ex,1,1);
  & ..

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

发帖回复