设为首页收藏本站|繁體中文 快速切换版块

 找回密码
 立即加入
搜索
查看: 1498|回复: 0

[讨论] 潮流运算牛拉matlab出错

[复制链接]

该用户从未签到

尚未签到

发表于 2015-6-6 20:23:30 | 显示全部楼层 |阅读模式

马上加入,结交更多好友,共享更多资料,让你轻松玩转电力研学社区!

您需要 登录 才可以下载或查看,没有账号?立即加入

×
%本程序的功能是用牛顿——拉夫逊法进行潮流计算
* @# Y2 |8 l- p% B1矩阵:1、支路首端号;2、末端号;3、支路阻抗;4、支路对地电纳
' G1 I# E( H$ A%         5、支路的变比;6、支路首端处于K侧为1,1侧为0( s" `  I9 U8 \% \9 {  e
% B2矩阵:1、该节点发电机功率;2、该节点负荷功率;3、节点电压初始值
0 C. }8 n) q8 s" [8 R! A%         4、PV节点电压V的给定值;5、节点所接的无功补偿设备的容量' _" g# f- I- \& O/ A0 ^
%         6、节点分类标号
; `- N  Z. Q0 M& N) p/ M! Q1 ^- gclear;, \& d% \# ]3 T) A
n=24;%input('请输入节点数:n=');  E* H( V% H* n& `$ }4 E' I; h
nl=38;%input('请输入支路数:nl=');
& }1 W) V( V/ j* j0 C8 ], Y2 Lisb=1;%input('请输入平衡母线节点号:isb=');4 I. r2 B6 d5 Q1 Q: t$ z
pr=0.00001;%input('请输入误差精度:pr=');
& n* D4 V/ G2 o. x! J3 o2 F3 tB1=[1 2 0.0026+0.0139i 0.4611i 1 0;
* i, B# ]2 Z$ q$ N+ D; j) w    1 3 0.0546+0.2212i 0.0572i 1 0;5 m8 z" K- v$ _+ e
    1 5 0.0218+0.0845i 0.0229i 1 0;  I  ]9 p5 I4 H) Y( V
    2 4 0.0328+0.1267i 0.0343i 1 0;0 x! j/ r4 p( J
    2 6 0.0497+0.192i 0.052i 1 0;+ i* \8 h! o: ~% T: {; t
    3 9 0.0308+0.119i 0.0322i 1 0;
4 b4 u; Z* l) R6 t9 t% V, w& U    3 24 0.0023+0.0839i 0 1 0;  Q1 z* t: T! B5 P: B( B6 u
    4 9 0.0268+0.1037i 0.0281i 1 0;
' u# m1 c4 E% a, m; |: N6 h2 j4 q    5 10 0.0228+0.0883i 0.0239i 1 0;* F6 z8 b* d7 e# Q* x- f0 \/ O
    6 10 0.0139+0.0605i 2.459i 1 0;
/ A; p8 W: t- n. i; ]( H0 L3 a1 [- D& d    7 8 0.0159+0.0614i 0.0166i 1 0;  p. ^2 X# }# Q7 j, s
    8 9 0.0427+0.1651i 0.0447i 1 0;6 u! c* h$ K4 H% L+ b2 `  B
    8 10 0.0427+0.1651i 0.0447i 1 0;0 g' V, u8 Q. b* ?. f/ Z
    9 11 0.0023+0.0839i 0 1 0;1 h( F. b" n8 U. p$ v* u
    9 12 0.0023+0.0839i 0 1 0;* N% @8 e8 q. {, V8 u5 i) G
    10 11 0.0023+0.0839i 0 1 0;, @  D- {  C! M6 L: V7 t( k0 i! V
    10 12 0.0023+0.0839i 0 1 0;
$ B. }5 _, \. H: N. {, h- S. ?    11 13 0.0061+0.0476i 0.0999i 1 0;
" @& N2 m) T& A    11 14 0.0054+0.0418i 0.0879i 1 0;
$ M- C& n# l. F, R2 v! U    12 13 0.0061+0.0476i 0.0999i 1 0;, D) g# e+ u$ F  X; C& ^7 p% P
    12 23 0.0124+0.0966i 0.203i 1 0;
0 U" g! I+ x8 z. Y7 J    13 23 0.0111+0.0865i 0.1818i 1 0;
) r* F, f/ u+ Q) |) K, D, f+ @- p4 C    14 16 0.005+0.0389i 0.0818i 1 0;
, _; }: t/ z% u) W    15 16 0.0022+0.0173i 0.0364i 1 0;
- r) R. U! W' e    15 21 0.0063+0.049i 0.103i 1 0;
* Y3 T5 b6 p# O( j/ n- x3 d    15 21 0.0063+0.049i 0.103i 1 0;
6 F) v$ `4 H& L+ I% @2 A6 k2 p    15 24 0.0067+0.0519i 0.1091i 1 0;
; f3 q; J$ m% n; i    16 17 0.0033+0.0259i 0.0545i 1 0;
& f. S* ?' @  L; s" J    16 19 0.003+0.0231i 0.0485i 1 0;. F% R% a$ H# r% G' Q( ~
    17 18 0.0018+0.0144i 0.0303i 1 0;
+ o& [" s1 Y8 Y( f+ z6 c    17 22 0.0135+0.1053i 0.2212i 1 0;
7 {2 i) Z, o5 M    18 21 0.0033+0.0259i 0.0545i 1 0;
9 T5 H9 }$ a# T9 t# D; l: e4 {2 a4 F    18 21 0.0033+0.0259i 0.0545i 1 0;$ Y8 M1 V  }- ]: ~/ Q' H+ d
    19 20 0.0051+0.0396i 0.0833i 1 0;$ d; ^1 X) S  q4 y2 E
    19 20 0.0051+0.0396i 0.0833i 1 0;
: ?4 j) S, L+ \  I; g( A    20 23 0.0028+0.0216i 0.0455i 1 0;1 ^. F( C2 M$ H0 m7 c
    20 23 0.0028+0.0216i 0.0455i 1 0;2 n$ F" b/ w) [' H4 Q1 }7 T& k1 {
    21 22 0.0087+0.0678i 0.1424i 1 0];%input('请输入由支路参数形成的矩阵: B1=');
5 D( I5 h/ M9 q) d5 bB2=[0.1788 1.0157 1 0 0 2;
) F( n4 r& \4 |( L# i7 n8 p1 K    0 0.9127 1 0 0 2;) P4 S$ n+ c% l
    0 1.6935 1 0 0 2;
( X1 y6 p/ h( d7 [    0 0.6958 1 0 0 2;
5 E. W1 v  q* f! x0 n1 V    0 0.6669 1 0 0 2;
9 m( G* h& v( B" @& B    0 1.2796 1 0 0 2;
0 b/ `7 _6 p; ~. K2 J& F    0 1.1748 1.05 1.05 0 1;
1 Q% Z: a1 W  y; @( G) t    0 1.6085 1 0 0 2;3 j; ~% J; b! h2 Z2 e) I1 N2 u! K
    0 1.6465 1 0 0 2;
# O1 Q2 q  ?; C, m5 g    0 1.8345 1 0 0 2;
( ~) J+ [- b$ H/ @    0 0 1 0 0 2;( `: ?& \+ Z6 z
    0 0 1 0 0 2;- K  V5 t" g. n- _/ O: e& C5 \+ c
    6.9529 2.4923 1 0 0 2;
. P; ]- B. T8 K5 B    0 1.8236 1 0 0 2;
3 ^( F- ?7 _. w0 P9 {& B    1.8235 2.9803 1 0 0 2;0 R+ _" n% r+ K+ g
    1.8235 0.9398 1 0 0 2;. C( ?2 c0 @8 T: s7 m* l$ D
    0 0 1 0 0 2;- T2 v/ N0 L/ p0 N
    4.7059 3.1322 1 0 0 2;
) P. F% v; _* Y5 I- u% o, h% I    0 1.7025 1 0 0 2;& R' [+ S2 h! l7 r  S
    0 1.2037 1 0 0 2;5 y& K# `) }" C
    4.7059 0 1 0 0 2;" q7 f, {& t( o( j% F
    0 0 1 0 0 2;
/ v& O2 z( c7 k; n. T9 n! r    7.7647 0 1 0 0 2;; F. S$ ^1 o) u+ x6 Q
    0 0 1 0 0 2];%input('请输入各节点参数形成的矩阵: B2=');9 N+ a1 ^" P* J
Y=zeros(n);e=zeros(1,n);f=zeros(1,n);V=zeros(1,n);sida=zeros(1,n);S1=zeros(nl);
/ m* `2 Z# u, w2 v8 Q: @%-------修改部分------------, M  W5 b/ x! |) s" z  `: j& _9 |
ym=0;
: v+ h2 n8 o3 E8 F! U) WSB=100;UB=220;
2 ?* ~5 B/ A0 R9 v& h%ym=input('您输入的参数是标么值?(若不是则输入一个不为零的数值)');5 B3 c$ I1 d* d6 K, W! r1 A
if ym~=0
! E  `4 R! ^5 v+ W& ~+ J! i    %SB=input('请输入功率基准值:SB=');
2 G+ O+ @1 E* \    %UB=input('请输入电压基准值:UB=');9 J9 f- h6 Y7 e" g
    YB=SB./UB./UB;
7 g4 z9 I9 ]  L- Y, l( |1 w    BB1=B1;
2 t3 l8 V% T  P& u    BB2=B2;
3 h2 d/ M$ Q5 J8 Y# i, J# x% y$ D    for i=1:nl
) p$ i4 E' m& h) G! J        B1(i,3)=B1(i,3)*YB;
, P$ d: v- L* U! T0 q- G- J* L        B1(i,4)=B1(i,4)./YB;
3 o  R* H2 N0 [2 _% N  N9 o& p    end( e6 q. q# ^& T6 d" t  l
    disp('B1矩阵B1=');5 ?5 K# E+ F7 H2 E. h+ t2 ?* e
    disp(B1)" \5 f, \- A) d, W6 q6 Z
    for i=1:n3 E4 b+ }! X# K8 Y
        B2(i,1)=B2(i,1)./SB;
- h, w. F/ Z. d+ T        B2(i,2)=B2(i,2)./SB;
+ u0 J! @1 {( ~! X  }        B2(i,3)=B2(i,3)./UB;. T( W. n9 b/ ^, `0 n( t9 m# p  o
        B2(i,4)=B2(i,4)./UB;3 s4 a9 i+ S; ?; ]4 `4 j
        B2(i,5)=B2(i,5)./SB;
1 u/ }3 U! {4 }( k  ^6 @    end; e) q/ ^) Z. F
    disp('B2矩阵B2=');8 l+ _* _/ `7 g2 ]1 e2 n& B
    disp(B2)
- S2 T3 G$ Y& ]# Oend
" q, c. `" B) }3 O, \% k  Q% % %---------------------------------------------------
1 ?& t% O: \) u# G; Xfor i=1:nl
2 N- E, N3 H5 E& i9 F( h  B) y%支路数
$ ?0 @" q* U5 R    if B1(i,6)==0
; r/ Y8 P. s- k8 V3 R' {%左节点处于低压侧
- M9 q- X5 s4 M' \) s1 j       p=B1(i,1);q=B1(i,2);
5 k. i2 Z; t& z2 G7 _8 F    else
: f% P9 C/ s& u5 Z, l( [       p=B1(i,2);q=B1(i,1);  o) }, z8 X( G
    end" `" d( [7 b/ U/ @
    Y(p,q)=Y(p,q)-1./(B1(i,3)*B1(i,5));
: u. |4 |4 k. V, [%非对角元. y& ]& @  {1 J& W  G9 o
    Y(q,p)=Y(p,q);
( j9 _3 b! ~% ]" N: L- V4 Q5 _    Y(q,q)=Y(q,q)+1./(B1(i,3)*B1(i,5)^2)+B1(i,4)./2;
- e5 I; Q! |- }) o! B%对角元K侧. t7 i2 {: ]* v' {! W1 P
    Y(p,p)=Y(p,p)+1./B1(i,3)+B1(i,4)./2
2 Z0 Q1 v* p0 _$ w# K9 n3 F$ }5 D%对角元1侧
" \" j4 K$ |: F& x( Hend
7 B' E% v' F/ E0 l+ R7 D%求导纳矩阵
9 X& o1 W0 |/ X/ }9 Gdisp('导纳矩阵 Y=');
6 N' C" x+ ^9 ?6 p4 X9 v0 odisp(Y)% R8 l( g: ?5 R; n
%----------------------------------------------------------
2 U7 N4 A7 R" q, y" ?+ zG=real(Y);B=imag(Y);& u; L) }& a) A0 B1 [
%分解出导纳阵的实部和虚部  
* G! ^& W. E- j5 ffor i=1:n
# n; a2 U2 A& b. A! ~%给定各节点初始电压的实部和虚部   
7 R# p" v( J0 {( i& k# s( Z4 |    e(i)=real(B2(i,3));
. U2 o: i" ^8 N" w$ P0 n    f(i)=imag(B2(i,3));6 c$ e% c9 D) t, H. I
    V(i)=B2(i,4);- Z+ _" B% f$ u  S: l
%PV节点电压给定模值
3 H# A! M( K; N: kend
, q( G7 R% U5 u3 U% [8 y" \for i=1:n
) z$ ]% `0 m+ c6 A; d' P%给定各节点注入功率   
" i9 {# h% p8 g    S(i)=B2(i,1)-B2(i,2);
. @6 _8 `8 F, M: Q%i节点注入功率SG-SL  3 w* i" v2 U5 R& |9 B
    B(i,i)=B(i,i)+B2(i,5);4 s- A5 X) X" h  Z# s/ e
%i节点无功补偿量  
' _  n7 B  z: fend
8 s, W0 b& t0 c7 l9 e. X$ Q%===================================================================
# u4 ^' \- z! ^  o% eP=real(S);Q=imag(S);
9 j/ o1 W2 d6 Q$ EICT1=0;IT2=1;N0=2*n;N=N0+1;a=0;" ^" L* ?8 @5 M% F& Z
while IT2~=06 k* X$ ^4 R9 b) t6 M. |( R4 C* q
      IT2=0;a=a+1;
, A9 a9 @. u6 R8 a, |9 ?      for i=1:n
, P4 p% f- X4 v( R, G3 H$ R          if i~=isb) W! s& s( `3 S2 G4 x9 r
%非平衡节点    . h3 j, h2 E0 c( m; r( ]
             C(i)=0;D(i)=0;/ j+ {) T! _. N( G5 J7 z; \
             for j1=1:n. Z! }4 f5 M" d/ e7 ~& v6 V
                 C(i)=C(i)+G(i,j1)*e(j1)-B(i,j1)*f(j1);%Σ(Gij*ej-Bij*fj)
. W- x" A& j. a" k& m                 D(i)=D(i)+G(i,j1)*f(j1)+B(i,j1)*e(j1);%Σ(Gij*fj+Bij*ej)
. o- b# i- Z4 d/ W! `" Y0 m% d             end
+ G" O! t6 f* N             P1=C(i)*e(i)+f(i)*D(i);%节点功率P计算eiΣ(Gij*ej-Bij*fj)+fiΣ(Gij*fj+Bij*ej). |; y2 @+ o2 \* X. q
             Q1=C(i)*f(i)-e(i)*D(i);%节点功率Q计算fiΣ(Gij*ej-Bij*fj)-eiΣ(Gij*fj+Bij*ej)
$ V1 m0 M0 h; V: [: \$ e$ \( X%求P',Q'   
! v. t( y" S  S. a6 y9 X4 j- E8 |9 `9 }             V2=e(i)^2+f(i)^2;
9 |: q% S# X1 X. F) R8 y%电压模平方
, e; e" G- E/ i# M! I3 Q$ _* e/ C     %========= 以下针对非PV节点来求取功率差及Jacobi矩阵元素 =========
' t' r7 o& g0 a2 a6 _5 B             if B2(i,6)~=3
' _. ^7 u1 S# P%非PV节点   
# e2 M3 A1 K( z" N6 |3 K# P1 K                DP=P(i)-P1;8 N8 ]; o; u- {0 z" j2 U' q7 v
%节点有功功率差   
$ q" E- {2 h, G' r                DQ=Q(i)-Q1;             %节点无功功率差  
1 T/ P* {& e( A/ v3 z5 m# P9 l%=============== 以上为除平衡节点外其它节点的功率计算 =================9 R( A5 C* H  W; ?
%================= 求取Jacobi矩阵 ===================
) K3 v7 g, Z2 W7 m2 k" Q                for j1=1:n
  x8 \" u$ v! @# q$ h8 h  g# i7 \                    if j1~=isb&j1~=i- G3 p0 @# J7 @
%非平衡节点&非对角元   
& X9 z& [: c  R' F( ^                       X1=-G(i,j1)*e(i)-B(i,j1)*f(i);
. e) v. J: _# @" q$ W% dP/de=-dQ/df      K7 _! L1 c3 J. m$ \
                       X2=B(i,j1)*e(i)-G(i,j1)*f(i);# ~, u2 @# z, x2 C1 [8 v/ j$ f
% dP/df=dQ/de   
  Z9 r; P: m& \, R5 q- \$ ~* O                       X3=X2;           % X2=dp/df  X3=dQ/de5 O: {% }% t& X$ X, y4 C, u
                       X4=-X1;          % X1=dP/de  X4=dQ/df* I- l3 k! e& g$ c* f
                        p=2*i-1;q=2*j1-1;& O3 y$ j0 R( g( D! \9 p, T
                        J(p,q)=X3;J(p,N)=DQ;m=p+1;! V' k  }/ C7 A6 G2 N! Q9 d( v
                        J(m,q)=X1;J(m,N)=DP;q=q+1;& [9 I! R/ X9 n4 [! l0 x
                        J(p,q)=X4;J(m,q)=X2;
3 X9 ~1 K- L) J4 m  P+ z                    elseif j1==i&j1~=isb  c% {/ ~) F( l& r
%非平衡节点&对角元   
. Y7 `' H, y/ n* x, o$ ~% m7 U                       X1=-C(i)-G(i,i)*e(i)-B(i,i)*f(i);% dP/de
2 ?; J- `3 r! l8 d                       X2=-D(i)+B(i,i)*e(i)-G(i,i)*f(i);% dP/df
# ]3 W! h" G2 w( D, f                       X3=D(i)+B(i,i)*e(i)-G(i,i)*f(i); % dQ/de
/ P$ ]# c' _6 \9 d; [                       X4=-C(i)+G(i,i)*e(i)+B(i,i)*f(i);% dQ/df; I* G8 ~) F$ n; E. S! E' @9 v
                       p=2*i-1;q=2*j1-1;J(p,q)=X3;J(p,N)=DQ;%扩展列△Q
- Y5 r. _4 N) I# z! O                       m=p+1;
' O  {) q* s. T* R: y9 T( F$ a                       J(m,q)=X1;q=q+1;J(p,q)=X4;J(m,N)=DP;%扩展列△P! F# J3 _5 t3 S
                       J(m,q)=X2;1 S# R2 r, ]2 R8 K7 ^
                    end! l# P4 q3 x# ?0 y
                end
! r( d9 U2 |& R2 m% q9 F             else1 i2 e( Y" m3 O/ Q" h$ R
    %=============== 下面是针对PV节点来求取Jacobi矩阵的元素 ===========3 _4 e# E, H! ]
                DP=P(i)-P1;9 }/ p+ z2 o3 ]3 Z7 g" F6 C$ C" e
% PV节点有功误差2 T5 ?+ D8 k, O+ d% Z
                DV=V(i)^2-V2;# O, Z3 K) Y, k0 [% ~3 \4 c
% PV节点电压误差    9 J% U. B, Y& h3 ]
                for j1=1:n5 z( C  B3 F3 `! E. e
                    if j1~=isb&j1~=i
: \+ H! m; v0 d%非平衡节点&非对角元    * w; D+ _! e# p+ L/ D+ m
                       X1=-G(i,j1)*e(i)-B(i,j1)*f(i);   % dP/de
4 v$ N7 x5 \9 z                       X2=B(i,j1)*e(i)-G(i,j1)*f(i);    % dP/df
. g% J. S: P) ]# v% J; \1 j7 P                       X5=0;X6=0;
! M+ E3 K; j9 {& i                       p=2*i-1;q=2*j1-1;J(p,q)=X5;J(p,N)=DV;. v, k: A9 z" G7 I
                       m=p+1;
3 Y0 k. n3 r( \  K' ]% `                       J(m,q)=X1;J(m,N)=DP;q=q+1;J(p,q)=X6;- x5 J! g$ q! `8 l# [+ ~
                       J(m,q)=X2;
# o2 s. Q1 v$ v' d5 f* Y% U                    elseif j1==i&j1~=isb
& }/ S# t" ~! b! @- W1 P%非平衡节点&对角元   
' M2 u- Y! M2 L! T. r                       X1=-C(i)-G(i,i)*e(i)-B(i,i)*f(i);% dP/de1 ]" f2 `  V& {2 ~+ Q
                       X2=-D(i)+B(i,i)*e(i)-G(i,i)*f(i);% dP/df
) T) H' j/ Y/ U' s  y# Z                       X5=-2*e(i);4 o5 `7 c# B( }$ Q
                       X6=-2*f(i);
" L1 _2 Q; T/ \5 z: v                       p=2*i-1;q=2*j1-1;J(p,q)=X5;J(p,N)=DV;
  L$ Q( \/ H3 `% F                       m=p+1;1 O0 D* N0 m+ N2 p5 Q2 o
                       J(m,q)=X1;J(m,N)=DP;q=q+1;J(p,q)=X6;
% D; S  Z6 P; b9 s) t                       J(m,q)=X2;+ h/ X. v( q9 h* P! c! ?
                    end% Y- [6 Z9 Q4 r' B% J
                end
$ ?1 U! l  b7 {3 E; }# V* z% g             end: d. M8 f1 `8 g  j9 x2 n
          end
, M  q; N8 C' R+ x* N) n; ~  [      end/ ?* f5 ?; K$ y* M9 J
%========= 以上为求雅可比矩阵的各个元素 =====================, m) _0 v1 f) e  f8 r5 k% f
      for k=3:N0                    % N0=2*n (从第三行开始,第一、二行是平衡节点)  D; L; w" s' [1 A* p" M$ k/ f1 W
          k1=k+1;N1=N;              % N=N0+1 即 N=2*n+1扩展列△P、△Q* e# a$ ~$ c! P8 \0 k
          for k2=k1:N1
! S! J' R! S3 {3 i% 扩展列△P、△Q    1 ?/ X) d7 G6 q1 i3 J5 [
              J(k,k2)=J(k,k2)./J(k,k);$ |+ y1 Y% |7 V% B% I
% 非对角元规格化    : k) ^4 r% ?% h; I# I+ [5 u
          end
9 R! L. W+ V( X; v9 V2 U1 V          J(k,k)=1;                 % 对角元规格化# h3 ^4 r, c. d' r% c" i! N
          if k~=3                   % 不是第三行: |  T( J3 {1 E
& q6 H! r' d! ]& s1 N4 x: h
%============================================================
, K1 Q4 ~' P9 b: F# \             k4=k-1;# o- t; h. e8 Q- U8 e
             for k3=3:k4
7 f2 l" o7 D4 W" o2 }: Y0 j: ~% 用k3行从第三行开始到当前行前的k4行消去' R  C- }6 J" V) s" H/ Z
                 for k2=k1:N1
+ o4 ~0 I1 P/ s5 w# h; E% k3行后各行下三角元素% ~' k0 H9 B; w# f& l* r4 ^( K
                  J(k3,k2)=J(k3,k2)-J(k3,k)*J(k,k2);%消去运算
. K3 Z" G* [, S3 T; M9 M, l                 end9 Z! k+ h; H( W' R3 n9 u
                 J(k3,k)=0;3 O2 D1 \! I- C9 Z9 s* X: N) u
             end
$ C4 ], V" f  I+ n! Q4 v- b/ C             if k==N0 " [& `5 v0 W6 U0 t% |6 e, b0 N) P
                 break;9 g( r/ g4 ~# i) \' G0 q3 X
             end
) o# \+ @% [1 Q+ W: T* @%==========================================
/ {; T3 ]8 ?$ p: n9 s% _; h             for k3=k1:N0) a/ y$ b8 _% V4 m
                for k2=k1:N1, K1 c, g; F: c) C) p  |
                    J(k3,k2)=J(k3,k2)-J(k3,k)*J(k,k2);%消去运算
, a# W7 U- F' A) ^                end
* t$ J  u- h. ^( d" D0 c                J(k3,k)=0;& |  z: T" A$ S+ f. W. m. S% A3 H0 z# {: F
             end
1 u9 s  x, l9 S6 U. u          else1 W" |& n4 i( Y  N: S
             for k3=k1:N0
- ~# A/ z$ V0 Q7 |2 R# x9 ~" r                 for k2=k1:N1
3 M) A8 x" Y$ o) I                     J(k3,k2)=J(k3,k2)-J(k3,k)*J(k,k2);%消去运算
* E0 |- G1 D& p- ~( ]$ N5 k                 end
8 m" D0 N: P# r) M3 z& X                 J(k3,k)=0;
. r; ~8 `9 }$ R) J' R* i             end7 _) S3 e' J: w8 u, q# X" I  q. x
          end
- H* \' {  U8 ^5 b3 G      end
) x) x' W) l. B2 P5 u      %====上面是用线性变换方式将Jacobi矩阵化成单位矩阵=====
+ k9 X  L" T/ \7 n2 e9 e      for k=3:2:N0-1& q: r; v% l# x" [, R6 A4 _
          L=(k+1)./2;+ i6 u2 M+ K4 F& g$ o& q  ?4 A
          e(L)=e(L)-J(k,N)    %修改节点电压实部
1 K3 I6 H- M  I          k1=k+1;
5 l9 d  x' \! l          f(L)=f(L)-J(k1,N)    %修改节点电压虚部
# K8 M9 I" D3 d% ^; J3 r# `3 G% u      end
  R$ q! B8 L5 I* {& u! s4 f: Q9 i      %------修改节点电压-----------
% ]3 f8 J3 d7 f2 |      for k=3:N0
1 C: b& }1 k* \# [, B* `          DET=abs(J(k,N));; J5 Q' ?4 W: p9 Y8 Q( V$ O4 ]9 z/ V
          if DET>=pr    %电压偏差量是否满足要求
7 i7 m1 A% q' y             IT2=IT2+1; %不满足要求的节点数加1) M3 |& R$ n/ j: y% s/ r  Z
          end3 ]  L5 k0 z, o+ j
      end
, f$ h4 l9 F( n      ICT2(a)=IT2;- g/ x# w' Z: I/ a$ O
      ICT1=ICT1+1;
& {5 m# s) B; t$ x, {  F) o8 u# w0 Wend
' m- w7 f) ^" F' ?%用高斯消去法解"w=-J*V"1 |( N: g5 W' G) D7 g4 w7 z
disp('迭代次数:');
3 I5 r5 W! C) E- _2 m, Sdisp(ICT1);! f) }1 O) A) C% a
disp('没有达到精度要求的个数:');" I& A5 B+ J# d- V
disp(ICT2);
  }' A) d3 g1 ^for k=1:n
