ehfield.h

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

#ifndef E_H
#define E_H
#include "antenna.h"
#include "vector.h"

typedef struct
{
	double r;
	double SINTH;
	double COSTH;
	double THINC;
	double PHINC;
}Observe;

typedef struct
{
	double EthM;
	double ErM;
	double HphM;
	doubleComplex Eth;
	doubleComplex Er;
	doubleComplex Hph;
}EH;
//doubleComplex K,ETA;
//double wavelength,w,XMU;

void Esetting(Mparameter P);									//将所求的媒质参数赋给场运算
////////////////////////////////////////////////////////////////////////
//
//                            目  标  点
//
////////////////////////////////////////////////////////////////////////
Observe verticalsetaimstation(doublexyzVector x);            //在直角坐标系中设置目标点的位置
Observe circlepolarsetaimstation(doubleVector r);  //在球坐标系中设置目标点的位置

////////////////////////////////////////////////////////////////////////
//
//                             矢  量  位(thinc)
//
////////////////////////////////////////////////////////////////////////
doubleComplex Athvector(double L);                     //辛普生法求矢量位函数
doubleComplex PthAvectorcell(double zi);                        //远去场矢量位被积函数
doubleComplex PthAvector(double L);                    //远去场辛普生法求矢量位函数
doubleComplex EthAvectorcell(double zi);                        //近场矢量位被积函数
doubleComplex EthAvector(double L);                    //近场辛普生法求矢量位函数

////////////////////////////////////////////////////////////////////////
//
//                             矢  量  位(r)
//
////////////////////////////////////////////////////////////////////////
doubleComplex Arvector(double L);                     //辛普生法求矢量位函数
doubleComplex PrAvectorcell(double zi);                        //远去场矢量位被积函数
doubleComplex PrAvector(double L);                    //远去场辛普生法求矢量位函数
doubleComplex ErAvectorcell(double zi);                        //近场矢量位被积函数
doubleComplex ErAvector(double L);                    //近场辛普生法求矢量位函数
EH EHivalue(Mparameter P,doublexyzVector Z);
EH EHivalue(Mparameter P,doubleVector Z);
EH EHisolve(double L);
/*
入口函数:verticalsetaimstation(double x,double y,double z);           
         circlepolarsetaimstation(double r,double THINC,double PHINC);
出口函数:Avector(Mparameter P,double a,double b);
		 Efield(doubleComplex Avect);
*/
#endif