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

咻咻永远

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

1. %程序名：qiantuihuitui_I_3.m
   0 Y& o& n6 Z4 p- g' l4 R  T
2. %功能：支路电流前推回推法求解潮流
   
3. clc
   # z3 \5 ~; {; }. H
4. clear all;
   
5. %--------------输入网络参数--------------
   
6. %1-支路编号，2-首节点，3-尾节点，4-自阻抗，5-尾节点复功率，6-支路性质（1-馈线段支路，2-变压器支路）,7-尾节点是否带负荷
   % p, O0 b5 h$ x3 H+ R2 H. X4 K8 C
7. DB=[1 1 2 0.000167+j*0.000208 0.42+j*0.31 1 1
   " {3 C7 X/ M8 k7 L! N
8.     2 2 3 0.000151+j*0.000188 6.15 1 0
   
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
    5 }) C7 T. `6 ?9 b
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
    ! [$ X/ U; N( B& D+ W: y
14.     8 4 9 0.000345+j*0.000430 0 1 0
    
15.     9 4 10 0.000287+j*0.000358 0 1 0
    ! G& ^1 y0 u, p- R! b6 q
16.     10  5 11 0.020563+j*0.321594 0 2 0
    4 Y% U) e$ a; p3 S
17.     11  6 12 0.020563+j*0.321594 0 2 0
    
18.     12  9 13 0.020563+j*0.321594 0 2 0
    
19.     13  10 14 0.020563+j*0.321594 0 2 0  
    
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
    ) L% g# S7 M9 g7 U/ C4 V1 t: j
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];
    : }' f! R: {% H' V% \
24. [n,m]=size(DB);  
    
25. B=[1 sin(2*pi/3)+j*cos(2*pi/3) sin(4*pi/3)+j*cos(4*pi/3)];
    ; A7 _+ u5 u0 \
26. C=[1.02 1 1.02 1 1 1 1.02 1.02 1 1 1 1 1 1 1 1 1 1];
    ; k+ f# |' K) e
27. U(:,1)=B(1,1)*C';
    
28. U(:,2)=B(1,2)*C';
    , U8 G8 a/ W3 H. z. i; c0 s9 _
29. U(:,3)=B(1,3)*C';
    # k9 E7 u  k4 f% a6 P# y
30. %-------------------求解潮流-----------------
    
31. for k=1:15
    
32. %   I(:,k)=((DB(:,5).')*(diag(1./U(:,k))))';
      F, c7 j3 M1 ~5 S1 L
33. for i=n:-1:1
    
34.   %如果尾节点带恒功率负荷，需计算节点上负荷注入的电流
    3 j: M0 x( Z$ B% I- Q0 U" ^* {+ B
35. if DB(i,7)==1
    
36. c=DB(i,5)/3;
    5 H- S1 M( C6 N  {( h9 S
37. d=c/U(i+1,(3*k-2));
    8 |  b3 r# P% d9 M8 }' r) X
38. IL(i,3*k-2)=conj(d);
    
39. d=c/U(i+1,3*k-1);
      t& u6 b6 {8 N  W3 X' X
40. IL(i,3*k-1)=conj(d);
    
41. d=c/U(i+1,3*k);
    
42. IL(i,3*k)=conj(d);
    
43. else
    
44. IL(i,3*k-2)=0;
    ! c6 j$ v" v# r
45. IL(i,3*k-1)=0;
    
46. IL(i,3*k)=0;
    
47. end
    6 A6 }% G( m4 U( O* X
48. %找出所有与尾节点相连的支路,计算进支电流（末端电流），存IKj
    
49. A=(find(DB(:,2)==(i+1)));
    , u" L0 P3 d( C7 H
50. if isempty(A)~=1
    7 {+ {9 p4 H! n; U7 q! R3 j. H
51. IKj(i,3*k-2)=IL(i,3*k-2)-sum(IKj(A,3*k-2));
    3 s* @$ y) f1 ?6 W
52. IKj(i,3*k-1)=IL(i,3*k-1)-sum(IKj(A,3*k-1));
    
53. IKj(i,3*k)=IL(i,3*k)-sum(IKj(A,3*k));
    & O3 N% Z4 |8 D# C* j; E& O: P3 h
54. else
    ' V* f7 r! _4 I1 U2 Z4 E+ Y+ f
55. IKj(i,3*k-2)=IL(i,3*k-2);
    
56. IKj(i,3*k-1)=IL(i,3*k-1);
    
57. IKj(i,3*k)=IL(i,3*k);
    
58. end
    ; o2 g/ n9 y/ S2 ]2 Z
59. %计算出支电流（始端电流），存IKi
    
60. a=DB(i,2);
    5 Y* w5 i3 v3 c( a5 L0 g7 O
61. b=DB(i,3);
    ( `& m% y( p' U2 }: G7 ?" A: B
62. Yi=1/DB(i,4)*eye(3);
    
63. if DB(A,6)==1
    # c4 y! A$ ?: ?, \: e# {9 ?- Y. G' X5 Z
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)]';
    2 W- P: C$ n7 E
