|
马上加入,结交更多好友,共享更多资料,让你轻松玩转电力研学社区!
您需要 登录 才可以下载或查看,没有账号?立即加入
×
%本程序的功能是用牛顿——拉夫逊法进行潮流计算1 G" u$ G% \0 B* {9 e7 F
% B1矩阵:1、支路首端号;2、末端号;3、支路阻抗;4、支路对地电纳# v4 p; n# ~. L6 z4 Q
% 5、支路的变比;6、支路首端处于K侧为1,1侧为0
" H% r% I; }2 X$ b% B2矩阵:1、该节点发电机功率;2、该节点负荷功率;3、节点电压初始值9 x; C- V( ~8 z
% 4、PV节点电压V的给定值;5、节点所接的无功补偿设备的容量+ K+ u# S8 L! W
% 6、节点分类标号
" t3 _9 R% Q, q2 W& @! ?9 zclear;& q2 g1 ~* D k* N. h9 @
n=24;%input('请输入节点数:n=');
/ b6 H, P. u, K# v/ Mnl=38;%input('请输入支路数:nl=');
4 o" ~% L' ~: A( p( @# kisb=1;%input('请输入平衡母线节点号:isb=');
9 t I0 { B3 {pr=0.00001;%input('请输入误差精度:pr=');- q7 Q' n+ ~/ ?8 G
B1=[1 2 0.0026+0.0139i 0.4611i 1 0;
' p Q; S1 w% Y4 f 1 3 0.0546+0.2212i 0.0572i 1 0;
! |, _- M8 [, q% H& V# E8 I 1 5 0.0218+0.0845i 0.0229i 1 0;
+ O, h) P5 [0 q$ A* k! Q 2 4 0.0328+0.1267i 0.0343i 1 0;
2 x7 a% y: J( O, C 2 6 0.0497+0.192i 0.052i 1 0;
9 n9 W5 Y4 z; ? ^ 3 9 0.0308+0.119i 0.0322i 1 0;- f" d8 w. j' U+ o' I8 P
3 24 0.0023+0.0839i 0 1 0;
0 d2 c+ v/ Q/ V. V- E( D 4 9 0.0268+0.1037i 0.0281i 1 0;# ?8 V7 \5 D( Y4 v5 x3 x
5 10 0.0228+0.0883i 0.0239i 1 0;6 g& y# l7 A, e
6 10 0.0139+0.0605i 2.459i 1 0;+ M3 N3 l# O: ?$ b
7 8 0.0159+0.0614i 0.0166i 1 0;
/ j) ^) \6 F+ k! a 8 9 0.0427+0.1651i 0.0447i 1 0;
& R0 L! O! R/ i' t8 n 8 10 0.0427+0.1651i 0.0447i 1 0;0 ^! T* C1 N$ o
9 11 0.0023+0.0839i 0 1 0; w5 S+ X: K2 W H9 [6 b
9 12 0.0023+0.0839i 0 1 0;; q; H" b' h x: S. f
10 11 0.0023+0.0839i 0 1 0;
7 q+ x! h4 L7 U2 ]( d 10 12 0.0023+0.0839i 0 1 0;/ }* s, H7 y, _$ _
11 13 0.0061+0.0476i 0.0999i 1 0;' ]% ^" v* O7 `$ [$ M
11 14 0.0054+0.0418i 0.0879i 1 0;! m# Q8 z. M% F$ a8 |" N
12 13 0.0061+0.0476i 0.0999i 1 0;: i8 R o. [" J- y
12 23 0.0124+0.0966i 0.203i 1 0;
: _, C3 ~6 R1 t 13 23 0.0111+0.0865i 0.1818i 1 0;" `* u: b5 L0 b# p) g
14 16 0.005+0.0389i 0.0818i 1 0;' W Y. G# j$ o9 j' m9 G7 B; V
15 16 0.0022+0.0173i 0.0364i 1 0;% L' l, g' Q: y Z
15 21 0.0063+0.049i 0.103i 1 0;8 Y) \9 B4 ~/ v. C
15 21 0.0063+0.049i 0.103i 1 0;
$ m6 K# F& N$ B# `! J/ r 15 24 0.0067+0.0519i 0.1091i 1 0;
% w. P# R7 I" j) f J 16 17 0.0033+0.0259i 0.0545i 1 0;
& N% E+ V4 k" |# [5 ~" ^ 16 19 0.003+0.0231i 0.0485i 1 0;+ x y: {+ O/ l# P* _; D! t
17 18 0.0018+0.0144i 0.0303i 1 0;5 f* v/ Q+ ~7 R& i( |% }$ ?3 Q
17 22 0.0135+0.1053i 0.2212i 1 0;
: |- ?2 C+ e+ j 18 21 0.0033+0.0259i 0.0545i 1 0;4 d2 W/ x" y3 ~( \
18 21 0.0033+0.0259i 0.0545i 1 0;
3 S7 K3 E3 ~/ H; Z* q' `! e 19 20 0.0051+0.0396i 0.0833i 1 0;+ X8 i3 n- z' n# V. q8 F
19 20 0.0051+0.0396i 0.0833i 1 0;
& C" }" s: G8 |4 l3 l1 h& G 20 23 0.0028+0.0216i 0.0455i 1 0;
( L# R* X' m! A+ n0 m5 }' ~ 20 23 0.0028+0.0216i 0.0455i 1 0;
- L+ c1 q' C( j: c' V3 ^& C! u 21 22 0.0087+0.0678i 0.1424i 1 0];%input('请输入由支路参数形成的矩阵: B1=');
1 L* k8 Z% Y4 C0 f. I, p- k" UB2=[0.1788 1.0157 1 0 0 2;; M p4 k6 }1 w8 |* D
0 0.9127 1 0 0 2;
3 z& L6 D' ~; K# N8 N 0 1.6935 1 0 0 2;
1 C! w6 L" g% R, M2 k 0 0.6958 1 0 0 2;
) p# U: u4 S: _+ I' V7 y0 X4 p 0 0.6669 1 0 0 2;! ?; z% ]. N. c4 l5 U
0 1.2796 1 0 0 2;
: [# {2 [3 | F 0 1.1748 1.05 1.05 0 1;% y) q& O% `) W9 G/ |1 ^
0 1.6085 1 0 0 2;
6 J# B6 o+ `7 ~( x6 L+ c o" X 0 1.6465 1 0 0 2;
% k, I0 S0 n/ s* e" ` 0 1.8345 1 0 0 2;
1 r$ }7 ~+ K8 g/ q: U! U; v0 J 0 0 1 0 0 2;
! i" g4 ^/ R r9 g. ~3 x' {- F 0 0 1 0 0 2;
0 v1 c+ O1 J# v! U8 y1 z4 z 6.9529 2.4923 1 0 0 2;: G) c6 F* N; t3 o; e1 C
0 1.8236 1 0 0 2;
, @8 m/ T" U. s n9 c8 J' X$ F0 c 1.8235 2.9803 1 0 0 2;
: k& O9 a+ N+ O$ \, W& S! @ 1.8235 0.9398 1 0 0 2;
6 D+ v# c' i6 s8 j. ] 0 0 1 0 0 2;- D( e- f7 g' j1 I1 d
4.7059 3.1322 1 0 0 2;, P" {7 Y$ E- V& j
0 1.7025 1 0 0 2;
% u9 A1 Z: i* L* W5 D$ R 0 1.2037 1 0 0 2;
$ y) ]! g" T7 d7 k# W 4.7059 0 1 0 0 2;5 D, P, X1 E# l. w1 g* o
0 0 1 0 0 2; B6 n/ C9 u# r3 |) I* B9 p4 k% C$ K
7.7647 0 1 0 0 2;) r, a8 l/ }" V J( @ |, x( Q/ R
0 0 1 0 0 2];%input('请输入各节点参数形成的矩阵: B2=');3 b: Q+ N1 u: h) h6 x& M) A
Y=zeros(n);e=zeros(1,n);f=zeros(1,n);V=zeros(1,n);sida=zeros(1,n);S1=zeros(nl);
6 y0 U2 f1 s3 k; J2 N%-------修改部分------------% V5 p( d8 C8 f! |
ym=0;2 K9 k" k( L0 A3 j3 D2 Y
SB=100;UB=220;
. T$ L/ y5 ?( j! v1 t- C* ^$ ^%ym=input('您输入的参数是标么值?(若不是则输入一个不为零的数值)');
3 U3 C" s( V' M- {+ [if ym~=0
+ [7 B6 [& x7 q8 O. O %SB=input('请输入功率基准值:SB=');! e5 U% W& \9 }
%UB=input('请输入电压基准值:UB=');6 g9 {7 H, |* l$ G
YB=SB./UB./UB;
, Q( u" V2 u- @4 T) p4 w' K m: y BB1=B1;
* T, n$ Y" x3 u BB2=B2;
& H* s4 u: [# } for i=1:nl
( U* q' s! N! I B1(i,3)=B1(i,3)*YB;
. ?) k! r5 L! [) D" l+ {/ } B1(i,4)=B1(i,4)./YB;
4 f3 M$ P* B6 k# F end
$ a+ d% P$ n7 k disp('B1矩阵B1=');. l1 S6 H/ S; f: ?( q
disp(B1)
" u4 c! r+ w7 L& t for i=1:n
# e' T. n) I0 C) _% H2 R& [ B2(i,1)=B2(i,1)./SB;
, H" W( f2 l; C* J7 O B2(i,2)=B2(i,2)./SB;+ I, r; H1 ~7 L$ G2 F
B2(i,3)=B2(i,3)./UB;
; p8 |& }! t* Q) t( w B2(i,4)=B2(i,4)./UB;
3 T. A7 |5 r% K+ } B2(i,5)=B2(i,5)./SB;" O1 T8 ~# `1 u$ ]" i
end
5 F# q k, b9 q disp('B2矩阵B2=');3 Z: b7 y2 ^4 q9 G" }
disp(B2)
_0 T3 \6 C, ?) c% F9 Bend* y7 N$ |6 q6 Q2 {% i% F
% % %---------------------------------------------------
) L8 b" }) n" V6 K% N4 dfor i=1:nl9 R/ T/ P- ]) I7 i* L* h- \. s3 q
%支路数4 |% m+ ~- @0 G; ?( f, e: C/ {* n* `
if B1(i,6)==0
" b# F$ X2 j; L5 I' F5 ~%左节点处于低压侧
G. O4 C- r) \8 d& F2 s p=B1(i,1);q=B1(i,2);
) P& h/ m. _8 m6 H& _1 d( { else) N4 ]( i' w& |; p6 W$ Z
p=B1(i,2);q=B1(i,1);- O9 W6 _' x/ W, x) S
end
0 d1 z7 U. Q/ h. w Y(p,q)=Y(p,q)-1./(B1(i,3)*B1(i,5));
+ }( Y K9 k% p* j%非对角元
; Z* x7 Q; i9 [, y0 M2 I Y(q,p)=Y(p,q);
# I, T) Q3 j* [$ z% }7 ?: y% e Y(q,q)=Y(q,q)+1./(B1(i,3)*B1(i,5)^2)+B1(i,4)./2;& }4 n2 C3 f: T, v) w" O
%对角元K侧- _9 I# I3 ]! j% z' d
Y(p,p)=Y(p,p)+1./B1(i,3)+B1(i,4)./24 ]% o& P+ h" d5 D/ y8 `
%对角元1侧
. A9 O( f! S" ~/ u/ V+ D$ H2 b, K8 ]/ Uend5 N( T# A }. I2 p
%求导纳矩阵
+ @, Q1 V" D. g! vdisp('导纳矩阵 Y=');) I7 ]/ ~5 h7 z5 H5 X: G& m
disp(Y)
/ v8 S' k- s! h" \! M8 J%----------------------------------------------------------3 j+ k1 W6 n3 x0 f2 w& J; z
G=real(Y);B=imag(Y);, v6 g- i7 O; n" |( R
%分解出导纳阵的实部和虚部 1 q$ r# K& x" |* v7 i" [
for i=1:n/ c1 k7 H4 i6 P! {: ~6 P
%给定各节点初始电压的实部和虚部 & e2 o2 }3 ~) e; y* E# n+ v
e(i)=real(B2(i,3));
: J. {# ~$ V& @% G2 L/ I. t f(i)=imag(B2(i,3));9 ~8 `+ V, S6 @1 a5 s+ @
V(i)=B2(i,4);
, O( W, A" H8 r: k: A9 @) n%PV节点电压给定模值
! E& R* Q( o8 wend: E( \2 S5 d, K
for i=1:n, I1 O/ [9 w1 ]5 C
%给定各节点注入功率 7 E3 A2 f. O+ d5 N
S(i)=B2(i,1)-B2(i,2);
! e6 G4 m/ }- o+ ~) ~6 L/ ^! W" d9 p%i节点注入功率SG-SL + S% T5 Y: F, E0 i- M, I
B(i,i)=B(i,i)+B2(i,5);
. |# v O+ j w( {! p%i节点无功补偿量 " {) R- ^8 M/ h: w: L
end% M, M7 K9 d- _1 b$ Y6 m( t! s
%===================================================================
i: o/ |# r1 i/ s+ A- w& _9 qP=real(S);Q=imag(S);3 v; \$ a4 G' F% Z$ }
ICT1=0;IT2=1;N0=2*n;N=N0+1;a=0;: O4 R" r' ?5 X; n+ y' I
while IT2~=0
9 P; M9 n" f/ z IT2=0;a=a+1;7 t$ D1 @" w3 l- q! a! Z* y
for i=1:n* u* {2 l3 Q; d; F- \2 R! _" r6 a. _
if i~=isb
) @8 h, d1 f1 V h) Q( k1 m# p4 x3 `%非平衡节点 % _. o4 B" o7 ^5 e! J, s8 [
C(i)=0;D(i)=0;+ Z2 [- ?; @ @; J" J
for j1=1:n
' S! s. _. \( R+ I: ?( r$ x C(i)=C(i)+G(i,j1)*e(j1)-B(i,j1)*f(j1);%Σ(Gij*ej-Bij*fj) H& ~2 I$ {& d( X( }7 S2 L
D(i)=D(i)+G(i,j1)*f(j1)+B(i,j1)*e(j1);%Σ(Gij*fj+Bij*ej)
1 _7 U( j- y8 f# K) [, N end
1 T( Z5 Z% C6 A P1=C(i)*e(i)+f(i)*D(i);%节点功率P计算eiΣ(Gij*ej-Bij*fj)+fiΣ(Gij*fj+Bij*ej)" L: Q( x4 G: s' i6 g t4 C
Q1=C(i)*f(i)-e(i)*D(i);%节点功率Q计算fiΣ(Gij*ej-Bij*fj)-eiΣ(Gij*fj+Bij*ej)
' s F, s! S7 h+ V%求P',Q' : m5 w. g, y/ @4 e r) x
V2=e(i)^2+f(i)^2;
; c6 V. N. `+ _+ p) _1 o7 m%电压模平方) C" R0 J) K* Z' c) M: `, I
%========= 以下针对非PV节点来求取功率差及Jacobi矩阵元素 =========0 ]# V% I+ H4 @
if B2(i,6)~=3
0 Y9 l+ e+ h% W7 R. ^%非PV节点
) ?0 C% v. d: x! U0 N DP=P(i)-P1;
! T; y# M; I2 `! B& o7 q+ q! K- l%节点有功功率差
7 u4 j9 ]$ Q, B: o7 O% u DQ=Q(i)-Q1; %节点无功功率差 ( ]! K g* F1 a# |3 Z/ e: B0 f8 q
%=============== 以上为除平衡节点外其它节点的功率计算 =================9 c, S5 P1 `) ?8 J$ O% {
%================= 求取Jacobi矩阵 ===================+ l2 j* q) g+ f7 N9 N( T
for j1=1:n9 _' W8 c0 v* }
if j1~=isb&j1~=i" S$ M9 H$ c3 t
%非平衡节点&非对角元 : @& T1 ^+ p4 v( g P+ D# s
X1=-G(i,j1)*e(i)-B(i,j1)*f(i);7 P$ x) }* J$ a( o8 w9 U9 ]0 `
% dP/de=-dQ/df
* z! \0 _1 g ~( j7 W" A! {1 f X2=B(i,j1)*e(i)-G(i,j1)*f(i);
# k! a& A% S) S: h9 k% dP/df=dQ/de
4 U6 v) _' x, }) q$ ? X3=X2; % X2=dp/df X3=dQ/de
9 ?; c" d! K0 g X4=-X1; % X1=dP/de X4=dQ/df( x& _: s/ B8 r- g7 g
p=2*i-1;q=2*j1-1;
( A. l: B" {7 m J(p,q)=X3;J(p,N)=DQ;m=p+1;
, h1 @# d' o$ g) W# @% B+ G) s2 }1 C J(m,q)=X1;J(m,N)=DP;q=q+1;9 h- ? T! s& [1 m H) i
J(p,q)=X4;J(m,q)=X2;0 ^4 e5 F1 X9 ~% a! t
elseif j1==i&j1~=isb
) h2 v! m' i0 B. f%非平衡节点&对角元 * L& H3 \% V, m. F
X1=-C(i)-G(i,i)*e(i)-B(i,i)*f(i);% dP/de
2 o' \) t& y% z3 F X2=-D(i)+B(i,i)*e(i)-G(i,i)*f(i);% dP/df8 y' a8 }3 r7 y# W
X3=D(i)+B(i,i)*e(i)-G(i,i)*f(i); % dQ/de0 A8 Z4 k5 E! G2 R6 Z- _
X4=-C(i)+G(i,i)*e(i)+B(i,i)*f(i);% dQ/df
- s# w1 ]3 p0 b( M( X p=2*i-1;q=2*j1-1;J(p,q)=X3;J(p,N)=DQ;%扩展列△Q% v) Y! e7 ]9 m. C1 ]1 C. U
m=p+1;
: J0 r- v( a- ^4 }, ?9 e J(m,q)=X1;q=q+1;J(p,q)=X4;J(m,N)=DP;%扩展列△P. v7 |" C7 z7 R+ o+ ~1 j
J(m,q)=X2;
0 i, H2 n0 ?* W% J* O' G0 G- M$ N end9 y: D4 G: Z0 v f z+ W; ^
end
$ c- N* B2 `. k9 e else1 f- \" E5 w9 b; F
%=============== 下面是针对PV节点来求取Jacobi矩阵的元素 ===========
5 i' e7 n+ Z5 ^! Z8 R' R3 ~5 ]) f& r DP=P(i)-P1;
! P9 f" a* P5 H# }: S% PV节点有功误差
$ m- L. V3 l- {( q DV=V(i)^2-V2;
2 I* L5 G6 M3 ~6 L' n( m% PV节点电压误差
) \: `& c' H7 R0 x for j1=1:n
6 d' U8 B' Z* f/ b7 p if j1~=isb&j1~=i: ^! g9 m) W" W# Z- E+ a
%非平衡节点&非对角元 1 T# f9 y' r$ h) B+ G: a& Z3 y
X1=-G(i,j1)*e(i)-B(i,j1)*f(i); % dP/de
# h [( @! j* C/ T: |% H9 I$ } X2=B(i,j1)*e(i)-G(i,j1)*f(i); % dP/df. F( a3 t: O9 v
X5=0;X6=0;6 Q; [4 p# K# ~6 Y% F% ?- r5 x0 t) O: {
p=2*i-1;q=2*j1-1;J(p,q)=X5;J(p,N)=DV;5 G! N# V/ g& _
m=p+1;
7 v% E2 U3 {3 J% p5 L* D J(m,q)=X1;J(m,N)=DP;q=q+1;J(p,q)=X6;
; p3 J2 l, ~' D$ g( U% j+ G5 Y% Y J(m,q)=X2;( @9 I e5 ]! L8 o
elseif j1==i&j1~=isb
$ X5 u5 T0 U. m2 k$ _: Y%非平衡节点&对角元
- g. s1 Y' D' j1 p2 T( n! j" o% c5 U( a X1=-C(i)-G(i,i)*e(i)-B(i,i)*f(i);% dP/de( h& ]) R1 N0 h3 o6 Z# ^3 o* F( V
X2=-D(i)+B(i,i)*e(i)-G(i,i)*f(i);% dP/df
6 f8 f# \0 H5 m) n$ B X5=-2*e(i);. p; b# V/ P5 G9 y3 m8 `
X6=-2*f(i);/ b( p* @. [ d2 J% H! b# a
p=2*i-1;q=2*j1-1;J(p,q)=X5;J(p,N)=DV;" i$ g$ W+ E$ E$ Q4 q! M
m=p+1;: e d1 ]- B; Q( H; c+ D, ?
J(m,q)=X1;J(m,N)=DP;q=q+1;J(p,q)=X6;, I! R( {1 {% s
J(m,q)=X2;
8 m/ w) n1 C, u( } end2 E4 u% J& f" E- L* z
end
7 i5 l2 ^; w; y4 V+ e" X- [3 I end! e0 A# f7 v N3 w
end0 J! G9 u( x- X0 m! b
end
2 x' F; p' t, n2 @2 x- U9 m7 d%========= 以上为求雅可比矩阵的各个元素 =====================3 X3 [3 V. I$ _# N( {" a# D9 W$ |
for k=3:N0 % N0=2*n (从第三行开始,第一、二行是平衡节点); p7 W9 L( c$ Y# ^
k1=k+1;N1=N; % N=N0+1 即 N=2*n+1扩展列△P、△Q
6 G# ^1 G& Q; I. K4 K' o$ {8 P9 b for k2=k1:N1
1 L* w# z' C2 B$ U% 扩展列△P、△Q 0 ~; u: e# w& p- d
J(k,k2)=J(k,k2)./J(k,k);& o$ {, C5 \3 P5 S1 s% {
% 非对角元规格化
3 D: _2 v! p; P" x6 Z: l* f6 ]( [ end8 p, H# r0 w5 U" R
J(k,k)=1; % 对角元规格化1 n% [- [+ M3 D! n3 `$ s
if k~=3 % 不是第三行% i& p0 Q" {$ [' W& {
9 k; @ n' [4 L T+ N9 ~%============================================================
3 \ `: J( N% R: G, o7 b5 p+ k k4=k-1;
' ]0 u( Y _* b* H0 g8 M4 |& q O for k3=3:k4$ e/ q [! o+ F" m0 ^
% 用k3行从第三行开始到当前行前的k4行消去9 g9 K/ `. m4 J1 Y4 u+ Q# D
for k2=k1:N1
) Y! F( x/ x& B: E% k3行后各行下三角元素' R) t# E( d3 Q& f. k' e
J(k3,k2)=J(k3,k2)-J(k3,k)*J(k,k2);%消去运算
8 N2 B& c2 A: V4 m end9 K( L0 a# _2 Q: t3 ^5 m5 s
J(k3,k)=0;
% d. r6 z: j$ ^4 @) E( K! V end7 U/ q4 V; E5 G: G% B
if k==N0 % P2 h0 }" y+ D C& h5 a
break;; x* e" w( ?, P6 M1 t' ^1 r
end3 s7 ]( I S7 Q% F. Z- P
%==========================================
: f* ]4 g5 E# K$ v for k3=k1:N0
' C/ ^% g. l3 \* x for k2=k1:N1
; ?5 f% V. a" L! A/ I J(k3,k2)=J(k3,k2)-J(k3,k)*J(k,k2);%消去运算- ?4 t/ H$ Y! ?- H, ^
end
6 E- J' a. S! @ J(k3,k)=0;
' e& O' z/ [0 G$ w9 h. z5 a end' ^8 \9 a1 s+ M$ B
else' q" p3 N$ \# c$ `' G
for k3=k1:N0
" L- g- Y. ]7 Y2 }2 e for k2=k1:N1
7 C. o# Q# X# n" c, V J(k3,k2)=J(k3,k2)-J(k3,k)*J(k,k2);%消去运算
^0 Y+ b) V s# y" {2 R0 i end9 U3 c' h3 t: P& E% q% H( b
J(k3,k)=0;/ D. T% a3 D8 z- z
end
/ V' p" G' E* k" d* v! J6 x+ G end& Z% h$ R+ n" I% l) G8 \, C; P
end
6 T4 K2 h) ]7 {7 z1 V m% L0 P %====上面是用线性变换方式将Jacobi矩阵化成单位矩阵=====& m4 V7 U6 X3 {3 }5 L% a H
for k=3:2:N0-1: V: H0 T6 V: |7 Y, k. ^3 v' w
L=(k+1)./2;3 `/ y# |/ z9 n* J5 }0 f# f
e(L)=e(L)-J(k,N) %修改节点电压实部
# n# l& ~) a9 Y) j k1=k+1;/ a ^! j/ W, t3 @+ o, }5 L; n1 `
f(L)=f(L)-J(k1,N) %修改节点电压虚部+ ^# ^* s/ q! a3 p% s: I
end
$ e7 h. @ U5 Y %------修改节点电压-----------7 ^3 I' v# C% t/ E* a- o& n
for k=3:N08 v0 v0 b) H' B- [2 o. i
DET=abs(J(k,N));
6 Z/ w7 `: h/ o- G) N9 l if DET>=pr %电压偏差量是否满足要求9 f8 @8 j% U, M; p( ]% d
IT2=IT2+1; %不满足要求的节点数加1
9 v- o* X% `& m" v" Z end
, t L* m) V' Y end
$ A3 o! O; N& r$ M4 |$ {8 { ICT2(a)=IT2;( D" f' d6 y/ D! i. \
ICT1=ICT1+1;& a* e' H1 M+ x) R1 w9 \9 E
end/ o, x C* T8 O2 R
%用高斯消去法解"w=-J*V"! Z. v- q" i2 B- B" L
disp('迭代次数:');
8 C. @& y$ c- o) j+ ?% c. I9 ^- Y- rdisp(ICT1);4 ~ p( u" M* ?4 Y/ l4 V
disp('没有达到精度要求的个数:');
* L. Q) I$ }; ?* odisp(ICT2);' |. w, i! G, x Z' b* s
for k=1:n
( V& V4 w2 R+ i% ]1 X) g% \7 ^ V(k)=sqrt(e(k)^2+f(k)^2);
+ Q2 F2 x, U9 F5 g6 R8 n6 X sida(k)=atan(f(k)./e(k))*180./pi;
) X' q0 i# W' \ m# C! x4 ~ E(k)=e(k)+f(k)*j;
$ y9 ]$ |8 _8 Y1 d* Mend
* n# E: _9 h9 v2 A* w%=============== 计算各输出量 ===========================; J" H; w) W2 J% t4 X9 v
disp('各节点的实际电压标幺值E为(节点号从小到大排列):');
) K3 D" K+ L, cdisp(E);) z' Z2 D2 ?% }2 L
EE=E*UB;/ {5 B$ x _2 J) v% p
disp(EE);) E: b( M/ f- h
disp('-----------------------------------------------------');
& |6 ^- I% o0 i. N4 L. tdisp('各节点的电压大小V为(节点号从小到大排列):');9 t% ~( ^2 g4 F7 h* `- h" c
disp(V);" {" j' n+ r8 c2 l R, Y
VV=V*UB;
/ e: l3 E4 d5 ^% [( o udisp(VV);
! q+ _8 D) X |% Y$ r5 jdisp('-----------------------------------------------------');
- V/ ~$ g) _4 E, q' ndisp('各节点的电压相角sida为(节点号从小到大排列):');
# A. ~& n; K/ t$ @disp(sida);
' ~, y$ q9 u u* E0 c" a6 Kfor p=1:n
; B4 j/ [" J( r% D" e6 O: T4 ^ C(p)=0;" D- G6 z/ m2 F0 a7 e) @
for q=1:n
' p- A1 u. r; g C(p)=C(p)+conj(Y(p,q))*conj(E(q));* E9 p" H& K; I% d9 q p5 X) H
end
' ]/ G4 s# l4 p* @) ?# r+ @ S(p)=E(p)*C(p);
- O. v) Q5 `2 ^2 cend5 f; p) `1 m4 Z% Y7 h1 z! {4 ]
disp('各节点的功率S为(节点号从小到大排列):');
a }% G# q4 Mdisp(S);
* [; F3 m! Z7 v- G7 I$ Udisp('~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~');& u0 [2 l/ `- z4 g
SS=S*SB;
% p4 S1 y# U6 p7 z0 D) U3 }disp(SS);+ d; v& ~7 t$ a
disp('-----------------------------------------------------');
7 q/ w, D, P8 c( Ddisp('各条支路的首端功率Si为(顺序同您输入B1时一致):');# j: \4 U' F* U. M/ P, U
for i=1:nl
2 U9 o+ J" [) \ |0 x3 l' u( M6 p8 W2 _ p=B1(i,1);q=B1(i,2);
- |) d* z" _4 {6 q3 S. b" d if B1(i,6)==0
1 Q1 P; Y* W9 _+ E) k9 x 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))));5 r8 w1 y' I! K8 e+ v2 N C& c4 a
Siz(i)=Si(p,q);: B# j9 [, ]# z+ H. L
else% N! c7 k4 U7 _4 ]% U% K5 T9 ^
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))));/ W" n, S4 O8 s1 Z$ J1 U' g7 R
Siz(i)=Si(p,q);
! c, A# `! Z- a! \ end, l2 ^" p X& A5 W/ C3 D1 N( T
disp(Si(p,q));6 i3 X* {; e9 E3 ]* x5 O
SSi(p,q)=Si(p,q)*SB;
2 F$ h, G5 D3 F, U ZF=['S(',num2str(p),',',num2str(q),')=',num2str(SSi(p,q))];
+ p: ~* _+ I: u7 h5 f) w0 { disp(ZF);+ Z0 W' a/ T7 `1 W7 d1 d* N
%disp(SSi(p,q));$ B ^+ {& G. p8 T
disp('-----------------------------------------------------');0 E* O" z( C! k: h
end
6 F6 x, N, z/ h" a6 s% t) xdisp('各条支路的末端功率Sj为(顺序同您输入B1时一致):');
" K! r2 Y" N( W! j6 x) D' Q, zfor i=1:nl8 d9 M7 `0 A2 r/ z. K \
p=B1(i,1);q=B1(i,2);/ k) k" H( F# L
if B1(i,6)==06 C8 W6 Z$ N. ?1 x
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 {+ f( ?5 j& U( W* F( X" `
Sjy(i)=Sj(q,p);
3 ]5 p# I4 b4 v# h. c' t _ else" f6 R8 C) u! }" g [' V9 _
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 s. e2 D- m2 A
Sjy(i)=Sj(q,p);
/ E# N W; U( X6 ` end0 G. M& r/ A* g7 R% x
disp(Sj(q,p));# i! R# Y8 z0 h1 t+ o
SSj(q,p)=Sj(q,p)*SB;
" w ^: N4 ~1 N& ^- n ZF=['S(',num2str(q),',',num2str(p),')=',num2str(SSj(q,p))];
+ e7 U: W0 ~$ }" [" L: ~7 ~/ B0 j disp(ZF);
; _# Y2 u, ~* w* E6 W3 M" U %disp(SSj(q,p));
) e$ X* X6 ^; K. [; a2 o' B disp('-----------------------------------------------------');/ {2 l2 M# I9 Z# D
end
0 V. Z" |$ O- Y" n. P m7 w2 Y5 @disp('各条支路的功率损耗DS为(顺序同您输入B1时一致):');" H# B( S2 Z% c8 e
for i=1:nl
; o* f! P/ G2 Y1 q* G: S p=B1(i,1);q=B1(i,2);
8 v& L1 L7 H0 _$ W DS(i)=Si(p,q)+Sj(q,p);; K7 J6 F. l0 O, m4 Z2 C9 c D& v
disp(DS(i));0 |3 J: t6 i: x! g& x; c" O
DDS(i)=DS(i)*SB;
4 @) K0 x8 Q& G3 ]4 H& w ZF=['DS(',num2str(p),',',num2str(q),')=',num2str(DDS(i))];
" q# W4 }! L: |! y, O" t$ @ disp(ZF);* Q2 R0 q1 c3 h# x; a
%disp(DDS(i));& L: t. R1 G% z6 ]# h% [
disp('-----------------------------------------------------');
. r; K7 U6 e4 G4 s2 u, z+ C! S( L/ aend& F k4 D# B5 u" m( l
figure(1);
& ^6 F; l; x" t+ b$ y, t* wsubplot(2,2,1);
' g) I' J- s* W) ~0 ]9 r# c% splot(V);
# U9 p& _/ d' b: x" z+ t+ n4 d- zxlabel('节点号');ylabel('电压标幺值');
7 w5 k( s# p9 @( _grid on;! ^+ e. O$ z& B! A& D- `
subplot(2,2,2);
% S1 D& l: [" lplot(sida);( @! I5 ]8 `# O
xlabel('节点号');ylabel('电压角度');
) E, ?! q3 |, I: }$ l2 zgrid on;! Q0 F1 K! ^% E" l
subplot(2,2,3);
3 y' C9 ^6 ^; h1 P; p! t8 x$ Ubar(S);
8 p, T: S$ i+ F1 {9 }( |ylabel('节点注入有功');
6 x' O! z+ G! O4 Fgrid on;
4 H* {4 ? [1 c1 e2 u( Asubplot(2,2,4); ]; {, B \0 M
bar(Siz);
`0 q' U# i$ w \ylabel('支路首端无功');
$ C5 N% ^5 H3 _0 Q: j* m. A4 ygrid on; |
|