IEEE14节点潮流程序

luojinwen

clear
2 V' |! h( ?' W0 l. ~2 Lr(14,10)=0.01938;x(14,10)=0.05917;r(10,11)=0.04699;x(10,11)=0.19797;r(10,1)=0.05811;x(10,1)=0.17632;r(14,2)=0.05403;x(14,2)=0.22304;r(10,2)=0.05695;x(10,2)=0.17388;
r(10,14)=0.01938;x(10,14)=0.05917;r(11,10)=0.04699;x(11,10)=0.19797;r(1,10)=0.05811;x(1,10)=0.17632;r(2,14)=0.05403;x(2,14)=0.22304;r(2,10)=0.05695;x(2,10)=0.17388;
r(11,1)=0.06701;x(11,1)=0.17103;r(1,2)=0.01335;x(1,2)=0.04211;r(2,12)=0.00000;x(2,12)=0.23488;r(1,3)=0.00000;x(1,3)=0.20452;r(3,13)=0.00000;x(3,13)=0.17615;
* D8 i4 H. k! z' lr(1,11)=0.06701;x(1,11)=0.17103;r(2,1)=0.01335;x(2,1)=0.04211;r(12,2)=0.00000;x(12,2)=0.23488;r(3,1)=0.00000;x(3,1)=0.20452;r(13,3)=0.00000;x(13,3)=0.17615;
r(1,4)=0.00000;x(1,4)=0.53894;r(3,4)=0.00000;x(3,4)=0.11001;r(4,5)=0.03181;x(4,5)=0.08450;r(12,6)=0.09498;x(12,6)=0.19890;r(12,7)=0.12291;x(12,7)=0.25581;
r(4,1)=0.00000;x(4,1)=0.53894;r(4,3)=0.00000;x(4,3)=0.11001;r(5,4)=0.03181;x(5,4)=0.08450;r(6,12)=0.09498;x(6,12)=0.19890;r(7,12)=0.12291;x(7,12)=0.25581;
$ s0 Y1 \) o: U5 a' J8 K5 D" br(12,8)=0.06615;x(12,8)=0.13027;r(4,9)=0.12711;x(4,9)=0.27038;r(5,6)=0.08205;x(5,6)=0.19207;r(7,8)=0.22092;x(7,8)=0.19988;r(8,9)=0.17093;x(8,9)=0.34802;
r(8,12)=0.06615;x(8,12)=0.13027;r(9,4)=0.12711;x(9,4)=0.27038;r(6,5)=0.08205;x(6,5)=0.19207;r(8,7)=0.22092;x(8,7)=0.19988;r(9,8)=0.17093;x(9,8)=0.34802;

& V- w( \6 B0 F% p3 vy(14,10)=0.02640;y(10,11)=0.02190;y(10,1)=0.01870;y(14,2)=0.02460;y(10,2)=0.01700;y(11,1)=0.01730;y(1,2)=0.00640;y(2,12)=-0.31063;y(1,3)=-0.10999;y(1,4)=-0.05936;
2 y. r5 k$ C1 f* ]* ]! yy(10,14)=0.02640;y(11,10)=0.02190;y(1,10)=0.01870;y(2,14)=0.02460;y(2,10)=0.01700;y(1,11)=0.01730;y(2,1)=0.00640;y(12,2)=0.28951;y(3,1)=0.10757;y(4,1)=0.05752;
# L- z. Q7 {" W3 L% ?- W. i
for m=1:14
    t(m)=0;
