|
|
马上加入,结交更多好友,共享更多资料,让你轻松玩转电力研学社区!
您需要 登录 才可以下载或查看,没有账号?立即加入
×
各位~~我根据下图的一个系统用Matlab编写了一段程序,是电流型前推回代法的三相配电系统潮流计算,但是现在结果不是我想要的,不收敛。请各位帮我看一下,提提建议。图片见附件了。还有个问题:已知的线路都是三相对称线路的阻抗值,那么线路的阻抗矩阵中的互阻抗应该如何计算,是什么样的形式,应用怎样的公式?这里我不是很清楚。- %程序名:qiantuihuitui_I_3.m 3 P( x1 E0 C( ~1 G* f6 J4 S
- %功能:支路电流前推回推法求解潮流 6 I& ]8 |" X& R$ V$ a, R2 E3 v
- clc 1 t2 \2 M$ k' [# B8 o. I
- clear all;
" u3 y1 q7 `2 V0 Q - %--------------输入网络参数-------------- / u; \& P8 Q( _
- %1-支路编号,2-首节点,3-尾节点,4-自阻抗,5-尾节点复功率,6-支路性质(1-馈线段支路,2-变压器支路),7-尾节点是否带负荷
) f& D7 f. f! ~" q5 v$ U/ P - DB=[1 1 2 0.000167+j*0.000208 0.42+j*0.31 1 1
0 M, ~8 E. q( c$ W - 2 2 3 0.000151+j*0.000188 6.15 1 0
+ r2 k, }6 K* ^ _( [; p- A5 I - 3 2 4 0.000066+j*0.000082 0.38+j*0.29 1 1
& ~1 k5 u. y3 f& E+ s9 d1 K5 Q - 4 2 5 0.000249+j*0.000310 0 1 0
: L# E- o% e7 N" a3 E; [ - 5 2 6 0.000172+j*0.000215 0 1 0' p0 O8 G. ?3 `( f# A6 N2 Z
- 6 4 7 0.000156+j*0.000195 6.06 1 01 w5 B0 ?$ u" S6 n8 g( S: o% `
- 7 4 8 0.000162+j*0.000202 6.04 1 0 ; G+ l7 K/ A+ v; ]& L; P
- 8 4 9 0.000345+j*0.000430 0 1 0
) o; \- j+ |& {" I& V/ l0 w - 9 4 10 0.000287+j*0.000358 0 1 0
7 v: g* u3 a3 Z* ~3 C* J - 10 5 11 0.020563+j*0.321594 0 2 0 5 {. T+ G$ N- g6 t- K8 z3 V
- 11 6 12 0.020563+j*0.321594 0 2 0
, g. p, v2 n6 ]/ m, s8 \8 `& K P - 12 9 13 0.020563+j*0.321594 0 2 0
, i* i4 u5 r% l! B# G" C/ v) Q - 13 10 14 0.020563+j*0.321594 0 2 0
+ w( C/ k! A& f- a& Y - 14 11 15 0.000237+j*0.000408 5.72+j*0.12 1 1 ! O% |1 J3 V3 W& Z: p
- 15 12 16 0.000237+j*0.000408 5.76+j*0.09 1 18 d' }4 L8 S" ?5 ^: O$ ?
- 16 13 17 0.000292+j*0.000502 5.86+j*0.11 1 15 ], k) x; a, X. N2 N
- 17 14 18 0.000274+j*0.000470 5.81+j*0.14 1 1]; 2 o+ y$ L" f4 q) N* `: L
- [n,m]=size(DB);
/ k$ l2 P6 x$ |1 u- W) o+ | - B=[1 sin(2*pi/3)+j*cos(2*pi/3) sin(4*pi/3)+j*cos(4*pi/3)];6 y1 ~/ q* y0 \$ y
- C=[1.02 1 1.02 1 1 1 1.02 1.02 1 1 1 1 1 1 1 1 1 1];
0 Z% j8 M! s( a! z - U(:,1)=B(1,1)*C';
, u" M2 h3 I& i& g3 ]8 C; r+ Y - U(:,2)=B(1,2)*C';! M8 J9 |# @' a: Y( N
- U(:,3)=B(1,3)*C';
; T) s3 }2 y: K' ]# m( K M( S2 | - %-------------------求解潮流-----------------
& J5 s; s9 [' q6 p- a - for k=1:15
( I! d5 z& D5 o- j- z - % I(:,k)=((DB(:,5).')*(diag(1./U(:,k))))'; $ u. u W8 I" u' d- k( s
- for i=n:-1:1
( V+ D0 r/ g; M w- s) k - %如果尾节点带恒功率负荷,需计算节点上负荷注入的电流
: m! [6 ^2 X0 D1 F3 D2 d/ p' o - if DB(i,7)==13 M) M( l, \5 o) R, T
- c=DB(i,5)/3;
. M0 Y/ E) E( f1 U - d=c/U(i+1,(3*k-2));+ p/ K1 m/ n% r) K# L
- IL(i,3*k-2)=conj(d);# b2 X7 u8 w) L: H s0 D5 B- v3 ?
- d=c/U(i+1,3*k-1);6 a. k$ a1 j x
- IL(i,3*k-1)=conj(d);$ H" [/ w& k9 |+ ?! Z4 L. C
- d=c/U(i+1,3*k);
2 R( N9 x5 S9 \: R - IL(i,3*k)=conj(d);
0 l5 @! P6 h/ G+ z9 G6 ]5 N4 [( g - else
( E& x) I2 \( C1 n, S S. @, i - IL(i,3*k-2)=0;
! d; E: n2 S, z( X2 \ - IL(i,3*k-1)=0;7 W# M+ y6 g+ a# _9 [6 S" f3 m
- IL(i,3*k)=0;
p$ w4 D! ~3 \! w - end3 e3 J4 t+ A% j Z
- %找出所有与尾节点相连的支路,计算进支电流(末端电流),存IKj
, G6 g& S# Z! ?9 n2 @+ ?% b - A=(find(DB(:,2)==(i+1)));
* b% I4 t6 q* T- }4 K0 ?$ [ - if isempty(A)~=12 @& s, }/ q# v K0 k
- IKj(i,3*k-2)=IL(i,3*k-2)-sum(IKj(A,3*k-2));, V( q8 C; u/ z7 X
- IKj(i,3*k-1)=IL(i,3*k-1)-sum(IKj(A,3*k-1));
- {# s6 ~5 k7 [3 ~1 u1 m* I - IKj(i,3*k)=IL(i,3*k)-sum(IKj(A,3*k));7 l% ]# H u& y- \! e
- else
' W- H; e% h5 S; u/ w' i7 V' A! s1 ` - IKj(i,3*k-2)=IL(i,3*k-2);
- L4 n$ S5 R# e6 P9 F# G( f - IKj(i,3*k-1)=IL(i,3*k-1);
1 r3 A. P. z! ]3 A - IKj(i,3*k)=IL(i,3*k);6 R4 E2 W5 F7 J) u* k
- end
' K: L+ o$ l. Z; Z6 Z0 J - %计算出支电流(始端电流),存IKi
2 S% T+ a8 o9 W1 P) N. D$ z - a=DB(i,2);
- `3 {- e! @ l H7 b - b=DB(i,3);
8 p) n2 e5 H+ {/ E0 n/ ^" R( |. F9 j - Yi=1/DB(i,4)*eye(3);
' _6 p3 G8 M* X8 T# v - if DB(A,6)==1$ U4 T+ ~" b: n, q
- F=0.5*Yi*[(U(a,3*k-2)+U(b,3*k-2)) (U(a,3*k-1)+U(b,3*k-1)) (U(a,3*k)+U(b,3*k))]'+[IKj(i,3*k-2) IKj(i,3*k-1) IKj(i,3*k)]';
N. H" J) j0 g, } V& b1 a - IKi(i,3*k-2)=F(1,1);
1 _& b J' e3 A0 l5 O - IKi(i,3*k-1)=F(2,1);
; C+ {3 g1 Q# ~6 [ - IKi(i,3*k)=F(3,1);/ G; ]+ T, l9 a+ Z) s9 e9 a1 t
- else
+ t0 S# X6 t' F* s! @/ c8 X - YT=1/real(DB(i,4))+j*(1/imag(DB(i,4)));2 s) |. I9 ?" h( [0 w
- YI=YT*eye(3);$ Q3 V) i% Z: H o! G: c/ x1 Z! ^
- YII=YT*eye(3);2 c4 H$ C2 m' J
- YIII=-YT*eye(3);; J2 ? l: V' P$ l, N- g9 j S
- D=inv(YIII)*(-[IKj(i,3*k-2) IKj(i,3*k-1) IKj(i,3*k)]'-YII*([U(b,3*k-2) U(b,3*k-1) U(b,3*k)]'));' l: h" m# j0 D: f9 R: t2 ~; o
- U(a,3*k+1)=D(1,1);# n1 O3 W# q% I6 T9 ^% \
- U(a,3*k+2)=D(2,1);
% N' O( ?2 F# f a& f4 Q7 W - U(a,3*k+3)=D(3,1);' F1 b& r7 \, \! v6 A
- E=YI*[U(a,3*k+1) U(a,3*k+2) U(a,3*k+3)]'-YIII*[U(b,3*k-2) U(b,3*k-1) U(b,3*k)]';
" l, X- P& i6 a: B* r- K" R8 O* V - IKi(i,3*k-2)=E(1,1);" c9 a' B# B. F4 E& v
- IKi(i,3*k-1)=E(2,1);9 m6 D" b7 R, m4 b) a
- IKi(i,3*k)=E(3,1);$ B O, K$ M6 v2 q9 z! X& E
- end
: ^* n! L% d& U" X1 z+ j% N0 N - end 9 b1 @: F0 ]7 o2 }& q; F ?
- %前推电压
+ ]% d+ z& L! b$ U* i# i4 D - for j=2:n( k! w$ x& X& N' O
- U(1,3*k-2)=1.02;+ f5 ^5 W; f7 K3 u/ J! r, C& C2 ?
- U(1,3*k-1)=1.02*(sin(2*pi/3)+j*cos(2*pi/3));
2 q- r$ @6 r% S - U(1,3*k)=1.02*(sin(4*pi/3)+j*cos(4*pi/3));
% Z$ N" z# u1 y - a=DB(j-1,2);
2 F' `( {0 u7 m6 b% w" ^, u - if DB(j-1,6)==11 \ K2 a9 T: P2 l$ v
- Yi=1/DB(j-1,4)*eye(3);
7 y8 C P# U2 `( K( @: @ - G1=[IKi(a,3*k-2) IKi(a,3*k-1) IKi(a,3*k)]';
0 `$ J: h# a) c7 M - G=[U(a,3*k-2) U(a,3*k-1) U(a,3*k)]'-DB(j-1,4)*eye(3)*(G1-0.5*Yi*[U(a,3*k-2) U(a,3*k-1) U(a,3*k)]');; d4 {. }; m: H* ?5 s* V( M
- U(j,3*k+1)=G(1,1); ) S9 m. C' ]% \+ B& q9 y! K
- U(j,3*k+2)=G(2,1);
9 R. u. _, k: P* \6 P. s - U(j,3*k+3)=G(3,1); % }& k4 `2 y8 w" x" W) [
- else& c5 J# D2 u+ y& s0 j" e: s% h# N
- YT=1/real(DB(i,4))+j*(1/imag(DB(i,4)));
; |" L* m2 @) m+ W: M; p: r6 r - YI=YT*eye(3); _: k0 X' a- }8 K' u
- YII=YT*eye(3);
, F& h$ }3 E( e n1 P2 u5 n0 t% D1 o- w - YIII=-YT*eye(3);
/ O! O2 ]! y4 T - H=inv(YIII)*([IKi(a,3*k-2) IKi(a,3*k-1) IKi(a,3*k)]'-YI*[U(a,3*k+1) U(a,3*k+2) U(a,3*k+3)]');
. U$ Z1 }+ D4 F) W - U(j,3*k+1)=H(1,1);% Q' F+ T4 S: q
- U(j,3*k+2)=H(2,1);5 x% d9 f" u: ^ X
- U(j,3*k+3)=H(3,1);
. g! M3 w$ o% t - end ; b% V- X& K* L' R5 @9 n
- end; Y5 Y5 c, N9 H ]$ v, G
- dU1=U(:,3*k+3)-U(:,3*k);
+ m5 a. t$ w$ G. E - dU2=U(:,3*k+2)-U(:,3*k-1);
) r1 j( P% |' H4 G# T7 t - dU3=U(:,3*k+1)-U(:,3*k-2);
V# |5 s: V; [7 W' z4 p - if (max(abs(dU1))<0.001&&max(abs(dU2))<0.001&&max(abs(dU2))<0.001) / G6 p9 U# x, F# O) a% r
- break;
7 {# r' Y) C. i; h, _, d1 X1 { - else ) [; i/ E9 @% I) q
- k=k+1; v- U; B) h7 d8 \$ {: B
- end * K; V, |( d0 m7 c3 H1 `/ K/ {
- end; J8 G( C+ m" G
- U, I1 w$ W$ k& G9 A! X# g) e/ B I A
- IL
8 k& p: C. |7 D6 _ - IKj
6 i2 ^( b, ~- t/ C% E$ j - IKi& ~! z% Q3 t8 o% P! Q8 V/ y1 \
- k
8 _) Y) H7 b* L" k
复制代码 |
|