电力系统潮流计算c语言程序,电力系统潮流计算C语言程序.doc

#include #include #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]; //系数矩阵1double B2[N-1][N-1]; //系数矩阵2double 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 ;cout << "请输入PV节点个数(r> r;m=n-r-1;cout << "PQ节点个数(n-r-1)为：" << "m=" << m << endl;cout << "请输入节点间电阻：\n";

for( i=0 ; i> R[j]; } }

cout << "\n请输入节点间电抗：\n"; for( i=0 ; i> X[j]; } }cout << "\n请输入节点对地容抗：\n"; for( i=0 ; i> XX[j]; } }cout << "\n请输入节点间变比：\n"; for( i=0 ; i