一维FDTD等离子体的Matlab,一维fdtd模拟程序 C语言

#include

#include

#include

#define NUM_of_ZAXIS 400

float gauss_pulse(float T,float t0,float spread);

int file_save(float* data,char* filename);

void main()

{

float ex[NUM_of_ZAXIS],hy[NUM_of_ZAXIS];

float obj_parameters[NUM_of_ZAXIS][4]; /*模型参数设置 */

float CA[NUM_of_ZAXIS],CB[NUM_of_ZAXIS],CP[NUM_of_ZAXIS],CQ[NUM_of_ZAXIS]; float ca,cb,cp,cq; /*真空时的参数*/

float var_ca,var_cb,var_cp,var_cq; /* */

float

var_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; /*边界吸收参数*/

float Epsilon,Mu,Pi,C; /*介电系数 Epsilon 0，磁导系数 Mu 0*/

float rel_epsz,rel_mu; /*相对介电系数，磁导系数*/

float e_sigma,h_sigma; /*电导率，磁导率*/

float dt,ddz;

float source,T;

int pos_driv_source;

int i,n,Nsteps;