这篇资料对于潮流计算很方便!标号,起始结点,终止结点,支路电阻参数,支路电抗参数,支路对地导纳参数
1,2,5,0.0,0.063,0.0,
2,5,9,0.019,0.072,0.075,
3,6,9,0.012,0.101,0.105,
4,3,6,0.0,0.059,0.0,
5,6,8,0.039,0.17,0.179,
6,4,8,0.017,0.092,0.079,
7,5,7,0.032,0.161,0.153,
8,4,7,0.01,0.085,0.088,
9,1,4,0.0,0.058,0.0,
潮流程序chaoliu.txt文档:
#include<stdio.h>
#include<math.h>
#define N 9 /*总结点数*/
#define M 6 /*PQ结点数*/
#define K 9 /*线路数*/
#define eps 1e-4
void guass(int n,int m,float c[],float b[][N],float x[]) /*高斯函数*/
{
float a[N][N],y[N];
int i,j,k;
for(i=0;i<m;i++)
{
y[i]=c[i];
for(j=0;j<m;j++) a[i][j]=b[i+n][j+n];
}
for(k=0;k<m-1;k++)
for(i=k+1;i<m;i++)
{
for(j=k+1;j<m;j++) a[i][j]=a[i][j]-a[k][j]*a[i][k]/a[k][k];
y[i]=y[i]-y[k]*a[i][k]/a[k][k];
}
for(i=m-1;i>=0;i--)
{
for(j=i+1;j<m;j++) y[i]-=a[i][j]*x[j];
x[i]=y[i]/a[i][i];
}
}
!