|
马上加入,结交更多好友,共享更多资料,让你轻松玩转电力研学社区!
您需要 登录 才可以下载或查看,没有账号?立即加入
×
%本程序的功能是用牛顿——拉夫逊法进行潮流计算
* @# 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; |
|