IEEE14节点潮流程序

luojinwen

clear
& J0 p( w! c6 Hr(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;
2 `' h7 B6 D" vr(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;
- D) S9 a, i+ f: C3 M: V$ A$ M; Zr(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;
8 u5 B! Q7 g: fr(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;
& Q1 ~; `' g/ ~& `7 X7 c# yr(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;
, D* o; A; O' J' {r(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;
. C3 ^  o4 w; ?2 O% S" `0 _, qr(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;

y(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;
: l" m- N( o0 V7 S8 Wy(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;

for m=1:14
, R. [/ V4 a8 ?* @! q    t(m)=0;
    for n=1:14
( {% \/ v% `7 ^6 \! `        t(m)=t(m)+y(m,n);
( E5 |- V5 z# C& f. O0 I9 Z2 L    end
    if m==4
        t(m)=t(m)+0.190;
: B- E/ w  e' q3 i: Y    end
end
' W& d( O  V  f8 Hfor m=1:14
5 \' s9 u% G3 S6 B7 f    d(m)=0;
    for n=1:14
        if m==n
            
# t9 L* m0 n) ?; J' t" q* }        elseif (r(m,n)==0)&(x(m,n)==0)
           d(m)=d(m)+0;
        else
4 ]0 T6 I+ m- z; e6 B# V           d(m)=d(m)+1/(r(m,n)+j*x(m,n));   
        end
    end
end
9 a: A- t* Y+ _& N+ ~for m=1:14
    for n=1:14
; C% D2 O- [: V$ f# Y        if m==n
           Y(m,n)=j*t(m)+d(m);
        elseif (r(m,n)==0)&(x(m,n)==0)
% g; R7 q$ g" t) ~8 [0 K# v) J+ i           Y(m,n)=0;
# k4 t1 W% q& d# e1 o        else
0 g' D) E; e4 G7 a           Y(m,n)=-1/(r(m,n)+j*x(m,n));
        end
    end
3 F3 L/ n: m% L1 e1 Gend
: ?+ h( l( p+ p8 I3 \4 `Y;
G=real(Y);B=imag(Y);
. K6 d% a% {/ Z8 h0 d% ~: \%给定节点电压初值和循环次数k=0
delt(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;
$ p4 q, _( l% P" wp(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
u(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     
         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)));
1 v1 H7 M" ?# v8 Q2 b0 n5 |qt(n)=u(m)*u(n)*(G(m,n)*sin(delt(m)-delt(n))-B(m,n)*cos(delt(m)-delt(n)));     
0 x& _5 `9 I6 v- A$ _& g             end
  Z4 J7 J2 d4 c( U6 l4 h. l$ l             pp(m)=p(m)-sum(pt);qq(m)=q(m)-sum(qt);
/ |$ {1 _# G7 Q$ J         else
             for n=1:N1+1
- Z+ Y' _* A* s, of(n)=u(m)*u(n)*(G(m,n)*cos(delt(m)-delt(n))+B(m,n)*sin(delt(m)-delt(n)));
             end
             pp(m)=p(m)-sum(f);qq(m)=0;
         end      
$ \1 k, F& G$ Z5 B. Y     end
( ?5 r: X; l2 ~6 t  g%计算雅克比矩阵各元素
for m=1:N1
7 ]5 a" Y! h9 d3 O9 H% G    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)));
% e* l8 A! M: k( t, d         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)));
' `( k# i  A0 A4 ]0 V: u  o' N# f         L0(n)=-u(m)*u(n)*(G(m,n)*sin(delt(m)-delt(n))-B(m,n)*cos(delt(m)-delt(n)));
    end
* E% {8 `6 b& N3 o* f: l    if m<=9
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)));
$ w: h6 N# x' O: rN(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)));
; R$ F9 e; A; [! W: E# L; cJ(m,m)=sum(j0)+u(m)*u(m)*(G(m,m)*cos(delt(m)-delt(m))+B(m,m)*sin(delt(m)-delt(m)));
L(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)));
9 T+ @8 X$ [3 \7 X& ~: i0 d    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)));
- k0 X- V4 S; z2 M0 |N(m,m)=0;
+ h7 g8 ^2 _/ J+ w8 m( ~J(m,m)=0;
/ p: b0 z# V/ m$ t' kL(m,m)=0;
    end
( X; T+ e( q8 H6 u7 Eend
for m=1:N1
    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);
end
' \" u" e& A" G/ ]2 D
6 W# O7 o2 g. f! s. xfor m=1:N1
' e4 r. G; h- d& @9 e    for n=1:N1   
        if     (m<=9)&(n<=9)&(m~=n)                       
$ }# X9 F+ |4 Q/ mH(m,n)= -u(m)*u(n)*(G(m,n)*sin(delt(m)-delt(n))-B(m,n)*cos(delt(m)-delt(n)));
- {% N# U: U9 l4 [2 H' ^J(m,n)=u(m)*u(n)*(G(m,n)*cos(delt(m)-delt(n))+B(m,n)*sin(delt(m)-delt(n)));
: q! H; C0 J+ |# g# U: D) ]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);
/ J# X& `4 N, }- I/ S6 hJJ(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))
- j! Y6 p+ s! A9 Z/ a1 vH(m,n)= -u(m)*u(n)*(G(m,n)*sin(delt(m)-delt(n))-B(m,n)*cos(delt(m)-delt(n)));
/ L6 P, Y6 V. q1 ~J(m,n)=u(m)*u(n)*(G(m,n)*cos(delt(m)-delt(n))+B(m,n)*sin(delt(m)-delt(n)));
3 A' P- D) z2 r% P" @: y) {+ JN(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;
* u3 u3 F4 [% [/ W3 v        elseif ((n==10)&(m<10))|((n==11)&(m<10))|((n==12)&(m<10))|((n==13)&(m<10))
$ t& D& p& A+ y' OH(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)));
4 `; @8 y! {1 E: n) Z( R- ]9 ^( 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)
H(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;
JJ(2*m,2*n-1)=0;JJ(2*m,2*n)=0;
        end
. i# e( E1 I; G+ A' O    end
# N. `) a! v) r1 _7 I4 y9 rend
) i5 l: f/ U. f( B% 对雅克比矩阵进行修正，使之变成非奇异阵，以便后续求逆计算
b=0;
for m=1:22
) g! q5 v2 Q$ N0 l& o& J( j    for n=1:22
) b' U' O& j+ K        if (m<=18)&(n<=18)
* o% y+ [$ p1 `* ^1 R7 y        A(m,n)=JJ(m,n);
: _7 G: A1 E3 b& L: ~        elseif (m>18)&(n<=18)
        A(m,n)=JJ(m+b,n);
; k6 V, Z+ p1 V! g- }% \        end   
    end
6 ]' o. J. P7 c9 y1 T3 w    if m>18
8 p' V# D  d3 g$ ^2 @* s        b=b+1;
    end
( }+ V6 z" g0 _' @6 i6 f+ Iend
0 S* I# x; f2 c- Z
# _* j' ^: ~2 ~! j' c1 ^b=0;
for n=1:22
    for m=1:22
        if n<=18
, [1 Y4 C- Y1 @: G+ S* L/ i        A(m,n)=A(m,n);
! o# {3 v' {& q. q        elseif (m<=18)&(n>18)
        A(m,n)=JJ(m,n+b);
) m) |* h& `$ @+ j( r% f        end
4 l/ u. F6 g8 P' d' ], w3 k    end
0 H' t% W0 r1 k- }    if n>18
: M7 {! m& H9 ?        b=b+1;
    end
end
0 w, Q9 Z4 d3 {$ H9 L0 |
# ?1 J  p7 }* h; V7 yb=0;
1 L9 O$ q) _! i) m1 i1 nfor m=1:22
& ?! e2 ]9 L8 Q3 W; R" X    for n=1:26
" c: T/ c5 K7 d+ X8 U. b- e3 r8 U* s        if (m>=19)&(n>=19)
" V/ G2 [5 A6 u        D(m,n)=JJ(m+b,n);
6 d/ s) \; U# L) H        end
    end
+ T' H  n% `2 O8 y; O/ ?    if m>18
        b=b+1;
+ s# D" M2 H0 u    end
2 W1 f3 c6 t) f7 u( |% i4 |: qend
. Z1 U, j9 `$ n+ _0 y4 K4 j
b=0;
( A- d& ?' j+ t) Vfor n=1:22
. r0 q7 c5 I/ R2 f0 R( h0 P, T    for m=1:22
, f0 S0 w5 F* O, `        if (m>=19)&(n>=19)
( r6 B9 A$ |6 y( d0 e4 D" _+ E& K        D(m,n)=D(m,n+b);
8 U0 V% U/ w! y: S, p" z. t        end
2 _- B* E7 J# Z3 t8 E: T0 y) g    end
" a2 t/ I1 w2 _7 i% \    if n>18
6 T' `+ O& Q/ T: V6 M        b=b+1;
    end
% A6 g- I7 q5 `- ]( W# R; dend
/ y1 D# _* a6 u6 J/ b
for m=1:22
    for n=1:22
        if (m>=19)&(n>=19)
7 e6 O- y9 p" j  Y  O# }/ ?        A(m,n)=D(m,n);
        else
+ g4 @' v% u8 L) |        A(m,n)=A(m,n);
: u$ G% ]) n& k) }' q" R) f        end
    end
  H3 H1 l& }& s- f6 `end
% 形成不平衡量的列矩阵
for m=1:N1
    PP(2*m-1)=pp(m);PP(2*m)=qq(m);
: V6 N, g0 C( c% b2 Q& P0 |" Oend
; y( `% Z# |) c0 _# IC=PP(1:1:18);
5 |$ @3 q* j) db=0;
for m=1:22
    if m<=18
    C=C;
  C0 E! |* q: v$ a" f    else
/ T! T8 k1 _, k* b, `- o% ~    C(m)=PP(m+b);
    b=b+1;
: W  D$ b0 |; T/ B% m    end
end
% 解修正方程式，得到修正量
uu=-inv(A)*C';
* r0 d( o. m- }: u8 i2 vprecision=max(abs(uu));
% 计算个节点电压新值，即修正后值
8 A# U' V, ], @1 u7 ^1 Efor n=1:N1
    if n<=9
" [9 [' p' e4 G6 L3 `    delt(n)=delt(n)+uu(2*n-1);
5 O6 O$ y/ p+ @$ P, U    u(n)=u(n)+uu(2*n);
    else
4 Z, F) n$ ^! ^6 L! d7 o    delt(n)=delt(n)+uu(n+9);
5 A. L- H# J; m8 {( o3 b    u(n)=u(n);
) ^9 F8 e( ^1 I% q    end
end
% 循环次数k加一，带入新值进入下一次迭代
k=k+1;
! V4 ]% d4 e+ r; Q! o( P8 y; \end
, V( x$ R+ {6 e7 E% R, {% 迭代完成后，求迭代次数、节点母线电压幅值与相角
/ i) ]% C0 B5 G' P( ?k-1,(delt*180/pi)',u'
% 计算各节点的功率和线路功率
for n=1:N1+1
! H' z( D7 Q% q    U(n)=u(n)*(cos(delt(n))+j*sin(delt(n)));
; P3 B* {6 Z: x& U9 l6 f3 _" r8 Wend
for m=1:N1+1
$ L: i$ d( ?5 G. I% `2 D    for n=1:N1+1
    I(n)=Y(m,n)*U(n);
    end
6 Q2 e% t3 A+ }  |% Z3 ?, R0 A    S(m)=U(m)*sum(conj(I));
, P6 c% m1 X5 ~9 i, f5 G! c6 Oend
S
: E8 v" G+ O, ^  O- S8 I1 ffor m=1:N1+1
4 s+ s: B: _% F# i- V! Z    for n=1:N1+1
        T(m,n)=U(m)*conj((U(m)-U(n))*(-Y(m,n)))+U(m)*conj(U(m)*j*y(m,n));
    end
+ |, F3 R) W  |! o- R  Vend
T

原随风

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

luojinwen

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

lotus_cao

楼主的分享精神不错。

wangwangmoon

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

heming007

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

color

牛人啊，佩服佩服

liulina_1013

这个程序matlab论坛里就有

姗姗来迟156

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