潮流计算的一个问题

zhuwan617

最近用MATLAB编写潮流程序的时候发现一个小问题：
0 g% A0 Y: w2 y/ [( O    一般在编写潮流时我们会设定功率的符号（负荷为流出母线，一般设为正值，发电机向母线注入功率，则功率数值设为负值），但是我发现改变一下符号（即负荷为负，发电机功率为正时）出来的结果与之前的计算结果不一致，请问各位大侠们如何解释这个问题！
    我查了下原因，两种不同的符号规定方法导致迭代修正量的符号发生了变化，然而对初值进行修正时还是采用：修正值=初值+修正量，这个公式，从而导致了收敛到另外一个平衡点。请问这样解释合理么？
0 h7 ?6 Z! j: T" L3 b& C! |) Q. F    谢谢各位大侠了哇~

木易山水

两种算法差别大么？我的看法是：潮流计算没有求出解析解，而是是数值解，lz说的两种方法是两种不同的过程，结果肯定不一样，这就牛拉法和高斯法之间的差别是一个道理。
只要是最终满足收敛条件就可以了，至于方向，和电路里面理解的差不多吧。

baobobo

雅克比矩阵的符号错了……

Daniel_smile

牛拉法潮流计算的数学模型是节点功率平衡方程，即
节点注入功率=节点流向与之相邻节点的功率，记为式（1）；牛拉法仅仅是求解这组非线性方程组的方法。
节点注入功率=发电机注入节点功率－该节点负荷功率 注：（功率都取正值，因为用的“－”号；对应你提到的负荷为负，发电机为正的情况（记为情况a））；当负荷为正、发电机为负（记为情况b）时，注意式（1）等号左侧与情况a相比是取了相反数，所以等号右侧的“节点流向与之相邻节点的功率”与情况a相比也要取相反数。
1 u/ D5 J3 l# z+ B6 \  P# \体现在牛拉法迭代的过程中，雅克比矩阵两张情况也要取相反数，而两种情况电压幅值与相角更新是相同的，不必修改。