交流社区 › 电力系统与电网技术 › 请教一个关于Matlab编写的IEEE14节点潮流计算的问题

迪拉克之海 发表于 2008-4-29 22:27:26

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

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

[问题已经解决]

[ 本帖最后由 迪拉克之海 于 2008-5-7 13:37 编辑 ]

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 发表于 2008-5-7 11:24:14

u =

Columns 1 through 11

    1.0516    1.0474    1.0550    1.0549    1.0497    1.0327    1.0581    1.0123    1.0138    1.0450    1.0100

Columns 12 through 14

    1.0900    1.0700    1.0600


Theta =

Columns 1 through 11

   -0.2599   -0.2628   -0.2578   -0.2632   -0.2645   -0.2794   -0.2323   -0.1790   -0.1520   -0.0872   -0.2229

Columns 12 through 14

   -0.2323   -0.2483         0


Sslack =

   2.3249 - 0.0250i


S =

Columns 1 through 5

      0             0.0507 + 0.0332i      0                  0                  0         
-0.0506 - 0.0329i      0            -0.0394 - 0.0251i      0                  0         
      0             0.0396 + 0.0254i      0                  0                  0         
      0                  0                  0                  0             0.0166 + 0.0087i
      0                  0                  0            -0.0165 - 0.0086i      0         
-0.0917 - 0.0280i      0                  0                  0            -0.0573 - 0.0220i
   0.2789 + 0.0672i      0                  0                  0                  0         
   0.1595 - 0.0673i      0                  0                  0                  0         
      0                  0                  0                  0                  0         
      0                  0                  0                  0                  0         
      0                  0                  0                  0                  0         
      0                  0                  0                  0                  0         
      0                  0             0.0752 + 0.0448i   0.0784 + 0.0262i   0.1785 + 0.0769i
      0                  0                  0                  0                  0         

Columns 6 through 10

   0.0928 + 0.0303i-0.2789 - 0.0591i-0.1595 + 0.0831i      0                  0         
      0                  0                  0                  0                  0         
      0                  0                  0                  0                  0         
      0                  0                  0                  0                  0         
   0.0579 + 0.0232i      0                  0                  0                  0         
      0                  0                  0                  0                  0         
      0                  0            -0.2789 + 0.2447i      0                  0         
      0             0.2789 - 0.2196i      0            -0.6058 + 0.1645i-0.5446 + 0.0169i
      0                  0             0.6109 - 0.1483i      0            -0.4065 - 0.0361i
      0                  0             0.5614 + 0.0342i   0.4157 + 0.0642i      0         
      0                  0            -0.2300 + 0.0825i      0            -0.7120 + 0.0395i
      0            -0.0000 + 0.1972i      0                  0                  0         
      0                  0                  0            -0.4441 + 0.2775i      0         
      0                  0                  0             0.7522 + 0.0930i   1.5728 - 0.1753i

Columns 11 through 14

      0                  0                  0                  0         
      0                  0                  0                  0         
      0                  0            -0.0746 - 0.0434i      0         
      0                  0            -0.0776 - 0.0247i      0         
      0                  0            -0.1764 - 0.0726i      0         
      0                  0                  0                  0         
      0             0.0000 - 0.1914i      0                  0         
   0.2339 - 0.0725i      0                  0                  0         
      0                  0             0.4441 - 0.2212i-0.7245 + 0.0210i
   0.7354 + 0.0592i      0                  0            -1.5296 + 0.3072i
      0                  0                  0                  0         
      0                  0                  0                  0         
      0                  0                  0                  0         
      0                  0                  0                  0         


PLoss =

    0.1349


Qg =

Columns 1 through 11

         0         0         0         0         0         0         0         0         0    0.5919    0.3119

Columns 12 through 13

    0.1972    0.5004


Pg =

Columns 1 through 11

         0         0         0         0         0         0         0         0         0    0.4000   -0.0000

Columns 12 through 13

   -0.0000   -0.0000


ans =

    0.1349


ans =

    0.8414

hhusunguoqiang 发表于 2008-6-18 11:01:15

matlab 角度计算单位都为弧度，是否换算一下就和真值差不多了！

s31xili 发表于 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

看着好复杂啊，不过我也将要编

查看完整版本: 请教一个关于Matlab编写的IEEE14节点潮流计算的问题