高等工程热力学作业.docVIP

高等工程热力学作业（编程） 第三章 实际气体状态方程 第四章 实际气体导出热力学性质与过程 题目： 一、用PR方程计算制冷剂R290、R600a和混合制冷剂R290/R600a：50/50wt%的PVT性质。 二、用PR方程计算制冷剂R290、R600a和混合制冷剂R290/R600a的导出热力学性质焓和熵。 源程序： 1、牛顿迭代法求Z function Z=newton(A,B,Z) err=1e-6; for n=0:1000 f=Z^3-(1-B)*Z^2+Z*(A-2*B-3*B^2)-(A*B-B^2-B^3); Z=Z-f/(3*Z^2-2*(1-B)*Z+(A-2*B-3*B^2)); if(abs(f)err) break end end %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 2、求a、b、Z、v等参数函数 function [v,Z,a,b,beta]=vv(p,T) R=8.31451; N1=[44.096 369.89 4.2512 0.1521]; N2=[58.122 407.81 3.6290 0.1840]; k1=0.37464+1.54226*N1(4)-0.26992*N1(4)^2; alph