nlf.h

计算电力潮流

/* 定义变量类型*/

typedef double PointData;
typedef struct{int number;int i_node;int j_node;PointData r; PointData x; PointData bk;}
BRANCH;  
typedef struct{int node;PointData p;PointData q; PointData v; PointData q_max; PointData q_min;}
GENRT;	
typedef struct{int node;PointData p;PointData q;}
LOAD;
typedef struct{int node;PointData Pg;PointData Qg;PointData Pl;PointData Ql;PointData v;int flag;}
NODE;
typedef struct{int node;PointData x;}//潮流计算输入数据
GROUND;

typedef struct{int node;PointData e;PointData f;}
VOLTAGE;
typedef struct{int node;PointData p;PointData q;PointData v2;PointData dp;PointData dq;
PointData dv2;}
PQV;                               
typedef struct{PointData aii;PointData bii;} //迭代修正向量
AB;

typedef struct{int addres; PointData Gii; PointData Bii;}
YY_MATRIX;
typedef struct{int number;PointData Gij;PointData Bij;}
Y_MATRIX;  
typedef struct{int addres;PointData data;} 
GENE_YYMATRIX;
typedef struct{int number;PointData data;}
GENE_YMATRIX;                      //矩阵数据变量

typedef struct{int node; int amount;}
DATA_II;                            //节点优化数据

/*..........................................................*/
#define EXIT {printf("......the space is over!\n");exit(0);}
#define SMALL 0.00001    // 溢出控制
#define BIG 100000.0
#define EPS 0.0000000001//迭代精度
/*............................................................*/
#include <stdio.h>
#include <math.h>
#include "stdlib.h"
//#include <iostream.h>

class nlf
{
  int iter_number;
  int N,NS,number_genrt,number_load,number_branch,number_pv,number_ground;
  int *pv,*ppv,*o_n,*n_o;
  FILE   *in,*out;
  BRANCH *branch;
  GENRT  *genrt;
  LOAD   *load;
  NODE *node;
  GROUND *ground;

  VOLTAGE *voltage;
  PQV  *dpqv;

  Y_MATRIX   *y_matrix;
  YY_MATRIX  *yy_matrix;
  GENE_YMATRIX  *gene_ymatrix;
  GENE_YYMATRIX *gene_yymatrix;
  AB         *ab;
  
  void readdata();// 潮流计算前期准备（包括节点优化和支路优化）
  void node_order();
  void re_order();
  void branch_rank();
  
  void form_ymatrix();//导纳阵形成
  void initia_pqv();
  void initia_dpqv();
  void gene_matrix();//修正矩阵
  void back_solve();//回代求解
  void amend_voltage();//修正电压
public:

  nlf();
  ~nlf();
  void BeginPowerflow();
  void iterative(); 
  void result();
  void free_point();
};