b 0.08664 RTc / Pc 0.08664 8.314 425 / 3.8 106 80.56 10 6
v=Z*R*T/P 用 MATLAB 编程实现，首先建立函数 myfun1 和函数 myfun2 function ff=myfun1(Z); T=460;P=1.52e+6; %R-K 方法 Tc=425;Pc=3.8e+6;w=0.193;R=8.314; b=0.08664*R*Tc/Pc; a=0.42747*R^2*Tc^2.5/Pc; v=Z*R*T/P; righteq=(R*T)/(v-b)-a/(T^0.5*v*(v+b)) ff=abs(P/righteq-1); function xishu=myfun2(Z) T=460;P=1.52e+6; %R-K 方法 Tc=425;Pc=3.8e+6;w=0.193;R=8.314; b=0.08664*R*Tc/Pc; a=0.42747*R^2*Tc^2.5/Pc; v=Z*R*T/P; xishu=exp(log(v/(v-b))+a*log(v/(v+b))/(T^0.5*b*R*T)+Z-1-log(Z)) f=P*xishu disp(['压缩因子：',num2str(Z)]); disp(['正丁烷的摩尔体积：',num2str(v)]); disp(['正丁烷的逸度系数：',num2str(xishu)]); disp(['正丁烷的逸度：',num2str(f)]);
例：求 C4 H10 在 460K，1.520MPa 下的逸度 方法 1. I.g 解： R-K 法：
Tc 425K
2. R-K
Pc 3.80MPa
3. 普/23.8 10 2.8955 10
a 0.42748R 2Tc 2 . /5 Pc 0.42748 8.314 2425
输入命令开始迭代 Z=0.5;ff=myfun1(Z); while ff>0.001 Z=Z+0.0001; if Z<1 ff=myfun1(Z) myfun2(Z) end end 结果为：压缩因子：0.8835 正丁烷的摩尔体积：0.002223 正丁烷的逸度系数：0.89431 正丁烷的逸度：1359345.5333 2．普遍化关联法 T=460;P=1.52e+6; %普遍化关联法 Tc=425;Pc=3.8e+6;w=0.193;R=8.314; Tr=T/Tc;Pr=P/Pc; BO=0.083-0.422/(Tr^1.6); B1=0.139-0.172/(Tr^4.2); B=BO+w*B1; xishu=exp(B*Pr/Tr); f=P*xishu; disp(['正丁烷的逸度系数：',num2str(xishu)]); disp(['正丁烷的逸度：',num2str(f)]); 正丁烷的逸度系数：0.89977 正丁烷的逸度：1367646.3748