plasma.h

一维等离子体的plrc—fdtd算法程序

#define CollisionalFreq 100.0e6
#define wb 1.0e8*2.0*pii//1e-25//3.0e11
double PlasmaFreq(int z)
{	return 5.0e7*2.0*pii;//sqrt(5.0)*17.83e9/2;sqrt((z-MZ/2+thick)/(2.0*thick));//wp.
}

double chi(int m,int z)
{
	complex<double> cha;
	double wp,vc;
	wp=PlasmaFreq(z);
	vc=CollisionalFreq;
	cha=wp*wp*dt*(vc+j*wb)/(vc*vc+wb*wb)-(wp*wp*(vc+j*wb))/((vc*vc+wb*wb)*(vc-j*wb))*(1.0-exp(-1.0*vc*dt+j*wb*dt))*exp(-m*vc*dt+j*double(m)*wb*dt);
    return real(cha);

}

double xi(int m,int z)
{
    complex<double> xx;
	double wp,vc;
	wp=PlasmaFreq(z);
	vc=CollisionalFreq;
	xx=wp*wp*(vc+j*wb)*dt/(vc*vc+wb*wb)/2.0-wp*wp*(vc+j*wb)/(vc*vc+wb*wb)/(vc-j*wb)/(vc-j*wb)/dt*(1.0-(1.0+(vc-j*wb)*dt)*exp(-1.0*(vc-j*wb)*dt))*exp(-double(m)*(vc-j*wb)*dt);
	//if(z==MZ/2)cout<<real(xx)<<endl;
	return  real(xx);
  	// return  0.0;
}

inline double dChi0(int z){	return chi(0,z)-chi(1,z);}
inline double dXi0(int z){	return xi(0,z)-xi(1,z);}