一维FDTD等离子体的Matlab,修正过的一维FDTD等离子体MATLAB代码(公式修正)

修正过的一维FDTD等离子体MATLAB代码

clear;

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%%%%%%%%%%%%%系统参数

TimeT=3000;%迭代次数

KE=2000;%网格数目

kc=400;%源的位置

kpstart=600;%等离子体开始位置

kpstop=1000;%等离子体终止位置

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%%%%%%%%%%物理参数

f=1.356e7;%电磁场频率

fpe=1e8;%等离子体频率

c0=3e8;%真空中波速

zdelta=1e-1;%网格大小

dt=zdelta/(2*c0);%时间间隔

u0=2e10;%碰撞频率

wpe=2*pi*fpe;%等离子体圆频率

epsz=1/(4*pi*9*10^9); % 真空介电常数

mu=1/(c0^2*epsz);%磁常数

ex_low_m1=0;

ex_low_m2=0;

ex_high_m1=0;

ex_high_m2=0;

a0=2*u0/dt+(2/dt)^2;

a1=-8/(dt)^2;

a2=-2*u0/dt+(2/dt)^2;

b0=wpe^2+2*u0/dt+(2/dt)^2;

b1=2*wpe^2-8/(dt)^2;

b2=wpe^2-2*u0/dt+(2/dt)^2;

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%%%%%%%%%%初始化电磁场

Ex=zeros(1,KE);

Ex_Pre=zeros(1,KE);

Hy=zeros(1,KE);

Hy1=zeros(1,KE);

Hy_Pre=zeros(1,KE);

Dx=zeros(1,KE);

Dx_Pre=zeros(1,KE);