mediu.cpp

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

#include "mediu.h"

//激励源
double wavelength;                         //波长
double f;                                  //脉冲的频率
double T;                                  //脉冲的周期
double tao;                                //脉冲持续的时间

double w,w0;                               //角频率

//媒质
double EPS;                                //天线所处媒质的介电常数	
double XMU;								   //天线所处媒质的磁导率
double WMG;                                //天线所处媒质的导电率

//介质参数
//double EPSr;                              //相对介电常数
//double XMUr;                              //相对磁导率
double B,V;                                 //复相对介电常数
doubleComplex K;                          //复数形式的传播常数
doubleComplex ETA;                       //复数形式的波阻抗倒数
doubleComplex EPSE;                       //复数形式的等效介电常
////////////////////////////////////////////////////////////////////////
//
//                             媒       质
//
////////////////////////////////////////////////////////////////////////
//设置天线所处媒质的介质参数
void setMedium(double EPSr,double XMUr,double b)  
{
	EPS=EPSr*EPS0;
	XMU=XMUr*XMU0;
	B=b;
}

////////////////////////////////////////////////////////////////////////
//
//                             激  励  源
//
////////////////////////////////////////////////////////////////////////
//设置脉冲波的参数
void setSource(double SourceT,double Sourcetao)			  
{
	T=SourceT;
	tao=Sourcetao;
}

////////////////////////////////////////////////////////////////////////
//
//                            参 数 设 置
//
////////////////////////////////////////////////////////////////////////
//设置介质等价参数、波的频率、电压V
Mparameter setparameter(double f)
{
	Mparameter P;
	w0=2*PI/T;
	w=2*PI*f;

	WMG=B*w*EPS;
	EPSE.re=EPS;
	EPSE.im=-EPS*B;
		
	K=doubleDMpy(w,sqrtC(doubleDMpy(XMU,EPSE)));	
	ETA=sqrtC(doubleDMpy(XMU,doubleCDiv(EPSE)));
	wavelength=2*PI/K.re;
	//	V=Amp*tao/T*sin(w*tao/2)/(w*tao/2);
	V=1;

	cout<<"Source\n";
	cout<<"w  "<<w<<"  T  "<<T<<"  tao  "<<tao<<"  V  "<<V<<endl<<endl;
	cout<<"Mediu\n";
	cout<<"wavelength  "<<wavelength<<endl;
	cout<<"EPSE  "<<EPSE.re<<"+i"<<EPSE.im<<endl;
	cout<<"XMU   "<<XMU<<endl;
	cout<<"WMG   "<<WMG<<endl;
	cout<<"K     "<<K.re<<"+i"<<K.im<<endl;
	cout<<"ETA   "<<ETA.re<<"+i"<<ETA.im<<endl<<endl;

	P.w=w;
	P.XMU = XMU;
	P.WMG = WMG;
	P.wavelength = wavelength; 
	P.K = K;
	P.ETA =ETA;
	P.EPSE = EPSE;
	P.V=V;

	return P;
}