|
马上加入,结交更多好友,共享更多资料,让你轻松玩转电力研学社区!
您需要 登录 才可以下载或查看,没有账号?立即加入
×
各位~~我根据下图的一个系统用MATLAB编写了一段程序,是电流型前推回代法的三相配电系统潮流计算,但是现在结果不是我想要的,不收敛。请各位帮我看一下,提提建议。图片见附件了。还有个问题:已知的线路都是三相对称线路的阻抗值,那么线路的阻抗矩阵中的互阻抗应该如何计算,是什么样的形式,应用怎样的公式?这里我不是很清楚。- %程序名:qiantuihuitui_I_3.m . X" O, ^0 M: I1 \9 Q9 K
- %功能:支路电流前推回推法求解潮流
6 Q+ P! A: i# V. a, h - clc
u) e+ v/ I1 G - clear all;
0 J9 f' K' s+ Y0 m% z4 F% {3 l - %--------------输入网络参数-------------- L, ^- y: I5 m k" l) ], Q
- %1-支路编号,2-首节点,3-尾节点,4-自阻抗,5-尾节点复功率,6-支路性质(1-馈线段支路,2-变压器支路),7-尾节点是否带负荷
' Z, C! ~$ v {- P- T - DB=[1 1 2 0.000167+j*0.000208 0.42+j*0.31 1 1% H" I5 ]3 y! x1 ~( g
- 2 2 3 0.000151+j*0.000188 6.15 1 0. y3 q' X3 m) f8 p/ z
- 3 2 4 0.000066+j*0.000082 0.38+j*0.29 1 1* ?9 |8 I4 h6 G. ~. G
- 4 2 5 0.000249+j*0.000310 0 1 0
- P7 l4 o6 S6 \# |) { - 5 2 6 0.000172+j*0.000215 0 1 07 P7 [' W% m5 B+ }3 @+ ~
- 6 4 7 0.000156+j*0.000195 6.06 1 08 ?+ }7 U, _3 ]3 z0 S$ J
- 7 4 8 0.000162+j*0.000202 6.04 1 0 & Q* L' C) S1 W* U7 E5 M* {9 i o
- 8 4 9 0.000345+j*0.000430 0 1 0
, R% D( h# w. r, i9 D4 k. x - 9 4 10 0.000287+j*0.000358 0 1 0
n* Z5 I& Y! v% R - 10 5 11 0.020563+j*0.321594 0 2 0
0 q3 |$ ^# T2 w1 ?. x1 `- U4 n' Q - 11 6 12 0.020563+j*0.321594 0 2 0- K) U* @2 L/ W4 d
- 12 9 13 0.020563+j*0.321594 0 2 0 5 c( |! H/ q0 E6 m
- 13 10 14 0.020563+j*0.321594 0 2 0
: y& Q1 L- J; v3 |# T - 14 11 15 0.000237+j*0.000408 5.72+j*0.12 1 1 8 ]6 l6 i' t8 w! b% S- T7 y; ?7 x
- 15 12 16 0.000237+j*0.000408 5.76+j*0.09 1 1' Q" w/ }4 a1 ?) w- l
- 16 13 17 0.000292+j*0.000502 5.86+j*0.11 1 14 p, }, I: \/ c' G: o: m ^& P
- 17 14 18 0.000274+j*0.000470 5.81+j*0.14 1 1];
) w/ P( D; w3 V' H; j! x - [n,m]=size(DB);
$ A& K. K5 A# K - B=[1 sin(2*pi/3)+j*cos(2*pi/3) sin(4*pi/3)+j*cos(4*pi/3)];
( Y" k$ Q9 K c. a; G - C=[1.02 1 1.02 1 1 1 1.02 1.02 1 1 1 1 1 1 1 1 1 1];
2 w3 r5 x) X( [# i9 Z - U(:,1)=B(1,1)*C';- V, T2 g* l) e% [/ x# n
- U(:,2)=B(1,2)*C';# A' }4 R+ c5 ^
- U(:,3)=B(1,3)*C';. N3 G8 ?& _" y
- %-------------------求解潮流-----------------
1 E9 I" u/ ?- d# I* ?3 k; q+ U, y - for k=1:15 " h1 Y9 }" R3 n8 a" U* H
- % I(:,k)=((DB(:,5).')*(diag(1./U(:,k))))'; ! M \! o3 l' E Y* p) K0 L1 K
- for i=n:-1:1
, g1 C4 V) Z& P% X7 G! Z# @( u. M - %如果尾节点带恒功率负荷,需计算节点上负荷注入的电流
' t& E9 L) v) @ - if DB(i,7)==1# B, x9 s i- z' ]
- c=DB(i,5)/3;% k& i0 A! V- d
- d=c/U(i+1,(3*k-2));1 ? @5 @9 @' y! f6 H
- IL(i,3*k-2)=conj(d);+ c" o/ b( C4 v4 P3 v
- d=c/U(i+1,3*k-1);; @' |6 ?5 O$ l4 ^6 \: v3 o
- IL(i,3*k-1)=conj(d);- E4 {6 g2 Q, J! s8 E6 A* _
- d=c/U(i+1,3*k);
9 A. Q4 C8 a6 ?! W- b - IL(i,3*k)=conj(d);! S; k, S7 N1 C/ V. d7 X8 P" s# t
- else
0 d6 b6 q4 f1 ^& w* \ - IL(i,3*k-2)=0;
" [ V ]5 @& d; q* Y( [) O - IL(i,3*k-1)=0;
' O6 f% _9 x" t; e - IL(i,3*k)=0;
, p: ]0 j6 W" n. Z% _ - end& R! h) ^0 x: n; ?; f' k7 B
- %找出所有与尾节点相连的支路,计算进支电流(末端电流),存IKj
+ u3 ^3 J @7 U" i. B - A=(find(DB(:,2)==(i+1)));% t0 A+ {3 R2 A& l# a
- if isempty(A)~=1
v; Z. A1 A1 U$ P5 q - IKj(i,3*k-2)=IL(i,3*k-2)-sum(IKj(A,3*k-2));
% L* s) h" a' ~) p4 Q - IKj(i,3*k-1)=IL(i,3*k-1)-sum(IKj(A,3*k-1));
/ q- L( @, q' F) D7 K. |1 ^) o - IKj(i,3*k)=IL(i,3*k)-sum(IKj(A,3*k));8 r! U7 y3 A. f m, M* o1 E4 E
- else
2 C0 Q6 j/ I: V2 N) O - IKj(i,3*k-2)=IL(i,3*k-2);7 V3 J% `8 H; l ~' S
- IKj(i,3*k-1)=IL(i,3*k-1);
9 }/ {# ~8 G3 o# Q: s) T - IKj(i,3*k)=IL(i,3*k);
$ _8 u# E6 ?: B* z$ V x; a- g7 b - end [6 u# A9 k/ o; e0 A
- %计算出支电流(始端电流),存IKi- t g( p- m& p+ K& [8 w* z
- a=DB(i,2); b* x# S; X0 A4 @% K
- b=DB(i,3);5 V; q* a" T- S# b ~
- Yi=1/DB(i,4)*eye(3);3 `5 G% u7 ?% T0 k4 |" e
- if DB(A,6)==1! J2 \3 w9 o3 v9 Q2 l/ m
- 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)]';6 c" l; w8 H# u5 p
- IKi(i,3*k-2)=F(1,1);
# J% E: ~. G$ C: W. f4 ~! v - IKi(i,3*k-1)=F(2,1);
# F, M' }3 Q1 u: o - IKi(i,3*k)=F(3,1);
; E8 s1 ]! m, I - else
3 D1 I, f1 u* ?& L - YT=1/real(DB(i,4))+j*(1/imag(DB(i,4)));( L) d- @+ Y5 C* W& w$ a0 F: ^
- YI=YT*eye(3);
' V \- U8 C' L6 L2 r1 W - YII=YT*eye(3);
- x. o. a0 J# @( m - YIII=-YT*eye(3);
K8 I) w# G5 u* S, `7 P2 V) U* B* P - 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)]'));
" Y; n( V- D5 S! `) n. d* x1 B - U(a,3*k+1)=D(1,1);% ^ g( M6 z( n) K2 g
- U(a,3*k+2)=D(2,1);
) h, ^: t' x- |. _ - U(a,3*k+3)=D(3,1);" P! V6 }/ t7 n& J) R% Z2 q
- 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)]';
" s1 j, G; a1 K* |% G - IKi(i,3*k-2)=E(1,1);+ P4 E: i( u* y' e3 \3 q
- IKi(i,3*k-1)=E(2,1);1 ?# o9 o- i5 b; G5 X8 O( R3 E
- IKi(i,3*k)=E(3,1);
7 H1 y( P9 }4 M: C: u - end
1 E! h5 g; G8 Y# `, S - end
% R# [4 }$ W+ ^* C! U - %前推电压 O% _' F& S1 N g
- for j=2:n
: x& z5 J% S' B: j1 j' ? W/ U c - U(1,3*k-2)=1.02;
S* j7 {" l7 s - U(1,3*k-1)=1.02*(sin(2*pi/3)+j*cos(2*pi/3));
9 G- U" n6 L4 h( t0 T - U(1,3*k)=1.02*(sin(4*pi/3)+j*cos(4*pi/3));" A& W* c0 d. l: h+ o- q; U4 C6 k; j
- a=DB(j-1,2);+ f1 I* S4 F. D
- if DB(j-1,6)==1
* ~* h; h) n1 `. M/ `1 J - Yi=1/DB(j-1,4)*eye(3);
1 N) U! v1 T( W; P( }$ |/ U - G1=[IKi(a,3*k-2) IKi(a,3*k-1) IKi(a,3*k)]';
/ \1 g5 L9 v; }0 d9 C - 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)]');
9 A3 n9 z r. r. ]8 A+ ~" j; Q - U(j,3*k+1)=G(1,1); ; }" }1 @! v' C& f
- U(j,3*k+2)=G(2,1);
! n+ m6 Q+ k" I3 j! ] - U(j,3*k+3)=G(3,1); 5 J' o" V1 p: t& O# \: _6 v
- else1 Y3 x, l4 G7 V b$ t
- YT=1/real(DB(i,4))+j*(1/imag(DB(i,4)));
9 L4 Z4 |5 S8 k" P - YI=YT*eye(3);
9 ^, D* m0 G/ x$ } - YII=YT*eye(3);6 ~3 s' W/ q5 X2 J/ u: Z
- YIII=-YT*eye(3);- T1 Y' M# k# i5 _4 `6 z7 U
- 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)]');
8 H1 w. l5 V8 _4 {& F( R2 k' M - U(j,3*k+1)=H(1,1);5 O, H# h* U6 h) h: ]. }
- U(j,3*k+2)=H(2,1);( ~+ w( o: b, J8 f# s
- U(j,3*k+3)=H(3,1);7 ?) \8 l- f& Y+ I/ t# a, m' a
- end & p9 \% F& B. E6 g* C$ P: [
- end
: I+ \7 D6 y* B+ S) g - dU1=U(:,3*k+3)-U(:,3*k); & ^6 B" w$ c3 ? Q! J/ J& c+ ]
- dU2=U(:,3*k+2)-U(:,3*k-1);
/ s$ u( e* A/ c7 W+ v3 s - dU3=U(:,3*k+1)-U(:,3*k-2);
9 s" I; j$ s# B/ Q6 P - if (max(abs(dU1))<0.001&&max(abs(dU2))<0.001&&max(abs(dU2))<0.001)
- b. b- \4 T m - break;
8 X8 b, Z, ]# H- ~) q! [! Y5 n+ | - else
6 }5 Y9 M$ D5 k8 \( l* S - k=k+1;
8 ?; M! O4 k% K$ O& | - end
+ H o& m2 Q5 e$ O8 j+ z - end
7 g$ F- c+ O+ c. C - U0 S: L- O6 K _
- IL3 }5 j+ x' x9 [9 w7 B( _) D9 {* M
- IKj. k5 T- c+ q e. x' H
- IKi
8 ~9 a9 U& l3 X u - k: d; P/ |# K* w- C6 d5 V6 z
复制代码 |
|