电力系统潮流计算C程序

  • 格式:doc
  • 大小:51.50 KB
  • 文档页数:16

电力系统潮流计算
注:这是一个基于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>
#include <stdio.h>
#define f1(i) (i-1)
/* 把习惯的一阶矩阵的下标转化为C语言数组下标*/
#define f2(i,j,n) ((i-1)*(n)+j-1)
/* 把习惯的二阶矩阵的下标转化为C语言数组下标*/
/****************************************************
* 本子程序根据所给的支路导纳及有关信息,形成结点 *
* 导纳矩阵,如打印参数K=1,则输出电导矩阵G和电纳矩B *。