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

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

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

  [复制链接]

该用户从未签到

尚未签到

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

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

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

×
我参考了一些书籍,用for循环写了个基于极坐标牛顿-拉夫逊法的IEEE14节点潮流计算代码,执行后发现迭代的最终结果与IEEE14给出的结果出入较大。。。- X; a$ Z5 J: A& s0 Q' B6 d
这是我的毕业设计的一部分,在计算结果不能够统一的问题上我卡了一星期了,非常着急,自己盯着自己的代码看又死活看不出来什么地方出了问题,因此发到这里来,想借助各位高手的力量、帮我看看我到底什么地方出错了,非常感谢!
5 S; b$ t8 X+ \% a6 A2 |# @9 M- Q7 K9 G8 |% X  {6 O* u/ [
[问题已经解决]& W3 C. o+ I2 W( z

& `9 B/ K* o/ Y- y7 M[ 本帖最后由 迪拉克之海 于 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 =
    % m# \0 S- {( i1 j; j$ t
    6 \9 B/ x: I/ I- R0 t  Columns 1 through 11
      Z7 n) R( S+ |4 s/ H* ^/ f+ y9 c# j2 Y
        1.0516    1.0474    1.0550    1.0549    1.0497    1.0327    1.0581    1.0123    1.0138    1.0450    1.0100
    ! K' m$ v/ n$ b
    2 Q, x$ z) m) Q4 i: d2 ]% f1 N* d  Columns 12 through 14
    8 L$ z3 S$ J; a4 o- _  j% @$ ~! _" E: V8 ]7 E* r
        1.0900    1.0700    1.06003 ^* X! ]* q: W

    9 R6 A4 z7 E$ ^5 ]6 z# h  s
    # L$ f7 h+ O( V& }Theta =
    8 ?  [, e% L7 m# m4 ?$ `
    , X0 V2 c# \2 ~' q# N  Columns 1 through 11
    , f  R) i  p) i2 h+ A* _/ ^. ^/ C9 {
       -0.2599   -0.2628   -0.2578   -0.2632   -0.2645   -0.2794   -0.2323   -0.1790   -0.1520   -0.0872   -0.2229+ g. ^  X/ B$ X4 x

    & m% g- @7 i; U  i3 f  Columns 12 through 141 j% g: d- i$ f% G$ c7 m

    6 g( m7 r* o- |; Q5 k   -0.2323   -0.2483         0
    + R8 h' J. N) q, O
    ; ~6 ?1 Q" s0 |5 P3 u9 f0 b+ J! Q( ^& H  Y, K9 A
    Sslack =
    8 v8 E3 L" B* K
    $ y. g8 o) G! z. r   2.3249 - 0.0250i( {2 O+ Z$ E" x
    5 j2 X* H/ |4 A& D# {1 o
    - Y+ F3 G+ j# ?( O
    S =
    & ~9 k% S' N+ S# I, }3 ]9 E  Z2 r3 r4 U* \6 r5 m
      Columns 1 through 5+ X4 d. i4 x9 e  \. o
    0 o( J" D$ _  o% s" u
            0             0.0507 + 0.0332i        0                  0                  0          5 U" d# u: `& }$ B
      -0.0506 - 0.0329i        0            -0.0394 - 0.0251i        0                  0          4 E' X. x9 x& I  k! W
            0             0.0396 + 0.0254i        0                  0                  0         
    8 [4 h/ \( p$ z6 |        0                  0                  0                  0             0.0166 + 0.0087i/ E; E" u& }8 \- z# R; V* o. [
            0                  0                  0            -0.0165 - 0.0086i        0          % n( W6 e9 H( n5 K1 [5 {
      -0.0917 - 0.0280i        0                  0                  0            -0.0573 - 0.0220i% a/ N2 t4 Y" _" y9 @8 I2 r$ e+ \4 w
       0.2789 + 0.0672i        0                  0                  0                  0         
    $ J- t! p; V2 i$ B6 o: N6 Y$ X   0.1595 - 0.0673i        0                  0                  0                  0         
    4 p' a3 G  r+ T        0                  0                  0                  0                  0         
    . F" \+ d( c' |* V  X* U        0                  0                  0                  0                  0          9 F, O* x# Q6 x$ u1 E& n: d
            0                  0                  0                  0                  0          # Z; [& F( Y3 P& ~9 t5 k+ i
            0                  0                  0                  0                  0         
    * i2 g. k6 x) D, ?        0                  0             0.0752 + 0.0448i   0.0784 + 0.0262i   0.1785 + 0.0769i
    : ^* i7 g. ]9 Q* F2 _# t        0                  0                  0                  0                  0         
    & o& I- |' q) \
    9 |: o; F; w% z, Q0 b& q1 S  Columns 6 through 10
    0 r5 u8 R% \8 M4 m7 `& B! Z
    , q- X  L5 ~! a1 j   0.0928 + 0.0303i  -0.2789 - 0.0591i  -0.1595 + 0.0831i        0                  0         
    9 |6 X4 w" b2 Q) Q1 L+ {# H        0                  0                  0                  0                  0         
    8 Q7 z1 z2 p) J# I- \        0                  0                  0                  0                  0          3 B% i' c" J8 ^) ^# ]% i# F
            0                  0                  0                  0                  0         
    & O: k% k9 t- R4 `' j. c* ?! _0 d   0.0579 + 0.0232i        0                  0                  0                  0          & C* ?# D% B/ ~7 _$ i
            0                  0                  0                  0                  0         
    2 z- {+ S$ E' r, f        0                  0            -0.2789 + 0.2447i        0                  0          : V5 F% ~. T! k
            0             0.2789 - 0.2196i        0            -0.6058 + 0.1645i  -0.5446 + 0.0169i9 }9 A' W& B# Z: A
            0                  0             0.6109 - 0.1483i        0            -0.4065 - 0.0361i5 R7 p4 S+ g8 \2 c/ Y5 r; {$ n6 |
            0                  0             0.5614 + 0.0342i   0.4157 + 0.0642i        0          ' K( x7 g; j7 V" u( n
            0                  0            -0.2300 + 0.0825i        0            -0.7120 + 0.0395i6 J( i' r! f3 e  Z; z8 P
            0            -0.0000 + 0.1972i        0                  0                  0          % q- y& k- i" k+ i
            0                  0                  0            -0.4441 + 0.2775i        0         
    2 {0 n5 ]' m0 |        0                  0                  0             0.7522 + 0.0930i   1.5728 - 0.1753i
    8 j. b" ]- a9 O! [
    % K! I+ h1 x+ E% v- U  Columns 11 through 14
    ; ^6 @* S+ x0 E5 Z+ i' D/ @' Q& U% B' ^) r
            0                  0                  0                  0         
    * G- Z) r% q( m  w9 x5 r% M        0                  0                  0                  0          ' f9 N, b0 h" G& |/ Y+ l/ B5 U& w
            0                  0            -0.0746 - 0.0434i        0          : `  i. Q5 q. Q
            0                  0            -0.0776 - 0.0247i        0         
    4 y4 e/ v1 f8 }; O        0                  0            -0.1764 - 0.0726i        0         
    & [/ Q- A4 m9 x5 }+ I- U4 P$ v        0                  0                  0                  0          9 x- U& Q4 k' s8 I1 D6 i5 T
            0             0.0000 - 0.1914i        0                  0          8 L1 P; I0 s+ v& T
       0.2339 - 0.0725i        0                  0                  0          ; _) q: F6 l3 d0 k
            0                  0             0.4441 - 0.2212i  -0.7245 + 0.0210i6 j) \! q* l* \: _2 o4 ]8 S
       0.7354 + 0.0592i        0                  0            -1.5296 + 0.3072i
    % b, Z! d. h3 I: C7 E" K7 `7 B* E        0                  0                  0                  0          2 E9 d/ S/ R# w2 q5 ~
            0                  0                  0                  0          1 Y/ `) s8 ?) G7 T/ x+ p
            0                  0                  0                  0          2 X% _/ o/ \6 Z9 D- L
            0                  0                  0                  0         
    8 D( b# x( m8 I+ M( l- X8 \/ e( q, C$ b. C# t* @& r( D

    7 G+ M8 J- [% ZPLoss =+ j+ X( i1 f4 @4 M3 [6 `# I

    % F- o& Q% q6 X, T0 f* M7 O    0.1349
    % r" g+ e$ d$ E/ ?6 F5 L+ Q  o
      ^4 x$ O1 U2 n/ F5 h, n7 E: g( u' x+ E3 k! h# W
    Qg =
    4 ]0 i! s# @3 S% Q& W
    , w1 M) h2 v( i3 w' ^& ?* n0 c  Columns 1 through 11
    . P2 A9 }9 {; `4 W  D2 R- n: W) U8 f8 f
             0         0         0         0         0         0         0         0         0    0.5919    0.3119
    ; m4 @9 O# x) n4 k& v# L* `
    + @& ~, r; E. Y5 h+ {. y' ~  Columns 12 through 13
    7 q4 |9 f6 b  \) Y8 S5 a7 P* g. }2 D1 ]3 d& o  }; m
        0.1972    0.5004% M, @3 @; k& O# B/ {% H- e5 ^% E
    ) |5 m2 I. d8 c" s/ Z
      ?8 j; i: ?/ x. P" x
    Pg =
    $ i% B$ v1 U% b4 t5 y$ |- f2 E/ S/ ~& I( |2 X
      Columns 1 through 11
    % F; M3 E9 M7 K4 [* J& b% X0 a- p3 j" ]1 ^+ {
             0         0         0         0         0         0         0         0         0    0.4000   -0.0000
    ! E+ c0 I) x. I/ D8 ^( w2 s0 O4 x7 Q" h) `( h1 Y, ?4 h$ o
      Columns 12 through 13
    2 X( X  d- S7 g! Y, w
      ^! x1 A! U. l   -0.0000   -0.0000
    ) f' @  y/ G) W2 i4 s+ k) s0 X8 B. R' E( `& G0 P* Y  }# j

    3 c9 Z! U7 |) |& dans =
    , ?/ ]4 R9 u6 S7 X) l" ?# L; b# {; E8 O. T. I# E7 H8 c
        0.1349
    8 G" V1 k: u! X- Q/ O; m3 j4 M8 w% Y1 h0 n  L" w4 f9 O# K

    , s5 x. I8 M# W# U1 i% ]& dans =
    - L; b1 R' v% D5 O/ b
    " U/ X# ]( @$ f; \    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-3-20 15:19

    Powered by Discuz! X3.5 Licensed

    © 2001-2025 Discuz! Team.

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