对分发,迭代法,牛顿法RK,SRK方程

1、对分法SRK方程#in clude#in clude#define R 8.314float fun（ float A, float B, float Z）float fz;fz=pow（ Z,3）-pow（乙 2）+（ABB*B）*ZA*B; return fz;main（）int i=l;float p0,t0,w,pl,tl,h,Z0,Zl,Z2;float a,b,m,A,B,t2,at,y,yl,y2;a=b=0.0;printf（请输入临界参数:n p0=）；scanf（”化&p0）;printf（请输入临界参数:n t0=）;scanf（”化&t0）;printf(请输入偏心因子

2、:n w=);scanfCf, &w);printfj请输入实际温度:n tl=)；scanf(“化&tl);printf(“请输入实际压力:n pl=n)；scanfCf,&pl);m=0.480+1.574*w-0.176*w*w;t2=tl/t0;at=(l+m*(l-sqrt(t2)*(l+m*(l-sqrt(t2);a=0.42748*R*R*pow(t0,2)/p0*at;b=0.08664*R*t0/p0;A=a*pl/(R*R*pow(tl,2);B=b*pl/(R*tl);Z1=O.O;Z2=3.O;yl=fun(A,B,Zl);y2=fun(A,B,Z2);Zl=%.5fp

3、rintf(第 0 次计算结果 yl=%.6fy2=%.6fZ2=%.5fn,yl,y2,Zl,Z2);doZ0=(Zl+Z2)/2;y=fun(A,B,Z0);if(yo)y2=y;Z2=Z0;if(V0.0001);printfCXn);printf(“共经过d次计算，求得压缩因子为Z=%.5f n,-i,Zl);对分法RK方程1#in clude#in clude#define R 8.314# define tO 126.2# define pO 3394000.0# define w 0.04float fun( float A, float B, float Z)float fz;

4、fz=pow(乙 3)pow(乙 2)+(ABB*B)*ZA*B; return fz;main()int i=l;float pl,tl,Z0,Zl,Z2;float a,b,A,B,t2,y,yl,y2;a=b=0.0;printf(氨的临界压力 p0=%.lfn,p0);printf(M的临界温度 tO=%.lfn/tO);printf(氨的偏心因子 w=%.2fn,w);printf(请输入实际温度:n tl=n)；scanf(n%f,&tl);printfCi#输入实际压力:n pl=)；scanf(n%f,&pl);a=0.42748*R*R*pow(t0,2)/p0;b=0.08

5、664*R*t0/p0;A=a*pl/(R*R*pow(tl,2);B=b*pl/(R*tl);Z1=O.O;Z2=3.O;yl=fun(A,B,Zl);y2=fun(A,B,Z2);Zl=%.5fprintf(第 0 次计算结果 yl=%.6fy2=%.6fZ2=%.5fn,yl,y2,Zl,Z2);doZ0=(Zl+Z2)/2;y=fun(A,B,ZO);if(y0)y2=y;Z2=Z0;if(V0.0001);printfCXn);printf(哄经过d次计算，求得压缩因子为Z=%.5f n,-i,Zl);对分法RK方程2#in clude #in clude#define R8.324

6、# define tO 126.2# define pO 3394000.0it define w 0.04float fun( float A, float B, float Z)float fz;fz=pow( Z,3)-pow(乙 2)+(ABB*B)*ZA*B; return fz;main()int i=l;float pl,tl,h,Z0,Zl,Z2;float a,b,m,A,B,t2,at,y,yl,y2;a=b=0.0;printf(”氨的临界压力 p0=%.lfn,p0); printf(M的临界温度 tO=%.lfnztO); printf(氨的偏心因子 w=%.2fn?

7、w);printf（“请输入实际温度:n tl=）；scanf(n%f,&tl);printf(请输入实际压力:n pl=H)；scanf(n%f,&pl);m=0.480+1.574*w-0.176*w*w;t2=tl/t0;at=(l+m*(l-sqrt(t2)*(l+m*(l-sqrt(t2);a=0.42748*R*R*pow(t0,2)/p0*at;b=0.08664*R*t0/p0;A=a*pl/(R*R*pow(tl,2);B=b*pl/(R*tl);Z1=O.O;Z2=3.O;yl=fun(A,B,Zl);y2=fun(A,B,Z2);Zl=%.5fprintf(第 0 次计算

8、结果 yl=%.6fy2=%.6fZ2=%.5fn,yl,y2,Zl,Z2);doZ0=(Zl+Z2)/2;y=fun(A,B,Z0);if(yo)y2=y;Z2=Z0;if(V0.0001);printfCXn);printf(“共经过d次计算，求得压缩因子为Z=%.5fn,-i,Zl);迭代法RK方程#in clude#in clude# define R 8314float fun( float A, float B, float h)float Z;Z=1.0/(l-h)-A/B*(h/(l+h);return Z;main()int i=l;float pO,tO,w，pl,tlhZ

9、O;float a,bABZa=b=0.0;printf(请输入临界参数:n p0= ”); scanf(%f”,&p0);printf(请输入临界参数:n t0=)； scanf(%f”,&tO);printf(请输入偏心因子:n w=)； scanf(%f”,&w);printf(请输入实际温度:n tl=H)； scanf(%f”,&tl);printf(请输入实际压力:n pl=); scanf(%f”,&pl);a=0.42748*R*R*pow(t0,2.5)/p0;b=0.08664*R*t0/p0;A=a*pl/(R*R*pow(tl,2.5);B=b*pl/(R*tl);pr

10、intf(请输入Z的迭代初值:n Z=“)；scanf(“f”,&Z);dozo=z；h=B/ZO;Z=fun(A,B,h);printfC第4 次迭代 %f %fn”,i,Z,h);i+；while(fabs(Z-Z0)0.0001);printf(” n“);printf(”共经过d次迭代,求得压缩因子为w=%f n”厂i,Z);迭代法SRK方程#in clude#in clude# define R 8.314float fun( float A, float B, float h)float Z;Z=1.0/(l-h)-A/B*(h/(l+h);return Z;main()int i

11、=l;float pO,tO，w，pl,tlhZO;float abmABZt2，ft;a=b=0.0;printf(请输入临界参数:n p0= “)； scanf(”f”,&p0);printf(请输入临界参数:n t0=)； scanf(”f”,&tO);printf(请输入偏心因子:n w=)； scanf(”f”,&w);printf(请输入实际温度:n tl=H); scanf(”f”,&tl);printf(请输入实际压力:n pl=); scanf(”f”,&pl);m=0.480+1.574*w-0.176*w*w;t2=tl/tO;ft=(l+m*(l-sqrt(t2)*(l

12、+m*(l-sqrt(t2); a=0.42748*R*R*pow(t0,2)/p0*ft;b=0.08664*R*t0/p0;A=a*pl/(R*R*pow(tl,2)； B=b*pl/(R*tl);printf(请输入Z的迭代初值:n Z=H)； scanf(%f”,&Z);h=B/Z;printf(第 0 次迭代 %.4f %.5fn,Z,h); dozo=z；Z=fun(A,B,h);h=B/Z;printf(M%d 次迭代 %.4f %.5fnHJzZ,h);i+；while(fabs(Z-Z0)0.0001);printf(n);printfC1共经过d次迭代，求得压缩因子为w=%

13、4fn：7Z);牛顿迭代法RK方程#in clude#in clude# define R 8.314# define tO 126.2# define pO 3394000.0# define w 0.04float fun( float A, float B, float Z)float fz;fz=pow(乙 3)pow( Z,2)+(A-B-B*B)*Z-A*B; return fz;float funl( float A, float B, float Z)float fzl;fzl=3*pow(Z,2)-2*Z+(A-B-B*B);return fzl;main()int i=l;

14、float plAl/ZOl;float abABt2“yl,y2;a=b=0.0;printfC氨的临界压力 pO=%.lfnH/PO); printff氨的临界温度 tO=%.lfn/tO); printf(”氨的偏 41、因子 w=%.2fn,w);printfC请输入实际温度:n tl=); scanf(“f”,&tl);printf(请输入实际压力:n pl=); scanf(“f”,&pl);printf(请输入牛顿迭代初值:n Zl=)； scanf(“f”, &Z1);a=0.42748*R*R*pow(t0,2.5)/p0; b=0.08664*R*t0/p0;A=a*pl/

15、(R*R*pow(tl,2.5); B=b*pl/(R*tl);ZO=Z1;doZ1=ZO;ZO 二 Zlfun(ABZl)/funl(ABZl);printf(”第d 次牛顿迭代 ZO=%f Zl=%fn,i/ZO,Zl)； i+；while( fabs(Zl-ZO)O.OOOl);printfC经过第d 次牛顿迭代/Z=%fn,/-i,ZO);牛顿迭代法SRK方程#in clude#in cludeit define R 8.314# define tO 126.2# define pO 3394000.0# define w 0.04float fun( float A, float B

16、, float Z)float fz;fz=pow(乙 3)pow( Z,2)+(A-B-B*B)*Z-A*B; return fz;float funl( float A, float Bz float Z)float fzl;fzl=3*pow(Z,2)-2*Z+(A-B-B*B); return fzl;main()int i=l;float plAl/ZOl;float abn%Abt2，ft“y：Ly2;a=b=0.0;printffM的临界压力 pO=%.lfn,pO); printffM的临界温度 tO=%.lfn,/tO); printf(M 的偏心因子 w=%.2fn,/w)

17、;printf(请输入实际温度:n tl=H); scanf(”f”,&tl);printf(请输入实际压力:n pl=); scanf(”f”,&pl);printf(请输入牛顿迭代初值:n Zl=)； scanf(”f”,&Z1);m=0.480+1.574*w-0.176*w*w; t2=tl/tO;拿 sqrt(t2)r(l+m未 sqrt(t2); a=0.42748* R*R*pow(tO/2.5)/pO*ft; b=0.08664*R*t0/p0;A=a*pl/(R*R*pow(tlz2-5);B=b*pl/(R*tl);ZO=Z1;doZl=Z0;Z0 二 Zlfun(ABZl

18、)/funl(ABZl);printf(N第4 次牛顿迭代 ZO=%f Z*%fn” jZ0,Zl);i+；while( fabs(Zl-Z0)0.0001);printf(n经过第d 次牛顿迭代/Z=%fn,/-i/ZO);公式法SRK方程#in clude#in clude# define R 8314# define tO 126.2# define pO 3394000.0# define w 0.04# define I (float) 1/3main()int i=l;float pbtlQUMZhg；float a，b，m,A，B，t2，at，q，p;a=b=0.0;printf

19、(氨的临界压力 pO=%.lfn,pO); printf(氨的临界温度 tO=%.lfnH,tO); printf(M 的偏心因子 w=%.2fn,w);printf(请输入实际温度:n tl=H); scanf(”f”,&tl);printf(请输入实际压力:n pl=); scanf(”f”,&pl);m=0.480+1.574*w-0.176*w*w; t2=tl/tO;at=(l+m*(l-sqrt(t2)*(l+m*(l-sqrt(t2); a=0.42748*R*R*pow(t0z2)/p0*at;b=0.08664*R*t0/p0;A=a*pl/(R*R*pow(tlz2)； B=b*pl/(R*tl);p=(3*(A-B-B*B)-l)/3; q=(9*(A-B-B*B)-2-27*A*B)/27;D=sqrt(q*q/4+p*p*p/27);h=D-q/2;g=-D-q/2;U=pow(hJ);V=pow(gJ);Z=U+V;printf(H解出 Z=%fn/Z);公式法SR方程#in clude#in clud