关于电力系统潮流计算牛顿拉夫逊算法的MATLAB程序中 形成节点导纳矩阵时有个语句理解不了，请大神指教

因为没完整的程序，我没看清楚，你的程序种a(k)的意义，我猜想你说的支路类型，可能是k与节点N的关系类型，一般是支路是n指向K，a(k)是1，否则是-1，这么一来，(a(k)^2=1。

后面支路潮流计算中，需要用到支路的方向。但导纳矩阵，不需要计算y(k)的方向。

关于潮流计算的基础理论，你看看：

1,王锡凡，现代电力系统分析，科学出版社。

2.张伯明，高等电力网络分析，清华大学出版社。

3.于尔铿，能量管理系统，科学出版社。

clc

G(1,1)=3.750B(1,1)=-11.250G(1,2)=-2.500B(1,2)=7.500

G(1,3)=0B(1,3)=0G(1,4)=-1.250B(1,4)=3.750

G(1,5)=0B(1,5)=0

G(2,1)=-2.500B(2,1)=7.500G(2,2)=10.834B(2,2)=-32.500G(2,3)=-1.667B(2,3)=5.000

G(2,4)=-1.667B(2,4)=5.000G(2,5)=-5.000B(2,5)=15.000

G(3,1)=0B(3,1)=0G(3,2)=-1.667B(3,2)=5.000G(3,3)=12.917B(3,3)=-38.75

G(3,4)=-10.000B(3,4)=30.000G(3,5)=-1.250B(3,5)=3.750

G(4,1)=-1.250B(4,1)=3.750G(4,2)=-1.667B(4,2)=5.000G(4,3)=-10.000B(4,3)=30.000

G(4,4)=12.917B(4,4)=-38.750G(4,5)=0B(4,5)=0

G(5,1)=0B(5,1)=0G(5,2)=-5.000B(5,2)=15.000G(5,3)=-1.250B(5,3)=3.750

G(5,4)=0B(5,4)=0G(5,5)=6.25B(5,5)=-18.750

Y=G+j*B

delt(1)=0delt(2)=0delt(3)=0delt(4)=0

u(1)=1.0u(2)=1.0u(3)=1.0u(4)=1.0

ps(1)=-0.60qs(1)=-0.10ps(2)=0.20qs(2)=0.20

ps(3)=-0.45qs(3)=-0.15ps(4)=-0.40qs(4)=-0.05

k=1precision=1

N1=4

while precision>0.00001

delt(5)=0u(5)=1.06

for m=1:N1

for n=1:N1+1

pt(n)=u(m)*u(n)*(G(m,n)*cos(delt(m)-delt(n))+B(m,n)*sin(delt(m)-delt(n)))

qt(n)=u(m)*u(n)*(G(m,n)*sin(delt(m)-delt(n))-B(m,n)*cos(delt(m)-delt(n)))

end

pi(m)=sum(pt)qi(m)=sum(qt)

dp(m)=ps(m)-pi(m)

dq(m)=qs(m)-qi(m)

end

for m=1:N1

for n=1:N1

if m==n

H(m,m)=-qi(m)-u(m)^2*B(m,m)N(m,m)=pi(m)+u(m)^2*G(m,m)

J(m,m)=pi(m)-u(m)^2*G(m,m)L(m,m)=qi(m)-u(m)^2*B(m,m)

JJ(2*m-1,2*m-1)=H(m,m)JJ(2*m-1,2*m)=N(m,m)

JJ(2*m,2*m-1)=J(m,m)JJ(2*m,2*m)=L(m,m)

else

H(m,n)=u(m)*u(n)*(G(m,n)*sin(delt(m)-delt(n))-B(m,n)*cos(delt(m)-delt(n)))

J(m,n)=-u(m)*u(n)*(G(m,n)*cos(delt(m)-delt(n))+B(m,n)*sin(delt(m)-delt(n)))

N(m,n)=-J(m,n)L(m,n)=H(m,n)

JJ(2*m-1,2*n-1)=H(m,n)JJ(2*m-1,2*n)=N(m,n)

JJ(2*m,2*n-1)=J(m,n)JJ(2*m,2*n)=L(m,n)

end

end

end

JJ

for m=1:N1

PP(2*m-1)=dp(m)

PP(2*m)=dq(m)

end

uu=inv(JJ)*PP'precision=max(abs(uu))

for n=1:N1

delt(n)=delt(n)+uu(2*n-1)

u(n)=u(n)+uu(2*n)*u(n)

end

k=k+1

end

for n=1:N1+1

U(n)=u(n)*(cos(delt(n))+j*sin(delt(n)))

end

for m=1:N1+1

I(m)=Y(5,m)*U(m)

end

S5=U(5)*sum(conj(I))

for m=1:N1+1

for n=1:N1+1

S(m,n)=U(m)*(conj(U(m))-conj(U(n)))*conj(-Y(m,n))%

end

end

end

S

%

用过的例子，若有标点问题，请调试后再用。

---