电力网潮流计算C++实现

基于C++的潮流计算，利用牛顿拉夫逊迭代法计算。

#include “stdafx.h“
#include “admittance_matrix.h“
#include “Complex.h“
#include 
#include 
#include 
using namespace std;

admittance_matrix::admittance_matrix（ListBus *bListBranch *l）
{
	float **G**B;
	N_bus=b->get_num（）;								
	N_branch=l->get_num（）;							
	int ij;
	NodeBranch *p;
	float y_ij[2];
	//动态创建二维数组
	G=new float *[N_bus];
	B=new float *[N_bus];

	for（i=0;i		G[i]=new float[N_bus];
		B[i]=new float[N_bus];
	}
	//初始化
	for（i=0;i		for（j=0;j			G[i][j]=0;
			B[i][j]=0;
		}
	//自导纳 互导纳同时生成
	p=l->head->link;
	while（p）{
		i=p->i;
		j=p->j;
		l->yij（ijy_ij）;
		B[i-1][i-1]+=p->B/2;
		B[j-1][j-1]+=p->B/2;
		G[i-1][i-1]+=y_ij[0];
		G[j-1][j-1]+=y_ij[0];
		B[i-1][i-1]+=y_ij[1];
		B[j-1][j-1]+=y_ij[1];
		if（y_ij[0]!=0）G[j-1][i-1]=G[i-1][j-1]=-y_ij[0];
		if（y_ij[1]!=0）B[j-1][i-1]=B[i-1][j-1]=-y_ij[1];
		if（p->Type==2）{									//非标准变比修正
			G[i-1][i-1]+=（1/（p->K*p->K） -1）*y_ij[0];
			B[i-1][i-1]+=（1/（p->K*p->K） -1）*y_ij[1];
			G[i-1][j-1]-=（（1/p->K）-1）*y_ij[0];
			G[j-1][i-1]=G[i-1][j-1];
			B[i-1][j-1]-=（（1/p->K）-1）*y_ij[1];
			B[j-1][i-1]=B[i-1][j-1];
		}
		p=p->link;
	}

	//稀疏处理
	int k=1;
	GII=new float [N_bus];
	BII=new float [N_bus];
	for（i=0;i	{			
		GII[i]=G[i][i];
		BII[i]=B[i][i];
	}

	table=new float *[N_branch+1];						//开辟N_branch+1个空间
	for（i=0;i	table[0][0]=table[0][1]=float（N_branch）;			//存储右上角非零元素个数，即不接地支路数
	table[0][2]=float（N_bus）;							//存储电网节点个数，即矩阵大小为N_bus*N_bus阶
	pos=new int [N_bus];					
	num=new int [N_bus];					
	
	pos[0]=1;
	for（i=0;i		num[i]=0;
		for（j=i+1;j			if（G[i][j]!=0||B[i][j]!=0）{
				num[i]++;
				table[k][2]=float（j）;
				table[k][0]=G[i][j];
				table[k][1]=B[i][j];
				k++;
			}
		}
		if（i	}

	for（i=0;i	{
		delete [] G[i];
		delete [] B[i];
	}
	delete []G;
	delete []B;
}

void admittance_matrix::sparse_print（）
{
	int i;
	cout<<“辅助向量pos:“<	for（i=0;i	cout<<‘\n‘<	cout<<“辅助向量num:“<	for（i=0;i	cout<<‘\n‘<	cout<<“对角元实部:“<	for（i=0;i	cout<<‘\n‘<	cout<<“对角元虚部:“<	for（i=0;i	cout<<‘\n‘<	cout<<“非对角元:“<	cout<	for（i=0;i		cout<			if（i==0） cout<<‘\t‘<<“//此行的三个值分别代表:非零电导、纳值的个数以及导纳矩阵的阶数“;
			cout<	}
}

void admittance_matrix::sparse_print_to_file（）{
	int i;
	ofstream out（“sparse_matrix.txt“）;
	out<<“辅助向量pos:“<	for（i=0;i	out<<‘\n‘<	out<<“辅助向量num:“<	for（i=0;i	out<<‘\n

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       5020  2009-05-07 22:05  admittance_matrix.cpp

     文件       7277  2009-05-08 20:05  bus_branch.cpp

     文件       9511  2009-05-08 18:52  Newton_Rufusion.cpp

     文件        233  2009-05-05 19:10  stdafx.h

     文件        826  2009-05-07 22:39  admittance_matrix.h

     文件       2217  2009-05-08 20:03  bus_branch.h

     文件       1470  2009-05-08 16:34  Newton_Rufusion.h

     文件        541  2009-05-08 20:30  关于输入文件bus.txt和branch.txt的说明.txt

     文件         67  2009-05-05 22:30  branch.txt

     文件       1935  2009-05-07 20:25  Complex.cpp

     文件       3002  2009-05-08 20:18  Project.cpp

     文件         73  2009-05-06 22:07  bus.txt

     文件        956  2009-05-07 20:25  Complex.h

     文件        212  2009-05-05 19:10  stdafx.cpp

----------- ---------  ---------- -----  ----

                33340                    14