设为首页收藏本站|繁體中文 快速切换版块

 找回密码
 立即加入
搜索
查看: 7629|回复: 20

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

 火 [复制链接]

该用户从未签到

尚未签到

发表于 2008-4-29 22:27:26 | 显示全部楼层 |阅读模式

马上加入,结交更多好友,共享更多资料,让你轻松玩转电力研学社区!

您需要 登录 才可以下载或查看,没有账号?立即加入

×
我参考了一些书籍,用for循环写了个基于极坐标牛顿-拉夫逊法的IEEE14节点潮流计算代码,执行后发现迭代的最终结果与IEEE14给出的结果出入较大。。。* E& n7 n2 P; f
这是我的毕业设计的一部分,在计算结果不能够统一的问题上我卡了一星期了,非常着急,自己盯着自己的代码看又死活看不出来什么地方出了问题,因此发到这里来,想借助各位高手的力量、帮我看看我到底什么地方出错了,非常感谢!+ N# n6 V; O7 V; B( u; d% `/ ?+ s

3 O- A( r0 @9 v  E7 \6 X[问题已经解决]$ O  a- v0 ?; ?

% q' P# _( w4 R$ \) A! X( n[ 本帖最后由 迪拉克之海 于 2008-5-7 13:37 编辑 ]
"真诚赞赏,手留余香"
还没有人打赏,支持一下
楼主热帖
帖文化:【文明发帖 和谐互动】 社区精神:【创新、交流、互助、共享】

该用户从未签到

尚未签到

发表于 2008-4-30 05:10:19 | 显示全部楼层
我也曾经用matlab做过用极坐标的牛拉法算潮流,遇到跟你一样的问题,建议你通过手算潮流,看你的第一次迭代是否正确,其实第一次迭代手算并不麻烦,因为很多值都是0,还有就是看你的程序是否是收敛在六次以内的。手算是调试程序最好的方法,不要怕麻烦嘛。
"真诚赞赏,手留余香"
还没有人打赏,支持一下
帖文化:【文明发帖 和谐互动】 社区精神:【创新、交流、互助、共享】

该用户从未签到

尚未签到

发表于 2008-4-30 10:11:44 | 显示全部楼层
仔细对照电力系统的书看一遍,看看是不是有公式输错,并联补偿支路没有考虑的情况,不过最后的结果确实与IEEE会有一些差别
"真诚赞赏,手留余香"
还没有人打赏,支持一下
帖文化:【文明发帖 和谐互动】 社区精神:【创新、交流、互助、共享】

该用户从未签到

尚未签到

发表于 2008-5-5 14:52:08 | 显示全部楼层
这个问题我也遇到了,迭代结果是无法计算得了,我也正有想法用手工算,不过现在手工算的话,是有点麻烦,所以我在准备想写一个只迭代一次,而且把每时刻计算的结果输出来对比一下咯
"真诚赞赏,手留余香"
还没有人打赏,支持一下
帖文化:【文明发帖 和谐互动】 社区精神:【创新、交流、互助、共享】
  • TA的每日心情
    开心
    2016-2-29 11:47
  • 签到天数: 1 天

    连续签到: 1 天

    [LV.1]初来乍到

    累计签到:2 天
    连续签到:1 天
    发表于 2008-5-7 11:24:14 | 显示全部楼层
    u =
    ! `/ w; j7 p. l. i
    3 l& k% V3 t' {; K) j# z9 S8 W$ r' E  Columns 1 through 11
    $ S2 j+ V9 v' f5 D6 Y- u# D( `; y4 z4 S% S. g6 }
        1.0516    1.0474    1.0550    1.0549    1.0497    1.0327    1.0581    1.0123    1.0138    1.0450    1.0100
      L3 A3 Q7 t5 Z0 ^& Q* C" v9 k6 x4 V  u  p! R" O' y3 {4 x, g7 A
      Columns 12 through 148 ]* Z8 Y7 J. ]& k" h% `, [

    3 I3 A7 @. \" D1 i# c    1.0900    1.0700    1.0600
    ( n% H7 G2 B' c4 u, c! h1 U' Y# ~* C. o7 Z

    - m. T6 t+ j1 k/ R5 V8 `  F- F2 lTheta =
    * f, w* X* z* `2 s. R
    6 y3 v4 E4 D' z1 l  Columns 1 through 11+ Z7 A, G8 J- B8 x* G2 n9 ]

    4 ?) g/ x7 v# O7 S   -0.2599   -0.2628   -0.2578   -0.2632   -0.2645   -0.2794   -0.2323   -0.1790   -0.1520   -0.0872   -0.2229
    ( L6 Q9 [! |. l% K
    ) `. }& w  b$ N) I+ G4 f! Y  Columns 12 through 14% i/ M/ v% F; Q& U- ^! K5 ]
    : T# T9 ~9 @; {" G
       -0.2323   -0.2483         0
    2 q( J7 S- w8 e: l- }  N' {  B: ~6 |; R9 L1 {$ i

    4 i3 Y2 c' o! cSslack =3 u: c0 @; R3 p1 D; |9 J6 R4 ?

    4 y. p. j0 n" [0 I; w, G7 I! S   2.3249 - 0.0250i
    8 ^; D# h9 R  H' B1 |) l9 a- u& w& P& p
    0 Y  d, f9 Y( K1 |5 {3 N' L
    S =" o# L; H' t8 _/ A, K5 L" v
    - \3 M& c6 X; Z* _5 \9 K
      Columns 1 through 5
    % V4 H8 @1 a9 N, Y$ w) J9 n; a; b0 C" V, |. x4 w  g8 T
            0             0.0507 + 0.0332i        0                  0                  0         
    4 s& @2 {% L3 b2 T9 l  -0.0506 - 0.0329i        0            -0.0394 - 0.0251i        0                  0          * V5 G2 I! R1 E# _7 p* J* T6 S* p! J
            0             0.0396 + 0.0254i        0                  0                  0          , m; V4 Z' n, p% ?. a
            0                  0                  0                  0             0.0166 + 0.0087i
    , p- S9 w8 |8 Q1 \9 o/ ~- ^( A        0                  0                  0            -0.0165 - 0.0086i        0          # k1 B, M; {& j, R( W
      -0.0917 - 0.0280i        0                  0                  0            -0.0573 - 0.0220i
    ! Q8 x: ^; i, j3 |: E4 X5 x" _, \   0.2789 + 0.0672i        0                  0                  0                  0          ( a1 R- N" @2 p+ x0 B  ?
       0.1595 - 0.0673i        0                  0                  0                  0          $ u3 Q1 k1 V* l& u' h+ E
            0                  0                  0                  0                  0         
    6 B- T; E- H  a3 E' v        0                  0                  0                  0                  0         
    : k" w5 A5 d2 F6 M        0                  0                  0                  0                  0         
    2 W7 _; \' R% d        0                  0                  0                  0                  0         
    9 G4 p% j7 g$ o8 B# r4 Y7 }, o        0                  0             0.0752 + 0.0448i   0.0784 + 0.0262i   0.1785 + 0.0769i# u) i/ \. P( V" a
            0                  0                  0                  0                  0          # G8 N1 X+ f' x9 {: i( u; w! C  s
    ! A( u- v% G5 `2 ~+ g3 m
      Columns 6 through 10
    + r6 P, E+ Y: C& f" ~" ]; E3 N: \; X
       0.0928 + 0.0303i  -0.2789 - 0.0591i  -0.1595 + 0.0831i        0                  0         
    % [* G( F8 M: f- |  i6 X        0                  0                  0                  0                  0         
    - k: b2 ]  v4 h        0                  0                  0                  0                  0          ) I3 V' s0 R6 ~, _1 \
            0                  0                  0                  0                  0          4 U' d- Q: e& g) M% C7 o) a5 E
       0.0579 + 0.0232i        0                  0                  0                  0          ! P8 ~4 `) O2 t% p! X
            0                  0                  0                  0                  0         
    1 q. u+ ^; j: P+ R5 u        0                  0            -0.2789 + 0.2447i        0                  0          0 u/ t( u3 m" n- X$ Q, ^% _8 i" j
            0             0.2789 - 0.2196i        0            -0.6058 + 0.1645i  -0.5446 + 0.0169i
    1 R8 w9 D" L* p4 y& P8 Z        0                  0             0.6109 - 0.1483i        0            -0.4065 - 0.0361i
    ' U2 C1 m) H! e" A        0                  0             0.5614 + 0.0342i   0.4157 + 0.0642i        0          0 g& R8 K6 D: y) ^8 x! J" @
            0                  0            -0.2300 + 0.0825i        0            -0.7120 + 0.0395i$ Y2 D& o5 O) m$ r# P
            0            -0.0000 + 0.1972i        0                  0                  0          , b$ F8 j$ \0 }* ~+ l4 L+ Z
            0                  0                  0            -0.4441 + 0.2775i        0          ( e) s. m$ _3 q! e
            0                  0                  0             0.7522 + 0.0930i   1.5728 - 0.1753i
    8 r' p" O3 [* y6 U
    . s$ J- L0 o& [9 B' |7 s& j  Columns 11 through 14
    / y3 ]: m, G# j8 P5 c9 m# ]6 t/ K3 _# V1 w  I. }6 y! U
            0                  0                  0                  0         
    0 i! G8 n9 S$ O1 t( z% b5 f* O        0                  0                  0                  0          4 F3 a+ h- [$ o
            0                  0            -0.0746 - 0.0434i        0          9 k1 A" N) P2 w& u% ]
            0                  0            -0.0776 - 0.0247i        0          ( }! u% J+ h/ [0 {
            0                  0            -0.1764 - 0.0726i        0          & N5 C1 h$ @" Y1 {  B
            0                  0                  0                  0          6 i$ N$ |$ \8 l7 s( ?, e
            0             0.0000 - 0.1914i        0                  0          5 j% z; ]! L6 @3 e( d$ v
       0.2339 - 0.0725i        0                  0                  0         
    8 a, L* a7 k1 [        0                  0             0.4441 - 0.2212i  -0.7245 + 0.0210i
    3 Q4 [' j+ W# U, G3 ^, l   0.7354 + 0.0592i        0                  0            -1.5296 + 0.3072i
      {/ g  m* N' S0 T4 ~! l# y- ?        0                  0                  0                  0         
    1 m0 ]' M: T* n' T7 Y; v        0                  0                  0                  0          . P( s) S7 M7 b" @% s5 M
            0                  0                  0                  0         
    9 \. P( T. O+ {        0                  0                  0                  0          + O% T: B% f, s" j

    * s" S+ U0 ~* Y: p7 T  O' G& _- c) o. T3 [0 ]$ I
    PLoss =
    ) C# d$ n' u- `0 U$ l2 @4 V) T
    % R! v2 i& ?- p" z3 ~    0.1349& b$ H" d) Y5 T+ x* ]& o1 K  r
    . f4 \5 c3 f, j- M! L

    ; r5 p0 H5 a: \Qg =! L( i* N3 K' X! K( R$ D
    8 R2 a, l8 y% V
      Columns 1 through 11
    / V# s4 l: A5 c( ?5 L2 C
    9 e, O5 G1 z( G         0         0         0         0         0         0         0         0         0    0.5919    0.3119
    / k& U# U! W0 T  N& w3 H' s! j( ]8 c2 _$ A& C
      Columns 12 through 13
    & Y- m( f" @/ ~0 z) W# p  g+ i' t8 o8 E8 b: ~2 P' Y7 ^! N' |- R
        0.1972    0.5004
    , O; z& S5 c. o4 Z% y0 K; o4 \' D  j/ s# U4 d- y# J7 K' x

    ! j: W" v5 p. g2 l: g& V. dPg =+ z6 H4 I" i" T7 p/ T( B

      {1 X2 V- q- n0 M' e  Columns 1 through 11- z. r: H: a( ~) Y- h' ?

    $ o5 [6 {- K+ u- E, r, ~( A         0         0         0         0         0         0         0         0         0    0.4000   -0.00009 y* w  z" f; y0 M5 g% t+ I

    7 F* c9 U- T! m1 p7 o. n  Columns 12 through 137 e1 ]: p* [6 m, g- A( Q! s0 J$ s

    ! ~+ J$ c6 F) W' w# x" X   -0.0000   -0.0000
    9 Z9 R, r! t0 p5 W1 L5 I" r+ k( f; w! g1 c, C+ d/ T, p
    # U7 W& w/ y# e5 U5 R: N$ R1 G
    ans =
    " |. s; v/ r% v" [' J0 j) Z9 }. [) B( E8 u5 d- f
        0.13496 U' G, @" X/ s% n9 J7 p$ r
    ( g  ?# G4 b' ?% z4 G# h

    % f* ?5 w* P$ L" ~( D! b* Rans =
    # z3 F9 J8 a+ k* E7 F7 ]. h$ x' r
    * A! R# V0 Z# O    0.8414
    "真诚赞赏,手留余香"
    还没有人打赏,支持一下
    帖文化:【文明发帖 和谐互动】 社区精神:【创新、交流、互助、共享】
  • TA的每日心情
    郁闷
    2018-7-12 14:08
  • 签到天数: 1 天

    连续签到: 1 天

    [LV.1]初来乍到

    累计签到:1 天
    连续签到:1 天
    发表于 2008-6-18 11:01:15 | 显示全部楼层
    matlab 角度计算单位都为弧度,是否换算一下就和真值差不多了!
    "真诚赞赏,手留余香"
    还没有人打赏,支持一下
    帖文化:【文明发帖 和谐互动】 社区精神:【创新、交流、互助、共享】
  • TA的每日心情
    郁闷
    2018-1-16 18:34
  • 签到天数: 15 天

    连续签到: 1 天

    [LV.4]偶尔看看III

    累计签到:15 天
    连续签到:1 天
    发表于 2008-6-25 16:51:23 | 显示全部楼层
    检查每一步计算以及第一次迭带的结果。正常的话3次就可以收敛了精度0.00000001。计算结果一概与IEEE14的参考结果在0.01以内。
    "真诚赞赏,手留余香"
    还没有人打赏,支持一下
    帖文化:【文明发帖 和谐互动】 社区精神:【创新、交流、互助、共享】

    该用户从未签到

    尚未签到

    发表于 2008-6-25 21:06:42 | 显示全部楼层
    matlab里面有个工具箱matpower是专门用于电力系统潮流和最优潮流计算的,里面都有现成的潮流程序,像你这种IEEE14节点的潮流计算很简单的,其实都不用自己编程,你把数据结构打进去就是了,或者你参考下他的程序把。matpower这个论坛里好像有人发过的,或者你跟我联系,我发个到你邮箱里,我的邮箱是wangxiaowan3280@sina.com
    "真诚赞赏,手留余香"
    还没有人打赏,支持一下
    帖文化:【文明发帖 和谐互动】 社区精神:【创新、交流、互助、共享】

    该用户从未签到

    尚未签到

    发表于 2009-3-26 13:00:27 | 显示全部楼层
    呵呵  看着好复杂啊
    "真诚赞赏,手留余香"
    还没有人打赏,支持一下
    帖文化:【文明发帖 和谐互动】 社区精神:【创新、交流、互助、共享】

    该用户从未签到

    尚未签到

    发表于 2009-5-8 22:14:04 | 显示全部楼层
    看着好复杂啊,不过我也将要编
    "真诚赞赏,手留余香"
    还没有人打赏,支持一下
    帖文化:【文明发帖 和谐互动】 社区精神:【创新、交流、互助、共享】
    您需要登录后才可以回帖 登录 | 立即加入

    本版积分规则

    招聘斑竹

    小黑屋|手机版|APP下载(beta)|Archiver|电力研学网 ( 赣ICP备12000811号-1|赣公网安备36040302000210号 )|网站地图

    GMT+8, 2026-10-20 11:42

    Powered by Discuz! X3.5 Licensed

    © 2001-2026 Discuz! Team.

    快速回复 返回顶部 返回列表