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

迪拉克之海

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

[问题已经解决]
" y" Y/ q# E7 o5 {5 w6 W& A
1 m7 e0 g  A  f9 h$ S[ 本帖最后由 迪拉克之海 于 2008-5-7 13:37 编辑 ]

xiao1806

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

flyingclouds

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

lgs_0003

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

yingfengli

u =

! W$ P* M* U3 G: c( {! L  Columns 1 through 11

/ a1 I1 I5 l/ j" E    1.0516    1.0474    1.0550    1.0549    1.0497    1.0327    1.0581    1.0123    1.0138    1.0450    1.0100

! a, Y4 C& \3 F  f4 s* {7 L8 J  Columns 12 through 14

' w, i& N" H9 b& _    1.0900    1.0700    1.0600


/ }$ ~4 M, O$ z  b3 f6 Z# W5 pTheta =
( q: a8 H  f* {: `
  Columns 1 through 11
: r6 I. n( N) {: f
   -0.2599   -0.2628   -0.2578   -0.2632   -0.2645   -0.2794   -0.2323   -0.1790   -0.1520   -0.0872   -0.2229

" c, Q$ x" {) H2 {  Columns 12 through 14
: w5 U: z; q4 Y4 ~+ X
   -0.2323   -0.2483         0
9 W) u* G4 l- \3 ?( k
& g' ?) ~' k, t3 d' z- S9 S8 U1 T
Sslack =
4 Z& P$ j' S9 Q$ ~9 g
* Q3 j( n  M! ~/ Q( j/ ]   2.3249 - 0.0250i
1 ]- n, E4 e5 b. c/ G
! v5 ~3 v  U$ e. I5 ?
  n  t( _( c6 x( JS =
2 P9 g. @9 M6 c# c: ?
! m+ x. u) Q3 |: B. j% r. {  Columns 1 through 5

        0             0.0507 + 0.0332i        0                  0                  0         
8 o+ O. F& V) ^% v  -0.0506 - 0.0329i        0            -0.0394 - 0.0251i        0                  0         
        0             0.0396 + 0.0254i        0                  0                  0         
) H# j; e& G& l) b5 R5 M        0                  0                  0                  0             0.0166 + 0.0087i
+ ^7 C" H5 Y) R/ `        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         
, x% C. r, s' }        0                  0                  0                  0                  0         
        0                  0                  0                  0                  0         
  g7 n7 z0 ^. k% i% c" D' z        0                  0                  0                  0                  0         
        0                  0                  0                  0                  0         
        0                  0             0.0752 + 0.0448i   0.0784 + 0.0262i   0.1785 + 0.0769i
        0                  0                  0                  0                  0         

$ G- b7 ~' A0 S; `  Columns 6 through 10

+ t1 c8 W) W: V/ Z9 l- |# @6 W; z   0.0928 + 0.0303i  -0.2789 - 0.0591i  -0.1595 + 0.0831i        0                  0         
# M& o+ i8 C" A  |) t        0                  0                  0                  0                  0         
        0                  0                  0                  0                  0         
        0                  0                  0                  0                  0         
/ s/ V  |5 M3 @/ A  [   0.0579 + 0.0232i        0                  0                  0                  0         
) f: d$ m6 I7 d* o; X+ {        0                  0                  0                  0                  0         
& @$ Y& t) p) s9 J& x9 y! Q        0                  0            -0.2789 + 0.2447i        0                  0         
        0             0.2789 - 0.2196i        0            -0.6058 + 0.1645i  -0.5446 + 0.0169i
4 c) h- o- R2 i% `        0                  0             0.6109 - 0.1483i        0            -0.4065 - 0.0361i
2 g  e8 N3 R( r$ A        0                  0             0.5614 + 0.0342i   0.4157 + 0.0642i        0         
        0                  0            -0.2300 + 0.0825i        0            -0.7120 + 0.0395i
        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
# @2 @5 j5 F3 j- K% g5 W
  Columns 11 through 14
2 S6 Z' ^4 r9 n/ g# m6 @& }
" o6 `; p: _+ i        0                  0                  0                  0         
        0                  0                  0                  0         
        0                  0            -0.0746 - 0.0434i        0         
        0                  0            -0.0776 - 0.0247i        0         
: q6 ^7 {) p$ N; r        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         
        0                  0             0.4441 - 0.2212i  -0.7245 + 0.0210i
   0.7354 + 0.0592i        0                  0            -1.5296 + 0.3072i
        0                  0                  0                  0         
        0                  0                  0                  0         
) p6 E- d" Q# d, ]; U        0                  0                  0                  0         
$ g$ j% J2 D9 T& [# [( E; V        0                  0                  0                  0         
+ z8 ]: c5 U4 B% u5 O2 x
2 `7 y* P/ Y2 M" b
PLoss =

    0.1349
  r/ Z4 U& o  t% t9 D

Qg =

" N( P* v5 F& H1 B2 f. G  Columns 1 through 11

         0         0         0         0         0         0         0         0         0    0.5919    0.3119

" s" V+ Y# k, @  Columns 12 through 13

    0.1972    0.5004
& u! m0 R. W5 k9 f: m4 z+ l
4 F% W% x- Z6 X7 v3 t: h! V
Pg =
& |8 k. C3 |  m
  Columns 1 through 11

         0         0         0         0         0         0         0         0         0    0.4000   -0.0000
+ K) }! f; t* R" w) W* o$ D6 `
  Columns 12 through 13
) H: L+ ^. j* Z0 h( A) a: R% U, f- T
   -0.0000   -0.0000
1 H% \- D' u  j. s' Y3 {7 T

ans =

    0.1349
5 t$ k2 B: t: L
5 E# u4 J) M7 A; z  j% o
ans =

# f, S% S! b  A$ D" n* d, n( \    0.8414

hhusunguoqiang

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

s31xili

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

wanru

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

yaoyaohsm

呵呵  看着好复杂啊

chenkang

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