一维fdtd模拟程序 C语言.doc

《一维fdtd模拟程序 C语言.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、#includeH-.8{8 #includej5>3Td. #includeww"HV;i PKM8MYvo #defineNUM_of_ZAXIS400nE3'm[) 2U;ImC1g floatgauss_pulse(floatT,floatt0,floatspread);IZ
kQmA= L:;~L/ intfile_save(float*data,char*filename);m:7$"oq

2、 ,!s;o6

3、*y

4、 sSMcF[]@2I WU-_Y^ voidmain()_8CE

5、<Cn {(hn@+hc     floatex[NUM_of_ZAXIS],hy[NUM_of_ZAXIS];O!D0hW4     floatobj_parameters[NUM_of_ZAXIS][4];  /*模型参数设置*/n~g,qEI;

6、oatca,cb,cp,cq;        /*真空时的参数*/;oQ*gd     floatvar_ca,var_cb,var_cp,var_cq;        /**/L)cy&"L

7、     floatvar_border,ex_low_m1,ex_low_m2,ex_high_m1,ex_high_m2,    ex_low_s1,ex_low_s2,ex_high_s1,ex_high_s2;    /*边界吸收参数*/y-TS?5Dr]     floatEpsilon,Mu,Pi,C; 

8、         /*介电系数Epsilon0，磁导系数Mu0*/#mNM5(o     floatrel_epsz,rel_mu;  /*相对介电系数，磁导系数*/]e'Ol$3U9=     floate_sigma,h_sigma;  /*电导率，磁导率*/.&AS-">Z     ?X1vU0c     floatdt,ddz;
%8)GuxG*     floatsource,T;4

9、ryt4B     intpos_driv_source;P

10、o_/BS     inti,n,Nsteps

11、;<

12、_/i/H     inte_low,e_high;ZqHh$QBD9 HG2i^y VuDS
jh 4+t9"SD     FILE*fp;#:{Bd8PS     4#$~gTc@     NGD?.^(G     /*初始化各个变量*/8}3dwr;-     Pi=3.14159;/c1=`OJ     Epsilon=8.85e-12;=@U~sl[     Mu=4*Pi*1e-7;NPws^ uFDJRQJ<     C=pow((float)

13、(1/(Epsilon*Mu)),(float)(0.5));SomA`y+ERn wu5]S)?*     ddz=0.01;1$Rua     dt=ddz/(2*C);eRkvNI     T=0;,E{z+:Es                   @^8tk3$Y                     //采用葛德彪书中一维FDTD公式中的参数hqW4.

14、&c     e_sigma=0;l
+BJh1^     h_sigma=0;z1tYD     var_ca=

15、0.5*e_sigma*dt/Epsilon;w=.w*?>     var_cb=dt/Epsilon;lP*_dt9     var_cp=0.5*h_sigma*dt/Mu;okfGd=&     var_cq=dt/Mu;:_;9&[H9ha hsIC5@s3     ca=(1-var_ca)/(1+var_ca);4.

16、-?qG     cb=var_cb/(1+var_ca);s`W
`w}     cp=(1-var_cp)/(1+var_cp);#9aB3C   

17、  cq=var_cq/(1+var_cp);07WZw1(; 4.)
hCb     var_border=(C*dt-ddz)/(C*dt+ddz);@/i;/$
 o<@b]ukl& f`w$KVZ1!w                   //边界处用来存数据的临时变量__HPwOCG7     ex_low