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

《一维等离子体FDTD的Matlab源代码两种方法.doc》由会员上传分享，免费在线阅读，更多相关内容在应用文档-天天文库。

1、一维等离子体FDTD的Matlab源代码(两种方法）%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%1=V-V< %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%*$g-:ILRuZ %%%%%%%%%%%%%%%%%%%%        1D     

2、         %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%oCz/HQoBk %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%k9L;!TH~1K %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

3、'D1xh~ %%%%%%%%%初始化QVgl(;lX clear;3^yK!-Wp( %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%utV_W& %%%%%%%%%%%%%%%%系统参数uwGc@xOgg, TimeT=3000;%迭代次数qIT@g"%}t KE=2000;%网格树木p4Z(^+Aa kc=450;%源的位置f3y=Wxk[ kpstart=500;%等离子体开始位置

4、2A:eI8^ kpstop=1000;%等离子

5、体终止位置'LDQgC*% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%A@#E@;lm %%%%%%%%%%%%%物理参数V
!~wj c0=3e8;%真空中波速3Jn;} zdelta=1e-9;%网格大小#GFr`o0$^ dt=zdelta/(2*c0);%时间间隔))Za&S*< f=900e12;%Gause脉冲的载频

6、iN~j&Be u0=57e12%碰撞频率m]&SNz= fpe=2000e12;%等离子体频率K(

7、}dl: wpe=2*pi*fpe;%等离子体圆频率m4Zk,1m.

8、 epsz=1/(4*pi*9*10^9);%真空介电常数$/]
,tSm mu=1/(c0^2*epsz);%磁常数;9#KeA_ ex_low_m1=0;yt2PU_), ex_low_m2=0;CvdN"k ex_high_m1=0;8FhdN ex_high_m2=0;v`r:=K a0=2*u0/dt+(2/dt)^2;'

9、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"
{` %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%&mvSiyKX %%%%%%%%%%%%%初始化电磁场?X;RLpEc

10、A 

11、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+iP6B Sx1=zeros(1,KE);(iGTACoF Sx2=zeros(1,KE);Wez5N Sx3=zeros(1,KE);U;I9bK8 Sx=zeros(1,KE);-.3w^D"l Dx=Ex;

12、nF/OPd Dx1=Ex;Qci]i)s$js Dx2=Ex;3N:D6
w-R Dx3=Ex;:i7;w%
B Ex1=Ex;cS+>J@L Ex2=Ex;WjjB