39 lines
1.2 KiB
Matlab
39 lines
1.2 KiB
Matlab
%%
|
||
%利用镜像法计算架空线场强
|
||
%依据 韦钢 李海峰 《同杆并架多回线下方的电场强度和感应电压》
|
||
%%
|
||
clc
|
||
clear
|
||
%坐标,单回
|
||
ConductorX=[-9,0,9];
|
||
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);];
|
||
%形成电位系数矩阵
|
||
%导线与镜像的距离
|
||
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));
|
||
%导线与导线的距离
|
||
conductor2conductorDistance=matCplxConductorXY-conj(matCplxConductorXY');
|
||
conductor2conductorDistance=conductor2conductorDistance-diag(diag(conductor2conductorDistance));
|
||
%计算电位系数
|
||
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;
|
||
%求电荷
|
||
QRI=P\Volt;
|
||
%求电场强度
|