请问用MATLAB编写潮流计算程序时，比如PQ法或者NR法，设置循环次数的意义在哪里呢？

在潮流计算程序中设这迭代次数上限是为了应对潮流发散的情况。

比如一个正确的潮流计算 一般用牛拉法是 迭代6次，PQ分解法是 迭代12次。若潮流出现发散情况，迭代到20次，dP和dQ仍然不能收敛到你要的精度（比如0.0001），此时我们判定潮流发散，若不设一个迭代上限，可想而知，程序会一直迭代下去，此时数值已完全失去意义。

潮流发散意味着计算失败，这可能由两个原因导致：1数据存在问题，初值不在收敛域以内；2程序本身问题，考虑不全面，有bug。

同时也可能产生两个现象，一种是 潮流直接发散，若干次迭代后，数值已经很离谱，此时数值肯定没有意义。还有一种是，雅克比矩阵接近病态，收敛极其缓慢（也可能是初值选取不当），此时适当的放大收敛次数限制，有一定可能得到收敛后的结果。

多编程序，多思考，希望对你有帮助！

1、60年代中期，基于导纳矩阵的牛顿—拉夫逊法。牛顿一拉夫逊法(简称牛顿法)是数学中解决非线性方程式的典型方法，有较好的收敛性。在解决电力系统潮流计算问题时，是以导纳矩阵为基础的，因此，只要我们能在迭代过程中尽可能保持方程式系数矩阵的稀疏性，就可以大大提高牛顿法潮流程序的效率。自从60年代中后期，在牛顿法中利用了最佳顺序消去法以后，牛顿法在收敛性、内存要求、速度方面都超过了阻抗法，成为60年代末期以后广泛采用的优秀方法。牛拉法的要点是把非线性方程式的求解过程变成反复地求解线性的修正方程式过程，即通常所称的逐次线性化过程。

2、70年代中期，PQ分解法。由于交流高压电网中输电线路等元件的R<<X，因此有功功率的变化主要决定于电压相位角的变化，而无功功率的变化则主要决定于电压模值的变化。这个特性反映在极坐标形式的牛顿法修正方程式的元素上，是N及J二个子块元素的数值相对于H、L二个子块的元素要小得多。

这个方法，根据电力系统的特点，抓住主要矛盾，对纯数学的牛顿法进行了改进，从而在内存容量及计算速度方面都大大向前迈进了一步。使一个32K内存容量的数字计算机可以计算1000个节点系统的潮流问题，此方法计算速度已能用于在线计算，作系统静态安全监视。目前，我国很多电力系统都采用了PQ分解法潮流程序。

3、在有些应用场合，对计算精度的要求不高，而对计算速度要求较高。如输电网规划初期，只需要考虑有功功率平衡的问题，而不需要考虑无功功率平衡和电压的问题，这时可以对潮流方程进行简化处理，用直流潮流进行计算。直流潮流方程是一个线性方程组，求解不需迭代，不存在收敛的问题；导纳矩阵是稀疏的，可利用稀疏技术进一步提高计算速度；当高压电网满足R<<X时时，计算误差通常在3%-10%之内，可满足对精度要求不高的场合。直流潮流不能计算节点的电压和无功功率潮流。

4、为满足不同的需求开发了各种潮流算法：动态潮流、保留非线性的潮流、最小化潮流计算法、自动调整潮流、最优潮流、交直流系统潮流、直流潮流、随机潮流、三相潮流，含有柔性元件的潮流，并行算法等。

---