实验数据及模型参数(1)ppt课件

1、第第1章章 实验数据及模型参数实验数据及模型参数拟合方法拟合方法n1.1 问题的提出 n1.2拟合的规范 n1.3单变量拟合和多变量拟合 n1.4解矛盾方程组 n1.5梯度法拟合参数 n1.6吸附等温曲线回归 1.1 问题的提出问题的提出n化工设计及化工模拟计算中，有大量的物性参数及各种设备参数。实验丈量得到的经常是一组离散数据序列xi ,yin图1-1所示为“噪声n图1-2所示为无法同时满足某特定的函数-202468101214161820050100150200Y X 024681005101520Y X 图1-1 含有噪声的数据图1-2 无法同时满足某特定函数的数据序列1.1 问题的提出

2、问题的提出n在化学化工中，许多模型也要利用数据拟合技术，求出最正确的模型和模型参数。n如在某一反响工程实验中，我们测得了如表1-1所示的实验数据：表1-11.1 问题的提出问题的提出确定在其他条件不变的情况下，转化率y和温度T的详细关系，现拟用两种模型去拟合实验数据，两种模型分别是：2111TcTbay2222)45(Tbacy (1-2) (1-3) 1.2 拟合的规范拟合的规范 向量Q与Y之间的误差或间隔有以下几种定义方法：1用各点误差绝对值的和表示2用各点误差按绝对值的最大值表示3用各点误差的平方和表示miiiyxR11)(iimiyxR)(max122212 )(Q(x)-YRyxRR

3、imii或(1-4) (1-5) (1-6) R称为均方误差1.2 拟合的规范拟合的规范n由于计算均方误差的最小值的原那么容易实现而被广泛采用。按均方误差到达极小构造拟合曲线的方法称为最小二乘法。同时还有许多种其他的方法构造拟合曲线，感兴趣的读者可参阅有关教材。本章主要讲述用最小二乘法构造拟合曲线。1.2 拟合的规范拟合的规范 实例实例n实验测得二甲醇DME的饱和蒸汽压和温度的关系如下表 ：序号温度 蒸气压 MPa1-23.70.1012-100.174300.2544100.3595200.4956300.6627400.880表1-2 DME饱和蒸气压和温度的关系由表1-2的数据观测可得，

4、DME的饱和蒸汽压和温度有正相关关系。1.2 拟合的规范拟合的规范 实例实例n假设以直线拟合p=a+bt,即拟合函数是一条直线。经过计算均方误差Q ( a , b )最小值而确定直线方程见图1-3 -30-20-10010203040500.00.20.40.60.81.0pt图1-3 DME饱和蒸汽压和温度之间的线性拟合2121)()(),(imiiimiipbtaptpbaQ拟合得到得直线方程为：tp0.01210.30324 相关系数R为0.97296，平均绝对偏向SD为0.05065。 (1-8) (1-7) 1.2 拟合的规范拟合的规范 实例实例假设采用二次拟合，经过计算下述均方误差

5、：21221012210)()(),(imiiimiiiptataaptpaaaQ拟合得二次方程为：2000150009570248450t.t.p1-91-10相关系数为R为0.99972，平均绝对偏向SD为0.0056。详细拟合曲线见图1-4 -30-20-10010203040500.00.20.40.60.81.0y=0.24845+0.00957 x+0.00015 x2压力, P(MPa)温度, t()图1-4 DME饱和蒸汽压和温度之间的二次拟合1.2 拟合的规范拟合的规范 实例实例 比较图1-3和图1-4以及各自的相关系数和平均绝对偏向可知：对于DME饱和蒸汽压和温度之间的关系

6、，在实验温度范围内用二次拟合曲线优于线性拟合。二次拟合曲线具有局限性，由图1-4察看可知，当温度低于-30时，饱和压力有升高的趋势，但在拟合的温度范围内，二次拟合的平均绝对偏向又小于一次拟合，故对物性数据进展拟合时，不仅要看在拟合条件下的拟合效果，还必需根据物性的详细性质，判别在拟合条件之外的物性变化趋势，以便使拟合公式在已做实验点数据之外运用。 1.3 单变量拟合和多变量拟合单变量拟合和多变量拟合n1.3.1单变量拟合n1.3.2 多变量的曲线拟合1.3.1 单变量拟合单变量拟合 线性拟线性拟合合 n给定一组数据xi,yi,i=1, 2 , , m ，做拟合直线p (x)=a + bx ,

7、均方误差为 ：2121)()(),(imiiimiiybxayxpbaQ1-11 Q (a , b)的极小值需满足： 0)(2),(1imiiybxaabaQ0)(2),(1iimiixybxabbaQ1.3.1 单变量拟合单变量拟合 线性拟线性拟合合n整理得到拟合曲线满足的方程：mimimiiiiimimiiiyxbxaxybxma111211)()()(miiimiimiimiimiiyxybaxxxm111211或 (1-12) 称式(1-12)为拟合曲线的法方程。1.3.1 单变量拟合单变量拟合 线性拟线性拟合合)()()(/()( 2112111211211121121112111m

8、imiiimimiimiiiimiimiimiiimiimiimiimiimiimiimiimiiimiimiixxmyxyxmbxxmyxxxyxxxmxyxxya可用消元法或克莱姆方法解出方程：1.3.1 单变量拟合单变量拟合 线性拟合线性拟合实例实例例1.1：下表为实验测得的某一物性和温度之间的关系数据，表中x为温度数据，y为物性数据。请用线性函数拟合温度和物性之间的关系。x131516212223252930313640y111011121213131214161713x42556062647072100130y1422142121241723341.3.1 单变量拟合单变量拟合 线性

9、拟合线性拟合实例实例n解：设拟合直线 ，并计算得下表： 编号xyxyx212345211315162122130956111011121234344143150176252264442018913121100121144144115661640将数据代入法方程组1-12中，得到：189133446164095695621ba 解方程得：a = 8.2084 , b = 0.1795 。拟合直线为：x . .p( x ) 1795020848 1.3.1 单变量拟合单变量拟合 二次拟合二次拟合函数函数n给定数据序列xi,yi,i=1, 2 , , m ,用二次多项式函数拟合这组数据。212210

10、12210)()(),(imiiimiiiyxaxaayxpaaaQ (1-13) 由数学知识可知，Q( a0 ,a1 ,a2 )的极小值满足：miiiiimiiiiimiiiixyxaxaaaQxyxaxaaaQyxaxaaaQ12221021221011221000)(20)(20)(21.3.1 单变量拟合单变量拟合 二次拟合二次拟合函数函数n整理上式得二次多项式函数拟合的满足条件方程：miiimiiimiimiimiimiimiimiimiimiimiiyxyxyaaaxxxxxxxxm121121014131213121121(1-14) 解此方程得到在均方误差最小意义下的拟合函数p

11、 ( x )。方程组1-14称为多项式拟合的法方程，法方程的系数矩阵是对称的。1.3.1 单变量拟合单变量拟合 二次拟合二次拟合函数函数n上面是二次拟合根本类型的求解方法，和一次拟合一样，二次拟合也可以有多种变型： 例如52310 x a x a ap ( x ) 套用上面的公式，我们可以得到关于求解此拟合函数的法方程 ：miiimiiimiimiimiimiimiimiimiimiimiiyxyxyaaaxxxxxxxxm1513121011018151816131513(1-15) 1.3.1 单变量拟合单变量拟合 二次拟合二次拟合函数函数n假设我们需求求解是下面的拟合函数：n参照上面的方

12、法，我们很容易得到求解该拟合函数的法方程： 5 . 1110)273(273lnxbxaaymiiimiiiimiimiimiimiimiimiimiimiimiiyxxyyaaaxxxxxxxxm15 . 1112101315 . 015 . 115 . 012115 . 11ln)273(273lnln)273()273()273()273()273(12731)273(27311.3.1 单变量拟合单变量拟合 二次拟合二次拟合实例实例n例1.2：请用二次多项式函数拟合下面这组数据。序号1234567x-3-2-10123y4230-1-2-51.3.1 单变量拟合单变量拟合 二次拟合二次

13、拟合实例实例n解：设 ，由计算得下表：2210 x ax a ap ( x ) 序号xyxyx2x2yx3x41234567-3-2-1012304230-1-251-12-4-30-1-4-15-39941014928368-30-1-8-45-7-27-8-1018270811610116811961.3.1 单变量拟合单变量拟合 二次拟合二次拟合实例实例n将上面数据代入式 (1-14) ，相应的法方程为：719602839028012807210210210aaaaaaaaa解方程得：a0 =0.66667 , a1 = -1.39286 , a2 = -0.13095 21309503

14、92861666670 x.x- . - .p (x ) 1.3.1 单变量拟合单变量拟合 二次拟合二次拟合实例实例n拟合曲线的均方误差：n结果见图 1-6。二次曲线的拟合程序可利用后面引见的单变量n次拟合程序。 -3-2-10123-6-4-2024y=0.66667-1.39286 x-0.13095 x2yY X 09524. 3)(712712iiiiiyxp图 1-6 拟合曲线与数据序列1.3.2 多变量的曲线拟合多变量的曲线拟合n实践在化工实验数据处置及模型参数拟合时，通常会碰到多变量的参数拟合问题。一个典型的例子是传热实验中努塞尔准数和雷诺及普兰德准数之间的拟合问题： 32ccP

15、rRecNu11-16 求出方程1-16中参数c1、c2、c3 这是一个有两个变量的参数拟合问题 1.3.2 多变量的曲线拟合多变量的曲线拟合为不失普通性，我们把它表达成以下方式：给定数据序列 用一次多项式函数拟合这组数据。设 ，作出拟合函数与数据序列的均方误差： m iyxxiii, ,3 ,2 1 ), , (2122110 x a x a ap ( x ) 212211012210)()(),(imiiimiiiyxaxaayxpaaaQ (1-17) 1.3.2 多变量的曲线拟合多变量的曲线拟合n由多元函数的极值原理，Q( a0 ,a1 ,a2 )的极小值满足：miiiiimiiiii

16、miiiixyxaxaaaQxyxaxaaaQyxaxaaaQ122211021122110112211000)(20)(20)(2整理得多变量一次多项式函数拟合的法方程：miiimiiimiimiimiiimiimiiimiimiimiimiiyxyxyaaaxxxxxxxxxxm12111210122121121211211112111-18 1.3.2 多变量的曲线拟合多变量的曲线拟合n经过求解方程1-18就可以得到多变量函数线性拟合时的参数。n我们可以经过对方程1-16两边同取对数，就可以得到以下线性方程： PrcReccNulnlnln)ln(3211-19 只需作如下变量代换：32

17、221110 ln lnln )ln(caPrxcaRexcaNuy并将实验数据代入法方程1-18就可以求出方程1-16中的系数。 1.3.2 多变量的曲线拟合多变量的曲线拟合 实实例例n例1.3： 根据某传热实验测得如下数据，请用方程1-16的方式拟合实验曲线。Nu1.1272.4162.2052.3121,4846.0387.325Re100 200 300 500 100 700 800Pr2410.3534n解：利用已给的VB程序，将数据依次输入，就可以得到方程1-16中的三个参数：1.3.2 多变量的曲线拟合多变量的曲线拟合 实例实例3 . 08 . 0023. 0321ccc那么1

18、-16式就变成了常见的光滑管传热方程：3 .08 .0023.0PrReNu n假设拟合方程的方式和方程1-16不同，那么需对上面提供的程序作适当修正。如对以下两个自变量的拟合函数： n其中n1和n2是知系数，我们可以将看作，看作，得到上面拟合函数的法方程： 1.3.2 多变量的曲线拟合多变量的曲线拟合 实例实例2221110nnx a x a ap ( x ) miinimiinimiiminimininiminimininiminiminiminiminiyxyxyaaaxxxxxxxxxxm12211112101222122111221221111211111221111-20 1.4

19、解矛盾方程组解矛盾方程组 n用最小二乘法求解线性矛盾方程的方法来构造拟合函数，并将其推行至恣意次和恣意多个变量的拟合函数。n给定数据序列xi,yi,i=1, 2 , , m ，做拟合直线p (x) = a0 + a1x ，假设要直线 p (x)过这些点，那么就有 p (xi ) = a0 + a1xi =yi, i=1, 2 , , m , 即 ：mmyxaayxaayxaa1022101110 mmyyyaaxxx211021 111矩阵方式： 1.4 解矛盾方程组解矛盾方程组n普通地，将含有n个未知量m个方程的线性方程组： mnmnmmnnnnyxaxaxayxaxaxayxaxaxa22

20、112222212111212111 mnmnmmnnyyyxxxaaaaaaaaa2121212222111211 矩阵方式 普通情况下，当方程数n多于变量数m，且m个方程之间线性不相关, 那么方程组无解，这时方程组称为矛盾方程组。 1.4 解矛盾方程组解矛盾方程组n方程组在普通意义下无解，也即无法找到n个变量同时满足m个方程。这种情况和拟合曲线无法同时满足一切的实验数据点相仿，故可以经过求解均方误差 极小意义下矛盾方程的解来获取拟合曲线。 n由数学的知识还将证明：方程组ATAX = AT b的解就是矛盾方程组AX = b 在最小二乘法意义下的解，这样我们只需经过求解ATAX = AT b就

21、可以得到矛盾方程的解，进而得到各种拟合曲线，为拟合曲线的求解提高了另一种方法。 22minbAX 1.4 解矛盾方程组解矛盾方程组n例如，拟合直线p (x ) = a0 +a1x的矛盾方程组ATAX = AT b的方式如下：mmmmyyyxxxaaxxxxxx2121102121 111111 111化简得到与式(1-12)一样的法方程：miiimiimiimiimiiyxyaaxxxm11101211 1.4 解矛盾方程组解矛盾方程组n对于n次多项式曲线拟合，要计算n Q ( a0 ,a1 , , an )n的极小问题。这与解矛盾方程组 ：miinini - yx+ a+ x + a a12

22、10) (mnmnmnnnnyxaxaayxaxaayxaxaa102221011110 mnyyyaaaA2110 或 与求 miinin110-yxaxaa12的极小问题是一回事。 1.4 解矛盾方程组解矛盾方程组在这里 nmmnnxxxxxxA111 2211故对离散数据xi,yi,i=1, 2 , , m ；所作的n次拟合曲线y= nini10 x+ a+ x+ aa，可经过解以下方程组求得： mTnTyyy A aaaAA2110 1-21 1.4 解矛盾方程组解矛盾方程组n假设拟合函数有n个自变量并进展一次拟合，那么其拟合函数为：n 1-22 miinimiinimiiimiiim

23、iinnminiminiminiminiminiminiminiminiminimiimiimiiminimiimiiyxyxyxyxyaaaaaxxxxxxxxxxxxxxxm1111211121012121111121111111131211121xnnnkkxaxaxaxaxaay1122110经过mmn次实验，丈量得到了m组 ),(, 1, 2, 1ininikiiixxxxxy的实数据，那么可得到上面n个自变量拟合函数的法方程1.4 解矛盾方程组解矛盾方程组n只需对法方程1-22稍加修正，就可以得到有n个自变量的恣意次方的拟合函数的法方程，经过法方程的求，就可以得到拟合函数中的各项系

24、数。 miiinmiiinmiiimiiimiinnmininmiiinmiiinmiinmiininmiinmiiinmiinmiinimiiimiimiimiinmiimiiyxyxyxyxyaaaaxxxxxxxxxxxxxxxxxxxxxm1,1, 112,1, 1111012,1, 2,1, 1,1,1, 112, 11, 1, 11, 11, 11, 2,.112, 11, 11,1, 21, 11-23 1.4 解矛盾方程组解矛盾方程组 实例实例n例 1.4：利用解矛盾方程的方法，用二次多项式函数拟合下面数据。 x -3 -2 -1 0 1 2 3 y 4 2 3 0 -1 -2

25、 -51.4 解矛盾方程组解矛盾方程组 实例实例n解：记二次拟合曲线为 ，构成法方程 2210)(x ax a a xf 72110 yyy A aaaAATnT71i471i371i271i371i271i71i271i2772222112722217217 111 111iiiiiiiiTxxxxxxxxxxxxxxxxxxxxAA19602802802807931421111001111421931941014932101231111111 1.4 解矛盾方程组解矛盾方程组 实例实例739132103249410149321012311111117127171iiiiiiiiTyxyxy

26、YA739119602802802807210aaa得到：解方程得到：a0 = 0.66667 , a1 = -1.39286 , a2 = -0.13095 213095039286166670 x.x - . -. f ( x ) 1.4 解矛盾方程组解矛盾方程组 实例实例n例 1.5：给出一组数据，见下表。用解矛盾方程的思绪将下面数据拟合成 的阅历公式。3 bx a f ( x ) x-3 -2 -1 2 4y14.3 8.3 4.7 8.3 22.71.4 解矛盾方程组解矛盾方程组 实例实例n解：列出法方程： YAbaAATT495436365511111111115165135133

27、5343332313534333231iiiiiiTxxxxxxxxxxxxxAA10623 .587 .223 . 87 . 43 . 83 .14648182711111YAT而：1.4 解矛盾方程组解矛盾方程组 实例实例n故法方程为：n解方程得：na = 10.675 , b = 0. n拟合曲线为：10623 .58495436365ba3137. 0678.10)(xxf1.5 梯度法拟合参数梯度法拟合参数n前面曾经提到函数拟合的目的是使拟合函数和实践丈量值之间的差的平方和为最小，也就求下面函数的最小值：min Q ( a0 ,a1 , , an ) miii - yAX P12)

28、(），（1-24 对于最小值问题,梯度法是用负梯度方向作为优化搜索方向。 1.5 梯度法拟合参数梯度法拟合参数n梯度是一个向量,假设们用向量变量U来表示一切的拟合系数a0 ,a1 , , an，用函数f(U)来替代Q ( a0 ,a1 , , an )，那么函数下降最快的方向为：nSk=- f(U) (1-25)n在梯度法中,新点由下式得到nUk+1=UK- k f(UK) (1-26) n 1.5 梯度法拟合参数梯度法拟合参数梯度法的计算步骤为：(1)选择初始点U0；(2)用数值法(或解析法)计算偏导数 ；(3)计算搜索方向向量： Sk= - ；(4)在Sk方向上作一维搜索，即求解单变量()

29、优化问题 f(Uk+ Sk) 由一维搜索的解k ，求出新点 iuf)(Uf)( Uk+1= Uk+kSk 1.5 梯度法拟合参数梯度法拟合参数n(5)作停顿搜索判别。假设不满足精度要求，前往步骤(2)，反复进展计算。梯度法停顿搜索的判据为：n这个算法的优点是迭代过程简单，要求的存贮也少，而且在远离极小点时，函数的下降还是比较快的。因此，常和其它方法结合，在计算的前期运用此法，当接近极小点时，再改用其它的算法，如共轭梯度法。 )(kUf1.5 梯度法拟合参数梯度法拟合参数n共轭梯度法的计算步骤为：n(1)选择初始点U0或其它方法计算得到的最后点；n(2)计算梯度g0= f(U0) ,以负梯度方向

30、作为初始搜索方向 S0=- g0n (3)在S0 方向上作一维搜索,得到新点U1；nU1= U0+ S0n (4)计算U1点的梯度g1=f(U1)。新的搜索方向S1 ,即共轭方向,为S0与g1的线性组合；nS1=- g1+ S0 00011ggggTT1.5 梯度法拟合参数梯度法拟合参数 对于k1，上式为 Sk+1=-gk+1+Sk 可以证明,由上式得到的方向Sk+1与Sk共轭。 对于多元函数,在n次搜索后(n为变量数)，令U0=Uk+1,然后回到第1步,重新计算共轭方向。 (5)作停顿搜索判据，假设满足，那么停顿搜索。否那么回到第2步，进展反复计算。kTkkTkgggg111.5 梯度法拟合

31、参数梯度法拟合参数 实例实例n例1-8 利用梯度法，用Antoine公式n拟合DEM饱和蒸气压和温度之间的关系。 n解：分析Antoine公式的方式，假设采用解矛盾方程法求解，在进展函数和变量变换后，仍需求进展对C的优化求解，而采用梯法，可直接优化求解，其优化函数为：)(CTBAep71)()()(iiCTBAPeUfi1.6 吸附等温曲线回归吸附等温曲线回归n1.6.1 吸附等温曲线的常见类型吸附等温曲线的常见类型 n1.6.2 几种常用的吸附等温曲线回归方法几种常用的吸附等温曲线回归方法n1.6.3 回归方法的比较回归方法的比较1.6.1 吸附等温曲线的常见类型吸附等温曲线的常见类型 n普通有物理吸附和化学吸附两种。n对于物理吸附而言,单位分量吸附剂吸附吸附质的多少(吸附量)是衡量吸附剂性能好坏的重要目的。n常见吸附等温曲线有以下五种类型，各种不同的类型阐明了n不同的吸附机理，以第一种为例，它是典型的单分子层吸附，其等温曲线的回归