牛顿－拉夫逊法潮流计算 基本步骤： 　　（1）形成节点导纳矩阵 　　（2）将各节点电压设初值U， 　　（3）将节点初值代入相关求式，求出修正方程式的常数项向量 　　（4）将节点电压初值代入求式，求出雅可比矩阵元素 　　（5）求解修正方程，求修正向量 　　（6）求取节点电压的新值 　　（7）检查是否收敛，如不收敛，则以各节点电压的新值作为初值自第3步重新开始进行狭义次迭代，否则转入下一步 　　（8）计算支路功率分布

IEEE 33节点系统（仅PQ节点）三相不平衡潮流的前推回代算法，MATLAB中直接运行出结果。可在本程序基础上做修改，进行进一步开发。-IEEE 33-node system (only the PQ node) three-phase unbalanced trend forward and backward substitution algorithm, results run directly MATLAB.