电力系统及其自动化--潮流计算--咻咻永远

咻咻永远 · 发表于 2013-5-21 16:12:22

马上加入，结交更多好友，共享更多资料，让你轻松玩转电力研学社区！

您需要登录才可以下载或查看，没有账号？立即加入

×

各位，我想请教一个问题，我现在做的船舶交流电力系统，想要仿真一个18节点破冰船舶电力系统，主要参考的文献为《船舶电力系统拓扑分析和潮流计算研究-周容华》，需要的是3相潮流计算，但是现在编写的程序不收敛，想请大家帮我看看问题主要出在哪里？源代码在下边。

%程序名：qiantuihuitui_I_3.m - ]" _% {$ ], y" P% @- E1 V
%功能：支路电流前推回推法求解潮流 # z% J' ~7 _' B$ E
clc
; [- ? i V( b
clear all; " n0 }, `/ a8 g1 \5 f
%--------------输入网络参数-------------- ( T2 a3 v% O J" e2 @5 a1 r3 P
%1-支路编号，2-首节点，3-尾节点，4-自阻抗，5-尾节点复功率，6-支路性质（1-馈线段支路，2-变压器支路）,7-尾节点是否带负荷
) M3 R2 ?$ j! J. Q6 j; { m
DB=[1 1 2 0.000167+j*0.000208 0.42+j*0.31 1 1' o( @3 [, o" O2 j3 ^' v) X
2 2 3 0.000151+j*0.000188 6.15 1 02 P1 Y7 P2 L' U- V F- ^0 D
3 2 4 0.000066+j*0.000082 0.38+j*0.29 1 1
5 h0 s# b8 A+ N3 R
4 2 5 0.000249+j*0.000310 0 1 0
" a& E2 _* ^6 l
5 2 6 0.000172+j*0.000215 0 1 0* j2 _1 s3 C/ h: e) q4 [" d
6 4 7 0.000156+j*0.000195 6.06 1 0
1 u* y1 Y" }; U# d) a5 D: @
7 4 8 0.000162+j*0.000202 6.04 1 0 : s! g0 }2 @, z4 f, ?* Y: F* B) Z& l
8 4 9 0.000345+j*0.000430 0 1 0 ' x' Q0 [/ U7 G
9 4 10 0.000287+j*0.000358 0 1 0) I; w/ {* O: D& P& D7 t) E" A
10 5 11 0.020563+j*0.321594 0 2 0 & V. t% n# d2 S2 g$ c
11 6 12 0.020563+j*0.321594 0 2 0
; {: K' O0 Y Z0 M
12 9 13 0.020563+j*0.321594 0 2 0 + |6 N0 t$ s0 h( p6 S/ ]' C
13 10 14 0.020563+j*0.321594 0 2 0 5 X# s" m+ d' F
14 11 15 0.000237+j*0.000408 5.72+j*0.12 1 1
2 I2 a$ R6 {4 J) p9 g. B/ m
15 12 16 0.000237+j*0.000408 5.76+j*0.09 1 1$ g; ^6 D I, h" B8 b
16 13 17 0.000292+j*0.000502 5.86+j*0.11 1 1
! {5 m1 A6 ?# [6 G f
17 14 18 0.000274+j*0.000470 5.81+j*0.14 1 1]; 4 v! Y: A& T4 Q) Y" X0 B6 A
[n,m]=size(DB); * b! o% E6 Z/ j+ w5 ?5 R# [
B=[1 sin(2*pi/3)+j*cos(2*pi/3) sin(4*pi/3)+j*cos(4*pi/3)];$ v, F1 p$ P) Y; e$ K
C=[1.02 1 1.02 1 1 1 1.02 1.02 1 1 1 1 1 1 1 1 1 1];
% x9 _6 G6 r2 o# p& T8 e' a6 X
U(:,1)=B(1,1)*C';
. `( j, D0 w) P9 Z
U(:,2)=B(1,2)*C';
, r) t* l) C" {( h. P: z
U(:,3)=B(1,3)*C';5 ]+ ^9 V4 B/ z! k3 M, ^& `
%-------------------求解潮流-----------------
0 ^) ]: M+ T, h* @0 [
for k=1:15
4 |& O6 ]6 o s9 b2 W0 {" u3 x- G
% I(:,k)=((DB(:,5).')*(diag(1./U(:,k))))'; : m8 e, t, t5 O; M% q; g' X# r W* ]) L
for i=n:-1:1 2 Z/ v8 ?7 g, g/ y# H
%如果尾节点带恒功率负荷，需计算节点上负荷注入的电流
3 F3 ~3 f1 E& [+ V, i
if DB(i,7)==11 Y# S2 d3 M2 V; L, R
c=DB(i,5)/3;
: B+ Y* h7 l$ u" ^) @2 x
d=c/U(i+1,(3*k-2));
8 u* {* h9 [/ W' C1 m. j
IL(i,3*k-2)=conj(d);- {2 j. i( t6 ~- U: ~
d=c/U(i+1,3*k-1);9 C; P- f4 m+ F1 J1 i
IL(i,3*k-1)=conj(d);# d) K8 U' Z$ X/ U1 c
d=c/U(i+1,3*k);8 D( X/ n' q2 Z
IL(i,3*k)=conj(d);2 X* ^5 u& I7 a @: u
else7 n" t' ^. ?# D+ I8 O
IL(i,3*k-2)=0;+ X$ s6 d9 h6 G% ?1 S! H! ?5 Z! }) c
IL(i,3*k-1)=0;
; B, y! K) E P
IL(i,3*k)=0;
0 A: {: h0 m! V
end1 |7 H& o; ]+ E3 |/ e
%找出所有与尾节点相连的支路,计算进支电流（末端电流），存IKj
4 G: j7 U6 d9 p0 i9 q
A=(find(DB(:,2)==(i+1)));
* h3 v& C' I: l# s# o
if isempty(A)~=1. d4 j7 [9 |8 b# {* j; H# r
IKj(i,3*k-2)=IL(i,3*k-2)-sum(IKj(A,3*k-2));
~! D" \. N( t- _) j8 U
IKj(i,3*k-1)=IL(i,3*k-1)-sum(IKj(A,3*k-1));
2 ?2 r5 v2 s' J6 f# A, _
IKj(i,3*k)=IL(i,3*k)-sum(IKj(A,3*k));
" L. |) c9 B/ P, g& q- ~+ t
else' ?- ~) g" a6 c) W* Q" ^
IKj(i,3*k-2)=IL(i,3*k-2);1 E. T: F L8 L% o
IKj(i,3*k-1)=IL(i,3*k-1);
6 C" Q0 j% [; _, {/ H( s7 s
IKj(i,3*k)=IL(i,3*k);7 _0 }+ s( Z6 |2 [6 x" W7 p
end& _6 {" `) T3 f+ g& z6 G6 {; g; q$ n
%计算出支电流（始端电流），存IKi
3 q+ `6 I: l i( g5 w) U6 G- [' w
a=DB(i,2);
& ^9 a4 u- N) U0 V7 W/ D
b=DB(i,3);+ L( Q" s/ `* ^" \* k
Yi=1/DB(i,4)*eye(3);/ B3 [' r' i& j/ k2 u3 l9 w& }
if DB(A,6)==1
4 O3 u7 a' t, P
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)]';
' | G4 z) O9 G6 A2 n
IKi(i,3*k-2)=F(1,1);
% T8 [9 m/ B2 H
IKi(i,3*k-1)=F(2,1);
5 y% D, U+ W# F
IKi(i,3*k)=F(3,1);2 `: X/ M; v; @0 a
else. W, ^$ A/ }, Q' G: W
YT=1/real(DB(i,4))+j*(1/imag(DB(i,4)));* s* f6 b0 t Q8 E# d2 |$ u# U, R5 W
YI=YT*eye(3);
& p3 ~ R/ T2 p+ o
YII=YT*eye(3);3 H2 T& F0 B# [) H
YIII=-YT*eye(3);
' ]9 o; D2 |9 t* c6 m
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)]'));3 }0 g% v y" x' d4 Z
U(a,3*k+1)=D(1,1);" t, |5 u$ H( D) Y$ ~
U(a,3*k+2)=D(2,1);4 i! Z. w, B% @
U(a,3*k+3)=D(3,1);
6 C. _* N) E$ Q# C9 o
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)]';% i+ p o) l3 ?& v
IKi(i,3*k-2)=E(1,1);
0 x5 o( k Y: L9 G2 E4 k4 ^
IKi(i,3*k-1)=E(2,1);2 A4 b* i4 U; Y7 s/ _
IKi(i,3*k)=E(3,1);
! E; a) b% K7 J$ l* e/ p2 Y' p- {
end
1 u: Y1 @& f9 F) @9 @! l1 y
end ) L( ] t. B5 p# m3 t
%前推电压
, {' _6 \5 T$ k8 d
for j=2:n
" f* P3 U2 b; J1 f% d9 C
U(1,3*k-2)=1.02;$ }6 I0 n* D8 y9 i
U(1,3*k-1)=1.02*(sin(2*pi/3)+j*cos(2*pi/3));
7 b8 i" r- }8 V4 d, }
U(1,3*k)=1.02*(sin(4*pi/3)+j*cos(4*pi/3));
' w, v+ O, z6 L3 F' O9 D
a=DB(j-1,2);2 l8 D7 z9 N; r2 E' s! T
if DB(j-1,6)==1
' \5 Y( [$ r# S
Yi=1/DB(j-1,4)*eye(3);
5 ^+ K& T# h1 I1 }4 P# Z
G1=[IKi(a,3*k-2) IKi(a,3*k-1) IKi(a,3*k)]';
8 N+ e2 v" n4 B$ Y* m2 C+ T
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)]');
1 |4 F0 Q' d) T+ H: u* j) E, l2 C0 M
U(j,3*k+1)=G(1,1);
8 o) t6 Y1 w( @) I/ b
U(j,3*k+2)=G(2,1);- q9 b2 u; C) ?& @5 @6 ^0 e' Z
U(j,3*k+3)=G(3,1);
( z C! m3 p7 G9 Y
else
" `, @- K1 ^# ^- S3 y; \
YT=1/real(DB(i,4))+j*(1/imag(DB(i,4)));
6 ^ u6 [/ O2 u. \
YI=YT*eye(3);" Q) y& r3 C0 b2 @+ U
YII=YT*eye(3);
$ ]1 @& J# d2 r4 V
YIII=-YT*eye(3);
+ P: i( S4 f. A9 g* f' 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)]');: |( D' B$ s4 a; p- h/ z1 q
U(j,3*k+1)=H(1,1);
# J5 x( X( I' N+ ^$ `
U(j,3*k+2)=H(2,1);1 P2 H1 J2 a5 c$ D2 }# B- }; O
U(j,3*k+3)=H(3,1);: s0 l% G1 @ r6 h" z
end , Z2 |/ Y! z8 h$ b& ^2 J$ `! I4 d
end0 M$ z# t- W& ^/ A7 X, h+ t
dU1=U(:,3*k+3)-U(:,3*k);
+ e) U" B1 O; I! I( i- p) z% c
dU2=U(:,3*k+2)-U(:,3*k-1); 5 y$ t! {% T2 l* G8 r1 S
dU3=U(:,3*k+1)-U(:,3*k-2); ( L# Z1 P: v# |& L) b
if (max(abs(dU1))<0.001&&max(abs(dU2))<0.001&&max(abs(dU2))<0.001)
1 p" |7 `. c5 n4 U
break; 4 j: G3 \* v. t
else
, q5 S* L$ Z, o. v+ L2 v) P
k=k+1; - r* b! i+ [+ }7 g( e, N
end 1 G3 J# ?8 H6 s* }
end
( q* g7 I) z: d3 y o! d
U: \/ M% P6 ?* w j; \! g
IL9 x& m$ J1 C5 J: Y8 S
IKj# _; r- W5 w/ U, v" _2 H1 W
IKi3 K1 { [, e0 d' d% _
k
! `% C8 p- @) e( L

复制代码

sun710504 · 发表于 2013-5-22 13:11:42

学习学习

账号		自动登录	找回密码
密码			立即加入

电力系统及其自动化--潮流计算--咻咻永远

马上加入，结交更多好友，共享更多资料，让你轻松玩转电力研学社区！

相关帖子

浏览过的版块