2 l+ X" ^2 t0 f: _# ~5 I    V(k)=sqrt(e(k)^2+f(k)^2);
4 s& t* ?# Y8 Y: `; w  C    sida(k)=atan(f(k)./e(k))*180./pi;/ H9 g8 B, J  t& o& H$ S- g- P8 O" F
    E(k)=e(k)+f(k)*j;0 l! D1 b7 }8 ~/ X- c  X
end
) T1 p! f1 M9 _3 ]* M8 q( ]%=============== 计算各输出量 ===========================; J$ v/ [, j5 F" L" S
disp('各节点的实际电压标幺值E为(节点号从小到大排列):');9 V$ m5 {# f, E# L
disp(E);
: \. Q  o  f: E1 HEE=E*UB;! W/ V# ^5 }4 j  r5 C' ]$ K
disp(EE);- u; k1 S1 y7 C3 F: L# O% |
disp('-----------------------------------------------------');
, ^6 Q/ I( C/ z9 i1 Z( M8 s$ ], Qdisp('各节点的电压大小V为(节点号从小到大排列):');
* e; a& m( R" Q9 y9 V+ Wdisp(V);
5 |  `# T8 f6 v# L  C/ u6 N+ VVV=V*UB;
  Y+ w  g$ W6 Hdisp(VV);
' V5 C* Q0 _* edisp('-----------------------------------------------------');' h& }3 }- n' A+ @
disp('各节点的电压相角sida为(节点号从小到大排列):');  g7 |, |9 V! b- p$ G
disp(sida);; Q6 b7 ?6 C, V9 Y/ ~
for p=1:n' j- z7 Y3 k' d1 n) l& c. x
    C(p)=0;
3 A5 _' F+ P; U    for q=1:n
7 l' f' W" b- @' D* x        C(p)=C(p)+conj(Y(p,q))*conj(E(q));1 x1 D6 m; Z7 h' H' o9 ^/ ?
    end: p: v' J  n9 {9 ?9 V
    S(p)=E(p)*C(p);" W1 |8 w) X7 W# C% e
end
/ x' ^! y# W1 |; }( L  u; v( ndisp('各节点的功率S为(节点号从小到大排列):');  ^, V: P- k) o* ]9 g6 F
disp(S);
# M4 K) q5 h# q( I: Ydisp('~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~');0 K0 t" o/ q9 p- \( x+ k- `1 A6 J
SS=S*SB;3 r  ^; E- u5 w- k
disp(SS);; S0 j+ f# Q  E! `" P  R
disp('-----------------------------------------------------');7 b% }+ b1 K, v8 w' y( ?4 r# V8 G
disp('各条支路的首端功率Si为(顺序同您输入B1时一致):');+ R4 G: P, j  @/ E
for i=1:nl) p3 F& k9 j  T# O9 t$ S
       p=B1(i,1);q=B1(i,2);# o" P6 r+ U: @/ f) \" G  k* N
       if B1(i,6)==0
5 B9 R7 M& [3 f5 a8 Z. E            Si(p,q)=E(p)*(conj(E(p))*conj(B1(i,4)./2)+(conj(E(p)*B1(i,5))-conj(E(q)))*conj(1./(B1(i,3)*B1(i,5))));. r* C  U% f$ ]* P% k
            Siz(i)=Si(p,q);- q/ t. g& u2 H  h
        else* W# b: h1 J) h( \0 ?' j
            Si(p,q)=E(p)*(conj(E(p))*conj(B1(i,4)./2)+(conj(E(p)./B1(i,5))-conj(E(q)))*conj(1./(B1(i,3)*B1(i,5))));
9 M/ h, c. y; B            Siz(i)=Si(p,q);; D! \& i# c! r3 p! I/ b
        end
6 {) Q" g/ f1 O- D- S- @8 t    disp(Si(p,q));
/ K3 ~8 B1 _% s& ^1 G$ o" X    SSi(p,q)=Si(p,q)*SB;
% W0 N/ ~) s. @( ~" K7 S' S    ZF=['S(',num2str(p),',',num2str(q),')=',num2str(SSi(p,q))];  |( e7 j8 N; k6 @3 B& e, a* x
    disp(ZF);
/ b0 {, w; T! r& v    %disp(SSi(p,q));. m7 a1 M# y! D# m
    disp('-----------------------------------------------------');0 `2 n9 ^0 `; n2 F5 Y
end8 V+ n* [9 k6 I$ j4 \2 @
disp('各条支路的末端功率Sj为(顺序同您输入B1时一致):');* \$ r4 b0 R  H
for i=1:nl
  P1 Z) {5 @  r3 r    p=B1(i,1);q=B1(i,2);0 N- a: n& u) O1 a
    if B1(i,6)==05 l+ S2 [( [4 f4 t9 M+ @* h7 s- I
        Sj(q,p)=E(q)*(conj(E(q))*conj(B1(i,4)./2)+(conj(E(q)./B1(i,5))-conj(E(p)))*conj(1./(B1(i,3)*B1(i,5))));. o4 T% K1 H6 I9 {# \! @  h0 d6 x6 m
        Sjy(i)=Sj(q,p);9 n1 K2 r- \8 v0 ~! m" \# r4 w
    else, s1 k6 I  i- s
    Sj(q,p)=E(q)*(conj(E(q))*conj(B1(i,4)./2)+(conj(E(q)*B1(i,5))-conj(E(p)))*conj(1./(B1(i,3)*B1(i,5))));
4 r4 ?! D6 ~+ ?/ _( C        Sjy(i)=Sj(q,p);
7 z, p. K9 G7 F2 Y7 y    end
+ l7 d  c* ?4 N" h+ @    disp(Sj(q,p));
, t% g- B3 l+ w, F. d    SSj(q,p)=Sj(q,p)*SB;
) p# L& j8 D1 h) g, I, w6 |; S* _6 u    ZF=['S(',num2str(q),',',num2str(p),')=',num2str(SSj(q,p))];
4 @9 |8 k/ W, i6 q, E    disp(ZF);4 v; d4 A- o* }( D# R
    %disp(SSj(q,p));! J; v0 B9 x. w2 ^, H, O& \+ M0 ~
    disp('-----------------------------------------------------');. D. U0 c/ [9 C% h9 D. ?
end1 }& c5 x6 M. m! Y% L6 C4 Q
disp('各条支路的功率损耗DS为(顺序同您输入B1时一致):');
1 n# H  v' Y$ F1 g8 T+ q5 lfor i=1:nl
/ r% E% _& g. F: P* g        p=B1(i,1);q=B1(i,2);; h' u1 j8 s, T2 W* H% B! y. p
    DS(i)=Si(p,q)+Sj(q,p);1 i! D$ M" p+ L, [- e  ?5 q5 T
    disp(DS(i));& a) n1 Z0 Z! s) l& c8 ~; I
    DDS(i)=DS(i)*SB;
  v7 a& w3 u$ m/ B  ^( Z    ZF=['DS(',num2str(p),',',num2str(q),')=',num2str(DDS(i))];% f4 x: @/ ]. ~+ J! B) M1 b
    disp(ZF);2 T, p* D5 p6 G6 f
    %disp(DDS(i));" v3 {! H1 C* ?8 i
    disp('-----------------------------------------------------');# ~, @/ i# ~( X  k6 w: L
end! Z. C' n# y4 W/ M; h4 p- w# q
figure(1);
$ Q- ^5 l3 g0 A' hsubplot(2,2,1);" X1 X- V5 }- I  g( i
plot(V);3 D0 k1 j! J0 B2 a! u/ S( p% J5 `
xlabel('节点号');ylabel('电压标幺值');
  r: K/ ?. t+ C" b1 O8 w" k$ P. Mgrid on;% U+ O  I" s  W+ L; n2 W" U/ ^
subplot(2,2,2);7 W. _" i# a6 q. a
plot(sida);
2 w. d+ ~6 u: T) b3 L9 Rxlabel('节点号');ylabel('电压角度');
' F: Y0 e$ J6 w  @  w, hgrid on;/ F1 M$ `7 z9 Y- j' i& v
subplot(2,2,3);- o2 e' l7 \  ?4 ^7 V
bar(S);, w( t* b$ t; r0 E* K5 K& g" Q
ylabel('节点注入有功');
2 k& Q' }, @; s$ V& T$ G" ngrid on;% d0 L( r2 k( G) E
subplot(2,2,4);
7 @( e) T5 l2 ?. s# Nbar(Siz);
3 D, b' t4 Q) Qylabel('支路首端无功');% k! f6 [7 `1 b
grid on;
"真诚赞赏,手留余香"
还没有人打赏,支持一下
楼主热帖
帖文化:【文明发帖 和谐互动】 社区精神:【创新、交流、互助、共享】
您需要登录后才可以回帖 登录 | 立即加入

本版积分规则

招聘斑竹

小黑屋|手机版|APP下载(beta)|Archiver|电力研学网 ( 赣ICP备12000811号-1|赣公网安备36040302000210号 )|网站地图

GMT+8, 2025-4-4 19:19

Powered by Discuz! X3.5 Licensed

© 2001-2025 Discuz! Team.

快速回复 返回顶部 返回列表