对分发迭代法牛顿法RKSRK方程docx.docx

对分发迭代法牛顿法RKSRK方程docx.docx

《对分发迭代法牛顿法RKSRK方程docx.docx》由会员分享，可在线阅读，更多相关《对分发迭代法牛顿法RKSRK方程docx.docx（14页珍藏版）》请在冰豆网上搜索。

对分发迭代法牛顿法RKSRK方程docx

对分法SRK方程

#include

#include

#defineR8.314

floatfun（floatA,floatB,floatZ）

{

floatfz;

fz=pow（Z,3）-pow（乙2）+（A・B・B*B）*Z・A*B;returnfz;

}

main（）

{

inti=l;

floatp0,t0,w,pl,tl,h,Z0,Zl,Z2;

floata,b,m,A,B,t2,at,y,yl,y2;

a=b=0.0;

printf（"请输入临界参数:

\np0="）；

scanf（”％化&p0）;

printf（"请输入临界参数:

\nt0="）;

scanf（”％化&t0）;

printf（"请输入偏心因子:

\nw="）;

scanfC^f",&w）;

printfj请输入实际温度:

\ntl="）；

scanf（“％化&tl）;

printf（“请输入实际压力:

\npl=n）；

scanfC^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=%.5f

printf（"第0次计算结果yl=%.6fy2=%.6f

Z2=%.5f\n",yl,y2,Zl,Z2）;

do

Z0=（Zl+Z2）/2;

y=fun（A,B,Z0）;

if（y>o）

{y2=y;Z2=Z0;}

if（V<0）

{yl=y;Zl=ZO;}

if（y=o）

{Zl=ZO;break;}

printf（"第％・24次计算结果y2%・6fy2=%.6fZl=%.5f

Z2=%.5f\n“jyby2,ZbZ2）;

i++；

}

while（（fabs（Zl-Z2））>0.0001）;

printfCXn"）;

printf（“共经过％d次计算，求得压缩因子为Z=%.5f\n",-i,Zl）;

对分法RK方程1

#include

#include

#defineR8.314

#definetO126.2

#definepO3394000.0

#definew0.04

floatfun（floatA,floatB,floatZ）

{

floatfz;

fz=pow（乙3）・pow（乙2）+（A・B・B*B）*Z・A*B;returnfz;

}

main（）

{

inti=l;

floatpl,tl,Z0,Zl,Z2;

floata,b,A,B,t2,y,yl,y2;

a=b=0.0;

printf（"氨的临界压力p0=%.lf\n",p0）;

printf（"M的临界温度tO=%.lf\n"/tO）;

printf（"氨的偏心因子w=%.2f\n",w）;

printf（"请输入实际温度:

\ntl=n）；

scanf（n%f",&tl）;

printfC'i#输入实际压力:

\npl="）；

scanf（n%f",&pl）;

a=0.42748*R*R*pow（t0,2）/p0;

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

printf（"第0次计算结果yl=%.6fy2=%.6f

Z2=%.5f\n",yl,y2,Zl,Z2）;

do

Z0=（Zl+Z2）/2;

y=fun（A,B,ZO）;

if（y>0）

{y2=y;Z2=Z0;}

if（V<0）

{yl=y;Zl=ZO;}

if（y=o）

{Zl=ZO;break;}

printf（"第％・201次计算结果y2%.6fy2=%.6fZl=%.5f

Z2=%.5f\n“jyby2,ZbZ2）;

i++；

}

while（（fabs（Zl-Z2））>0.0001）;

printfCXn"）;

printf（哄经过％d次计算，求得压缩因子为Z=%.5f\n",-i,Zl）;

}

对分法RK方程2

#include#include

#defineR8.324

#definetO126.2

#definepO3394000.0

itdefinew0.04

floatfun（floatA,floatB,floatZ）

{

floatfz;

fz=pow（Z,3）-pow（乙2）+（A・B・B*B）*Z・A*B;returnfz;

}

main（）

