基于MATLAB的架空输电线的等值电路参数计算

500kV输电线路长600km，采用三分裂导线3×LGJQ-400，分裂间距为400mm，三相水平排列，相间距离为11m，LGJQ-400导线的计算半径r=13.6mm。试计算输电线π型等值电路的参数：
（1）不计线路参数的分布特性；
（2）近似计及分布特性；
（3）精确特性计及分布特性。

%计算半径mm：
r=13.6
%分裂间距mm：
d=400

%线路长km：
l=600

%相间距离m：
D=11

%首先计算该线路的单位长度电路参数：
%单根导线的自几何均距mm：
Ds=0.9*r
Ds=12.24

%分裂导线的自几何均距mm：
Dsb=(Ds*d^2)^(1/3)
Dsb=125.11mm

%相间几何均距m：
Deq=1.26D
Deq=13.86

%每相导线组的等值半径mm：
req=(r*d^2)^(1/3)
req=129.58
%单位长度的直流电阻Ω/km:
s=3*d
r0=ρ/s
r0=0.0262
%单位长度的电抗Ω/km：
x0=0.1445*log10(Deq/Dsb)
x0=0.295
%单位长度的电纳S/km：
b0=7.58e-6/(log10(Deq/req))
b0=3.735e-6
%(1)不计线路的分布特性参数
%单位：Ω
R=r0*l
R=15.72

%单位：Ω
X=x0*l
x=177

%单位：S
B=b0*l
B=2.241e-3
（2）近似计及分布特性（利用修正参数）
kr=1-1/3*b0*l^2
kx
4000
=1-1/6*(x0*b0-r0^2*b0/x0)*l^2
kb=1+1/12*x0*b0*l^2

%单位：Ω
R=kr*r0*l
R=13.64

%单位：Ω
X=kx*x0*l
X=165.39

%单位：S
B=kb*b0*l
B=2.316e-3

%精确计及分布参数
Zc=((r0+x0*j)/(b0*j))^(1/2)
Zc=281.315-12.468j
gama=((r0+x0*j)*(b0*j))^(1/2)
gama=4.6567e-5+1.0507e-3j

Z=Zc*sinh(gama8*l)
Z=13.703+165.61j

Y=Z*(cosh(gama*l)-1)/Zc/sinh(gama*l)
Y=6.267e-6+2.3179e-3j