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

咻咻永远

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

1. %程序名：qiantuihuitui_I_3.m
   
2. %功能：支路电流前推回推法求解潮流
   
3. clc
   " n! v$ K; h( T) v, b6 w: Q! P: t  W
4. clear all;
   ( X; p. [. z* ?3 V9 ^4 p
5. %--------------输入网络参数--------------
   8 R( Y* Z5 s+ [3 l! Y
6. %1-支路编号，2-首节点，3-尾节点，4-自阻抗，5-尾节点复功率，6-支路性质（1-馈线段支路，2-变压器支路）,7-尾节点是否带负荷
   * K2 y8 u+ {, B' B5 X" F& o
7. DB=[1 1 2 0.000167+j*0.000208 0.42+j*0.31 1 1
   
8.     2 2 3 0.000151+j*0.000188 6.15 1 0
   4 o8 a, A, |' P* f) t: O; n1 C  z8 [
9.     3 2 4 0.000066+j*0.000082 0.38+j*0.29 1 1
   
10.     4 2 5 0.000249+j*0.000310 0 1 0
    
11.     5 2 6 0.000172+j*0.000215 0 1 0
    
12.     6 4 7 0.000156+j*0.000195 6.06 1 0
    
13.     7 4 8 0.000162+j*0.000202 6.04 1 0
    ; P4 \# e. K1 E  O3 O
14.     8 4 9 0.000345+j*0.000430 0 1 0
    
15.     9 4 10 0.000287+j*0.000358 0 1 0
    
16.     10  5 11 0.020563+j*0.321594 0 2 0
    5 z# Q- ]5 \/ y4 }
17.     11  6 12 0.020563+j*0.321594 0 2 0
      w, N! B' L& ~; o
18.     12  9 13 0.020563+j*0.321594 0 2 0
    
19.     13  10 14 0.020563+j*0.321594 0 2 0  
    ( ~! |: [# D! E& Q+ Z  ~3 k- D
20.     14  11 15 0.000237+j*0.000408 5.72+j*0.12 1 1
    
21.     15  12 16 0.000237+j*0.000408 5.76+j*0.09 1 1
    
22.     16  13 17 0.000292+j*0.000502 5.86+j*0.11 1 1
    
23.     17  14 18 0.000274+j*0.000470 5.81+j*0.14 1 1];
    * A# U5 x4 m) ~: D1 l( V) j
24. [n,m]=size(DB);  
    3 T9 m) P! ~  i. n, `2 s3 o4 }# ^7 E' {
25. B=[1 sin(2*pi/3)+j*cos(2*pi/3) sin(4*pi/3)+j*cos(4*pi/3)];
    
26. C=[1.02 1 1.02 1 1 1 1.02 1.02 1 1 1 1 1 1 1 1 1 1];
    
27. U(:,1)=B(1,1)*C';
    
28. U(:,2)=B(1,2)*C';
    0 ^- N' A* ]# V* h6 D! k% i
29. U(:,3)=B(1,3)*C';
    # G) }! C  `& D+ |& u' w
30. %-------------------求解潮流-----------------
    
31. for k=1:15
    0 E: ]1 f# d; I* d, [; t
32. %   I(:,k)=((DB(:,5).')*(diag(1./U(:,k))))';
    
33. for i=n:-1:1
    6 A7 [! G5 b. @9 \
34.   %如果尾节点带恒功率负荷，需计算节点上负荷注入的电流
    
35. if DB(i,7)==1
    / s% v+ k* w0 ~6 P! Y2 ~
36. c=DB(i,5)/3;
    & h" d; ^! J5 l2 k. m+ ~
37. d=c/U(i+1,(3*k-2));
    ' ]1 U! w/ z' ]/ g: _. q. k6 v
38. IL(i,3*k-2)=conj(d);
    7 o7 _6 P( R; u
39. d=c/U(i+1,3*k-1);
    
40. IL(i,3*k-1)=conj(d);
    
41. d=c/U(i+1,3*k);
    + {" b8 l& e. u1 k" R  \
42. IL(i,3*k)=conj(d);
    ; s: N* i1 @! d; K6 q" c+ h
43. else
    
44. IL(i,3*k-2)=0;
    
45. IL(i,3*k-1)=0;
    4 Y# |2 j1 d) h
46. IL(i,3*k)=0;
    
47. end
    
48. %找出所有与尾节点相连的支路,计算进支电流（末端电流），存IKj
    
49. A=(find(DB(:,2)==(i+1)));
    
50. if isempty(A)~=1
    
51. IKj(i,3*k-2)=IL(i,3*k-2)-sum(IKj(A,3*k-2));
    
52. IKj(i,3*k-1)=IL(i,3*k-1)-sum(IKj(A,3*k-1));
    ) n6 b3 n3 A9 M" H
53. IKj(i,3*k)=IL(i,3*k)-sum(IKj(A,3*k));
    
54. else
    
55. IKj(i,3*k-2)=IL(i,3*k-2);
    4 h4 q* J$ G7 ?% T6 {3 J  M" b9 ]
56. IKj(i,3*k-1)=IL(i,3*k-1);
    . \/ s$ T  j/ V' N' {
57. IKj(i,3*k)=IL(i,3*k);
    7 r$ p& J" X! ]* L+ P  g% E
58. end
    
59. %计算出支电流（始端电流），存IKi
    ) l  n$ H/ J: ?  e& B! ^
60. a=DB(i,2);
    2 _7 a; G1 k7 \/ v
61. b=DB(i,3);
    
62. Yi=1/DB(i,4)*eye(3);
    8 J0 q( i" S' @+ P( d! c
63. if DB(A,6)==1
    + @  [* H3 |7 a3 x- H
64. 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)]';
    
65. IKi(i,3*k-2)=F(1,1);
    " [2 W3 M# F7 [8 E+ {& M& @* a# R
66. IKi(i,3*k-1)=F(2,1);
    - Z, G3 @8 K5 H9 [  N7 B1 ?
67. IKi(i,3*k)=F(3,1);
    
68. else
    
69. YT=1/real(DB(i,4))+j*(1/imag(DB(i,4)));
    
70. YI=YT*eye(3);
    
71. YII=YT*eye(3);
    
72. YIII=-YT*eye(3);
    / J( K/ n& m4 q' q9 \% `6 r
73. 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)]'));
    2 f$ o4 d0 R+ P* |! W
74. U(a,3*k+1)=D(1,1);
    
75. U(a,3*k+2)=D(2,1);
    
76. U(a,3*k+3)=D(3,1);
    
77. 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)]';
    $ v5 M! I+ }) `+ _0 t, Q5 f( n* ~
78. IKi(i,3*k-2)=E(1,1);
    
79. IKi(i,3*k-1)=E(2,1);
    
80. IKi(i,3*k)=E(3,1);
    
81. end
    
82. end
    
83. %前推电压
    
84. for j=2:n
    4 u5 J$ j: k+ q8 l. g  @
85. U(1,3*k-2)=1.02;
    ( T" F, n" d  d. w$ A
86. U(1,3*k-1)=1.02*(sin(2*pi/3)+j*cos(2*pi/3));
    
87. U(1,3*k)=1.02*(sin(4*pi/3)+j*cos(4*pi/3));
    6 `0 T0 j$ o' i# J: b& {' j
88. a=DB(j-1,2);
    
89. if DB(j-1,6)==1
    
90. Yi=1/DB(j-1,4)*eye(3);
    
91. G1=[IKi(a,3*k-2) IKi(a,3*k-1) IKi(a,3*k)]';
    - p! E) t1 b& \8 R- ^; |
92. 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)]');
    
93. U(j,3*k+1)=G(1,1);
    % T) H/ g: N8 n: p0 Y  V3 a
94. U(j,3*k+2)=G(2,1);
    ( _* [# n' J; z! {9 |. \) w
95. U(j,3*k+3)=G(3,1);
    
96. else
    8 d) ?0 J  q5 `: @9 r3 @
97. YT=1/real(DB(i,4))+j*(1/imag(DB(i,4)));
    
98. YI=YT*eye(3);
    0 I& U' l5 E" a9 {
99. YII=YT*eye(3);
    
100. YIII=-YT*eye(3);
     
101. 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)]');
     
102. U(j,3*k+1)=H(1,1);
     
103. U(j,3*k+2)=H(2,1);
     4 K+ |. J  G3 c% F9 r
104. U(j,3*k+3)=H(3,1);
     $ ^7 K: J0 N( m0 Y
105. end
     
106. end
     
107. dU1=U(:,3*k+3)-U(:,3*k);
     
108. dU2=U(:,3*k+2)-U(:,3*k-1);
     5 {& x* ?! ^" d+ u& h8 J0 {. p* g9 t
109. dU3=U(:,3*k+1)-U(:,3*k-2);
     , D5 Y' p4 ~$ {. W7 m; [
110. if (max(abs(dU1))<0.001&&max(abs(dU2))<0.001&&max(abs(dU2))<0.001)
     ; E* l& O4 f+ ]' S; P
111. break;
     " O+ }2 Q% Q+ X2 L' O4 G# e
112. else            
     " e" _0 A& M4 f
113. k=k+1;
     % I0 O* J5 G6 [. e% o
114. end
     
115. end
     / _4 u# O. Z1 d4 q: F
116. U
     
117. IL
     % t" c9 ~, |' f  g
118. IKj
     
119. IKi
     / `" b8 A1 [+ c
120. k
     

复制代码

sun710504

学习学习