potential/Field.m

%%
%<EFBFBD><EFBFBD><EFBFBD>þ<EFBFBD><EFBFBD>񷨼<EFBFBD><EFBFBD><EFBFBD><EFBFBD>ܿ<EFBFBD><EFBFBD>߳<EFBFBD>ǿ
%<EFBFBD><EFBFBD><EFBFBD><EFBFBD> Τ<EFBFBD><EFBFBD> <EFBFBD><EFBFBD><EFBFBD> <EFBFBD><EFBFBD>ͬ<EFBFBD>˲<EFBFBD><EFBFBD>ܶ<EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD>·<EFBFBD><EFBFBD>ĵ糡ǿ<EFBFBD>Ⱥ͸<EFBFBD>Ӧ<EFBFBD><EFBFBD>ѹ<EFBFBD><EFBFBD> 
%%
clc
clear
%<EFBFBD><EFBFBD><EFBFBD><EFBFBD>,<EFBFBD><EFBFBD><EFBFBD><EFBFBD>
ConductorX=[-9,0,9];
ConductorY=[14.5,14.5,14.5];
%<EFBFBD><EFBFBD><EFBFBD><EFBFBD>H<EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD>߾<EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD>ĸ߶<EFBFBD>
H=diag(ConductorY);
r=0.03/2*eye(length(ConductorX));%<EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD>Լ<EFBFBD><EFBFBD>ξ<EFBFBD><EFBFBD><EFBFBD>
MirrorConductorX=[-9,0,9];
MirrorConductorY=-[14.5,14.5,14.5];
%<EFBFBD><EFBFBD>ѹ<EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD>ң<EFBFBD>
%<EFBFBD><EFBFBD><EFBFBD><EFBFBD>
Volt=[525; 525*exp(1j*240*pi/180); 525*exp(1j*120*pi/180);];
%<EFBFBD>γɵ<EFBFBD>λϵ<EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD>
%<EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD>뾵<EFBFBD><EFBFBD><EFBFBD>ľ<EFBFBD><EFBFBD><EFBFBD>
cplxConductorXY=ConductorX+1j*ConductorY;
cplxMirrorConductorXY=MirrorConductorX+1j*MirrorConductorY;
matCplxConductorXY=repmat(cplxConductorXY,length(ConductorX),1);
matCplxMirrorConductorXY=repmat(conj(cplxMirrorConductorXY'),1,length(ConductorX));
conductor2MirrorDistance=abs(matCplxConductorXY-matCplxMirrorConductorXY);
conductor2MirrorDistance=conductor2MirrorDistance-diag(diag(conductor2MirrorDistance));
%<EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD>뵼<EFBFBD>ߵľ<EFBFBD><EFBFBD><EFBFBD>
conductor2conductorDistance=matCplxConductorXY-conj(matCplxConductorXY');
conductor2conductorDistance=conductor2conductorDistance-diag(diag(conductor2conductorDistance));
%<EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD>λϵ<EFBFBD><EFBFBD>
eslong=1;
P1=1/pi/eslong*log(2*H./r);
P1(isnan(P1))=0;
P2=1/pi/eslong*log(conductor2MirrorDistance./conductor2conductorDistance);
P2(isnan(P2))=0;
P=P1+P2;
%<EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD>
QRI=P\Volt;
%<EFBFBD><EFBFBD><EFBFBD>糡ǿ<EFBFBD><EFBFBD>