65. IKi(i,3*k-2)=F(1,1);
    $ X1 {/ H% @8 f. e$ Y: a! {" W
66. IKi(i,3*k-1)=F(2,1);
    
67. IKi(i,3*k)=F(3,1);
    $ o- w! A1 r! R9 W% y' p
68. else
    
69. YT=1/real(DB(i,4))+j*(1/imag(DB(i,4)));
    ' N. [+ R1 g$ y8 Y1 T# `6 b4 i
70. YI=YT*eye(3);
    
71. YII=YT*eye(3);
    $ R4 e/ Z& m% `1 |9 m/ x
72. YIII=-YT*eye(3);
    ( F, R8 p) p  z  o' F
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)]'));
    
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)]';
    
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);
    + k5 X7 @2 x) Z- w: @/ D
81. end
    
82. end
    % r3 s) V# p; X' e, r& K
83. %前推电压
    
84. for j=2:n
    
85. U(1,3*k-2)=1.02;
    / e3 K6 C% P1 w  E
86. U(1,3*k-1)=1.02*(sin(2*pi/3)+j*cos(2*pi/3));
    * R2 q& ?0 O7 r; g8 S: P
87. U(1,3*k)=1.02*(sin(4*pi/3)+j*cos(4*pi/3));
    
88. a=DB(j-1,2);
    : ~9 \, K4 T) v2 ]$ d7 n/ y
89. if DB(j-1,6)==1
    ' _$ N4 n6 c1 y! _1 w
90. Yi=1/DB(j-1,4)*eye(3);
    5 I8 M! t- ^% L1 F# k4 u9 p
91. G1=[IKi(a,3*k-2) IKi(a,3*k-1) IKi(a,3*k)]';
    
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);
    3 A! M0 i) H/ b
94. U(j,3*k+2)=G(2,1);
    8 i+ D* m6 F2 t1 }! I
95. U(j,3*k+3)=G(3,1);
    8 u8 Q1 x* s4 i* e: h% \
96. else
    6 `' j' C) ]& c) j- o
97. YT=1/real(DB(i,4))+j*(1/imag(DB(i,4)));
    4 u- R' M( f7 e1 W1 {; n% u
98. YI=YT*eye(3);
    
99. YII=YT*eye(3);
    
100. YIII=-YT*eye(3);
     0 C, h8 |0 D. n. T3 V$ d
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)]');
     , V3 H+ v- z0 v6 I
102. U(j,3*k+1)=H(1,1);
     
103. U(j,3*k+2)=H(2,1);
     
104. U(j,3*k+3)=H(3,1);
     
105. end
     
106. end
     # I/ ~3 O$ V  f! S/ s* z. v
107. dU1=U(:,3*k+3)-U(:,3*k);
     
108. dU2=U(:,3*k+2)-U(:,3*k-1);
     ) t6 h, t/ q" m  `' a+ l# a
109. dU3=U(:,3*k+1)-U(:,3*k-2);
     
110. if (max(abs(dU1))<0.001&&max(abs(dU2))<0.001&&max(abs(dU2))<0.001)
     
111. break;
     
112. else            
     
113. k=k+1;
     + ^& x  z8 E+ m7 J
114. end
     ' i: ?& W& O$ r# D0 ^9 g
115. end
     / [; V3 ]2 ]0 _5 M, _" g: K
116. U
     
117. IL
     
118. IKj
     * I6 b4 |; ^2 u' K% Z- [
119. IKi
     1 V. Z( \( ?5 L" a2 d: ~* i1 @. M
120. k
     

复制代码

sun710504

学习学习