facat
/
potential
1
0
Fork 0
Code Issues Pull Requests Releases Wiki Activity

按《同杆并架多回线下方的电场强度和感应电压》 做，想和“设计手册”例题2-1比较，做了一半没有继续下去。 

Signed-off-by: dmy@lab <dmy@lab.com>
This commit is contained in:
dmy@lab 2015-09-18 13:16:16 +08:00
parent 0e7f2af5d1
commit 2ad26b0a31
1 changed files with 28 additions and 41 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

69
Field.m
View File

@ -4,48 +4,35 @@
%% %%
clc clc
clear clear
%电压（复数，至上到下，先左后右） %坐标,单回
%双回 ConductorX=[-9,0,9];
Volt=[220; 220*exp(1j*240*pi/180); 220*exp(1j*120*pi/180);220; 220*exp(1j*240*pi/180); 220*exp(1j*120*pi/180) ]; ConductorY=[14.5,14.5,14.5];
%计算H，导线距地面的高度
H=diag(ConductorY);
r=0.03/2*eye(length(ConductorX));%导线自几何均距
MirrorConductorX=[-9,0,9];
MirrorConductorY=-[14.5,14.5,14.5];
%电压（复数，先左后右）
%单回
Volt=[525; 525*exp(1j*240*pi/180); 525*exp(1j*120*pi/180);];
%形成电位系数矩阵 %形成电位系数矩阵
H=diag([28.5, 21, 27.5, 28.5, 21, 27.5 ]); %导线与镜像的距离
r=diag([0.2,0.2,0.2,0.2,0.2,0.2]); cplxConductorXY=ConductorX+1j*ConductorY;
Dcomma=[0,58.5,55,62.42,58.98,55.5; cplxMirrorConductorXY=MirrorConductorX+1j*MirrorConductorY;
58.5,0,51.5,58.98,55.5,52.1; matCplxConductorXY=repmat(cplxConductorXY,length(ConductorX),1);
55, 51.5 0,55.5,52.1 48.7; matCplxMirrorConductorXY=repmat(conj(cplxMirrorConductorXY'),1,length(ConductorX));
62.42,58.98,55.55,0,58.5,55; conductor2MirrorDistance=abs(matCplxConductorXY-matCplxMirrorConductorXY);
58.98,55.55,52.13,58.5,0,51.5; conductor2MirrorDistance=conductor2MirrorDistance-diag(diag(conductor2MirrorDistance));
55.55,52.13,48.73,55,51.5,0]; %导线与导线的距离
D=[ 0 3.51 7.02 7.2 8.27 10.4; conductor2conductorDistance=matCplxConductorXY-conj(matCplxConductorXY');
3.51 0 3.51 8.27 7.8 8.8 conductor2conductorDistance=conductor2conductorDistance-diag(diag(conductor2conductorDistance));
7.02 3.51 0 10.4 8.8 8.4 %计算电位系数
7.2 8.27 10.48 0 3.5 7.02
8.27 7.8 8.8 3.51 0 3.51
10.4 8.8 8.4 7.02 3.51 0];
eslong=1; eslong=1;
P=1/pi/eslong*log(2*H./r); P1=1/pi/eslong*log(2*H./r);
P(isnan(P))=0; P1(isnan(P1))=0;
P2=1/pi/eslong*log(Dcomma./D); P2=1/pi/eslong*log(conductor2MirrorDistance./conductor2conductorDistance);
P2(isnan(P2))=0; P2(isnan(P2))=0;
P=P+P2; P=P1+P2;
%求电荷 %求电荷
QRI=inv(P2)*Volt; QRI=P\Volt;
%导线坐标 %求电场强度
ABCx=[-3.6; -3.9; -4.2; 3.6; 3.9; 4.2];
ABCy=[31; 27.5; 24; 31; 27.5; 24];
%求地面场强
plotY=zeros(60,1);
plotX=zeros(60,1);
for I=-30:30
x=I;
y=0;
ERy=sum( (ABCy-y)./( (ABCy-y).^2+(ABCx-x).^2 ).*real(QRI)./2/pi/eslong );
EIy=sum( (ABCy-y)./( (ABCy-y).^2+(ABCx-x).^2 ).*imag(QRI)./2/pi/eslong );
ERx=sum( (ABCx-x)./( (ABCy-y).^2+(ABCx-x).^2 ).*real(QRI)./2/pi/eslong );
EIx=sum( (ABCx-x)./( (ABCy-y).^2+(ABCx-x).^2 ).*imag(QRI)./2/pi/eslong );
E=sqrt(ERy.^2+EIy.^2+ERx.^2+EIx.^2);
plotY(I+31)=E;
plotX(I+31)=x;
end
plot(plotX,plotY);