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

 找回密码
 立即加入
搜索
查看: 30789|回复: 195

自己编写的Matlab-牛顿拉夫逊法潮流计算程序(短学期课程)

   火... [复制链接]
  • TA的每日心情
    开心
    2017-8-29 08:13
  • 签到天数: 2 天

    连续签到: 1 天

    [LV.1]初来乍到

    累计签到:2 天
    连续签到:1 天
    发表于 2008-8-29 21:42:53 | 显示全部楼层 |阅读模式

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

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

    ×
    自己编写的MATLAB-牛顿拉夫逊法潮流计算程序(我的短学期课程设计)
    ! f/ @( ?- t/ p' o# A4 k! o会在后面的帖子中进行详细的解释,附带上自己搜集的Matlab潮流计算论文8 h3 Q  Y) @  e" g1 S6 S4 q- }

    8 X; @# t0 M: b% o) n在My_NRflow.rar 中已经包括了测试文件和详细的正确答案1 L4 }- G2 R; |: o0 S% F

    / E, }- L* w( @, h, |[ 本帖最后由 zjuee 于 2008-8-29 21:57 编辑 ]-%-96814-%--%-96822-%--%-97248-%-

    My_NRflow.rar

    1.15 MB, 下载次数: 697, 下载积分: 威望 -2 点, 学分 -5 点

    程序+论文

    典型测试系统数据.rar

    163.92 KB, 下载次数: 340, 下载积分: 威望 -2 点, 学分 -5 点

    测试数据

    评分

    参与人数 1威望 +2 收起 理由
    qiuwei_002 + 2 支持原创,感谢分享!

    查看全部评分

    "真诚赞赏,手留余香"
    还没有人打赏,支持一下
    楼主热帖
    帖文化:【文明发帖 和谐互动】 社区精神:【创新、交流、互助、共享】
  • TA的每日心情
    郁闷
    2018-7-20 10:04
  • 签到天数: 4 天

    连续签到: 4 天

    [LV.2]偶尔看看I

    累计签到:4 天
    连续签到:4 天
    发表于 2018-7-17 23:04:29 | 显示全部楼层
    楼主,可以给我私法一下吗,非常需要,万分感谢!!!!只要测试数据就好了
    "真诚赞赏,手留余香"
    还没有人打赏,支持一下
    帖文化:【文明发帖 和谐互动】 社区精神:【创新、交流、互助、共享】
    回复 推荐 踩下

    使用道具 举报

  • TA的每日心情
    开心
    2019-5-12 22:02
  • 签到天数: 166 天

    连续签到: 1 天

    [LV.7]常住居民III

    累计签到:388 天
    连续签到:5 天
    发表于 2023-9-13 13:28:52 | 显示全部楼层
    支持原创,感谢分享
    "真诚赞赏,手留余香"
    还没有人打赏,支持一下
    帖文化:【文明发帖 和谐互动】 社区精神:【创新、交流、互助、共享】
    回复 推荐 踩下

    使用道具 举报

  • TA的每日心情
    开心
    2017-8-29 08:13
  • 签到天数: 2 天

    连续签到: 1 天

    [LV.1]初来乍到

    累计签到:2 天
    连续签到:1 天
     楼主| 发表于 2008-8-29 21:46:54 | 显示全部楼层

    各个子文件介绍!!!!

    1.powerflow(主程序)
    1 n- e4 l. H5 Z( }/ [2.y(节点导纳矩阵)' s) B9 A' `' a, @8 y% a8 X
    3. run_NR (利用NR法计算潮流)
    2 X7 Q' w" i- f4. openfile (打开文件)% X- v, y" O3 u
    5. printY (显示Y矩阵)
    7 f( |1 G( k/ {( U+ ^$ s7 O6. newnode (节点重新编号)
    2 p) ^# ]1 a* k$ v7. Renode (还原节点)7 O0 }. ], l; V& b
    8. form_jac (求取j矩阵)$ C1 n, `, I- C$ g3 S
    9. dPQ (求取有功无功不平衡量)# ?/ K% a# {" `2 P! g
    10. PowerResult(计算节点电压相角)
    ( q& T# q! Y4 |9 a11. PowerLosses(计算支路性息)
    # H! m, M& H' y: b+ T12. printData(显示支路性息)
    "真诚赞赏,手留余香"
    还没有人打赏,支持一下
    帖文化:【文明发帖 和谐互动】 社区精神:【创新、交流、互助、共享】
  • TA的每日心情
    开心
    2017-8-29 08:13
  • 签到天数: 2 天

    连续签到: 1 天

    [LV.1]初来乍到

    累计签到:2 天
    连续签到:1 天
     楼主| 发表于 2008-8-29 21:49:30 | 显示全部楼层

    运行方式!!!

    1 打开My_NRflow程序包
    % g7 ?4 E8 F+ B# e, M2 将matalb路径改为My_NRflow程序包所在地址
    / A* c- u1 B+ @4 ~* X4 a2 ^3 打开powerflow(主程序)
    ( \4 W& G8 p  G# s4 运行并选择输入文件 共有6个text2 y6 x+ ]0 J7 u, z: _
    5 在My_NRflow程序包所在地址上一层文件夹内有结果文件——output.dat1 h" b' U7 k. d( r
    6 对比运行结果与正确答案
    "真诚赞赏,手留余香"
    还没有人打赏,支持一下
    帖文化:【文明发帖 和谐互动】 社区精神:【创新、交流、互助、共享】
  • TA的每日心情
    开心
    2017-8-29 08:13
  • 签到天数: 2 天

    连续签到: 1 天

    [LV.1]初来乍到

    累计签到:2 天
    连续签到:1 天
     楼主| 发表于 2008-8-29 21:51:58 | 显示全部楼层

    母线(节点)数据

    格式:节点编号, 节点电压,节点相角,节点注入有功,节点注入无功,节点类型
    ) H* s$ x/ I- I* Q8 a/ Q 7 H& y3 P0 n( V* Q3 `7 m6 @
    说明:其中节点类型取值为1、2、3,分别对应PQ节点、PV节点和平衡节点;节点注入功率如果是负荷,取负值;否则取正值;节点编号的原则是先PQ节点,然后PV节点,最后平衡节点。' [. n& x; L) ], B5 d4 `
    对于课本201页上的例题3-3,首先形成一个数据文件test.m,内容如下:
    $ l& N: i/ l9 ]5 c3 ^$ U/ p0 i; H/ e; [
    节点编号, 节点电压,节点相角,节点注入有功,节点注入无功,节点类型! i* y6 q4 D" Y0 a: ?8 `* M
    bus=[
    4 H) h( \( C( }) @+ D     1   1.00  0.00   -1.60  -0.80  1;
    $ _. G9 Q3 a0 r, W- T; n9 b     2   1.00  0.00   -2.00  -1.00  1;
    - c+ u3 q3 n0 D$ e$ _5 g, q     3   1.00  0.00   -3.70  -1.30  1;
    ' c% d  k' u* E     4   1.05  0.00    5.00   0.00   2;
    9 c$ {7 S* V* _- o$ _3 N     5   1.05  0.00    0.00   0.00   3];
    "真诚赞赏,手留余香"
    还没有人打赏,支持一下
    帖文化:【文明发帖 和谐互动】 社区精神:【创新、交流、互助、共享】
  • TA的每日心情
    开心
    2017-8-29 08:13
  • 签到天数: 2 天

    连续签到: 1 天

    [LV.1]初来乍到

    累计签到:2 天
    连续签到:1 天
     楼主| 发表于 2008-8-29 21:53:04 | 显示全部楼层

    线 路 数 据 格 式

    对于线路(包括对地支路和变压器),原始数据格式定义如下:
    ) X5 S8 c7 o* L& W! I2 G( ]7 v* M9 P1 i9 {, y, P5 x4 D' S
    节点I, 节点J, 线路电阻R,线路电抗X,线路电导G,线路电纳B,变压器变比K
    , ~. U0 {# W! K6 T其中,K=0表示普通线路,K>0表示变压器线路的非标准变比; 8 [  \2 e7 |8 r
    line = [
    7 F' g9 q2 _0 u& C( n' ]   1   2   0.04    0.25     0.0    0.25    0;
    1 Y5 v* |5 m: B2 u/ n   1   3   0.10    0.35     0.0    0.0      0;
    . K$ h, L, e3 M1 A   2   3   0.08    0.30     0.0    0.25    0;, \* v; @4 o# F% p, t" t) P4 F
       5   3   0.00    0.03     0.0    0.0     1.05;
    + Q0 @+ m# V4 `  `& [   4   2   0.00    0.015   0.0    0.0     1.05];
    "真诚赞赏,手留余香"
    还没有人打赏,支持一下
    帖文化:【文明发帖 和谐互动】 社区精神:【创新、交流、互助、共享】
  • TA的每日心情
    开心
    2017-8-29 08:13
  • 签到天数: 2 天

    连续签到: 1 天

    [LV.1]初来乍到

    累计签到:2 天
    连续签到:1 天
     楼主| 发表于 2008-8-29 22:00:40 | 显示全部楼层

    节点重新编号的问题

    在电力系统的潮流计算中,总是假设节点编号依据PQ节点、PV节点、平衡节点的次序进行编号的。但是,实际电力系统的节点编号具有一定的随意性,并不总是按照这个原则形成原始数据的。此外,有些时候需要改变节点类型,如改变原来的平衡节点3为PV节点,将原来的PV节点4变为平衡节点等。因此,当原始节点数据不按照PQ节点、PV节点、平衡节点的次序进行编号时,需要我们自己写一个节点重新编号子程序。 ) C7 p& \% T( G
        我们老师给了相应的程序,其实挺简单的,就是排一下矩阵,呵呵 就是 newnode (节点重新编号)程序
    "真诚赞赏,手留余香"
    还没有人打赏,支持一下
    帖文化:【文明发帖 和谐互动】 社区精神:【创新、交流、互助、共享】
  • TA的每日心情
    开心
    2019-11-15 20:29
  • 签到天数: 11 天

    连续签到: 3 天

    [LV.3]偶尔看看II

    累计签到:11 天
    连续签到:3 天
    发表于 2008-8-29 22:35:51 | 显示全部楼层
    支持原创的东西% [9 `, X, h# u  u6 L8 s
    应该给楼主加分呵呵
    "真诚赞赏,手留余香"
    还没有人打赏,支持一下
    帖文化:【文明发帖 和谐互动】 社区精神:【创新、交流、互助、共享】
  • TA的每日心情
    开心
    2017-8-29 08:13
  • 签到天数: 2 天

    连续签到: 1 天

    [LV.1]初来乍到

    累计签到:2 天
    连续签到:1 天
     楼主| 发表于 2008-8-29 23:04:35 | 显示全部楼层
    "真诚赞赏,手留余香"
    还没有人打赏,支持一下
    帖文化:【文明发帖 和谐互动】 社区精神:【创新、交流、互助、共享】
  • TA的每日心情
    开心
    2019-10-7 20:14
  • 签到天数: 2 天

    连续签到: 1 天

    [LV.1]初来乍到

    累计签到:2 天
    连续签到:1 天
    发表于 2008-8-29 23:31:27 | 显示全部楼层
    看了一下,写的很好,而且解释的很全面。因为自己也写过,看到了就会有些亲切的感觉,而且也找出了很多自己写的不足的地方。多谢楼主的分享。
    "真诚赞赏,手留余香"
    还没有人打赏,支持一下
    帖文化:【文明发帖 和谐互动】 社区精神:【创新、交流、互助、共享】

    该用户从未签到

    尚未签到

    发表于 2008-9-1 09:15:19 | 显示全部楼层

    谢谢分享

    辛苦楼主啦,谢谢共享
    "真诚赞赏,手留余香"
    还没有人打赏,支持一下
    帖文化:【文明发帖 和谐互动】 社区精神:【创新、交流、互助、共享】
    您需要登录后才可以回帖 登录 | 立即加入

    本版积分规则

    招聘斑竹

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

    GMT+8, 2024-5-4 15:43

    Powered by Discuz! X3.5 Licensed

    © 2001-2024 Discuz! Team.

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