电力系统潮流计算C.docVIP

电力系统潮流计算 注：这是一个基于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支