请教一个关于Matlab编写的IEEE14节点潮流计算的问题

迪拉克之海

我参考了一些书籍，用for循环写了个基于极坐标牛顿-拉夫逊法的IEEE14节点潮流计算代码，执行后发现迭代的最终结果与IEEE14给出的结果出入较大。。。
  i! K, `- u8 ?4 U这是我的毕业设计的一部分，在计算结果不能够统一的问题上我卡了一星期了，非常着急，自己盯着自己的代码看又死活看不出来什么地方出了问题，因此发到这里来，想借助各位高手的力量、帮我看看我到底什么地方出错了，非常感谢！

[问题已经解决]

. M. U. o. _. C: G[ 本帖最后由 迪拉克之海 于 2008-5-7 13:37 编辑 ]

xiao1806

我也曾经用matlab做过用极坐标的牛拉法算潮流，遇到跟你一样的问题，建议你通过手算潮流，看你的第一次迭代是否正确，其实第一次迭代手算并不麻烦，因为很多值都是0，还有就是看你的程序是否是收敛在六次以内的。手算是调试程序最好的方法，不要怕麻烦嘛。

flyingclouds

仔细对照电力系统的书看一遍，看看是不是有公式输错，并联补偿支路没有考虑的情况，不过最后的结果确实与IEEE会有一些差别

lgs_0003

这个问题我也遇到了，迭代结果是无法计算得了，我也正有想法用手工算，不过现在手工算的话，是有点麻烦，所以我在准备想写一个只迭代一次，而且把每时刻计算的结果输出来对比一下咯

yingfengli

u =
' ?# x* V8 C9 l, C2 p
  Columns 1 through 11
4 J& R. L8 K( `( t4 Y9 h; x: t$ C$ u( F
    1.0516    1.0474    1.0550    1.0549    1.0497    1.0327    1.0581    1.0123    1.0138    1.0450    1.0100
/ z' Q9 R9 y) y
7 e$ j! c8 J' C  i9 g' y2 @5 i  Columns 12 through 14

    1.0900    1.0700    1.0600
9 u8 k0 }  _$ p  |
; ?- g9 k% C' n- o( Q% s, m
6 u+ C; P" D) ^$ t: |Theta =
9 e0 N1 Q1 F$ C8 T! {
  Columns 1 through 11
! L/ S4 n' `  w& K+ B
" H6 x  ?; m( T! G' ^! F. S7 x$ Q9 H. U   -0.2599   -0.2628   -0.2578   -0.2632   -0.2645   -0.2794   -0.2323   -0.1790   -0.1520   -0.0872   -0.2229

  Columns 12 through 14

   -0.2323   -0.2483         0


Sslack =

2 c% i  q! ^6 @5 B+ r- Q! d+ U   2.3249 - 0.0250i

8 X4 O) ?, B* t9 X7 k& Z
S =

# n5 |8 k8 F% N2 K/ o9 P  Columns 1 through 5

        0             0.0507 + 0.0332i        0                  0                  0         
  -0.0506 - 0.0329i        0            -0.0394 - 0.0251i        0                  0         
% V) |& a/ a8 E        0             0.0396 + 0.0254i        0                  0                  0         
        0                  0                  0                  0             0.0166 + 0.0087i
        0                  0                  0            -0.0165 - 0.0086i        0         
  -0.0917 - 0.0280i        0                  0                  0            -0.0573 - 0.0220i
   0.2789 + 0.0672i        0                  0                  0                  0         
! S# ^' s) A3 w+ c- m2 B   0.1595 - 0.0673i        0                  0                  0                  0         
        0                  0                  0                  0                  0         
" M2 i/ d9 ?# x& ?' G        0                  0                  0                  0                  0         
( X9 r9 q0 P* b& `        0                  0                  0                  0                  0         
( ^9 T% q5 p! X0 J        0                  0                  0                  0                  0         
4 p# D: ~2 J8 d) Y) i        0                  0             0.0752 + 0.0448i   0.0784 + 0.0262i   0.1785 + 0.0769i
        0                  0                  0                  0                  0         
, a; U2 O- Z; y; q
  Columns 6 through 10

   0.0928 + 0.0303i  -0.2789 - 0.0591i  -0.1595 + 0.0831i        0                  0         
3 n( q! y/ A4 T  L; ?$ l        0                  0                  0                  0                  0         
        0                  0                  0                  0                  0         
- O3 @! e" k; ~6 U: X" H5 m        0                  0                  0                  0                  0         
% Q. ^/ @* F( I! u6 P# @0 r   0.0579 + 0.0232i        0                  0                  0                  0         
        0                  0                  0                  0                  0         
        0                  0            -0.2789 + 0.2447i        0                  0         
        0             0.2789 - 0.2196i        0            -0.6058 + 0.1645i  -0.5446 + 0.0169i
: K1 p. d+ _2 l# B6 `3 o& I) l        0                  0             0.6109 - 0.1483i        0            -0.4065 - 0.0361i
2 W* H' t* p! Z7 I9 s        0                  0             0.5614 + 0.0342i   0.4157 + 0.0642i        0         
4 U9 f+ x/ e" y        0                  0            -0.2300 + 0.0825i        0            -0.7120 + 0.0395i
; j7 T6 o/ k; g        0            -0.0000 + 0.1972i        0                  0                  0         
        0                  0                  0            -0.4441 + 0.2775i        0         
        0                  0                  0             0.7522 + 0.0930i   1.5728 - 0.1753i

