www.pudn.com > loopantenna.rar > point.m, change:2002-04-17,size:1097b


function[EField, HField]=... 
    point1(Point,eta_,K,DipoleMoment,DipoleCenter) 
 
%POINT Radiated/scattered field at a point of a dipole array  
%   or a single dipole. Gives exact near- and far-fields. Outputs 
%   individual contribution of each dipole. 
% 
%   Observation point                   Point(1:3)          
%   Array of dipole moments             DipoleMoment(1:3,1:EdgesTotal)  
%   Array of dipole centers             DipoleCenter(1:3,1:EdgesTotal) 
%   E-field at the observation point    E(1;3,1:EdgesTotal) 
%   H-field at the observation point    H(1;3,1:EdgesTotal) 
% 
%   Copyright 2002 AEMM. Revision 2002/03/11  
%   Chapter 3 
 
C=4*pi; 
ConstantH=K/C; 
ConstantE=eta_/C; 
     
m=DipoleMoment; 
c=DipoleCenter; 
r       =repmat(Point,[1 length(c)])-c(1:3,:); 
PointRM =repmat(sqrt(sum(r.*r)),[3 1]); 
EXP     =exp(-K*PointRM); 
PointRM2=PointRM.^2; 
C=1./PointRM2.*(1+1./(K*PointRM)); 
D=repmat(sum(r.*m),[3 1])./PointRM2; 
M=D.*r; 
HField=ConstantH*cross(m,r).*C.*EXP; 
EField=ConstantE*((M-m).*(K./PointRM+C)+2*M.*C).*EXP;