def.h

二维FDTD算法

#ifndef DEFINITION_H
#define DEFINITION_H

/****************范围设置******************/
//总迭代范围
const int DOMAINX=60;  
const int DOMAINY=60;

//点源位置(中心)
const int IC=DOMAINX/2;
const int JC=DOMAINY/2;

/*************电磁参数**************/
#define PI 3.1415926535897932384626433832795

const double V=3.0*pow(10.0,8.0);         //波速
const double MIU0=4.0*PI*pow(10.0,-7.0);  //真空磁导率
const double EPSILON0=1.0/(4.0*PI*9.0*pow(10.0,9.0));  //真空电导率
const double Z0=sqrt(MIU0/EPSILON0);      //真空波阻抗

const double CFL=0.5;          //*迭代常数(c*deltat/deltax)
const double Coefficient1=CFL;  //总系数1

/******************源*******************/
const int	 T=1500;						    //迭代时间
const double FREQUENCY=3.0*pow(10.0,9.0);       //频率
const double WAVELENGTH=V/FREQUENCY;            //波长
const double DELTAX=0.005;  //0.05;				//空间步
const double DELTAT=DELTAX*CFL/V;	//时间步
const double TIME_period=1/(FREQUENCY*DELTAT);

//脉冲源设置
const double T0=25;  //2.0*pow(10.0,-9.0); //*中心
const double SPREAD=8.0;				   //*宽度

//PML厚度
const int NPML=10;

//PML测试
double maxA=0.;
double maxB=0.;
double EzA[1000];        //测试向量A
double EzB[1000];        //测试向量B

//场外部变量(x,y,z表示极化);
static double Ez[DOMAINX][DOMAINY];
static double Hx[DOMAINX][DOMAINY];
static double Hy[DOMAINX][DOMAINY];
#endif