{

inti=l;

floatpl,tl,h,Z0,Zl,Z2;

floata,b,m,A,B,t2,at,y,yl,y2;

a=b=0.0;

printf（”氨的临界压力p0=%.lf\n",p0）;printf（"M的临界温度tO=%.lf\n"ztO）;printf（"氨的偏心因子w=%.2f\n"?

w）;

printf（“请输入实际温度:

\ntl="）；

scanf（n%f",&tl）;

printf（"请输入实际压力:

\npl=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=%.5f

printf（"第0次计算结果yl=%.6fy2=%.6f

Z2=%.5f\n",yl,y2,Zl,Z2）;

do

{

Z0=（Zl+Z2）/2;

y=fun（A,B,Z0）;

if（y>o）

{y2=y;Z2=Z0;}

if（V<0）

{yl=y;Zl=ZO;}

if（y=O）

{Zl=ZO;break;}

printf（"第％・24次计算结果y2%・6fy2=%.6fZl=%.5f

Z2=%.5f\n“jyby2,Z：

bZ2）;

i卄；

}

while（（fabs（Zl-Z2））>0.0001）;

printfCXn"）;

printf（“共经过％d次计算，求得压缩因子为Z=%.5f\n,,,-i,Zl）;

}

迭代法RK方程

#include

#include

#defineR8314

floatfun（floatA,floatB,floath）

{

floatZ;

Z=1.0/（l-h）-A/B*（h/（l+h））;

returnZ;

main（）

inti=l;

floatpO,tO,w，pl,tlhZO;

floata,bABZ

a=b=0.0;

printf（"请输入临界参数:

\np0=”）;scanf（"%f”,&p0）;

printf（"请输入临界参数:

\nt0="）；scanf（"%f”,&tO）;

printf（"请输入偏心因子:

\nw="）；scanf（"%f”,&w）;

printf（"请输入实际温度:

\ntl=H）；scanf（"%f”,&tl）;

printf（"请输入实际压力:

\npl="）;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）;

printf（"请输入Z的迭代初值:

\nZ=“）；

scanf（“％f”,&Z）;

do

{

zo=z；

h=B/ZO;

Z=fun（A,B,h）;

printfC'第％4次迭代%f%f\n”,i,Z,h）;

i++；

}

while（（fabs（Z-Z0））>0.0001）;

printf（”\n“）;

printf（”共经过％d次迭代,求得压缩因子为w=%f\n”厂i,Z）;

迭代法SRK方程

#include

#include

#defineR8.314

floatfun（floatA,floatB,floath）

{

floatZ;

Z=1.0/（l-h）-A/B*（h/（l+h））;

returnZ;

}

main（）

{

inti=l;

floatpO,tO，w，pl,tlhZO;

floatabmABZt2，ft;

a=b=0.0;

printf（"请输入临界参数:

\np0=“）；scanf（”％f”,&p0）;

printf（"请输入临界参数:

\nt0="）；scanf（”％f”,&tO）;

printf（"请输入偏心因子:

\nw="）；scanf（”％f”,&w）;

printf（"请输入实际温度:

\ntl=H）;scanf（”％f”,&tl）;

printf（"请输入实际压力:

\npl="）;scanf（”％f”,&pl）;

m=0.480+1.574*w-0.176*w*w;

t2=tl/tO;

ft=（l+m*（l-sqrt（t2）））*（l+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的迭代初值:

\nZ=H）；scanf（"%f”,&Z）;

h=B/Z;

printf（"第0次迭代%.4f%.5f\n",Z,h）;do

zo=z；

Z=fun（A,B,h）;

h=B/Z;

printf（M^%d次迭代%.4f%.5f\nHJzZ,h）;

i++；

}

while（（fabs（Z-Z0））>0.0001）;

printf（"\n"）;

printfC1共经过％d次迭代，求得压缩因子为w=%・4f\n：

7Z）;

}

牛顿迭代法RK方程

#include

#include

#defineR8.314

#definetO126.2

#definepO3394000.0

#definew0.04

floatfun（floatA,floatB,floatZ）

{

floatfz;

fz=pow（乙3）・pow（Z,2）+（A-B-B*B）*Z-A*B;returnfz;

}

floatfunl（floatA,floatB,floatZ）

{

floatfzl;

fzl=3*pow（Z,2）-2*Z+（A-B-B*B）;

returnfzl;

}

main（）

