一维等离子体FDTD的Matlab源代码(两种方法).docVIP

一维等离子体FDTD的Matlab源代码(两种方法） %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 1=V-V ?%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% *$g-:ILRuZ ?%%%%%%%%%%%%%%%%%%%%????????1D??????????????%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% oCz/HQoBk ?%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% k9L;!TH~1K ?%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% D1xh~ ?%%%%%%%%%初始化 Q\Vgl(;lX ?clear; 3^yK!-Wp( ?%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% utV_ W ?%%%%%%%%%%%%%%%%系统参数 uwGc@xOgg, ?TimeT=3000;%迭代次数 qIT@g%}t ?KE=2000;%网格树木 p4Z(^+Aa ?kc=450;%源的位置 f3y=Wxk[ ?kpstart=500;%等离子体开始位置 |2A:eI8 ^ ?kpstop=1000;%等离子体终止位置 LDQgC*% ?%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% A@#E@ ;lm ?%%%%%%%%%%%%%物理参数 V 
!~wj ?c0=3e8;%真空中波速 3 Jn ;} ?zdelta=1e-9;%网格大小 #GFr`o0$^ ?dt=zdelta/(2*c0);%时间间隔 ) )ZaS* ?f=900e12;%Gause脉冲的载频 C*hokqqP ?d=3e-15%脉冲底座宽度 .e-#
yET ?t0=2.25/f;%脉冲中心时间 uXiN~j Be ?u0=57e12%碰撞频率 m]SNz= ?fpe=2000e12;%等离子体频率 K (|}dl: ?wpe=2*pi*fpe;%等离子体圆频率 m4Zk\,1m.| ?epsz=1/(4*pi*9*10^9); % 真空介电常数 $/ ]
,tSm ?mu=1/(c0^2*epsz);%磁常数 ;9#KeA _ ?ex_low_m1=0; yt2PU_), ?ex_low_m2=0; CvdNk ?ex_high_m1=0; 8  FhdN ?ex_high_m2=0; v` r:=K ?a0=2*u0/dt+(2/dt)^2; hf8ZEW9 ?a1=-8/(dt)^2; +H2Qk4XFB ?a2=-2*u0/dt+(2/dt)^2; 2t,zLwBdnJ ?b0=wpe^2+2*u0/dt+(2/dt)^2; *lb$E]=! ?b1=2*wpe^2-8/(dt)^2; F:ELPs4 ?b2=wpe^2-2*u0/dt+(2/dt)^2; Vw\
{` ?%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% m vSiyKX ?%%%%%%%%%%%%%初始化电磁场 ?X;RLpEc|A ?Ex=zeros(1,KE); %XTI-B/K ?Ex_Pre=zeros(1,KE); rM l@3hP ?Hy=zeros(1,KE); i6N,jFU ?Hy_Pre=zeros(1,KE); E} .^kc[(4 ?Dx=zeros(1,KE); {XHh8_ ^ ?Dx_Pre=zeros(1,KE); K+iP 6B ?Sx1=zeros(1,KE); (iGTACoF ?Sx2=zeros(1,KE); We z 5N ?Sx3=zeros(1,KE