c++电力系统潮流程序

#include      /*使用输入/输出函数*/
#include       /*使用setw函数*/
#include          /*使用数学函数*/
#include 
#include 
#include 
#include 

void main（）
{
	double **Zarray;					
	double **G;
	double **B;
cout<<“电研104 武星华 电力系统潮流计算“<<“\n“<										/*定义数组Zarray[i][j]为网络基本参数矩阵*/
										 
										 
	int ij;
    int N=0;
	int L=0;							/*定义N-节点个数，L-支路个数*/

	cout<<“请输入节点个数N:“<	cin>>N;

	cout<<“请输入支路个数L:“<	cin>>L;
cout<<“\n“<
 
	if （（Zarray = new double*[L]） ==NULL）
	{
 		cout << “Could not allocate. “;

 		exit（1）;
 	}
	for（j=0; j<6; j++）
 	{
 		if （（Zarray[j]= new double[6]） == NULL ）
		{
			cout << “Could not allocate. “;
			exit（1）;
		}
 	}	
	
	
    for（i=0; i 	{
		for（j=0; j<6; j++）
		{
			Zarray[i][j]=0;
		}
	}								/*将 Zarray[i][j]分配空间并且初始化*/

//////////////////////////////////////////////////////////////////    /*读文件形成节点导纳距阵*/
	i=0;
	int a1=0 ;
	int	a2=0 ;
	double a3=0 ;
	double a4=0 ;
	double a5=0 ;
	double a6=0 ;
 	ifstream fin（“网络基本参数矩阵.txt“）;
	while （!fin.eof（））
 	{
		fin >> a1 >> a2 >> a3 >> a4 >> a5 >> a6 ;
		
		Zarray[i][0]=a1;
		Zarray[i][1]=a2;
		Zarray[i][2]=a3;
		Zarray[i][3]=a4;
		Zarray[i][4]=a5;
		Zarray[i][5]=a6;

		i++;
		
 	}
    cout << endl ;

	cout << “网络基本参数矩阵:“ << endl;

	for（ i=0; i	{
		for （j=0; j<6; j++）
		{	
			cout<< setw（10） << Zarray[i][j];
		}
		cout << endl;
	}
	cout<<“\n“<	

	if （（G= new double*[N]） ==NULL）
 	{
 	
 		cout << “Could not allocate. “;
 		exit（1）;
 	}
 	for（j=0; j 	{
 		if （（G[j]= new double[N]） ==NULL）
 		{
 			cout << “Could not allocate. “;
 			exit（1）;
 		}
 	}

    for （i=0; i	{
		for （j=0; j		{
			G[i][j]=0;
		}
	}

	if （（B= new double*[N]） ==NULL）
 	{
 	
 		cout << “Could not allocate. “;
 		exit（1）;
 	}
 	for（j=0; j 	{
 		if （（B[j]= new double[N]） ==NULL）
 		{
 			cout << “Could not allocate. “;
 			exit（1）;
 		}
 	}

    for （i=0; i	{
		for （j=0; j		{
			B[i][j]=0;
		}
	}
	int JK;
	double rkxb;
	for（i=0;i	{
		J=int （Zarray[i][0]）;
		K=int （Zarray[i][1]）;
		r=Zarray[i][2];
		x=Zarray[i][3];
		b=Zarray[i][4];
		k=Zarray[i][5];

		if（k!=0）
		{
			G[J-1][J-1]=G[J-1][J-1];
			B[J-1][J-1]=B[J-1][J-1]-1/（x*k*k）;

			G[K-1][K-1]=G[K-1][K-1];
			B[K-1][K-1]=B[K-1][K-1]-1/x;

			G[J-1][K-1]=G[J-1][K-1];
			B[J-1][K-1]=B[J-1][K-1]+1/（x*k）;
			

			G[K-1][J-1]=G[K-1][J-1];
			B[K-1][J-1]=B[K-1][J-1]+1/（x*k）;

		}
		else
		{
			G[J-1][J-1]=G[J-1][J-1]+r/（r*r+x*x）;
			B[J-1][J-1]=B[J-1][J-1]+b-x/（r*r+x*x）;

			G[K-1][K-1]=G[K-1][K-1]+r/（r*r+x*x）;
			B[K-1][K-1]=B[K-1][K-1]+b-x/（r*r+x*x）;

			G[J-1][K-1]=G[J-1][K-1]-r/（r*r+x*x）;
			B[J-1][K-1]=B[J-1][K-1]+x/（r*r+x*x）;
			

			G[K-1][J-1]=G[K-1][J-1]-r/（r*r+x*x）;
			B