{

inti=l;

floatplAl/ZO^l;

floatabABt2“yl,y2;

a=b=0.0;

printfC'氨的临界压力pO=%.lf\nH/PO）;printff"氨的临界温度tO=%.lf\n"/tO）;printf（”氨的偏41、因子w=%.2f\n,,,w）;

printfC'请输入实际温度:

\ntl="）;scanf（“％f”,&tl）;

printf（"请输入实际压力:

\npl="）;scanf（“％f”,&pl）;

printf（"请输入牛顿迭代初值:

\nZl="）；scanf（“％f”,&Z1）;

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）;

ZO=Z1;

do

{

Z1=ZO;

ZO二Zl・fun（ABZl）/funl（ABZl）;

printf（”第％d次牛顿迭代ZO=%fZl=%f\n,,,i/ZO,Zl）；i++；

}

while（fabs（Zl-ZO）>O.OOOl）;

printfC'经过第％d次牛顿迭代/Z=%f\n,,/-i,ZO）;

}

牛顿迭代法SRK方程

#include

#include

itdefineR8.314

#definetO126.2

#definepO3394000.0

#definew0.04

floatfun（floatA,floatB,floatZ）

floatfz;

fz=pow（乙3）・pow（Z,2）+（A-B-B*B）*Z-A*B;returnfz;

}

floatfunl（floatA,floatBzfloatZ）

{

floatfzl;

fzl=3*pow（Z,2）-2*Z+（A-B-B*B）;returnfzl;

}

main（）

{

inti=l;

floatplAl/ZO^l;

floatabn%Abt2，ft“y：

Ly2;

a=b=0.0;

printff'M的临界压力pO=%.lf\n",pO）;printff'M的临界温度tO=%.lf\n,,/tO）;printf（"M的偏心因子w=%.2f\n,,/w）;

printf（"请输入实际温度:

\ntl=H）;scanf（”％f”,&tl）;

printf（"请输入实际压力:

\npl="）;scanf（”％f”,&pl）;

printf（"请输入牛顿迭代初值:

\nZl="）；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;

do

{

Zl=Z0;

Z0二Zl・fun（ABZl）/funl（ABZl）;

printf（N第％4次牛顿迭代ZO=%fZ*%f\n”jZ0,Zl）;

i++；

}

while（fabs（Zl-Z0）>0.0001）;

printf（n经过第％d次牛顿迭代/Z=%f\n,,/-i/ZO）;

}

公式法SRK方程

#include

#include

#defineR8314

#definetO126.2

#definepO3394000.0

#definew0.04

#defineI（float）1/3

main（）

{

inti=l;

floatpbtlQUMZhg；

floata，b，m,A，B，t2，at，q，p;

a=b=0.0;

printf（"氨的临界压力pO=%.lf\n",pO）;printf（"氨的临界温度tO=%.lf\nH,tO）;printf（"M的偏心因子w=%.2f\n",w）;

printf（"请输入实际温度:

\ntl=H）;scanf（”％f”,&tl）;

printf（"请输入实际压力:

\npl="）;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=%f\n"/Z）;

}

公式法SR方程

#include

#include

#defineR8314

#definetO126.2

#definepO3394000.0

#definew0.04

#defineI（float）1/3

main（）

{

inti=l;

floatpbtlQUMZhg；

floata,b,m,A,p;

a=b=0.0;

printf（"氨的临界压力pO=%.lf\n",pO）;printf（"氨的临界温度tO=%.lf\n,,/tO）;

printf（”氨的偏41、因子w=%.2f\nHzw）;

printff"请输入实际温度:

\ntl=H）;scanf（”％f”,&tl）;

printf（'"请输入实际压力:

\npl="）;

scanf（“％化&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）;

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（g,l）;

Z=U+V;

printfC*解出Z=%f\n,,/Z）;

}

三次方程求根公式

设一元三次方程厂。

在复数集中的根是A-2,X3,那么勺-（-手+（（#）2+（铲+（铲，

七=S*升J炸）2+尋+T号馬~+（铲，

-1士向-1-阪

其中曲=一，勺=2°

3ac-2?

2b^-9abc+•27a^d

这里P・PL，孑而