7 [; [' P$ q& n    for n=1:14
3 X6 p' I, S0 H! [4 _        t(m)=t(m)+y(m,n);
    end
- V5 `/ _& v, i5 M    if m==4
; k: P% h1 V' e9 _* u. N. c        t(m)=t(m)+0.190;
    end
end
7 z$ o2 |9 @. m% afor m=1:14
    d(m)=0;
! g" f7 J( o. A! E: C0 u" M% |    for n=1:14
( R1 X6 M$ _2 ]/ _  i, f$ v% ~        if m==n
            
0 N9 d6 q# ?  ~, r        elseif (r(m,n)==0)&(x(m,n)==0)
           d(m)=d(m)+0;
/ A- W. w. Q( M. ?4 V6 j5 g! H3 y        else
           d(m)=d(m)+1/(r(m,n)+j*x(m,n));   
) N8 Q) y+ d& U7 s& @6 m4 @) j# S        end
    end
/ b! ?3 ]; c9 a8 h3 Z  Mend
" k+ u8 j6 O3 ~$ Q+ n' R; \2 efor m=1:14
    for n=1:14
        if m==n
           Y(m,n)=j*t(m)+d(m);
        elseif (r(m,n)==0)&(x(m,n)==0)
& i- K" n9 S4 o* C* D& S           Y(m,n)=0;
        else
           Y(m,n)=-1/(r(m,n)+j*x(m,n));
        end
8 w% V& S+ m* Q0 a0 w9 x% x" [" Q+ o' u- b    end
2 X) j* L& T" M1 T. c) fend
Y;
( f8 _# o( O* p; |5 YG=real(Y);B=imag(Y);
%给定节点电压初值和循环次数k=0
$ u$ C# J* D: {$ e2 ]7 fdelt(10)=0;delt(11)=0;u(1)=1.0;delt(1)=0;u(2)=1.0;delt(2)=0;delt(12)=0;u(3)=1.0;delt(3)=0;delt(13)=0;u(4)=1.0;delt(4)=0;u(5)=1.0;delt(5)=0;u(6)=1.0;delt(6)=0;u(7)=1.0;delt(7)=0;u(8)=1.0;delt(8)=0;u(9)=1.0;delt(9)=0;
1 @$ ]5 }) B& l" f/ r) ap(10)=0.183;p(11)=-0.942;p(1)=-0.478;q(1)=0.039;p(2)=-0.076;q(2)=-0.016;p(12)=-0.112;p(3)=0;q(3)=0;p(13)=0;p(4)=-0.295;p(5)=-0.09;q(5)=-0.058;p(6)=-0.035;q(6)=-0.018;p(7)=-0.061;q(7)=-0.016;p(8)=-0.135;q(8)=-0.058;p(9)=-0.149;q(9)=-0.05;
k=0;precision=1;N1=13;
%求节点功率的不平衡量�6�2P和�6�2Q
while precision>0.00001
2 S8 H( F5 o# wu(10)=1.0450;u(11)=1.01;u(12)=1.07;u(13)=1.09;u(14)=1.06;delt(14)=0;q(4)=-0.166+(u(4))^2*0.190;
     for m=1:N1     
6 x, e, Q/ c$ y         if m<=9
             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)));     
) W) l5 @) D( i8 K2 O% M             end
( _& E9 J7 M) r, W             pp(m)=p(m)-sum(pt);qq(m)=q(m)-sum(qt);
         else
             for n=1:N1+1
f(n)=u(m)*u(n)*(G(m,n)*cos(delt(m)-delt(n))+B(m,n)*sin(delt(m)-delt(n)));
1 H$ i  f/ _8 _4 U! T9 ]# n             end
             pp(m)=p(m)-sum(f);qq(m)=0;
6 \8 e- V$ ?, Z6 @% O/ o# ]         end      
     end
%计算雅克比矩阵各元素
( S( W7 o  A$ A# J  dfor m=1:N1
    for n=1:N1+1
         h0(n)=u(m)*u(n)*(G(m,n)*sin(delt(m)-delt(n))-B(m,n)*cos(delt(m)-delt(n)));
( Y2 Y3 t8 Z* N: o6 f" y3 h- ]         n0(n)=-u(m)*u(n)*(G(m,n)*cos(delt(m)-delt(n))+B(m,n)*sin(delt(m)-delt(n)));
         j0(n)=-u(m)*u(n)*(G(m,n)*cos(delt(m)-delt(n))+B(m,n)*sin(delt(m)-delt(n)));
) S% ]3 O! S6 b) ^( S) k9 s         L0(n)=-u(m)*u(n)*(G(m,n)*sin(delt(m)-delt(n))-B(m,n)*cos(delt(m)-delt(n)));
    end
    if m<=9
. F9 R" c+ C7 m# G3 s% ?0 aH(m,m)=sum(h0)- u(m)*u(m)*(G(m,m)*sin(delt(m)-delt(m))-B(m,m)*cos(delt(m)-delt(m)));
/ p9 X2 E) H" N% {( d& |N(m,m)=sum(n0)-2*u(m)^2*G(m,m)+u(m)*u(m)*(G(m,m)*cos(delt(m)-delt(m))+B(m,m)*sin(delt(m)-delt(m)));
J(m,m)=sum(j0)+u(m)*u(m)*(G(m,m)*cos(delt(m)-delt(m))+B(m,m)*sin(delt(m)-delt(m)));
; Q* y  X& s5 jL(m,m)=sum(L0)+2*u(m)^2*B(m,m)+u(m)*u(m)*(G(m,m)*sin(delt(m)-delt(m))-B(m,m)*cos(delt(m)-delt(m)));
( R* U% H# Z& v    else
H(m,m)=sum(h0)- u(m)*u(m)*(G(m,m)*sin(delt(m)-delt(m))-B(m,m)*cos(delt(m)-delt(m)));
$ I* j/ k1 G& z8 ^# T9 O5 @% ?N(m,m)=0;
1 t& x& d2 C4 BJ(m,m)=0;
$ ~6 d7 I) Q0 b) Y6 }: N0 b, U+ WL(m,m)=0;
0 W, V+ `+ s' Z: D4 I: ^: l" v    end
4 L1 B4 z  r4 `end
for m=1:N1
    JJ(2*m-1,2*m-1)=H(m,m);JJ(2*m-1,2*m)=N(m,m);
4 ?( r& A" g+ H# j    JJ(2*m,2*m-1)=J(m,m);JJ(2*m,2*m)=L(m,m);
end
  V7 L% o5 B7 ~, C: D& ]" G
for m=1:N1
* h, w3 }1 E' C9 z    for n=1:N1   
        if     (m<=9)&(n<=9)&(m~=n)                       
H(m,n)= -u(m)*u(n)*(G(m,n)*sin(delt(m)-delt(n))-B(m,n)*cos(delt(m)-delt(n)));
+ p8 D1 _. e4 q+ S, D& @J(m,n)=u(m)*u(n)*(G(m,n)*cos(delt(m)-delt(n))+B(m,n)*sin(delt(m)-delt(n)));
8 e4 X- g8 d% N0 a. f+ vN(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);
( I' g7 f' }, |( J* F$ zJJ(2*m,2*n-1)=J(m,n);JJ(2*m,2*n)=L(m,n);
        elseif ((m==10)&(n<10))|((m==11)&(n<10))|((m==12)&(n<10))|((m==13)&(n<10))
H(m,n)= -u(m)*u(n)*(G(m,n)*sin(delt(m)-delt(n))-B(m,n)*cos(delt(m)-delt(n)));
- I2 s6 W. \( u: D8 z; q' O& mJ(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);
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)=0;JJ(2*m,2*n)=0;
' H! g  X. M$ {2 r        elseif ((n==10)&(m<10))|((n==11)&(m<10))|((n==12)&(m<10))|((n==13)&(m<10))
) u9 K5 E( ?$ t. u; LH(m,n)= -u(m)*u(n)*(G(m,n)*sin(delt(m)-delt(n))-B(m,n)*cos(delt(m)-delt(n)));
9 O' Y1 r; ]5 e, i3 ?# n6 A+ SJ(m,n)=u(m)*u(n)*(G(m,n)*cos(delt(m)-delt(n))+B(m,n)*sin(delt(m)-delt(n)));
( R, |" h6 F9 c7 TJJ(2*m-1,2*n-1)=H(m,n);JJ(2*m-1,2*n)=0;
JJ(2*m,2*n-1)=J(m,n);JJ(2*m,2*n)=0;
        elseif (m>=10)&(n>=10)&(m~=n)
6 B& p* o$ N) [/ bH(m,n)= -u(m)*u(n)*(G(m,n)*sin(delt(m)-delt(n))-B(m,n)*cos(delt(m)-delt(n)));
JJ(2*m-1,2*n-1)=H(m,n);JJ(2*m-1,2*n)=0;
2 Y" {3 A) B# qJJ(2*m,2*n-1)=0;JJ(2*m,2*n)=0;
) R! P+ e8 B3 t        end
    end
end
% 对雅克比矩阵进行修正，使之变成非奇异阵，以便后续求逆计算
b=0;
' G4 ^5 \+ g6 V) H' }for m=1:22
    for n=1:22
        if (m<=18)&(n<=18)
0 G- v* k1 p7 ^, C# J' d. n        A(m,n)=JJ(m,n);
        elseif (m>18)&(n<=18)
        A(m,n)=JJ(m+b,n);
& _& L* x+ y% O# Y2 P        end   
    end
    if m>18
1 X, L( z% ]& c        b=b+1;
    end
! `0 T% E6 g: D3 t+ Cend
1 b: a, a& S1 x- s( t
2 n& B0 f! Q* ?: bb=0;
for n=1:22
6 ^" D$ I* P5 a/ d% @    for m=1:22
1 p5 _# _+ j& K0 m+ Z        if n<=18
        A(m,n)=A(m,n);
( e- U& e' S& B. f0 ^6 @7 @        elseif (m<=18)&(n>18)
. ]$ T# z# R" }( q- R2 C8 k        A(m,n)=JJ(m,n+b);
        end
8 D  U) e- X4 J: F    end
+ q7 f% J- ?) t0 U    if n>18
& [* V9 l5 t$ ]2 V1 F        b=b+1;
    end
% C3 u9 \# h* @+ A  t. v& @end

b=0;
for m=1:22
  Y' r% B% h( D. Z+ A    for n=1:26
        if (m>=19)&(n>=19)
        D(m,n)=JJ(m+b,n);
% }/ u2 ]4 R1 J0 E: l9 f        end
; n. A( M# K/ Y+ g  j" R5 K    end
& e$ B. H/ T% Z    if m>18
& U! J) m* n1 q3 k        b=b+1;
& C6 p. K7 s3 @9 n( K( r2 a/ H* p1 {9 X    end
end

b=0;
1 Z3 H! }- ]  j0 q, e3 Xfor n=1:22
    for m=1:22
, w) Z9 S  x6 U( n- r# R) @1 y1 o        if (m>=19)&(n>=19)
        D(m,n)=D(m,n+b);
        end
: T9 ^+ }! ^& p1 R) L" ]    end
    if n>18
