第1章实验数据及模型参数分析.pptVIP

第1章 实验数据及模型参数拟合方法 1.1 问题的提出 1.2拟合的标准 1.3单变量拟合和多变量拟合 1.4解矛盾方程组 1.5梯度法拟合参数 1.6吸附等温曲线回归 1.1 问题的提出 化工设计及化工模拟计算中，有大量的物性参数及各种设备参数。实验测量得到的常常是一组离散数据序列（xi ,yi） 图1-1所示为“噪声” 图1-2所示为无法同时满足某特定的函数 1.1 问题的提出 在化学化工中，许多模型也要利用数据拟合技术，求出最佳的模型和模型参数。 如在某一反应工程实验中，我们测得了如表1-1所示的实验数据： 1.1 问题的提出 确定在其他条件不变的情况下，转化率y和温度T的具体关系，现拟用两种模型去拟合实验数据，两种模型分别是： 1.2 拟合的标准 向量Q与Y之间的误差或距离有以下几种定义方法： （1）用各点误差绝对值的和表示 （2）用各点误差按绝对值的最大值表示 （3）用各点误差的平方和表示 1.2 拟合的标准 由于计算均方误差的最小值的原则容易实现而被广泛采用。按均方误差达到极小构造拟合曲线的方法称为最小二乘法。同时还有许多种其他的方法构造拟合曲线，感兴趣的读者可参阅有关教材。本章主要讲述用最小二乘法构造拟合曲线。 1.2 拟合的标准 实例 实验测得二甲醇（DME）的饱和蒸汽压和温度的关系如下表 ： 1.2 拟合的标准 实例 如果以直线拟合p=a+bt,即拟合函数是一条直线。通过计算均方误差Q ( a , b )最小值而确定直线方程（见图1-3） 1.2 拟合的标准 实例 如果采用二次拟合，通过计算下述均方误差： 1.2 拟合的标准 实例 比较图1-3和图1-4以及各自的相关系数和平均绝对偏差可知： 对于DME饱和蒸汽压和温度之间的关系，在实验温度范围内用二次拟合曲线优于线性拟合。 二次拟合曲线具有局限性，由图1-4观察可知，当温度低于-30℃时，饱和压力有升高的趋势，但在拟合的温度范围内，二次拟合的平均绝对偏差又小于一次拟合，故对物性数据进行拟合时，不仅要看在拟合条件下的拟合效果，还必须根据物性的具体性质，判断在拟合条件之外的物性变化趋势，以便使拟合公式在已做实验点数据之外应用。 1.3.1单变量拟合 1.3.2 多变量的曲线拟合 1.3.1 单变量拟合 线性拟合 给定一组数据（xi,yi）,i=1, 2 , …, m ，做拟合直线p (x)=a + bx , 均方误差为 ： 1.3.1 单变量拟合 线性拟合 整理得到拟合曲线满足的方程： 1.3.1 单变量拟合 线性拟合 1.3.1 单变量拟合 线性拟合实例 例1.1：下表为实验测得的某一物性和温度之间的关系数据，表中x为温度数据，y为物性数据。请用线性函数拟合温度和物性之间的关系。 1.3.1 单变量拟合 线性拟合实例 解：设拟合直线 ，并计算得下表： 1.3.1 单变量拟合 二次拟合函数 给定数据序列（xi,yi）,i=1, 2 , …, m ,用二次多项式函数拟合这组数据。 1.3.1 单变量拟合 二次拟合函数 整理上式得二次多项式函数拟合的满足条件方程： 1.3.1 单变量拟合 二次拟合函数 上面是二次拟合基本类型的求解方法，和一次拟合一样，二次拟合也可以有多种变型： 1.3.1 单变量拟合 二次拟合函数 如果我们需要求解是下面的拟合函数： 1.3.1 单变量拟合 二次拟合实例 例1.2：请用二次多项式函数拟合下面这组数据。 1.3.1 单变量拟合 二次拟合实例 解：设 ，由计算得下表： 1.3.1 单变量拟合 二次拟合实例 将上面数据代入式 (1-14) ，相应的法方程为： 1.3.1 单变量拟合 二次拟合实例 拟合曲线的均方误差： 结果见图 1-6。二次曲线的拟合程序可利用后面介绍的单变量n次拟合程序。 1.3.2 多变量的曲线拟合 实际在化工实验数据处理及模型参数拟合时，通常会碰到多变量的参数拟合问题。一个典型的例子是传热实验中努塞尔准数和雷诺及普兰德准数之间的拟合问题： 1.3.2 多变量的曲线拟合 为不失一般性，我们把它表达成以下形式： 给定数据序列 用一次多项式函数拟合这组数据。 设