matlab-物料衡算.docVIP

《工程应用数学》作业5 材研1401 642080502011 黄惠洪 题目：1.用不同温度下固体铅的定压热容数据，求298K时的绝对熵。 2. 第一题：实验数据的数值积分问题 用不同温度下固体铅的定压热容数据，求298K时的绝对熵。数据见下表： 由热力学方法可知，绝对熵与定压热容有如下关系： 本题可用多种方法来解决，如Simpson数值积分法，或采用的形式，或先用数据点拟合方程得回归参数，然后对回归方程进行积分（回归方程可用，，也可在此基础上探索更好的回归方程形式）。 1．方法与思路 采用变型公式：求解给定区间数值积分问题。 2. 过程与结果 1）积分区间的确定： 根据网上教程，matlab对于(-∞, +∞)的定义是x=-inf: inf，但是正常情况下是无需这样定义的，因系统默认变量的取值是从负无穷到正无穷的。 2 ) 编写程序Entropy_Int1.m并保存： 这里将最小二乘B样条拟合法和三次样条插值函数计算作比较，命令行程序如下： function Entropy_Int clear all, clc % 读入数据 T = [5. 10. 15. 20. 25. 30. 50. 70. 100. 150. 200. 250. 298.]; Cp = [0.305 2.8 7.0 10.8 14.1 16.5 21.4 23.3 24.5 25.4 25.8 26.2 26.5]; lnT=log(T); % 用最小二乘B样条拟合法和三次样条插值函数计算 lnTi = linspace(lnT(1),lnT(end),299); knots = 3; K = 3; % 三次B样条 sp = spap2(knots,K,lnT,Cp); cs = csapi(lnT,Cp); % 生成三次样条插值函数cs % 绘制浓度拟合曲线 Ti = linspace(lnT(1),lnT(end),299); Cp_B = fnval(sp,Ti); Cp_c = fnval(cs,Ti); plot(lnT,Cp,ro,Ti,Cp_B,b-,Ti,Cp_c,Bl-.) xlabel(LnT) ylabel(Cp) legend(实验值,B样条拟合,立方样条拟合) % 进行数值积分 pp1= fnint(sp); s1= fnval(pp1,lnT) pp2= fnint(cs); s2= fnval(pp2,lnT) endEntropy_Int，结果如下。 s1 = Columns 1 through 11 0 1.1069 3.3382 5.9047 8.5613 11.2145 20.9711 28.7027 37.3706 47.4362 54.7432 Columns 12 through 13 60.5284 65.1646 s2 = Columns 1 through 11 0 0.6604 2.5965 5.1400 7.9193 10.7146 20.5046 28.0510 36.5924 46.7244 54.0893 Columns 12 through 13 59.8897 64.5187 3.结果分析 1）查阅CRC化学和物理手册，可得到固体铅298K时的绝对熵为64.8，与计算结果之间的差别很小。△s1=(s1-s)/s=0.5%，△s2=(s2-s)/s=0.4%。 2）对比B样条拟合和立方样条拟合，可以发现，立方样条拟合和实验值的匹配度更高。最后拟合出的结果也是立方样条更接近绝对熵的理论值。但是B样条拟合可以适当的曾大结点数knots来提高其拟合程度。 3）其中，T=0K时的实验数据由于等于-∞，如果采用-inf带入求解会出现如下警告： Warning: All data points with NaN or Inf as their site will be ignored. 说明-∞将被系统忽略。因此直接将此点舍去来让程序可以正常运行。便可计算出结果。 准确的汽液平衡数据可以用来计算液相活度系数和Gibbs过量自由能。对下表的双组分汽液平衡数据，计算液相活度系数，并对结果进行热力学的一致性检验，假定气相是理想的，即 ，其中 为纯组分的饱和蒸汽压。数据来源：S. Weissman and S. E. Wood, Vapor-Liquid Equilibrium of Benzene-2,2,4-trimethylpentane Mixtures, J. Chem. Phys., vo