牛顿拉夫逊法潮流计算matlab程序

牛顿拉夫逊法潮流计算程序

matlab

%电力系统的潮流计算,以下程序参考文献 《电力系统毕业设计》中国水利电力出版社 %(该文献用极坐标下的牛顿——拉夫逊方法实现，在此为了与课本一致做了修改)

%为了计算方便将原来的下标做以下修改: S2 S3 S4 S5 U2 U3 U4 U5 改为S1 S2 S3 S4 U1 U2 %U3 U4 ,即原题的平衡点1就变为现在的平衡点5

%1.形成节点导纳矩阵,

yb55=6.250-18.750j;yb51=-5.000+15.000j;yb52=-1.250+3.750j;yb53=0.000-0.000j;yb54=0.000-0.000j;

yb15=-5.000+15.000j;yb11=10.834-32.500j;yb12=-1.667+5.000j;yb13=-1.667+5.000j;yb14=-2.500+7.500j;

yb25=-1.250+3.750j;yb21=-1.667+5.000j;yb22=1

2.917-38.750j;yb23=-10.000+30.000j;yb24=0.000-0.000j;

yb35=0.000-0.000j;yb31=-1.667+5.000j;yb32=-10.000+30.000j;yb33=12.917-38.750j;yb34=-1.250+3.750j;

yb45=0.000-0.000j;yb41=-2.500+7.500j;yb42=0.000-0.000j;yb43=-1.250+3.750j;yb44=3.750-11.250j;

YB=[yb11 yb12 yb13 yb14 yb15; yb21 yb22 yb23 yb24 yb25 ;yb31 yb32 yb33 yb34 yb35; yb41 yb42 yb43 yb44 yb45 ;yb51 yb52 yb53 yb54 yb55];

%计算各节点功率的不平衡量设U=E+jF ；Y=G+Bj；

E(1)=1.00;E(2)=1.00;E(3)=1.00;E(4)=1.00; F(1)=0;F(2)=0;F(3)=0;F(4)=0; G=real(YB);B=imag(YB); %设S=P+Bj；

S(1)=0.20+0.20i;S(2)=-0.45-0.15i;S(3)=-0.40-0.05i;S(4)=-0.60-0.10i;

P=real(S);Q=imag(S);

k=0;precision=1; N1=4;

while precision > 0.00000001 E(5)=1.06;F(5)=0;

for m=1:N1 for n=1:N1+1

%计算Pi,Qi,设Pi=Pt;Qi=Qt

Pt(n)=(E(m)*(G(m,n)*E(n)-B(m,n)*F(n))+F(m)*(G(m,n)*F(n)+B(m,n)*E(n)));

Qt(n)=(F(m)*(G(m,n)*E(n)-B(m,n)*F(n))-E(m)*(G(m,n)*F(n)+B(m,n)*E(n)));

end

%设P,Q的改变量为dP,dQ

dP(m)=P(m)-sum(Pt); dQ(m)=Q(m)-sum(Qt); end

for m=1:N1 for n=1:N1+1

%计算Hij Nij Jij Lij

H(m,n)=-B(m,n)*E(m)+G(m,n)*F(m);

N(m,n)=G(m,n)*E(m)+B(m,n)*F(m);

J(m,n)=-B(m,n)*F(m)-G(m,n)*E(m);

L(m,n)=G(m,n)*F(m)-B(m,n)*E(m); end end