电力系统潮流计算程序
#include <math.h>
#include <iostream.h>
#define N 10
//#define n 5
void main()
{
double t; //中间变量
int n; //统计节点个数
int i,j,k; //循环变量
int r; //PV 节点个数
//int r=0;
int m; //PQ 节点个数
double R[N][N];//节点间的电阻
double X[N][N];//节点间的电抗
double g1[N][N]={0}; //节点间的电导
double b1[N][N]={0}; //节点间的电纳
double G[N][N]={0};
double B[N][N]={0};
double RR[N][N]={0};
double XX[N][N]={0};
double KK[N][N]={0};
double gg1[N][N]={0}; //节点间对地的电导
double bb1[N][N]={0}; //节点间对地的电纳
double P[N]={0}; //节点有功初值,平衡节点除外
double Q[N]={0}; //节点无功初值,平衡节点除外
double B1[N-1][N-1]; //系数矩阵 1
double B2[N-1][N-1]; //系数矩阵 2
double A1[N-1][N-1]; //系数矩阵 1 的逆矩阵
double A2[N-1][N-1]; //系数矩阵 2 的逆矩阵
double a[N][2*N]; //求逆矩阵时用的矩阵
double w=0,s=0,l,m1; //求逆矩阵时用的变量
double U[N]; //节点电压幅值
double rad[N]; //节点电压幅角
double cp[N]; //有功不平衡量
double cq[N]; //无功不平衡量
double ec; //容许误差
cout << "请将平衡节点设为零,接下来设 PQ 节点,最后是 PV 节点。";
cout << "请输入网络节点数:";
cin >> n ;