|
|
马上加入,结交更多好友,共享更多资料,让你轻松玩转电力研学社区!
您需要 登录 才可以下载或查看,没有账号?立即加入
×
各位~~我根据下图的一个系统用Matlab编写了一段程序,是电流型前推回代法的三相配电系统潮流计算,但是现在结果不是我想要的,不收敛。请各位帮我看一下,提提建议。图片见附件了。还有个问题:已知的线路都是三相对称线路的阻抗值,那么线路的阻抗矩阵中的互阻抗应该如何计算,是什么样的形式,应用怎样的公式?这里我不是很清楚。- %程序名:qiantuihuitui_I_3.m
" i5 M8 X0 z7 Y5 _) L: A6 D9 i! V - %功能:支路电流前推回推法求解潮流
6 s3 N( l, b' F, Y4 l - clc ' r" r) a' ~- a! y
- clear all;
g3 ?4 ?- B1 Q; h - %--------------输入网络参数--------------
% l$ z% y. N+ j# s, M1 q - %1-支路编号,2-首节点,3-尾节点,4-自阻抗,5-尾节点复功率,6-支路性质(1-馈线段支路,2-变压器支路),7-尾节点是否带负荷# w# L$ J, k2 f
- DB=[1 1 2 0.000167+j*0.000208 0.42+j*0.31 1 1
& x% A$ y# w$ H0 Q+ `5 y - 2 2 3 0.000151+j*0.000188 6.15 1 00 N2 z6 m( J4 @. c3 B8 |& z
- 3 2 4 0.000066+j*0.000082 0.38+j*0.29 1 1
# o6 b p' G" m' `: ] - 4 2 5 0.000249+j*0.000310 0 1 0
0 y0 z2 [+ v9 G) D$ O" T6 } - 5 2 6 0.000172+j*0.000215 0 1 0
/ V9 H# L# K/ z; p$ T6 ~ - 6 4 7 0.000156+j*0.000195 6.06 1 0" y* L* |+ S3 W, T& S0 U, U
- 7 4 8 0.000162+j*0.000202 6.04 1 0 # ^8 G8 G. ]5 F! @9 i
- 8 4 9 0.000345+j*0.000430 0 1 0 * R, n+ l. v% _9 {; S
- 9 4 10 0.000287+j*0.000358 0 1 05 J: ~% u+ p/ v
- 10 5 11 0.020563+j*0.321594 0 2 0 + [" C! R8 Y6 ]- c
- 11 6 12 0.020563+j*0.321594 0 2 0: O K; M9 i$ x8 P0 E4 C
- 12 9 13 0.020563+j*0.321594 0 2 0 . m9 o$ [* f% K
- 13 10 14 0.020563+j*0.321594 0 2 0
W( P0 ^* _% q: W3 z8 a - 14 11 15 0.000237+j*0.000408 5.72+j*0.12 1 1 - x4 N9 m+ K, z5 a( i1 n) {/ o
- 15 12 16 0.000237+j*0.000408 5.76+j*0.09 1 1
# q& J/ v& P9 d+ m2 d - 16 13 17 0.000292+j*0.000502 5.86+j*0.11 1 1- b! d6 V: {5 X+ c# Z! G
- 17 14 18 0.000274+j*0.000470 5.81+j*0.14 1 1];
6 W1 g; f# q; C - [n,m]=size(DB);
8 }2 n' e8 g3 K5 n5 j: n/ Q& t: [ - B=[1 sin(2*pi/3)+j*cos(2*pi/3) sin(4*pi/3)+j*cos(4*pi/3)];$ h( w3 q( ]) p, J5 P/ t6 P7 U
- C=[1.02 1 1.02 1 1 1 1.02 1.02 1 1 1 1 1 1 1 1 1 1];
" `; u4 c) u+ Y - U(:,1)=B(1,1)*C';
# N5 p; b4 x v" S - U(:,2)=B(1,2)*C';
" j5 b R7 Y1 f; y* [" X - U(:,3)=B(1,3)*C';
8 @- ], n" }5 }. ] - %-------------------求解潮流----------------- 0 _+ |6 ]4 w0 j( a2 y8 Z% O
- for k=1:15 ; w Z ^, G: c- y: h
- % I(:,k)=((DB(:,5).')*(diag(1./U(:,k))))'; - w/ F% c- r* P! ~% V4 X6 v/ |
- for i=n:-1:1 ! ]7 G9 `0 g. v6 n+ U
- %如果尾节点带恒功率负荷,需计算节点上负荷注入的电流
1 n. N' F4 H# Y S0 d1 r7 A( O - if DB(i,7)==19 A1 p( Y5 z) _6 X/ V1 r( X. |
- c=DB(i,5)/3;/ h* U0 W0 @4 b- m. @3 p: f: {4 D
- d=c/U(i+1,(3*k-2));
Z' ~+ `; f( r& c - IL(i,3*k-2)=conj(d);5 e1 {. R+ c4 a8 W
- d=c/U(i+1,3*k-1);
- R$ a1 m# U$ b$ d3 ?( M% T7 L* O( N - IL(i,3*k-1)=conj(d);
S9 u" l1 Y" ~: d5 Y# K, u9 m6 P' v - d=c/U(i+1,3*k);) a& v! y; y$ u7 }6 U. `7 ~* _7 l2 B* A
- IL(i,3*k)=conj(d);* w- d: {9 H2 K5 n$ k
- else+ }9 t* n! }* A' B9 ?3 \ I
- IL(i,3*k-2)=0;, _) v+ C* r1 E
- IL(i,3*k-1)=0;
4 `* T/ u4 N& {# ?7 C: i& q2 @0 | - IL(i,3*k)=0;
2 M) d7 ^ b" w1 S1 w2 m3 e - end
. w4 l( J8 V+ ^3 F1 b, ~6 N - %找出所有与尾节点相连的支路,计算进支电流(末端电流),存IKj. a( D* [8 k# ]
- A=(find(DB(:,2)==(i+1)));. v4 c# b! Z- J6 ]+ m1 T) v* N
- if isempty(A)~=17 B- ]. u% K/ R" b2 r* s( s' X& G' W
- IKj(i,3*k-2)=IL(i,3*k-2)-sum(IKj(A,3*k-2));6 J2 P2 w; p b0 S: J! z
- IKj(i,3*k-1)=IL(i,3*k-1)-sum(IKj(A,3*k-1));' k |2 B3 D9 A3 W& ?& X
- IKj(i,3*k)=IL(i,3*k)-sum(IKj(A,3*k));
, m& U& D' `+ T. U. ] - else
. U5 T: @: c% h; ` - IKj(i,3*k-2)=IL(i,3*k-2); r" d. g5 J; O1 G: {# a0 Q. }
- IKj(i,3*k-1)=IL(i,3*k-1);" P' ]) H0 f7 U8 b9 \
- IKj(i,3*k)=IL(i,3*k);, I i" B& }0 ? v7 @( a
- end4 E7 I7 ^, s# J; g& F4 u* S
- %计算出支电流(始端电流),存IKi; i0 }, h1 J6 { _ I: U5 w
- a=DB(i,2);0 R% l3 K9 y; C# H1 ?8 |% i
- b=DB(i,3);
/ |7 k9 N' C- ? W5 W - Yi=1/DB(i,4)*eye(3);8 C% {+ I0 |4 Y! i! T- J0 g
- if DB(A,6)==1( l. w" x% D( ?: n G: g
- 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)]';+ V; x: W R. `9 ~, J
- IKi(i,3*k-2)=F(1,1);
/ b- T. W% U- c& f; f+ U+ T* ?' R - IKi(i,3*k-1)=F(2,1);
; W( o* G4 B8 n* l& M2 e - IKi(i,3*k)=F(3,1);
( f& o2 {" B5 s9 N4 t, I - else
3 x2 u0 o$ g8 S; i - YT=1/real(DB(i,4))+j*(1/imag(DB(i,4)));6 z7 a) Z# s4 T& ]. ?9 J8 ]
- YI=YT*eye(3);* W+ p4 H1 m& r# \5 \
- YII=YT*eye(3);
$ \- C1 {' p* u" D( D+ } - YIII=-YT*eye(3);# J) B1 C& t8 j7 i
- 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)]'));- l5 c P9 ?7 P* P- V
- U(a,3*k+1)=D(1,1);3 l5 D/ Z" m( s. ]7 a
- U(a,3*k+2)=D(2,1);
, [3 J! f" M& q - U(a,3*k+3)=D(3,1);
8 {8 f7 {! }5 V |0 j9 Y4 Y - 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)]';2 K# {6 L, a# n$ b2 _! F* W" g' a) P
- IKi(i,3*k-2)=E(1,1);% J: n4 P4 o( O$ S7 e$ [& g# H* T
- IKi(i,3*k-1)=E(2,1);
# ]6 B' e3 O. u. l' C E - IKi(i,3*k)=E(3,1);
! {- V2 n: k' G' Q1 \; p - end
2 n4 { l, [. @( U0 Y1 B% L - end 7 {. G5 j7 b( I9 L7 W2 V: N! v
- %前推电压 ' \3 [7 t9 T* B
- for j=2:n
. r3 J: d! v7 u/ d* h! { - U(1,3*k-2)=1.02;& x) M/ U3 u5 S
- U(1,3*k-1)=1.02*(sin(2*pi/3)+j*cos(2*pi/3));! ^$ }* d: q( C0 y2 `
- U(1,3*k)=1.02*(sin(4*pi/3)+j*cos(4*pi/3));
% q- y; M( }9 C2 W" I) u% [ - a=DB(j-1,2);
3 Z$ z7 C: F* Q3 ~ - if DB(j-1,6)==1# a; x1 ]/ Y# u/ W0 \
- Yi=1/DB(j-1,4)*eye(3);; I- L8 v0 W, z+ l! ^/ X/ n
- G1=[IKi(a,3*k-2) IKi(a,3*k-1) IKi(a,3*k)]';
; i5 x- p* l, `: W+ E8 f - 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)]');
5 w4 l$ A$ I7 W0 c/ \2 w - U(j,3*k+1)=G(1,1);
; A3 d7 `$ ^/ }: Y$ {% Y - U(j,3*k+2)=G(2,1);
7 l. A7 D! j) T. ]+ E - U(j,3*k+3)=G(3,1);
+ g9 C$ x; i. a) a" s3 L - else
" m; [! X/ z: N- t+ h9 b' [ - YT=1/real(DB(i,4))+j*(1/imag(DB(i,4)));' b, y' b3 A! G* y# S% z' D
- YI=YT*eye(3);# U. ~" x" p" G' ]7 A
- YII=YT*eye(3);
& j& C% q5 I6 y u3 u8 \7 b' W - YIII=-YT*eye(3);# j2 K0 Q6 |! f9 R
- 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)]');3 _2 M! z8 N. _8 e1 b5 c
- U(j,3*k+1)=H(1,1);6 z# |7 _; `8 [7 a* y
- U(j,3*k+2)=H(2,1);* }: A& a" I, b/ q" {* J
- U(j,3*k+3)=H(3,1);
4 h/ T/ ]( }; b7 a( M - end
2 g3 s9 @7 U" h( M, e - end) H; a5 t- ~* w4 S; U; A. c
- dU1=U(:,3*k+3)-U(:,3*k);
4 ]. b( R; u3 f1 l - dU2=U(:,3*k+2)-U(:,3*k-1); , k' g; M; G: A5 D8 ^
- dU3=U(:,3*k+1)-U(:,3*k-2); 0 i/ a% I! N0 [
- if (max(abs(dU1))<0.001&&max(abs(dU2))<0.001&&max(abs(dU2))<0.001)
T( A6 p/ [7 M9 y* L4 _ - break; ' }2 b* t- z3 A, t' U
- else
% M" {9 g6 I z+ i/ ^: |/ ~ - k=k+1;
b5 W* ^: A7 z9 C- R! Q - end
+ z1 @+ [7 W g" h3 j - end
: ?5 n9 b5 h( c" V8 a- x - U6 W' ~) d: L3 c
- IL
( g# v3 ~$ A4 c/ A5 k% `' \' K% G - IKj4 b, h. r, V! S/ g8 c- R
- IKi
) ?# Z% C" B5 k: c9 F - k
( |( z: ~1 f2 D- I
复制代码 |
|