华北电力大学 电力系统稳态潮流上机计算程序结果

{
unsigned int i,j,k;
//消元过程
for(i=0;i<Dimension;i++)
{
//规格化过程:(每次消元前,先规格化,以便以下各行消元时,消元系数直接取待消
//列元素的值即可,也便于回代过程,而运算量并不增加)
for( j = i+1; j < Dimension; j++ )
D[y][y].diandao+=lineg;//尾节点自导纳中电导
D[y][x].dianna=D[x][y].dianna;//首节点与尾节点之间导纳与尾节点与首节点之间导纳相等
D[y][x].diandao=D[x][y].diandao;
}
if(gLineData[i].Type==2)//对线路为变压器的类型 对变压器进行π型等值 计算方法与导线型相同
{
if( FactorMatrix[j][i] != 0 ) //如果第j行第i列元素本就是0,则不需本列对应的消元过程
{
for( k = i + 1; k < Dimension; k++ ) //当FactorMatrix[i][k]=0,a[j][k]值不变,可省去运算
{
D[x][y].dianna+=(-1)*lineb/gLineData[i].Conductance;
D[x][y].diandao+=(-1)*lineg/gLineData[i].Conductance;
D[x][x].diandao+=lineg/gLineData[i].Conductance+(1-gLineData[i].Conductance)*lineg/(gLineData[i].Conductance*gLineData[i].Conductance);
{
FactorMatrix[i][j] = FactorMatrix[i][j] / FactorMatrix[i][i];
}
ConstVector[i] = ConstVector[i]/FactorMatrix[i][i];
for( j = i+1; j < Dimension; j++ ) //消去第i列(从i+1行到最后一行)
};
struct daona D[NODE_TOTAL_NUM][NODE_TOTAL_NUM];//定义一个二维数组存放导纳元素并对定义的数组赋0
int j;
for (i=0;i<NODE_TOTAL_NUM;i++)
for (j=0;j<NODE_TOTAL_NUM;j++)
{
D[i][j].dianna=0;
printf("\n节点数据\n");
for(i=0;i<NODE_TOTAL_NUM;i++)
{
fscanf(fp,"%d %d %f %f",&gNodeData[i].Index,
&gNodeData[i].Type,
&gNodeData[i].FirstInput,
float Resistance; //resistance of the line
float Reactance; //reactance of the line
float Conductance; //Line: conductance B/2,transformer:change ratio
lineb=(-1)*gLineData[i].Reactance/tab; //第i条线路电纳
if(gLineData[i].Type==1)//对线路为导线的类型
{
D[x][y].dianna+=(-1)*lineb;//首尾节点之间互导纳中的电纳
&gLineData[i].Reactance,
&gLineData[i].Conductance);
printf("%u %u %u %u %u %f %f %f\n",gLineData[i].Index,
gLineData[i].Type,
{
unsigned int Index; //node index
unsigned int Type; //node type: PQ:1,PV:2,balance point:0
float FirstInput; //PQ or PV:active power,balance point:V
D[x][y].diandao+=(-1)*lineg;//首尾节点之间互导纳中的电导
D[x][x].diandao+=lineg;//首节点自导纳中电导
D[x][x].dianna+=lineb+gLineData[i].Conductance;//首节点自导纳中电纳
D[y][y].dianna+=lineb+gLineData[i].Conductance;//尾节点自导纳中电纳
ConstVector[i-1] = ConstVector[i-1] - FactorMatrix[i-1][j] * ConstVector[j]; //a[i][k]=0时可以不算
}
}
return;
}
int main(int argc, char* argv[])
float SecondInput; //PQ:reactive power,PV:V,balance point:angle
};
struct NodeData gNodeData[NODE_TOTAL_NUM]; // define a global variable of node
unsigned int Status; //line status: on:1,off:0
unsigned int SrcNode; //source node index of this line
unsigned int DestNode;//destination node index of this line
{
printf("Hello World!\n");
//第一阶段 输入节点和线路数据
FILE* fp; //file structure
int i;
fp=fopen(NODE_DATA_FILENAME,"rb"); //open and read the node data输入节点数据
#define NODE_DATA_FILENAME "E:\\loadflow\\node.txt"
#define LINE_TOTAL_NUM 9
struct LineData
{
unsigned int Index; //line index
unsigned int Type; //line type: Line:1,tranformer:2,line connected to ground:3
}
}
//回代过程
for( i = Dimension-1; i > 0; i-- ) //Dimension-2:最后一个变量可直接获得,从n-1个变量求起
{
for(j = Dimension-1; j > i-1 ; j-- )
{
if( FactorMatrix[i-1][j] != 0 )
&gLineData[i].Type,
&gLineData[i].Status,
&gLineData[i].SrcNode,
&gLineData[i].DestNode,
&gLineData[i].Resistance,
};
struct LineData gLineData[LINE_TOTAL_NUM];
#define LINE_DATA_FILENAME "E:\\loadflow\\line.txt"
void SolveEquation(unsigned int Dimension,float FactorMatrix[14][14], float ConstVector[14]) //求解方程的函数 Ax=B
&gNodeData[i].SecondInput);
printf("%d %d %f %f\n",gNodeData[i].Index,
gNodeData[i].Type,
gNodeData[i].FirstInput,
gNodeData[i].SecondInput);
y=gLineData[i].DestNode-1;
tab=pow(gLineData[i].Resistance,2)+pow(gLineData[i].Reactance,2);//第i条线路电阻与电抗平方之和
lineg=gLineData[i].Resistance/tab; //第i条线路电导
gLineData[i].Status,
gLineData[i].SrcNode,
gLineData[i].DestNode,
gLineData[i].Resistance,
gLineData[i].Reactance,
}