电力系统潮流计算注：这是一个基于 N－R 法的潮流计算通用程序，仅提供了子程序，需要做些处理才能成为一个可运行的计算程序！此程序非我原创，仅与大家共享!!！/***＊*＊＊＊＊*＊**＊＊＊**＊*＊＊＊**＊*＊＊＊＊*＊***＊**＊＊＊****＊＊＊*＊*****＊＊***＊＊＊*＊** 这里提供的是电力系统潮流计算机解法的五个子程序,采纳的方法是 *＊ Newton_Raphson 法。 ** 程序中所用的变量说明如下： ＊＊ N：网络节点总数. M：网络的 PQ 节点数。 *＊ L：网络的支路总数. N0：雅可比矩阵的行数. ＊* N1：N0+1 K：打印开关.K=1，则打印；否则,不打印。** K1：子程序 PLSC 中推断输入电压的形式。K1=1,则为极座标形式。否则** 为直角坐标形式. ＊* D：有功及无功功率误差的最大值. ＊＊ G(I，J):Ybus 的电导元素（实部）。 *＊ B(I，J)：Ybus 的电纳元素（虚部）。 *＊ G1(I） :第 I 支路的串联电导。 B1(I）:第 I 支路的串联电纳. *＊ C1（I) :第 I 支路的 pie 型对称接地电纳. ＊＊ C(I,J)：第 I 节点 J 支路不对称接地电纳. ＊* CO（I) :第 I 节点的接地电纳. ＊＊ S1(I) ：第 I 节点的起始节点号. E1(I)：第 I 节点的终止节点号. *＊ P（I) ：第 I 节点的注入有功功率. Q(I）:第 I 节点的注入无功功率.** P0（I） ：第 I 节点有功功率误差. Q0（I):第 I 节点无功功率误差. ** V0（I) :第 I 节点(PV 节点)的电压误差(平方误差）. ＊＊ V（I） ：第 I 节点的电压误差幅值。 ＊* E（I) ：第 I 节点的电压的实部。 F(I):第 I 节点的电压的虚部. ＊＊ JM（I,J）：Jacoby 矩阵的第 I 行 J 列元素。 ＊＊ A(I，J)：修正方程的增广矩阵，三角化矩阵的第 I行 J 列元素,运算结 ** 束后 A 矩阵的最后一列存放修正的解。 ** P1（I) ：第 I 支路由 S1（I）节点注入的有功功率. ＊* Q1(I) :第 I 支路由 S1（I)节点注入的无功功率. ＊* P2(I） ：第 I 支路由 E1(I)节点注入的有功功率. ＊＊ Q2(I） :第 I 支路由 E1(I）节点注入的无功功率。 ＊* P3（I） :第 I 支路的有功功率损耗。 ** Q3（I) ：第 I 支路的无功功率损耗。 ** ANGLE（I）：第 I 节点电压的角度。 ＊**＊***＊*＊＊＊＊＊*＊＊**＊＊＊＊*＊***＊****＊*＊**＊**＊*＊*＊**＊＊＊＊＊**＊＊*＊＊***＊＊*＊*/＃include 〈math.h>#inc...