/ }+ ?: D* N4 }) \! h  Columns 11 through 14
. f5 I- g5 u! t$ a7 d" e0 _
) @/ N# N% w$ Z# w* U        0                  0                  0                  0         
4 _# C8 P6 ^! ]. C; E' n        0                  0                  0                  0         
1 N2 T$ O! H$ G4 j- B& j$ k        0                  0            -0.0746 - 0.0434i        0         
        0                  0            -0.0776 - 0.0247i        0         
        0                  0            -0.1764 - 0.0726i        0         
        0                  0                  0                  0         
        0             0.0000 - 0.1914i        0                  0         
   0.2339 - 0.0725i        0                  0                  0         
/ a* S3 [1 v/ L* k/ w# h; _        0                  0             0.4441 - 0.2212i  -0.7245 + 0.0210i
3 P8 W4 @  [) n3 E   0.7354 + 0.0592i        0                  0            -1.5296 + 0.3072i
5 b# ]7 A, d  }, t3 c        0                  0                  0                  0         
        0                  0                  0                  0         
8 U! K7 x$ w- s* T2 p2 W. P3 R) Z5 j# Y        0                  0                  0                  0         
        0                  0                  0                  0         
+ M/ Z% H) \1 I" \& q
, A! y6 n# X% s
& j  V) |- e0 ]" q' gPLoss =
2 v; O; ]+ B8 w/ [4 N
( v* J2 V  Y+ ?) e    0.1349

' u* L( e7 Q* b0 X' j3 g" I
! H5 h6 V! v1 u9 L$ _: j. LQg =
4 r/ g) x: S- Z" y
' V4 c" W3 o3 i8 S9 a; l$ N  Columns 1 through 11
' K+ `' m2 e& Y$ X7 O3 y
         0         0         0         0         0         0         0         0         0    0.5919    0.3119

  Columns 12 through 13
( D$ [* O/ @. H; W; d! N- R
    0.1972    0.5004
" \1 u; _1 D) h; a+ e3 g  b

Pg =

0 w6 E7 _: {0 h/ n8 @  Columns 1 through 11
8 I1 [1 g3 U: y# y( W
# [& u& ]2 u0 ]- g+ u         0         0         0         0         0         0         0         0         0    0.4000   -0.0000
# b# i- _! x7 c- Z9 K
  Columns 12 through 13
( C& _1 E6 h% G3 ?; U
! N7 C* G( z- h  F( @   -0.0000   -0.0000

& d" B/ m& n8 L$ Z0 r( m6 J: e  d
8 u: \. |' S+ C6 \% Ians =

1 i, K7 f% B* U8 O+ c! C    0.1349
0 e- W& G: w* s! ?' R' ^/ m! R0 R
7 b7 b3 `; `1 B6 {
+ I5 e! Q! A; z& R8 Y* D- C: gans =

    0.8414

hhusunguoqiang

matlab 角度计算单位都为弧度，是否换算一下就和真值差不多了！

s31xili

检查每一步计算以及第一次迭带的结果。正常的话3次就可以收敛了精度0.00000001。计算结果一概与IEEE14的参考结果在0.01以内。

wanru

matlab里面有个工具箱matpower是专门用于电力系统潮流和最优潮流计算的，里面都有现成的潮流程序，像你这种IEEE14节点的潮流计算很简单的，其实都不用自己编程，你把数据结构打进去就是了,或者你参考下他的程序把。matpower这个论坛里好像有人发过的，或者你跟我联系，我发个到你邮箱里，我的邮箱是wangxiaowan3280@sina.com

yaoyaohsm

呵呵  看着好复杂啊

chenkang

看着好复杂啊，不过我也将要编