fdtd_1d_hzey.h

2DFDTD程序

// FDTD_1D_HzEy.h: interface for the FDTD_1D_HzEy class.
//
//one dimensional FDTD algorithm with total field - scatterd field formulation
//without boundary conditions, the incident field is determined analytically
//
//////////////////////////////////////////////////////////////////////

#if !defined(AFX_FDTD_1D_HZEY_H__5A7A1439_641E_495C_A4A7_7F45E3C24B49__INCLUDED_)
#define AFX_FDTD_1D_HZEY_H__5A7A1439_641E_495C_A4A7_7F45E3C24B49__INCLUDED_

#if _MSC_VER > 1000
#pragma once
#endif // _MSC_VER > 1000

class FDTD_1D_HzEy  
{
protected:
	double *Ey_1D, *Hz_1D;

private:
	//the parametrs of the FDTD algorithm
	double pi;
	double *K_E1_a, *K_E1_b, *K_E2_a, *K_E2_b;

	int nx, n_PML; //the length of the computational space [cells]
	
	int source_type; //the type of the excitation
	double H0, t0, tw, omega, phi; //the parameters of the excitation
	double dx, dt;
	double eps_0, mu_0, eps_r, mu_r; //the speeed of the light in the media
	int i;

public:
	FDTD_1D_HzEy();
	virtual ~FDTD_1D_HzEy();

	BOOL Init_Main_Param_1D(int n_x, int n_pml, double Eps_r, double Mu_r, 
		                    double d_t, double d_x);
	
	void Set_PML_Param_1D();
	
	void calc_Hz_1D(double t);
	void calc_Ey_1D();
	
	void Init_Gauss_1D(double H_0, double t_0, double t_w);
	void Init_Sinus_1D(double H_0, double om, double Phi);
	void Init_Wave_Packet_1D(double H_0, double t_0, double t_w, double om, double Phi);

	void Get_Data_1D(double *&X, double *&Y);
	
};

#endif // !defined(AFX_FDTD_1D_HZEY_H__5A7A1439_641E_495C_A4A7_7F45E3C24B49__INCLUDED_)