#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];
	//動(dòng)態(tài)創(chuàng)建二維數(shù)組
	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;
		}
	//自導(dǎo)納?互導(dǎo)納同時(shí)生成
	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）{									//非標(biāo)準(zhǔn)變比修正
			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];						//開(kāi)辟N_branch+1個(gè)空間
	for（i=0;i	table[0][0]=table[0][1]=float（N_branch）;			//存儲(chǔ)右上角非零元素個(gè)數(shù)，即不接地支路數(shù)
	table[0][2]=float（N_bus）;							//存儲(chǔ)電網(wǎng)節(jié)點(diǎn)個(gè)數(shù)，即矩陣大小為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<<“對(duì)角元實(shí)部:“<	for（i=0;i	cout<<‘\n‘<	cout<<“對(duì)角元虛部:“<	for（i=0;i	cout<<‘\n‘<	cout<<“非對(duì)角元:“<	cout<	for（i=0;i		cout<			if（i==0）?cout<<‘\t‘<<“//此行的三個(gè)值分別代表:非零電導(dǎo)、納值的個(gè)數(shù)以及導(dǎo)納矩陣的階數(shù)“;
			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

?屬性????????????大小?????日期????時(shí)間???名稱(chēng)
-----------?---------??----------?-----??----

?????文件???????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??關(guān)于輸入文件bus.txt和branch.txt的說(shuō)明.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