fdtd_1d_ezhy.h

2维时域有限差分模拟

//////////////////////////////////////////////////////////////////////
// fdtd_1D.h: interface for the fdtd_1D class.
//////////////////////////////////////////////////////////////////////
#pragma once

class CFdtd_1D_EzHy  
{
private:
	double *Ez, *Hy;
	int nx; //Dimension of the simulation space
	double E0, omega, phi, t0, tw, alfa;
	int Source_Type;

	double pi, eps_0, mu_0;
	int jel_source; // 1 - Gauss, 2 - Sin, 3 - SinGauss

	//parameters of the boundary condition and materials
	double K_E, K_H, K_MUR_1, K_MUR_2;
	double Ez_n_1[2];
	double Ez_n[2];
	double Ez_nxm2_n;

public:
	CFdtd_1D_EzHy();
	virtual ~CFdtd_1D_EzHy();
	
	bool Init(int n_x, double dt, double dx, double eps_r, double mu_r);
	
	void Init_Sin_Source(double E_0, double om, double Phi);
	void Init_Gauss_Source(double E_0, double t_0, double t_w);
	void Init_GaussSin(double E_0, double om, double Phi, double t_0, double t_w);

	void calc_Ez(double time);
	void calc_Hy();

	void Get_Data(double *&ez, double *&hy);
};