2 a  `4 `8 I" p        b=b+1;
  o; F5 ^* E) \  G# N$ @( o8 d    end
end

$ O( @, z0 ?* Z/ q) pfor m=1:22
5 }" R& D$ p7 [8 ^8 `: m4 W$ Y    for n=1:22
0 u& h+ i$ ]$ q( L6 c2 l        if (m>=19)&(n>=19)
        A(m,n)=D(m,n);
        else
  y4 {' v4 Q2 V        A(m,n)=A(m,n);
7 H" U3 M+ ]3 I9 j4 G: D. ~0 ]9 `        end
5 i5 [$ N/ q1 @! B    end
end
% 形成不平衡量的列矩阵
  ^% ]$ B2 u" X! U3 pfor m=1:N1
    PP(2*m-1)=pp(m);PP(2*m)=qq(m);
end
5 T2 B# H) E& n( d$ f- hC=PP(1:1:18);
b=0;
for m=1:22
    if m<=18
  t  P8 e- a; ]+ o    C=C;
* E  L, w- N" G0 e# D, Z    else
    C(m)=PP(m+b);
0 r/ g+ X) s" w4 V4 [) j8 X7 l    b=b+1;
    end
end
, q7 K5 u% H9 j2 f' L% 解修正方程式，得到修正量
uu=-inv(A)*C';
precision=max(abs(uu));
: E0 L: l4 E% O& H+ K( O% 计算个节点电压新值，即修正后值
( B* w0 M; x6 K# tfor n=1:N1
    if n<=9
    delt(n)=delt(n)+uu(2*n-1);
    u(n)=u(n)+uu(2*n);
    else
  b; J! _. w4 A) g  f" F6 Y    delt(n)=delt(n)+uu(n+9);
/ ~5 `* D: ?$ Y4 s% C+ c    u(n)=u(n);
1 @3 N: G% r  v: b7 k    end
6 ~8 I% v! X- l6 ^; Jend
% 循环次数k加一，带入新值进入下一次迭代
  t2 x9 o$ y1 j' G0 h7 xk=k+1;
, [- L3 J; t; Xend
1 I+ p7 K; m! P* {# |% 迭代完成后，求迭代次数、节点母线电压幅值与相角
k-1,(delt*180/pi)',u'
9 V  H0 w4 T  u- T% 计算各节点的功率和线路功率
for n=1:N1+1
    U(n)=u(n)*(cos(delt(n))+j*sin(delt(n)));
3 j8 K6 s( R- p8 \% Hend
for m=1:N1+1
    for n=1:N1+1
    I(n)=Y(m,n)*U(n);
    end
    S(m)=U(m)*sum(conj(I));
end
S
for m=1:N1+1
; l# q4 W" s& H    for n=1:N1+1
9 X: @" `: t: x. @! D        T(m,n)=U(m)*conj((U(m)-U(n))*(-Y(m,n)))+U(m)*conj(U(m)*j*y(m,n));
    end
+ N$ _6 X6 `3 _end
T

原随风

多谢楼主分享。。。。。。。。。

luojinwen

这个好像算的结果不是很准确，但能学习下思想

lotus_cao

楼主的分享精神不错。

wangwangmoon

现在我们编程的缺点是不能用矩阵的思想全局的看待系统问题，会用很多for循环去实现功能

heming007

太复杂了，看不懂、有14节点编的算节点导纳矩阵的程序么

color

牛人啊，佩服佩服

liulina_1013

这个程序matlab论坛里就有

姗姗来迟156

多谢楼主，可以借鉴来学习啦