马上加入,结交更多好友,共享更多资料,让你轻松玩转电力研学社区!
您需要 登录 才可以下载或查看,没有账号?立即加入
×
clear %清除工作空间原有的变量
5 {4 K1 |) ^! G8 { L7 t& t, ~; Nclc
2 a2 Q0 W% |3 [" y3 Ja=textread('IEEE14.dat'); %读数据文件* Q9 Y/ G: x0 [! T, ~
n=a(1,1); %节点数
8 t0 D2 s7 E5 K' d5 Nline=a(1,2); %支路数: @7 L1 I: Z* H! G+ A8 \% s
blance=a(1,3); %平衡节点
3 c, V6 }; Y G- ASB=a(1,4); %基准功率
) ?( y. ]/ h8 t6 f3 u% ^kmax=a(1,5); %最大迭代次数) i6 y8 @- |7 ^: C, I6 r
ab=a(2,1); %计算精度
$ M9 C7 H3 ]- Pt=find(a(:,1)==0); 0 E( I9 j$ T, }1 M
xianlu=a([3:t(1)-1],:); %形成线路参数矩阵$ P( h8 q b' b( P
jiedi=a([t(1)+1:t(2)-1],:); %接地支路参数矩阵
0 c8 B* g8 d: }% l3 Hbianya=a([t(2)+1:t(3)-1],:); %形成变压器支路参数矩阵
: Z; b+ Q5 P* Y* z' V* }& Bjiedian=a([t(3)+1:t(4)-1],:);%节点功率数据矩阵3 J7 M# H& S q) F' ?
pv=a([t(4)+1:t(5)-1],:); %pv节点数据矩阵
7 s/ T1 }4 c5 n4 x& b: m%读线路参数
0 }- ~3 E- _5 v" [- y0 T3 Nlinei=xianlu(:,2)';
' B3 W/ J% I- h8 S) X3 g, elinej=xianlu(:,3)';
. M9 u1 B6 J- ~3 kr=xianlu(:,4)';& G. d- d+ Q( p: h5 }3 w8 r
x=xianlu(:,5)';
6 ~( f2 v6 F4 l: j1 L- a; A" F! pb0=xianlu(:,6)';
4 J g9 ~7 d5 _6 V& E1 {
# k& N( y0 b: @g=r./(r.*r+x.*x);
/ S6 t' d0 o- H/ ]b=-x./(r.*r+x.*x);
- Y8 ?, W% G% o: vG=sparse(linei,linej,-g,n,n)+sparse(linej,linei,-g,n,n)+sparse(linei,linei,g,n,n)+sparse(linej,linej,g,n,n);
2 O# Y0 J! ]: {/ E$ mB=sparse(linei,linej,-b,n,n)+sparse(linej,linei,-b,n,n)+sparse(linei,linei,b+b0,n,n)+sparse(linej,linej,b+b0,n,n);
' A! r' T. C; _* h7 _7 a/ W; I% 读变压器参数+ G: S0 [0 q/ l* W/ k4 G5 D) {
linei=bianya(:,2)';, p" f9 e% `0 J) a
linej=bianya(:,3)';5 j; f) G* T6 e' z
r=bianya(:,4)';7 W& a; |6 K2 g* e K* h6 }
x=bianya(:,5)';0 I( K$ ~0 \+ m4 n, i5 Q! ^
k0=bianya(:,6)';& Y" U+ n% f" L& D% y
g=r./(r.*r+x.*x);' d. P6 g4 _# z3 V
b=-x./(r.*r+x.*x);
: k9 D& v0 r8 Z6 f) YG=G+sparse(linei,linej,-g./k0,n,n)+sparse(linej,linei,-g./k0,n,n)+sparse(linei,linei,g./k0./k0,n,n)+sparse(linej,linej,g,n,n);
2 K' E6 t6 A; Z+ O- zB=B+sparse(linei,linej,-b./k0,n,n)+sparse(linej,linei,-b./k0,n,n)+sparse(linei,linei,b./k0./k0,n,n)+sparse(linej,linej,b,n,n);
3 n5 c8 {5 G' c' a% 读接地支路参数8 `0 T0 \7 w( c; G, \2 t
jiedidian=jiedi(:,1)';
* I4 e* G* N4 {: X% tbranchib=jiedi(:,2)';
' V, u. s/ p( F2 rB=B+sparse(jiedidian,jiedidian,branchib,n,n);
$ N1 T7 X' S5 l u( P& \# b% R%节点导纳矩阵
0 d: D z5 d+ `; w' S- H, VY=sparse(G+i*B);" W5 d# J0 p- n0 t% z1 t9 m3 R
% 读节点功率数据
4 i2 W! n! y/ ?+ r. [pg=jiedian(:,2)/SB;( Z7 s9 t) ~5 g* i
pd=jiedian(:,4)/SB;
/ {: e- B: }$ l' p& g, }" D/ hqg=jiedian(:,3)/SB;
* Y& c! u! r' sqd=jiedian(:,5)/SB;4 b l& @; j6 v% ^) m, @$ e+ V m
* Q" N" `, v' F- m( Y. Y4 Q
p=zeros(n,1);
6 O' `& l% }: f$ ?1 ]1 lq=zeros(n,1);" v4 B$ B* j x9 ?
p(jiedian(:,1))=pg-pd;3 I- ~) y1 _5 L0 O8 g9 j, e* r
q(jiedian(:,1))=qg-qd;3 g p) L# w4 t' }1 k8 r4 K" V* e" ]
% 读pv节点数据
9 ~2 p/ g' z$ w+ V- Y6 I$ tpvjiedian=pv(:,1);' O8 |; Y9 N0 L, b& y
vpv=pv(:,2);$ |2 i, }! d Q+ W F
v=ones(n,1); %电压初始值
, @, ~* M9 G5 I3 z1 e- dv(pvjiedian)=vpv; %pv节点(及平衡节点)电压已知
: Q% ~1 | @2 S' ]' ldelt=zeros(n,1); %相角初始值" V* X0 e5 z1 P/ Z. B! B
deltp=zeros(n,1); % 失配有功功率初值" b/ }$ O5 \$ N3 G+ r; L! w8 `5 `) @
deltq=zeros(n,1); % 失配无功功率初值- }' {7 ?8 b' E$ u
3 @. \1 x( _8 L
%迭代求潮流计算% f% T# N! f) \7 i* x* E( A8 t
for diedai=1:10 %置最大循环次数
% ?3 a. i9 i: y2 c clear i
: K* v1 G6 g9 p: j& q % 失配功率方程
; H, n' s+ n J% c5 r7 W- ? V=v.*cos(delt)+i*v.*sin(delt);, F0 h+ b$ v. i' q8 q
a2=conj(Y*V).*V;
) M9 Z1 W' j, g% n1 _9 x deltp=p-real(a2);
# z$ U" U' G$ O4 S$ ]. J1 T deltq=q-imag(a2);
) ]" S# t5 _: C& A: u7 F0 h
- q( A3 V6 ]; X+ `0 t' ` % 求雅可比矩阵8 Z t0 E9 e' `. k5 V! c$ z
for i=1:n
* A% h+ v! s6 C R* p: `: H+ Z for j=1:n$ x: K8 Q4 }7 a9 u/ L! r
if i~=j- _( q& o/ c# ?; x
H(i,j)=-v(i)*v(j)*(G(i,j)*sin(delt(i)-delt(j))-B(i,j)*cos(delt(i)-delt(j)));
3 ~9 r) t4 K3 M else
6 z; j3 X9 \1 ^% ]; ^3 I H(i,i)=v(i)*v(i)*B(i,i)+q(i)-deltq(i);
% |( K! s$ I- @. n. R) |- [7 I& V' V end
6 X0 l$ B# W* f o end
t5 L4 {! I: {9 W end
( g( Y+ Z- G3 x. e6 X) b2 R * o# u/ Q$ Y% @
for i=1:n
# h' ~% B4 `+ R# ^) B/ o for j=1:n
" a) A' I' D4 V2 J* Y if i~=j% r: V& s+ x. z. C
N(i,j)=-v(i)*(G(i,j)*cos(delt(i)-delt(j))+B(i,j)*sin(delt(i)-delt(j)));
0 M) v& u M/ W& L X else
) L3 m# R# r2 b& ~' ^; p) a% k N(i,i)=-v(i)*G(i,i)-(p(i)-deltp(i))/v(i);
; E1 o2 ]8 M9 a/ R1 @( P end: ~; M. E- O; H, Q# F l
end
; Y5 O' Q' U) T7 }5 s: v6 N end * G1 f3 l) T5 ]2 ?( [
2 p$ F: I- w2 q. ^2 ^. M( c/ B
for i=1:n2 L1 ?; U+ V( q" I
for j=1:n
& W9 N, X0 L8 F! c% H F+ A if i~=j
) R8 a) [4 x0 O J(i,j)=v(i)*v(j)*(G(i,j)*cos(delt(i)-delt(j))+B(i,j)*sin(delt(i)-delt(j)));- o, \9 C( P1 k* i
else
5 f5 M6 f$ t/ }! \) s9 D1 R8 W7 Z J(i,i)=v(i)*v(i)*G(i,i)-(p(i)-deltp(i));* Z: Y9 r. U$ M8 _: U: h5 g% m' Q
end. @# H7 R% t# y) [2 e% ^
end0 d y3 i# I* D# z; C7 C
end 1 E" ^ d$ |+ C8 D- P- x) y* N
4 ]! M$ D8 L R7 g) K9 I
for i=1:n
+ {' G0 L% S( p for j=1:n
G" [- o& R) k2 F) e2 ~ s if i~=j2 Q+ ^: W! K6 I; K
L(i,j)=-v(i)*(G(i,j)*sin(delt(i)-delt(j))-B(i,j)*cos(delt(i)-delt(j)));
( Y4 g1 L+ K& l% k% b else
0 z) r- P E$ S3 d& @) y' P# h L(i,i)=v(i)*B(i,i)-(q(i)-deltq(i))/v(i);
9 h7 s0 O" @& H" ?+ U end
+ P/ d8 c- Q% u) T. Y: |# F. K! W' N end
y" H3 Y: ?! w6 Y end
# P- T( I: K2 \5 S. U %修正雅克比矩阵+ T' q5 D7 n7 _# |
Ya=[H N;J L];
9 V9 v5 K* f; g) H Ya([blance blance+n],:)=0;
3 D% |4 Z0 C" e# J% w! d, t" U Ya(:,[blance blance+n])=0;
- X: d2 g/ T4 r$ R2 E Ya(blance,blance)=1; %平衡节点
! z3 L' O! [ D( r4 J9 A1 g Ya(blance+n,blance+n)=1;: F8 T+ t' q+ }. l. ? {
Ya(pvjiedian+n,:)=0;% {% X7 A) r3 p/ c: ^* h/ i
Ya(:,pvjiedian+n)=0;( [2 Q4 C5 [" p+ k
for h=1:length(pvjiedian)
) t4 @ x; v$ U/ g Ya(pvjiedian(h)+n,pvjiedian(h)+n)=1;
" {0 T& A7 m' j. I _ end
4 B! X, V/ h7 { / Z8 J P% l6 i+ N
% & y6 x3 ]1 C/ o) L
deltp(blance)=0; %平衡节点deltp=0' m( H$ f, W& d& t+ f" k. @3 h
deltq(blance)=0; %平衡节点deltq=0
7 `8 s* H5 q8 i* h: k; q1 g4 n 3 s6 y" ~! k+ _3 V5 P
deltq(pvjiedian)=0; %pv节点deltq=0
' d5 ], `1 G( Z5 _7 Z$ |3 q8 |2 I deltpq=[deltp;deltq];' x* B* F# f% u
if max(abs(deltpq))<ab % 如果结果满足精度
2 m# v }6 R, J4 [4 x break % 退出循环 : v% E0 t0 j& {: Z7 |* M
end : @: {4 Y) c( a% D7 t, {
R=Ya\(-deltpq);
/ f4 r0 _8 \6 Z2 Q3 j( ] delt=delt+R([1:n]);
) N) V- _6 C& P deltt=delt*180/pi
9 ~) D/ y8 j* n, G4 Z6 M* Y1 Q v=v+R([n+1:2*n])/ y n& U1 ~) G& `/ O$ g
end
1 H. Z) x; d9 e* z/ c7 a1 d4 } S5 |. l* _( U I
' S F$ j$ r' i) z X$ I9 Z6 J9 r/ i; r数据格式1 B' R+ P I! I* s: C
| 14 | 20 | 1 | 100 | 20 | 0.1 | | | 1.00E-06 | 1 | | | | | | | 1 | 1 | 2 | 0.01938 | 0.05917 | 0.0264 | | | 2 | 1 | 5 | 0.05403 | 0.22304 | 0.0246 | | | 3 | 2 | 3 | 0.04699 | 0.19797 | 0.0219 | | | 4 | 2 | 4 | 0.05811 | 0.17632 | 0.017 | | | 5 | 2 | 5 | 0.05695 | 0.17388 | 0.0173 | | | 6 | 3 | 4 | 0.06701 | 0.17103 | 0.0064 | | | 7 | 4 | 5 | 0.01335 | 0.04211 | 0 | | | 11 | 6 | 11 | 0.09498 | 0.1989 | 0 | | | 12 | 6 | 12 | 0.12291 | 0.25581 | 0 | | | 13 | 6 | 13 | 0.06615 | 0.13027 | 0 | | | 14 | 7 | 8 | 0 | 0.17615 | 0 | | | 15 | 7 | 9 | 0 | 0.11001 | 0 | | | 16 | 9 | 10 | 0.03181 | 0.0845 | 0 | | | 19 | 12 | 13 | 0.22092 | 0.19988 | 0 | | | 20 | 13 | 14 | 0.17093 | 0.34802 | 0 | | | 4 | 14 | 9 | 0.12711 | 0.27038 | 0 | | | 5 | 10 | 11 | 0.08205 | 0.19207 | 0 | | | 0 | | | | | | | | 9 | 0.19 | | | | | | | 0 | | | | | | | | 1 | 4 | 7 | 0 | 0.20912 | 0.978 | 0.9 | | 2 | 4 | 9 | 0 | 0.55618 | 0.9697 @4 I5 v8 A3 G- c3 D% a1 v
0.9 | 1.1 | | 3 | 5 | 6 | 0 | 0.25202 | 0.932 | 0.9 | | 0 | | | | | | | | 1 | 60 | 0 | 0 | 0 | | | | 2 | 65 | 42.4 | 21.7 | 12.7 | | | | 3 | 0 | 23.39 | 94.2 | 19 | | | | 4 | 0 | 0 | 47.8 | -3.9 | | | | 5 | 0 | 0 | 7.6 | 1.6 | | | | 6 | 85 | 12.24 | 11.2 | 7.5 | | | | 7 | 0 | 0 | 0 | 0 | | | | 8 | 0 | 17.36 | 0 | 0 | | | | 9 | 0 | 0 | 29.5 | 16.6 | | | | 10 | 0 | 0 | 9 | 5.8 | | | | 11 | 0 | 0 | 3.5 | 1.8 | | | | 12 | 0 | 0 | 6.1 | 1.6 | | | | 13 | 0 | 0 | 13.5 | 5.8 | | | | 14 | 0 | 0 | 14.9 | 5 | | | | 0 | | | | | | | | 1 | 1.06 | -40 | 50 | | | | | 2 | 1.045 | -40 | 50 | | | | | 3 | 1.01 | 0 | 40 | | | | | 6 | 1.07 | -30 | 40 | | | | | 8 | 1.09 | -30 | 45 | | | | | 0 | | | | | | | | 1 | 105 | 2.45 | 0.005 | 50 | 200 | | | 2 | 44.4 | 3.51 | 0.005 | 20 | 100 | | | 6 | 40.6 | 3.89 | 0.005 | 20 | 100 | | | 0 | | | | | | | | 0 | | | | | | | | 0 | | | | | | | | | | | | | |
|