工程数学作业第十五次方健.docx

姓名：方健 学号：652081701073 专业：化学工程 班级：化工1503 问题：利用Minitab研究二硫化碳饱和蒸汽压模型回归问题。 温度/C蒸汽压/mmHg温度/C蒸汽压/mmHg温度/C蒸汽压/mmHg温度/C蒸汽压/mmHg-701.6-3026.210198.050955.6-603.5-2046.520297.5601170.4-507.1-1078.830432.7701558.0-4014.00127.340616.7二硫化碳的基本性质：临界温度为273.05，临界压力是72.87atm。 分析： 对于实验数据进行分析，在得到模型参数之前，对于数据自变量之间的线性相关性进行检验，是发现回归模型应用的可靠性和准确性受限制的有效方法，因为自变量之间的线性相关性，使得无法区分它们对应变量的作用，回归模型参数时会遇到几乎是奇异的数据矩阵，这样的模型参数具有很大的不确定性。 对于回归模型的选择，回归模型的选择是数据模型化过程中经常遇到的问题，了解掌握一些基于数据统计分析原理得到的方法和原则对回归模型的选择会有很大的帮助。对于预测物质的饱和蒸汽压，其中有三参数的安托尼方程模型拟合。 采用三参数的Antoine（安托尼）方程进行回归分析 设定A，B，C三个参数的初始值为1 方法一 算法 Gauss-Newton 最大迭代次数 200 公差 0.00001 得到拟合曲线如下： 得到了残差图如下： 工作表的计算结果汇总如下： 方法 算法 Gauss-Newton 最大迭代次数 200 公差 0.00001 参数的起始值 参数 值 A 1 B 1 C 1 方程 ln(P) = 8.19011 - 11033.7 / (T/K - 246669) 参数估计 标准误 参数 估计 估计值 A 8 * B -11034 * C -246669 * ln(P) = A + B / (T/K + C) 失拟 没有仿行。 Minitab 无法根据纯误差进行失拟检验。 汇总 迭代 4 最终 SSE 87.2015 DFE 12 MSE 7.26679 S 2.69570 拟合线：ln(P)与T/K ln(P)残差图 结果分析与讨论： 从上面的回归结果可以知道，二硫化碳的饱和蒸汽压随着温度的曲线如图的拟合曲线所示，拟合的方程为ln(P) = 8.19011 - 11033.7 / (T/K - 246669)可以知道A=8.19011，B=-11033.7，C=-2466669，可以看到残差近似服从正态分布，对于二硫化碳的饱和蒸汽压随着温度变化都可以在图中看出来。