求助：使用matlab求潮流计算出现的问题

lgs_0003

由于本人现在正在研究MATLAB进行电力系统的仿真与计算，参考了孟祥萍编著的《电力系统分析》中牛顿法求潮流计算的方法。我自己建立了一个电力系统数学模型，根据该书的参考程序，没有什么大的问题，根据该书上的数据，计算出来的结果与该书的结果一致。问题就出现在我自己建立的数学模型无法进行迭代，如果误差设置成50，它就一次迭代完成了，而设置25，就处于一个无穷迭代了。我经过调试，发现该我输入的数据，迭代得到的最小误差居然是在31.x。经过昨晚一个晚上的修改数据，还是没有办法实现完成迭代，是否有人对这个问题研究过呢？这个迭代方法跟很多书本上的计算机算法一样的，是符合潮流计算的。在数据上，是有什么要求的吗？我设置的参数不可能和该书上的程序一样的，而且我也参考了王锡凡的《现代电力系统分析》里面的潮流计算方法，跟孟的书一样。
附上我的参数：
请输入节点数：n=7
8 m* q# H( y0 k) ]' M2 L请输入支路数：nl=7
请输入平衡母线节点号：isb=1
请输入误差精度：pr=10
  {& d3 L: A! X; @: R0 w) h请输入由支路参数形成的矩阵：B1=[1 4 7.017+52.045i 0 1.05 0;2 5 0.615i 0 1.05 0;3 6 0.216i 0 1.05 0;4 5 4.21+30.859i 0 1 0;4 7 4.21+30.859i 0 1 0;5 6 5.614+41.1453i 0 1 0;6 7 7.017+51.43i 0 1 0]
请输入由节点参数形成的矩阵：B2=[0 0 1.05 1.05 0 1;0.4 0 1.05 1.05 0 3;0.5 0 1.05 1.05 0 3;0 0 1 0 0 2;0 0 1 0 0 2;0 0.3+0.1i 1 0 0 2;0 0.9+0.3i 1 0 0 2]
请输入由节点号及其对地阻抗形成的矩阵：X=[1 0;2 0;3 0;4 0;5 0;6 0;7 0]
$ e/ v) J. P6 {" ?$ u" R节点1，2，3接发电机，而且节点3还带负荷，节点7为无穷大系统。

nanrui

这是本人本科毕设时用的程序，原数据亦是《现代电力系统分析》上一个3机5节点数据。计算是成功的。

lgs_0003

你说的那个模型参数是可以计算啊，我的程序也是可以运行，结果和他计算的参考结果一样，我是说我的问题出现在参数的设定上，而不是程序上。明确地说就是线路和节点的参数是否符合实际的运行方式。呵呵，不过还是非常感谢你的回复。

lixiaohui

由于做毕业设计需要用matlab编写的PQ分解法潮流计算程序，而自己的水平有限，所以只能求助论坛里的XDJM了。在这了先谢谢大家了。

lskai

你的功率设置上有问题
# A% i- v7 b0 b: @9 C1 @  i都要除100

zhangli000

学习了

Pikazhu

回复 2# nanrui
/ T9 w5 f9 v" w% q

   谢谢~~找了好久这样的小程序。