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

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

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

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

    连续签到: 1 天

    [LV.1]初来乍到

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

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

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

    ×
    自己编写的MATLAB-牛顿拉夫逊法潮流计算程序(我的短学期课程设计)
    ( y: I8 w, r) U会在后面的帖子中进行详细的解释,附带上自己搜集的Matlab潮流计算论文) D- R2 y; q, u6 U6 z# O

    + |' G3 ~6 L% z  p2 e- T在My_NRflow.rar 中已经包括了测试文件和详细的正确答案
    ! H1 e$ s( G/ e
    ) U" c' M* \3 @[ 本帖最后由 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

    累计签到:393 天
    连续签到:2 天
    发表于 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(主程序); E  b* ]6 s/ {, L. F
    2.y(节点导纳矩阵)
    # ~3 P8 ]# l& _, Z# K3. run_NR (利用NR法计算潮流)
    & m7 B) r6 A2 _, ?4. openfile (打开文件)
    5 r# g( N# o# v# N5. printY (显示Y矩阵). g9 I8 @5 Q, I1 H
    6. newnode (节点重新编号)
    3 Z+ U% H( P+ X" E7. Renode (还原节点)' r: s4 T% v9 Y3 ~- f1 T8 {
    8. form_jac (求取j矩阵)
    0 t( _$ K0 h" D+ w6 D9. dPQ (求取有功无功不平衡量)! ~, }0 _/ R" S7 D1 d  F  [
    10. PowerResult(计算节点电压相角)  Z  P) m& i# r0 q& |& H
    11. PowerLosses(计算支路性息): w+ \; a# o, q) S) ?1 k! ]
    12. printData(显示支路性息)
    "真诚赞赏,手留余香"
    还没有人打赏,支持一下
    帖文化:【文明发帖 和谐互动】 社区精神:【创新、交流、互助、共享】
  • TA的每日心情
    开心
    2017-8-29 08:13
  • 签到天数: 2 天

    连续签到: 1 天

    [LV.1]初来乍到

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

    运行方式!!!

    1 打开My_NRflow程序包! E$ v; T2 D6 }: {- ?9 W
    2 将matalb路径改为My_NRflow程序包所在地址
    1 G2 \5 m4 W6 U" E7 U3 打开powerflow(主程序)
    + a; u& u, Z2 U; t8 ?" K7 B' y4 运行并选择输入文件 共有6个text' V; {9 j# f. \" N" B# i
    5 在My_NRflow程序包所在地址上一层文件夹内有结果文件——output.dat
    0 D, H" j! |$ U+ B6 ^6 对比运行结果与正确答案
    "真诚赞赏,手留余香"
    还没有人打赏,支持一下
    帖文化:【文明发帖 和谐互动】 社区精神:【创新、交流、互助、共享】
  • TA的每日心情
    开心
    2017-8-29 08:13
  • 签到天数: 2 天

    连续签到: 1 天

    [LV.1]初来乍到

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

    母线(节点)数据

    格式:节点编号, 节点电压,节点相角,节点注入有功,节点注入无功,节点类型
    $ e7 R4 w1 X( T6 w / @1 B# H( q6 C8 a# H0 T# Z
    说明:其中节点类型取值为1、2、3,分别对应PQ节点、PV节点和平衡节点;节点注入功率如果是负荷,取负值;否则取正值;节点编号的原则是先PQ节点,然后PV节点,最后平衡节点。! r6 Q8 }$ A5 m
    对于课本201页上的例题3-3,首先形成一个数据文件test.m,内容如下:
    ) I5 m7 s% B. D  E. G! c# ^! U$ v* L9 ?7 ^( ?; A% C* y2 j" f3 t
    节点编号, 节点电压,节点相角,节点注入有功,节点注入无功,节点类型% X0 d1 y( H4 w! T' R
    bus=[& c2 g, U* n3 h
         1   1.00  0.00   -1.60  -0.80  1;
    2 T/ p5 j; c" {1 y1 s( W2 y     2   1.00  0.00   -2.00  -1.00  1;
    9 ?, ]# s+ u* ~     3   1.00  0.00   -3.70  -1.30  1;/ P- H6 L' x. j; c/ u4 ^0 D
         4   1.05  0.00    5.00   0.00   2;
    $ P& C* U* k0 \7 I! S( d/ A     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 | 显示全部楼层

    线 路 数 据 格 式

    对于线路(包括对地支路和变压器),原始数据格式定义如下:
    # z& h& b8 P* S+ k" M" Y- \
    ; m6 a% o% Y# ]# v+ @  Q" Q节点I, 节点J, 线路电阻R,线路电抗X,线路电导G,线路电纳B,变压器变比K! l: c1 J. y& k. U! e; f
    其中,K=0表示普通线路,K>0表示变压器线路的非标准变比;
    * W2 K6 O6 }6 _2 r! f4 t. Y, jline = [
    0 h* w; e/ ^2 q$ H7 ?  ^   1   2   0.04    0.25     0.0    0.25    0;
    % b  W7 k5 o' I7 U! M8 k   1   3   0.10    0.35     0.0    0.0      0;
    0 j! r( X& {0 i. l   2   3   0.08    0.30     0.0    0.25    0;
    : N7 I) y8 X5 F% x! f( `8 }6 W$ W   5   3   0.00    0.03     0.0    0.0     1.05;
    3 d) X* k- @  t6 ^   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节点、平衡节点的次序进行编号时,需要我们自己写一个节点重新编号子程序。
    # {' {9 t( h5 d2 t  z    我们老师给了相应的程序,其实挺简单的,就是排一下矩阵,呵呵 就是 newnode (节点重新编号)程序
    "真诚赞赏,手留余香"
    还没有人打赏,支持一下
    帖文化:【文明发帖 和谐互动】 社区精神:【创新、交流、互助、共享】
  • TA的每日心情
    开心
    2019-11-15 20:29
  • 签到天数: 11 天

    连续签到: 3 天

    [LV.3]偶尔看看II

    累计签到:11 天
    连续签到:3 天
    发表于 2008-8-29 22:35:51 | 显示全部楼层
    支持原创的东西
    " E+ o0 L/ i) e! m: g: ^+ y3 Q应该给楼主加分呵呵
    "真诚赞赏,手留余香"
    还没有人打赏,支持一下
    帖文化:【文明发帖 和谐互动】 社区精神:【创新、交流、互助、共享】
  • 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, 2026-3-18 10:14

    Powered by Discuz! X3.5 Licensed

    © 2001-2025 Discuz! Team.

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