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

1、牛顿拉夫逊法潮流计算m atlab程序作者:日期：%电力系统的潮流计算，以下程序参考文献电力系统毕业设计中国水利电力出版社% 该文献用极坐标下的牛顿一一拉夫逊方法实现，在此为了与课本一致做了修改)%为了计算方便将原来的下标做以下修改：S2 S3 S4 S5 U2U3 U4 U5改为S1 S2S3 S4 U1 U 2%U3 U4，即原题的平衡点1就变为现在的平衡点5%1.形成节点导纳矩阵，y b55=6. 2 50 - 18.750j;y b 51=- 5 . 0 00+1 5 .000 j ;y b 52=- 1.2 5 0+3.750j;yb53 = 0 .000-0. 0 00j； y b

2、54= 0.0 00 -0.000j；y b 1 5 =-5.0 0 0+15. 0 0 0j;yb11=10.834- 3 2. 5 00j; yb12=- 1 .6 6 7+ 5.000j; y b 1 3 = 1.667+5. 0 00j; yb14=-2.5 0 0 +7.5 0 0 j;yb25=-1.250+ 3.7 5 0 j ;yb 2 1=- 1 . 667+ 5 . 0 00j;yb22= 1 2.917-38.750j; yb23 = 10. 000 + 30.0 0 0j ;yb24= 0 .0 00-0. 0 0 0j;yb35=0.000-0. 0 0 0 j;y

3、b 3 1 = - 1 . 66 7 +5.000j;yb3 2 =-1 0 .000+3 0 .000j;yb 3 3=12.9 1 7-38. 75 0 j;yb34= - 1. 25 0 +3.7 5 0 j;yb4 5 =0. 0 0 0-0.000j ; yb41=- 2.5 0 0+7.5 0 0j;yb42=0 . 0 00- 0 . 0 0 0j;yb4 3 =-1.250+ 3 .7 5 0 j； y b 44= 3 . 7 5 0-11.250j;YB =yb11 y b1 2 y b1 3 yb 14 ybi 5; yb21 y b 22 yb 23 yb 24 yb25

4、 ;yb3 1 yb 3 2 yb 3 3 yb34 yb35 ; y b 41yb42 yb43 yb4 4 yb45 ; y b 5 1 yb5 2 yb5 3 y b 54 y b 5 5 ;%计算各节点功率的不平衡量设U=E+jF ;Y= G+Bj;E(l)=1. 0 0 ；E(2 )=1.0 0 ;E( 3 )= 1 .0 0 ；E(4 )=1 . 00;F(1) = 0;F(2)=0; F (3)= 0 ; F(4 )= 0 ;G= r eal(YB);B= i ma g (Y B);%设 S= P+Bj;S (1)= 0 .20+0.20i;S (2)=-0. 450. 15i;

5、 S (3)=-0 . 40-0. 0 5 i; S (4 )=-0. 6 0- 0 .10i;P=r e al ( S );Q= i m a g ( S);k = 0; p r e c is i o n= 1 ;N1=4;whil e pre cision > 0.00 0 000 0 1E (5)=1.06;F (5)=0 ;for m = 1 :N1for n = 1 ： N1+1%计算P 旧1设i=Pt;Q i =QtP t( n )=(E( m )*( G (m,n)*E(n ) - B (m,n)*F(n) ) +F (m )*(G(m,n) *F (n)+B(m,n) *E

6、 ( 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);en d%设，Q的改变量为dP,dQd P (m)=P (m) sum(P t); dQ (m)=Q (m)-sum (Q t);en dfor m =1:N1for n=1 : N1+ 1%计算 Hi j Nij J i j LijH (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 ,

7、n)*F(m)-G(m, n ) * E(m);L(m,n)= G (m,n) *F (m)-B (m,n) *E(m ); end endfor m= 1 :N1 for n= 1:N1 + 1Bi(n)= G (m,n) *F(n) + B(m,n)*E(n);A i (n)=G( m ,n) * E -B(m , n)*F(n ); end%计算 Hi i ,Ni i ,J i i ,Lii,由公式 4-44 b 左侧公式实现，sum(A i), s um(Bi)用 于实现公式中的sigerma从j至U n的求和；H(m,m尸sum(Bi)-(B( m,m)*E(m) +G(m,m)*F

8、(m) )+2*G(m,m) * F(m);N(m,m)=s u m( A i) -(G ( m,m)*E(m)-B( m , m)*F(m) ) + 2*G (m,m)*E (m);J(m ,m )=-2* B (m , m )*F(m ) + s um(Ai)-(G(m,m)*E(m ) -B(m,m) * F ( m);L (m,m) = 2*B(m,m) * E(m) ( s um(Bi) (B ( m,m) * E ( m)+G(m,m)* F(m);end%设雅可比矩阵为JJ,以下语句用来实现雅可比矩阵中对角线上元素H N J L的排列for m=1:N1JJ(2*m-1, 2 *

9、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 );end%以下语句用于实现雅可比矩阵非对角线上元素的排列for m=1： N1fo r n =1 : N1 if m= = nelseH(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)

10、 *E (m);J J( 2 *m-1 , 2*n-1) = H (m, n );J J ( 2 * m-1,2 * n尸N (m, n);JJ(2*m,2 * n-1)=J(m, n );J J(2*m,2 * n)=L (m ,n);en d ende nd%设由P,Q的改变量组成的8X1矩阵为PQ,由E,F的改变量组成的 8X1矩阵为dU for m=1:N 1PQ (2 *m-1)=dP( m) ;PQ(2*m尸 dQ(m);e ndd U = i nv(JJ)*PQ'precisi o n = m ax (ab s (dU);for n=1:N1F(n)=F( n ) +dU

11、(2* n -1);E(n) = E (n)+dU( 2 * n );e ndfor n = 1 : N1+1U (n)=E (n)+ (F ( n ) *j ;endk = k+ 1 ;k-1 , d U= d U', PQ,Uend%计算S (5)，也就是题目中的S1,即平衡节点功率 for m=1 : N1+1I (m)= Y B( 5 ,m) *U(m);e ndS(5) = U(5) * s u m (c o n j (I)%设网络总损耗为Ss计算输电效率 eff i c i enc yfor m = 1: N 1+ 1S0(m )= S (m); P(m尸r e al(S(m); e ndS s=su m(S0)efficien