C++油气集输程序

石油人在北京

#include
#include
int main()
{
int i,temp=1,flag=0;
double sum=0,sumC=0,Mpj=0;
double z[8],w[8]={0};
double C[8]= {0,25.9,6.6,12.3,12.4,12.4,16.5,913.9};
int M[8]={0,16,30,44,58,72,86,472};
for(i=1;i<=7;i++)
sum+=C[i];sumC=sum;
for(i=1;i<=7;i++)
{
w[i]=C[i]/sum;Mpj+=w[i]/M[i];
}
Mpj=1/Mpj;
for(i=1;i<=7;i++)
{
z[i]=(Mpj*w[i])/M[i];
}
double sumx=0,Mid,sumz=0;
double y[8],sumy=0;
double WL[8]={0},sumWL=0,Wg[8]={0},sumWg=0;
double x[8]={0};
double K1[8]={0,100,17,5.2,2.4,0.8,0.23,0.018};
double K2[8]={0,255,43,13,4.95,1.65,0.6,0.0375};
for(flag=0;flag<2;flag++)
{
double L1=0,V2=1;
sumx=0;
while(fabs(sumx-1)>1e-6)
{
Mid=(L1+V2)/2;sumx=0;
for(i=1;i<=7;i++)
{
x[i]=z[i]/(Mid+(1-Mid)*K1[i]);sumx+=x[i];
}
if(sumx-1>0)
L1=Mid;
else
V2=Mid;

}
printf("在第%d级分离后，组成为: L=%f, V=%f\n\n",flag+1,Mid,1-Mid);

sumy=0;
for(i=1;i<=7;i++)
{
y[i]=K1[i]*x[i];
sumy+=y[i];
}
double Pg;
sum=0;
for(i=1;i<7;i++)
sum+=y[i]*M[i];
sum=sum+y[7]*114;
Pg=sum/22.4;
printf("在第%d级分离后，气相密度为: %f \n\n",flag+1,Pg);
temp--;
sumWL=0;
sumWg=0;
for(i=1;i<=7;i++)
{
WL[i]=(x[i]*Mid*C[i])/z[i];
sumWL+=WL[i];
Wg[i]=C[i]-WL[i];
sumWg+=Wg[i];
}

printf(" Ci Xi Yi Zi WLi Wgi\n");
for(i=1;i<=6;i++)
{
printf("\n");
printf(" C%d %f %f %f %f %f %f\n",i,C[i],x[i],y[i],z[i],WL[i],Wg[i]);
}
printf("\n");
printf(" C7+%f %f %f %f %f %f\n",C[7],x[7],y[7],z[7],WL[7],Wg[7]);
printf("\n");
printf("合计%f %f %f %f %f %f\n\n",sumC,sumx,sumy,sumz,sumWL,sumWg);

if(flag==0)
{
for(i=1;i<=7;i++)
{
K1[i]=K2[i];z[i]=x[i];sumz+=z[i];
}
}
}
sumy=0;
}

chenh1964

正需要，非常感谢楼主！！

血麒麟

谢谢楼主分享！

zhang.limin

zccrg001

谢谢分享，下载学习