现代化工计算技术 曲线拟合 1_3 曲线拟合.ppt

* 单元1 数据处理 单元1 数据处理 * 单元1 数据处理 * 任务3 曲线拟合 曲线拟合 关联函数的选择和线性化 线性最小二乘法 采用Excel曲线拟合 采用自编软件曲线拟合 *曲线拟合效果分析 * 单元1 数据处理 * 任务3 曲线拟合 ☆将离散的数据描述成数学表达式的方法称曲线拟合，或者说经验建模。 一、关联函数的选择和线性化 实测数据关联成数学模型的方法一般有以下几种: ①具有一定的理论依据,可直接根据机理选择关联函数的形式。 ②尚无任何理论依据,但已有一些经验公式可选择。 ☆ 热容的经验公式(附表1-6):cp=a+bT+cT2+dT3 ☆反应动力学方程: ☆反应速率常数方程(Arrhenius方程): * 单元1 数据处理 * 任务3 曲线拟合 ☆将离散的数据描述成数学表达式的方法称曲线拟合，或者说经验建模。 一、关联函数的选择和线性化 实测数据关联成数学模型的方法一般有以下几种: ①具有一定的理论依据,可直接根据机理选择关联函数的形式。 ②尚无任何理论依据,但已有一些经验公式可选择。 ③没有任何经验可循的情况 对于此类情况，通常只能将实验数据画出图形与已知函数图形进行比较，选择图形接近的函数形式作拟合模型。 ☆对于一些非线性模型，应事先将其变换成线性形式，即线性化处理 ; ☆表1-15列出了化工中常用的几种函数类型及线性化的方法。此表中所列均为单变量问题，经线性化处理后的线性模型均可统一用式(1-4)表示,即: Y=A+BX （1-4） * 单元1 数据处理 * 任务3 曲线拟合 ☆表1-15 常用函数线性化方法 * 单元1 数据处理 * 任务3 曲线拟合 一、关联函数的选择和线性化 ☆常用函数线性化举例 幂函数:y=axb 两边取对数:lgy=lga+blgx 令: Y=lgy , X=lgx ,A=lga , B=b 则 Y=A+BX 指数函数:y=aebx 两边取对数:lny=lna+bx 令: Y=lny, X=x,A=lna , B=b 则: Y=A+BX 返回 * 单元1 数据处理 * 任务3 曲线拟合 二、线性最小二乘法 (一)一元线性模型中待定系数的确定方法 ☆弹簧在弹性限度内荷重与长度之间的关系遵循胡克定律 y ?=a+bx (1-8) ☆实验测定弹簧在弹性限度内荷重与长度之间的数据 荷重xi(kg) 长度yi(cm) 0 30.00 2 31.25 4 32.58 6 33.71 8 35.01 10 36.20 12 37.31 14 38.79 16 40.04 由实验数据作图可得一组直线,如图1-26 ☆要解决的问题:如何使拟合的函数式(1-8)计算值与实验值一致! * 单元1 数据处理 * 任务3 曲线拟合 二、线性最小二乘法 (一)一元线性模型中待定系数的确定方法 ☆弹簧在弹性限度内荷重与长度之间的关系遵循胡克定律 y ?=a+bx (1-8) 荷重xi(kg) 0 2 4 6 8 10 12 14 16 长度yi(cm) 30.00 31.25 32.58 33.71 35.01 36.20 37.31 38.79 40.04 由实验数据作图可得一组直线,如图1-26 ☆要解决的问题:如何使拟合的函数式(1-8)计算值与实验值一致!即对任一点i而言,使:?i =yi－ yi ?= yi－( a+bxi)称“残差”最小! ☆最佳效果:每点?i =0,不可能, ??i =0? 可能出现“+”,“－”相抵,解决措施? ??i?=0,但数学上不太好处理,解决方法? ☆最小二乘法:就是采用拟合函数计算值与原实验值之差(残差)的平方和最小来确定拟合函数中各参数的方法称~! * 单元1 数据处理 * 任务3 曲线拟合 二、线性最小二乘法 思考:一元函数如何求其极值 如:现有一长度为Lm的绳子,请设计一种最佳的围法使其绕一固定的墙壁围得最大的面积？ 如：L=8m时，Smax=8m2 思考：上述面积函数的极值属于什么类型极值问题？ 一元函数求极值问题!方法:求导！ * 单元1 数据处理 * 任务3 曲线拟合 二、线性最小二乘法 (一)一元线性模型中待定系数的确定方法 二元函数求极值! * 单元1 数据处理 * 任务3 曲线拟合 二、线性最小二乘法 (一)一元线性模型中待定系数的确定方法 称xi的离差,平方和 称yi的离差,平方和 两离差乘积 (二)一元线性模型中待定系数的计算举例(例1-6表1-16数据) * 单元1 数据处理 * 任务3 曲线拟合 二、线性最小二乘法 (二)一元线性模型中待定系数的计算举例(表1-10数据) * 单元1 数据处理 *