对分发迭代法牛顿法RKSRK方程docx文档格式.docx

1、t0）;请输入偏心因子:n w=scanfCf, &w）;printfj请输入实际温度:n tl=scanf（“化&tl）;printf（“请输入实际压力:n pl=n）；,&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=%.5f第

2、 0 次计算结果 yl=%.6f y2=%.6fZ2=%.5fn,yl,y2,Zl,Z2）;doZ0=（Zl+Z2）/2;y=fun（A,B,Z0）;if（yo）y2=y;Z2=Z0;if（V0.0001）;printfCXnprintf（“共经过d次计算，求得压缩因子为Z=%.5f n,-i,Zl）;对分法RK方程1#define tO 126.2#define pO 3394000.0#define w 0.04fz=pow（乙 3）pow（乙 2）+（ABB*B）*ZA*B;float pl,tl,Z0,Zl,Z2;float a,b,A,B,t2,y,yl,y2;氨的临界压力 p0=%

3、.lfn,p0）;M的临界温度 tO=%.lfn/tO）;氨的偏心因子 w=%.2fn,w）;请输入实际温度:n tl=n）；scanf（n%fprintfCi#输入实际压力:n pl=a=0.42748*R*R*pow（t0,2）/p0;y=fun（A,B,ZO）;第201 次计算结果 y2%.6f y2=%.6f Zl=%.5fprintf（哄经过d次计算，求得压缩因子为Z=%.5f n对分法RK方程2 #in cludeprintf（” n“）;printf（”共经过d次迭代,求得压缩因子为w=%f n”厂i,Z）;迭代法SRK方程# define R 8.314float pO,tO，

4、w，pl,tlhZO;float abmABZt2，ft;n p0= “）； scanf（”f”,&n tl=H）;t2=tl/tO;ft=（l+m*（l-sqrt（t2）*（l+m*（l-sqrt（t2）; a=0.42748*R*R*pow（t0,2）/p0*ft;A=a*pl/（R*R*pow（tl,2）； B=b*pl/（R*tl）;n Z=H）；h=B/Z;第 0 次迭代 %.4f %.5fn,Z,h）; doprintf（M%d 次迭代 %.4f %.5fnHJzZ,h）;nprintfC1共经过d次迭代，求得压缩因子为w=%4fn：7Z）;牛顿迭代法RK方程fz=pow（乙 3）

5、pow（ Z,2）+（A-B-B*B）*Z-A*B;float funl（ float A, float B, float Z）float fzl;fzl=3*pow（Z,2）-2*Z+（A-B-B*B）;return fzl;float plAl/ZOl;float abABt2“yl,y2;氨的临界压力 pO=%.lfnH/PO）; printff氨的临界温度 tO=%.lfn printf（”氨的偏 41、因子 w=%.2fn,w）; scanf（“f”,&请输入牛顿迭代初值:n Zl= scanf（“f”, &Z1）; b=0.08664*R*t0/p0;ZO=Z1;Z1=ZO;ZO

6、二 Zlfun（ABZl）/funl（ABZl）;printf（”第d 次牛顿迭代 ZO=%f Zl=%fn,i/ZO,Zl）； i+；while（ fabs（Zl-ZO）O.OOOl）;经过第d 次牛顿迭代/Z=%fn,/-i,ZO）;牛顿迭代法SRK方程it define R 8.314float funl（ float A, float Bz float Z） return fzl;float abn%Abt2，ft“y：Ly2;printffM的临界压力 pO=%.lfn,pO）; printffM的临界温度 tO=%.lfn,/tO）;M 的偏心因子 w=%.2fn,/w）; t2=

7、tl/tO;拿 sqrt（t2）r（l+m未 sqrt（t2）; a=0.42748* R*R*pow（tO/2.5）/pO*ft;A=a*pl/（R*R*pow（tlz2-5）;Zl=Z0;Z0 二 Zlfun（ABZl）/funl（ABZl）;printf（N第4 次牛顿迭代 ZO=%f Z*%fn” jZ0,Zl）;while（ fabs（Zl-Z0）printf（n经过第d 次牛顿迭代/Z=%fn,/-i/ZO）;公式法SRK方程#define R 8314#define I （float） 1/3float pbtlQUMZhg；float a，b，m,A，B，t2，at，q，p;氨

8、的临界压力 pO=%.lfn氨的临界温度 tO=%.lfnH,tO）;M 的偏心因子 w=%.2fn a=0.42748*R*R*pow（t0z2）/p0*at;A=a*pl/（R*R*pow（tlz2）；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方程float a,b,m,A,p;氨的临界温度 tO=%.lfn,/tO）;printf（”氨的偏41、因子 w=%.2fnHzw）;printffprintf（a=0.42748*R*R*pow（t0/2.5）/p0;A=a*pl/（R*R*pow（tl/2.5）;q=（9*（A-B-B*B）-2-27*A*B）/27;V=pow（g,l）;printfC*解出 Z=%fn,/Z）;三次方程求根公式设一元三次方程厂。在复数集中的根是A-2, X3,那么 勺-（-手+（#）2+（铲 + （铲，七=S *升J炸）2 +尋+ T号馬+（铲，-1士向 -1-阪其中曲=一，勺=2 3ac - 2? 2b - 9abc + 27ad这里PPL，孑 而