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

迪拉克之海

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

* k3 Y* \0 N. O9 U* h) c& o[问题已经解决]

' N: y; I6 k8 U! |1 @[ 本帖最后由 迪拉克之海 于 2008-5-7 13:37 编辑 ]

xiao1806

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

flyingclouds

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

lgs_0003

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

yingfengli

u =
4 J1 E- [" p( j* ?2 @. Z2 q$ ^
# C! }  }" B: |% R8 g% q  Columns 1 through 11

    1.0516    1.0474    1.0550    1.0549    1.0497    1.0327    1.0581    1.0123    1.0138    1.0450    1.0100
4 y: R" Y+ a. I) N9 m
8 h+ R4 p+ T* |7 n: s) H  Columns 12 through 14
& H4 F" _* p( e* N+ v& ^2 A1 @
    1.0900    1.0700    1.0600

9 l* m7 d  J' H: T' D& `9 N% w
' `) b+ A$ Q, |% L' {Theta =

" T/ U& F% m* J9 G! ]/ E  Columns 1 through 11

0 W% v& W" b" P2 z   -0.2599   -0.2628   -0.2578   -0.2632   -0.2645   -0.2794   -0.2323   -0.1790   -0.1520   -0.0872   -0.2229

: X! ]$ v$ v( |: ?( A" T# B8 D  Columns 12 through 14
8 t3 w' N4 D6 ~
+ C+ ^, @/ n- ^5 U' ~. G$ l   -0.2323   -0.2483         0
3 J+ v* W, A5 E

8 N% U% N0 y0 _. e( B7 r2 RSslack =
( `) H# [, j) k1 V
   2.3249 - 0.0250i


! [- |% d4 _$ ^S =
  e+ ]) S' _) B
2 U5 C) x+ L; [/ }- b/ c+ O, c  Columns 1 through 5
+ d# C5 L" t8 X0 \8 |8 B; e  k
( L- ]& N, n1 ^4 c# B6 C: I        0             0.0507 + 0.0332i        0                  0                  0         
  -0.0506 - 0.0329i        0            -0.0394 - 0.0251i        0                  0         
        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         
   0.1595 - 0.0673i        0                  0                  0                  0         
        0                  0                  0                  0                  0         
/ [+ ^' n; n5 Z% P4 V. a        0                  0                  0                  0                  0         
5 F" h4 a7 h6 t& e( |2 \% C7 j        0                  0                  0                  0                  0         
        0                  0                  0                  0                  0         
2 B6 Z  _: u4 O, u/ u4 C' d- ^% N3 g& s6 c        0                  0             0.0752 + 0.0448i   0.0784 + 0.0262i   0.1785 + 0.0769i
        0                  0                  0                  0                  0         

  Columns 6 through 10

   0.0928 + 0.0303i  -0.2789 - 0.0591i  -0.1595 + 0.0831i        0                  0         
7 p) u. h. m2 S$ H4 s+ S! b        0                  0                  0                  0                  0         
        0                  0                  0                  0                  0         
        0                  0                  0                  0                  0         
& l, R" N% `8 s   0.0579 + 0.0232i        0                  0                  0                  0         
        0                  0                  0                  0                  0         
        0                  0            -0.2789 + 0.2447i        0                  0         
2 \  j/ B3 ]% K: b, [' L: v) t. p        0             0.2789 - 0.2196i        0            -0.6058 + 0.1645i  -0.5446 + 0.0169i
        0                  0             0.6109 - 0.1483i        0            -0.4065 - 0.0361i
        0                  0             0.5614 + 0.0342i   0.4157 + 0.0642i        0         
0 U, v) e' @% T3 N5 U        0                  0            -0.2300 + 0.0825i        0            -0.7120 + 0.0395i
8 b: m& H1 p( y        0            -0.0000 + 0.1972i        0                  0                  0         
        0                  0                  0            -0.4441 + 0.2775i        0         
& n$ R: c3 M, _( i& T        0                  0                  0             0.7522 + 0.0930i   1.5728 - 0.1753i
. N1 x# i" X! E. w' a/ b
  Columns 11 through 14
% G  ~( v) t# M) B$ M
) N4 K. Q4 y+ h" S# j$ Q, g        0                  0                  0                  0         
        0                  0                  0                  0         
  d& j, c( k1 L! V/ Z1 m1 h+ Z( ^        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         
% x. T7 h5 U6 m, j; q        0             0.0000 - 0.1914i        0                  0         
1 R! B3 K+ w/ N, t& k7 _   0.2339 - 0.0725i        0                  0                  0         
        0                  0             0.4441 - 0.2212i  -0.7245 + 0.0210i
   0.7354 + 0.0592i        0                  0            -1.5296 + 0.3072i
7 O8 B) `8 s7 T6 P; l8 E        0                  0                  0                  0         
        0                  0                  0                  0         
1 L' ]! \/ |2 Y$ R  f        0                  0                  0                  0         
5 x/ O2 q; w4 @& R: w6 Q! Q        0                  0                  0                  0         


PLoss =
& r6 j6 w* s7 j
    0.1349


Qg =
/ K, |% k$ [5 }5 G3 i
  Columns 1 through 11

         0         0         0         0         0         0         0         0         0    0.5919    0.3119

  Columns 12 through 13
) E1 n. w+ C9 C3 x% J6 E1 J* E
4 z, v# t& F  }! n8 }    0.1972    0.5004
: H1 e4 t+ D: m, l; b
' ]8 E/ }, ^% ~' X5 S" i3 B
' T" v, d8 h0 z( X* iPg =

  Columns 1 through 11
2 c7 W" I! z- {/ m0 x
         0         0         0         0         0         0         0         0         0    0.4000   -0.0000

  Columns 12 through 13
: D" y  }' ~& |& n
   -0.0000   -0.0000
# S$ J/ L1 b6 w- z" L; O
0 t* q% Y$ `8 H3 {$ \$ A* _! B) D
ans =

    0.1349

  t$ h. W7 Y8 [7 s2 l) Z- G
ans =

    0.8414

hhusunguoqiang

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

s31xili

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

wanru

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

yaoyaohsm

呵呵  看着好复杂啊

chenkang

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