#include?
#include
#include
#define?NUMBER??50?

void?main（）
{
	double?LAMD（int?xint?yint?zfloat?hdouble?q）;
	
	FILE?*fp1*fp2*fp3;	
????int?C[NUMBER][5]A[NUMBER][NUMBER]M[NUMBER][NUMBER];
	int?flagdc=0;
	int?zrijlnmwTkp;//m?為節(jié)點(diǎn)數(shù)?n為管段數(shù)??k為已知壓力節(jié)點(diǎn)號?p為接點(diǎn)壓力
????float?vmdQ[NUMBER][NUMBER]jz;
	?//??黏度??密度??節(jié)點(diǎn)流量???局部阻力
?	double?ldsumq[NUMBER][NUMBER]G[NUMBER][NUMBER]H[NUMBER][NUMBER]N[NUMBER][NUMBER]F[NUMBER][NUMBER];
????//摩阻系數(shù)????管段流量???????????阻抗矩陣???????????A*G??????????????導(dǎo)納矩陣????????????????系數(shù)矩陣??
	double?x[NUMBER]?tempt[NUMBER][NUMBER]W[NUMBER]DP[NUMBER]wc[NUMBER][NUMBER];
???	//????壓力矩陣?????????管段流量暫存數(shù)組???速度數(shù)組???管段壓降
	fp1=fopen（“a1.dat““r“）;
	
	fscanf（fp1“%d?%d?%d?%d?%d?%f?%f?%d?%f“&m&n&w&k&p&v&md&T&jz）;
?????for（i=0;i	fscanf（fp1“%f?“&Q[i][0]）;

	?/*for（i=0;i		for（j=0;j			wc[i][j]=0;*/??//?誤差數(shù)組
?
	for（i=0;i	{		
		for（j=0;j<5;j++）
		fscanf（fp1“%d“&C[i][j]）;//?起點(diǎn)??終點(diǎn)??管長??管徑??管材
	
	}
	fclose（fp1）;
	
	for（i=0;i		for（j=0;j			A[i][j]=0;
		for（i=0;i		{?
			A[C[i][0]-1][i]=1;
			
			A[C[i][1]-1][i]=-1;
		}

/*for（i=0;i		{	for（j=0;j	??????printf（“%d“A[i][j]）;
		??printf（“\n“）;
		}*/
			
		
		
		for（i=0;i			for（j=0;j				G[i][j]=0;	
		
			for?（i=0;i	
		q[i][0]=（3.1416*C[i][3]*C[i][3]*w）*3600/（4*pow（106））;	
				
		???	/*for（i=0;i			{?
			????printf（“%lf“q[i][0]）;
				printf（“\n“）;?
			}*/
			
				for?（i=0;i			?t[i][0]=0;
			
				
again:		for（i=0;i???
	while（fabs（（q[i][0]-t[i][0]）/q[i][0]）>0.01）	
			
	{
	??????for?（i=0;i		???t[i][0]=q[i][0];

			for?（i=0;i			{	
		????????ld=LAMD（?C[i][3]?C[i][4]wvt[i][0]）;
			???G[i][i]=1/（（6.26*（1e+7））*ld*（1+jz）*t[i][0]*md*T*C[i][2]/（pow（C[i][3]5）*273））;	

	???????????}	
		for?（i=0;i			{
			?for（l=0;l			?{
				{?sum=0;
				for（j=0;j					sum=sum+A[i][j]*G[j][l];
				}
				H[i][l]=sum;
			?}
			
			}
	?			
	??	/*for?（i=0;i		{
			sum=0;
			for（j=0;j				sum=sum+A[i][j]*t[j][0];
			
			Q[i][0]=sum;
		}*/
	

??/*	for（j=0;j	{
		printf（“%lf“-
			Q[j][0]）;
		??printf（“\n“）;
		
	}*/
		for（i=0;i			for（j=0;j				M[i][j]=A[j][i];
			
	??/*	for（i=0;i			{
				for（j=0;j			
	????????	printf（“%5d“M[i][j]）;
	???????	????printf（“\n“）;
			}??*/

			
			for?（i=0;i			{
				for（l=0;l				{
					{?sum=0;
					for（j=0;j						sum=sum+H[i][j]*M[j][l];
					}
					N[i][l]=sum;
				}
				
			}
		?
????	?	for（i=0;i			{
				for（j=0;j			
	????????	printf（“%10.3lf“N[i][j]）;
	???????	????printf（“\n“）;
			}


	//改寫導(dǎo)納矩陣//	?
??????for（i=0;i				N[i][k-1]=-Q[i][0]-N[i][k-1]*p;
	?
			/*for（i=0;i			{
				for（j=0;j			
	????????	printf（“%lf“N[i][

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

?????文件???????7401??2009-06-06?19:31??低壓燃?xì)夤芫W(wǎng)計(jì)算\Cpp1.cpp

?????文件????????134??2009-04-20?14:46??低壓燃?xì)夤芫W(wǎng)計(jì)算\a1.dat

?????目錄??????????0??2009-06-16?17:06??低壓燃?xì)夤芫W(wǎng)計(jì)算

-----------?---------??----------?-----??----

?????????????????7535????????????????????3