实验数据的数值积分问题.docVIP

实验数据的数值积分问题 用不同温度下固体铅的定压热容数据，求298K时的绝对熵。数据见下表： 由热力学方法可知，绝对熵与定压热容有如下关系： 采用MATLAB来计算，这里将最小二乘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]; CpT = Cp./T; T = [0 T]; CpT = [0 CpT]; % 用最小二乘B样条拟合法和三次样条插值函数计算 knots = 3; K = 3; % 三次B样条 sp = spap2(knots,K,T,CpT); cs = csapi(T,CpT); % 生成三次样条插值函数cs % 绘制浓度拟合曲线 Ti = linspace(T(1),T(end),299); CpT_B = fnval(sp,Ti); CpT_c = fnval(cs,Ti); plot(T,CpT,ro,Ti,CpT_B,b-,Ti,CpT_c,Bl-.) xlabel(T) ylabel(Cp/T) legend(实验值,B样条拟合,立方样条拟合) % 进行数值积分 pp = fnint(cs); s = fnval(pp,T) s = 0 0.0433 0.8733 2.7820 5.3352 8.1104 10.9071 20.7236 28.2599 36.8085 46.9207 54.2919 60.0880 64.7233 在解决本问题时，需首先要考虑积分下限，即时的情形，因此时的为0/0。为此，可考虑先求出接近0K时的数值。可有多种方法此数值，如通过图的样条插值曲线读出，或采用回归拟合以及插值计算的方法。这里采用上述方法得到的数值是，。由于本题的温度数据间隔不等，直接采用样条插值积分是一较好的方法。运用IMSL数学库中的样条插值CSINT和依据样条插值系数的积分CSITG，可以求得绝对熵，下面给出调用的主程序和计算结果。 C To obtain the absolute entropy of solid lead at 298 K from C the experimental data of the heat capacity at constant pressure CP of C solid lead at various temperatures up to and including 298 K C by the cubic spline interpolant and integration from IMSL library. C the absolute entropy in CRC Handbook of Chemistry and Physics is 64.8 J K-1 mol-1. INTEGER NDATA, I, NINTV, NOUT PARAMETER (NDATA=14) REAL A, B, BREAK(NDATA), CSCOEF(4,NDATA), CSITG, ERROR, EXACT, F, FDATA(NDATA), FI, FLOAT, VALUE, X, XDATA(NDATA) INTRINSIC FLOAT EXTERNAL CSINT, CSITG, UMACH DATA (XDATA(I),I=1,NDATA)/0.05,5.,10.,15.,20.,25.,30.,50., 70.,100.,150.,200.,250.,298./ DATA (FDATA(I),I=1,NDATA)/0.003,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/ C Compute cubic spline interpolant DO 10 I=1,NDATA FDATA(I)=FDATA(I)/XDATA(I) 10 CONTINUE C cubic spline interpolant CALL CSINT (NDATA, XDATA, FDATA, BREAK, CSCOEF) C cubic spline integra