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

 找回密码
 立即加入
搜索

电力研学论坛»交流社区 › 电气工程研究与技术特色版区 › 电力系统与电网技术 › 请教一个关于Matlab编写的IEEE14节点潮流计算的问题 ...
123下一页
返回列表 发新帖
查看: 7631|回复: 20

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

 火 [复制链接]
迪拉克之海

该用户从未签到

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

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

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

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

$ }* v% l9 d7 C2 m+ P[问题已经解决]* `/ Q% [) n' P" Q; \
: G8 N. t# H% [$ ?8 o
[ 本帖最后由 迪拉克之海 于 2008-5-7 13:37 编辑 ]
潮流计算的节点, ieee33节点潮流计算, 潮流计算平衡节点, ieee57节点潮流计算
"真诚赞赏，手留余香"
打赏
还没有人打赏，支持一下
楼主热帖

相关帖子

帖文化：【文明发帖 和谐互动】 社区精神：【创新、交流、互助、共享】
参与“宣传中心”推广活动！让更多的人参与进来，让更多的人可以帮助您！+学分+金币
回复

使用道具 举报

xiao1806

该用户从未签到

尚未签到
发表于 2008-4-30 05:10:19 | 显示全部楼层
我也曾经用matlab做过用极坐标的牛拉法算潮流，遇到跟你一样的问题，建议你通过手算潮流，看你的第一次迭代是否正确，其实第一次迭代手算并不麻烦，因为很多值都是0，还有就是看你的程序是否是收敛在六次以内的。手算是调试程序最好的方法，不要怕麻烦嘛。
"真诚赞赏，手留余香"
打赏
还没有人打赏，支持一下
帖文化：【文明发帖 和谐互动】 社区精神：【创新、交流、互助、共享】
参与“宣传中心”推广活动！让更多的人参与进来，让更多的人可以帮助您！+学分+金币
回复 推荐 踩下

使用道具 举报

flyingclouds

该用户从未签到

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

使用道具 举报

lgs_0003

该用户从未签到

尚未签到
发表于 2008-5-5 14:52:08 | 显示全部楼层
这个问题我也遇到了，迭代结果是无法计算得了，我也正有想法用手工算，不过现在手工算的话，是有点麻烦，所以我在准备想写一个只迭代一次，而且把每时刻计算的结果输出来对比一下咯
"真诚赞赏，手留余香"
打赏
还没有人打赏，支持一下
帖文化：【文明发帖 和谐互动】 社区精神：【创新、交流、互助、共享】
参与“宣传中心”推广活动！让更多的人参与进来，让更多的人可以帮助您！+学分+金币
回复 推荐 踩下

使用道具 举报

yingfengli
  • TA的每日心情
    开心
    2016-2-29 11:47

    • 签到天数: 1 天

    连续签到: 1 天

    [LV.1]初来乍到

    累计签到：2 天
    连续签到：1 天
    发表于 2008-5-7 11:24:14 | 显示全部楼层
    u =! K. n- Z/ A8 l. O( M! G) ]8 e
    7 a" ^; Q+ j* o: ?
      Columns 1 through 11
    3 F6 [, k$ O) O  O& w( S+ z- T  O. ]! |9 O
        1.0516    1.0474    1.0550    1.0549    1.0497    1.0327    1.0581    1.0123    1.0138    1.0450    1.01001 C. m6 w0 {) G: Z

    : v, q, q5 s, `0 R  c5 a  Columns 12 through 14, \7 {2 P& l9 E
    / r3 F  k& c: E. l& L
        1.0900    1.0700    1.0600" [9 V$ v$ x2 f! u: H
    4 p4 i: z8 a/ Y( i! n

    ; G8 n- F5 Y# n: V% L! N8 v( o2 lTheta =
    . d5 T# }2 H! H$ r
    9 b& J. H( P- S8 ~  Columns 1 through 11
    1 R; n, u& p, V: v, I( K& u1 C5 W$ K6 \
       -0.2599   -0.2628   -0.2578   -0.2632   -0.2645   -0.2794   -0.2323   -0.1790   -0.1520   -0.0872   -0.22293 {- a5 c+ B) M

    8 f& i7 A  y6 P" y0 a4 f  Columns 12 through 148 k0 o$ L" M* o& O* L' C8 T* D2 `
    9 o" Z% I1 l" s, \9 l9 m
       -0.2323   -0.2483         03 ?0 K: j9 E. x" h$ t: T/ U6 h

    / W' j% y- r# Z6 J: Q; E7 Q6 m1 `5 x$ f# l
    Sslack =
    % z0 ^/ Q3 n% q9 o
    & n: z) ~! W9 T& Y   2.3249 - 0.0250i  o# w9 B) q. I/ V
    8 P# ~6 H. ~# t% y5 ^3 L9 U3 }; Q0 I

    ' b9 N# o, O8 }5 N" Y* sS =  F2 x0 E6 A, }

    , u& Q3 i! s7 {" z% [8 c+ G2 f" g  Columns 1 through 5! h1 z  Q) }# d

    / h5 M& r8 w, y( P4 {9 i: W: c        0             0.0507 + 0.0332i        0                  0                  0          ( C# S( X# X, @/ B. }9 P
      -0.0506 - 0.0329i        0            -0.0394 - 0.0251i        0                  0         
    - `3 z1 u4 x- e% [; v2 `" n, A        0             0.0396 + 0.0254i        0                  0                  0          & a) O% ~0 x+ W8 |
            0                  0                  0                  0             0.0166 + 0.0087i
    : a; _' i; g( ~3 m& P1 r0 p        0                  0                  0            -0.0165 - 0.0086i        0          + C7 u7 \) J9 _7 Z- L/ e+ \
      -0.0917 - 0.0280i        0                  0                  0            -0.0573 - 0.0220i
    2 T' @, [. v) C/ X4 m# k% e   0.2789 + 0.0672i        0                  0                  0                  0            D5 v3 K7 r% R$ D
       0.1595 - 0.0673i        0                  0                  0                  0          * E3 B/ u+ ?6 \$ G( A* c' C2 \
            0                  0                  0                  0                  0          ' v$ {0 c8 m! g
            0                  0                  0                  0                  0          " j7 |' G: E/ `; A! H$ a
            0                  0                  0                  0                  0          1 L2 Z# d+ c+ L) Q. S- P+ j" M
            0                  0                  0                  0                  0         
    ( V6 t; G: K- J0 f& _$ y3 M8 D; \        0                  0             0.0752 + 0.0448i   0.0784 + 0.0262i   0.1785 + 0.0769i
    " ~  l1 M" r$ |( g. R, |        0                  0                  0                  0                  0          5 q; L: D3 L/ t- P. t8 J
    % r) @8 B& @5 X0 y0 ?
      Columns 6 through 10* V2 d5 {' b; |0 ^( ?. U- Y
    - @/ g, X/ f* L5 l- F" a
       0.0928 + 0.0303i  -0.2789 - 0.0591i  -0.1595 + 0.0831i        0                  0          $ x+ V, i, N; {( _
            0                  0                  0                  0                  0          / s1 l1 ]  K# V& X1 K4 e' g  ?
            0                  0                  0                  0                  0          , r0 Q* C2 s- c, @- U  Q
            0                  0                  0                  0                  0          + s5 D( q0 b4 c( }- g- n3 O$ ?& i
       0.0579 + 0.0232i        0                  0                  0                  0         
    7 C+ ~; i' y/ y6 K1 \! a2 X5 G( Z+ O        0                  0                  0                  0                  0          $ H2 v3 W7 G) ^9 U6 b
            0                  0            -0.2789 + 0.2447i        0                  0          / g7 j) T1 s2 x1 N! J0 l
            0             0.2789 - 0.2196i        0            -0.6058 + 0.1645i  -0.5446 + 0.0169i
    # r. Z6 B8 D: v5 }$ F: S9 i9 D        0                  0             0.6109 - 0.1483i        0            -0.4065 - 0.0361i
    % m! L+ ]* }9 j! B        0                  0             0.5614 + 0.0342i   0.4157 + 0.0642i        0         
    " ~+ Q* D8 g0 S8 {7 u2 R        0                  0            -0.2300 + 0.0825i        0            -0.7120 + 0.0395i
    ! e# c% W; m! @4 R. H        0            -0.0000 + 0.1972i        0                  0                  0          1 m2 I( T% P. d+ z' v2 M$ i
            0                  0                  0            -0.4441 + 0.2775i        0          - l* \0 n1 p% F$ d
            0                  0                  0             0.7522 + 0.0930i   1.5728 - 0.1753i" O' k/ ]+ u# a* w' U" z, j3 L
    ; i6 w& W1 K4 d+ y5 e/ e
      Columns 11 through 14
    7 h- ~& _7 D# m- P
    . R/ n! x% t! v) ?3 Y2 N# {. r        0                  0                  0                  0          , A  `- k% _2 B8 }
            0                  0                  0                  0          ) R+ Z6 d: S; s9 {
            0                  0            -0.0746 - 0.0434i        0         
    , u* ]+ M2 C1 h8 L8 \" v        0                  0            -0.0776 - 0.0247i        0          1 @! _2 k; ^0 \- p- C
            0                  0            -0.1764 - 0.0726i        0          4 x* E; W0 B. Z' x
            0                  0                  0                  0          # X% i6 }2 J9 q  c
            0             0.0000 - 0.1914i        0                  0          : }* L/ f' S4 C
       0.2339 - 0.0725i        0                  0                  0         
    0 `4 I( t7 `0 w  u, b        0                  0             0.4441 - 0.2212i  -0.7245 + 0.0210i
    ' l0 }, g. w; v( {0 j9 \+ m   0.7354 + 0.0592i        0                  0            -1.5296 + 0.3072i
    1 d/ u/ f5 V0 m. ^        0                  0                  0                  0         
    9 ~; O* u8 K9 l2 @2 e        0                  0                  0                  0         
    / c3 ^) |0 `6 b1 _4 r        0                  0                  0                  0          / h" Z! q( I7 ?6 j! |
            0                  0                  0                  0          ' V& b3 F! {3 F
    ' N% y' S* k! H, ~7 e5 \( R9 q5 e4 s

    5 D# u+ Z7 c+ V' s0 L. g" mPLoss =8 _- W: q% O3 I# d5 d. E
    " }) n5 Q5 y. M4 e5 ^9 ~% J/ ^
        0.1349$ }, r1 M$ J- ?( v2 t4 X
    . [( D4 r- K! N9 t  A& ~
    : O0 q2 R2 h8 X. l4 I5 W5 G2 e
    Qg =
    * Z& c# ~" _4 A6 j& ^2 X! R# l0 _# z* y8 K
      Columns 1 through 11
    1 ?) M9 _) ^8 e: _/ }5 f5 {% F6 U0 Y3 f* E
             0         0         0         0         0         0         0         0         0    0.5919    0.3119' n, z# k1 |" I! c- T
    ' K: x$ V$ J# M1 W- z
      Columns 12 through 13# G! c* ]0 i5 T2 V' A( D

    6 l$ y  L& |- O5 j& F' \    0.1972    0.50049 S; Y+ {3 p( g# G+ k* ]
    1 H* F4 m1 e" A5 [/ e$ p  E$ z

    8 q$ c+ d& [; ~1 F$ g/ PPg =
    ! N! J# {0 _; P! x/ C8 K$ u" F# T5 H) o* e7 T2 u; g
      Columns 1 through 117 r" F. k. f1 v0 i- C

    ( x! l8 I& U  P1 o- v         0         0         0         0         0         0         0         0         0    0.4000   -0.0000
    0 ]* Z& W+ M1 @6 s) m! V; v+ ]. d8 d! N
      Columns 12 through 13
    # C$ l) B- w  p$ {' o: ~  Y7 r6 u! Y& P
       -0.0000   -0.0000$ ~* [( ?% O( C& g* f! t

    8 M6 n1 O  n+ U/ W* J. X1 _, c1 B6 M* T& Y1 h+ U2 @* k
    ans =
    + P8 H# C' _( \
    8 |$ V8 I: h# U: s    0.1349
    ) J2 U) S2 |( G4 z. G3 P$ [
    & j6 _9 z/ {! a8 S9 d; v% ~( {" Y! v% V* N( l( d9 `! E* j' r' i- X
    ans =% g+ r! W1 o8 m

    6 ?% {4 N6 K$ d- Z8 ~, |; K    0.8414
    "真诚赞赏，手留余香"
    打赏
    还没有人打赏，支持一下
    帖文化：【文明发帖 和谐互动】 社区精神：【创新、交流、互助、共享】
    参与“宣传中心”推广活动！让更多的人参与进来，让更多的人可以帮助您！+学分+金币
    回复 推荐 踩下

    使用道具 举报

    hhusunguoqiang
  • TA的每日心情
    郁闷
    2018-7-12 14:08

    • 签到天数: 1 天

    连续签到: 1 天

    [LV.1]初来乍到

    累计签到：1 天
    连续签到：1 天
    发表于 2008-6-18 11:01:15 | 显示全部楼层
    matlab 角度计算单位都为弧度，是否换算一下就和真值差不多了！
    "真诚赞赏，手留余香"
    打赏
    还没有人打赏，支持一下
    帖文化：【文明发帖 和谐互动】 社区精神：【创新、交流、互助、共享】
    小贴士：资料共享版块上传附件“专家分”值+5，“学分”+1
    回复 推荐 踩下

    使用道具 举报

    s31xili
  • 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以内。
    "真诚赞赏，手留余香"
    打赏
    还没有人打赏，支持一下
    帖文化：【文明发帖 和谐互动】 社区精神：【创新、交流、互助、共享】
    参与“宣传中心”推广活动！让更多的人参与进来，让更多的人可以帮助您！+学分+金币
    回复 推荐 踩下

    使用道具 举报

    wanru

    该用户从未签到

    尚未签到
    发表于 2008-6-25 21:06:42 | 显示全部楼层
    matlab里面有个工具箱matpower是专门用于电力系统潮流和最优潮流计算的，里面都有现成的潮流程序，像你这种IEEE14节点的潮流计算很简单的，其实都不用自己编程，你把数据结构打进去就是了,或者你参考下他的程序把。matpower这个论坛里好像有人发过的，或者你跟我联系，我发个到你邮箱里，我的邮箱是wangxiaowan3280@sina.com
    "真诚赞赏，手留余香"
    打赏
    还没有人打赏，支持一下
    帖文化：【文明发帖 和谐互动】 社区精神：【创新、交流、互助、共享】
    参与“宣传中心”推广活动！让更多的人参与进来，让更多的人可以帮助您！+学分+金币
    回复 推荐 踩下

    使用道具 举报

    yaoyaohsm

    该用户从未签到

    尚未签到
    发表于 2009-3-26 13:00:27 | 显示全部楼层
    呵呵  看着好复杂啊
    "真诚赞赏，手留余香"
    打赏
    还没有人打赏，支持一下
    帖文化：【文明发帖 和谐互动】 社区精神：【创新、交流、互助、共享】
    参与“宣传中心”推广活动！让更多的人参与进来，让更多的人可以帮助您！+学分+金币
    回复 推荐 踩下

    使用道具 举报

    chenkang

    该用户从未签到

    尚未签到
    发表于 2009-5-8 22:14:04 | 显示全部楼层
    看着好复杂啊，不过我也将要编
    "真诚赞赏，手留余香"
    打赏
    还没有人打赏，支持一下
    帖文化：【文明发帖 和谐互动】 社区精神：【创新、交流、互助、共享】
    小贴士：资料共享版块上传附件“专家分”值+5，“学分”+1
    回复 推荐 踩下

    使用道具 举报

    123下一页
    返回列表 发新帖
    高级模式
    B Color Image Link Quote Code Smilies
    您需要登录后才可以回帖 登录 | 立即加入

    本版积分规则

    招聘斑竹

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

    GMT+8, 2026-10-20 13:21

    Powered by Discuz! X3.5 Licensed

    © 2001-2026 Discuz! Team.

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