main.cpp

计算了在有耗介质中半波天线的近场和原场的电场和磁场分量

#include "EHfield.h"
void main()
{
	FILE *fp;
	fp=fopen("result.txt","w");

	//激励源的参数
	double	f=1.14e8;
	double	T=1e-4;
	double	tao=T*(2e-5);
	double	Amp=1;
	
	//相对介质参数
	double	EPSr=78;
	double	XMUr=1;
	double	B=0.036;
	
	//天线参数	
	double lcoffie=0.375;								//天线长度与波长的系数比
	double radiu=0.00318;
	double EthM;
		
	Mparameter M;
	EH P;
		
//	double EM[19],HM[19];
	doubleComplexMatrix *E,*H;
	doubleComplex Avect;							//在场点处的位矢量

	E=doubleComplexMatrixAlloc(19,1);
	H=doubleComplexMatrixAlloc(19,1);

	//设置媒质和天线
	setMedium(EPSr,XMUr,B);
	setSource(T,tao);
	M=setparameter(f);									//设置天线、介质、激励源的参数
	setAntenna(M,lcoffie,radiu);					    
	
	Antsolvecoffie();									//求天线上电流的分布
	PrintIline(fp);
	
/*	doubleVector Z;
	Z.R=1.125*M.wavelength;
	Z.thinc=PI/2;
	Z.phinc=0;
	P=EHivalue(M,Z);
	EthM=P.EthM;
	//设置场点的位置
	for(int i=0;i<37;i++)
	{
		Z.thinc=i*PI/18;
		P=EHivalue(M,Z);
	//	fprintf(fp,"ehinc  %d  度  模值：E[%d]=%f\n",i*10,i,P.EthM/EthM);
		fprintf(fp,"%d,%f;\n",i*10,P.EthM/EthM);
		cout<<"ehinc "<<right<<setfill(' ')<<setw(3)<<i*10<<"度  模值：H["<<i<<"]＝"<<P.EthM/EthM<<endl;
	}
	cout<<endl;*/
}