化工计算方法化工计算方法插值法.pptVIP

3 插 值 法 插 值 法 目的：运用计算机方便地求解函数值和进行微积分等计算 思路：寻找列表函数的近似解析表达式 (不知道 f(x) 的函数表达形式) 插 值 法 3.2 拉格朗日(Lagrange)插值3.2.1 线性插值（2个节点） 由两点式可看出， L1(x) 是由两个线性函数 线性插值图示 3.2.2 二次插值（三点插值或抛物线插值） 二次插值（三点插值或抛物线插值） 利用三个二次插值基函数可组合成二次插值多项式 3.2.3 n次插值 3.2.3 n次插值 3.2.3 n次插值 用11点构造10次多项式插值的龙格现象 某些点插值结果误差很大，函数两端震荡加剧 在节点很多的场合，通常不宜采用高次插值 分段的低阶插值往往效果更好 例题 3-1 已知函数表 MATLAB有可供直接调用的插值函数。以下参考程序调用一元插值函数interp1 例题 3-2 实验测得某物质在20℃下，其水溶液浓度c (重量％) 3.3 二元插值 函数值随两个自变量变化的情形： 化工原理中摩擦系数既随雷诺数Re也随相对粗糙度变化； 二维传热时，某点的温度同时随轴向和径向位置变化; 二元函数以表格形式表示时称为二元列表函数。二元列表函数的插值原理和方法与上述一元情况类似； 将一元线性插值推广到二元的情况就是二元线性插值； 一元三点插值推广到二元的情况就是二元三点插值； 假设两个自变量的插值结点分别为xi 和yj，其对应的函数值为zij，在x和y方向各取n个点，构成二元插值基函数，然后进行线性组合得到二元n点拉格朗日插值公式为 例题 3-4 已知函数表 例题 3-3 精馏计算举例 精馏苯-甲苯混合物，已知进料F=8570 kg，苯的浓度（质量分率），回流比R＝3，试分别计算进料状态参数q=1.23和q=0时的理论塔板数和各板液气相组成。苯-甲苯平衡数据（摩尔分率）如下表 例题 3-3 精馏计算举例 3. 操作线方程（精馏段，提馏段） 精馏段 提馏段 式中提馏段液体量和气体量分别为 例题 3-3 精馏计塔计算步骤 精馏塔计算程序框图 3.4 样条插值 3.4.1 三次样条插值函数 3.4.2 三次样条插值函数S(x)的构建 三次样条插值函数S(x)的表达式 Mj的求解 Mj的求解 Mj的求解 本章要求 * * 什么叫插值？ 31.164 19.923 12.34 7.3766 4.2474 2.3346 1.2262 0.6082 饱和蒸汽压kN/m2 292.99 70 251.12 60 209.30 50 167.51 40 125.69 30 83.90 20 42.04 10 0 0 焓kJ/kg 温度 oC 如何查水在27oC、32.7oC的饱和蒸汽压和焓? 水的物理性质 函数关系：函数值和自变量的关系以表格给出，称列表函数 列表函数的特点： ① 自变量与函数值一一对应； ② 函数值有很可靠的精确度； ③ 自变量与函数间的解析表达式可能不清楚，或者解析表达式非常复杂不便于计算（如为无穷级数等）； ④ 没有直接给出未列出点的函数值，不便于进行微分和积分以及计算机计算。 ＃ 代数插值问题描述（怎样进行插值） 列表函数例 在区间上连续 设函数 已知函数 f (x) 在各点(x0 , x1 …)上的值 方法：建立一个次数不超过n的代数多项式 P(x) 来近似 f(x) 31.164 19.923 12.34 7.3766 4.2474 2.3346 1.2262 0.6087 饱和蒸汽压f(x) x7= 70 x6= 60 x5= 50 x4= 40 x3= 30 x2= 20 x1= 10 x0= 0 温度 x 用代数多项式P(x)近似列表函数 y 插值方式： 全节点插值——用全部节点构造插值多项式（通过全部节点） 分段插值——只用部分节点构造插值多项式（只通过部分节点） 本节讨论 一元n 次拉格朗日插值(n +1个点) 线性插值（2个节点） 抛物线插值（3个节点） ＃ 条件：满足 Pn(xi)=yi (i=0,1,…n) Pn(x）称为函数 f(x) 的插值多项式 点 x0, x1, … xn 叫做插值节点 [a，b]为插值区间。 y7= 31.164 y6= 19.923 y5= 12.34 y4= 7.3766 y3= 4.2474 y2= 2.3346 y1= 1.2262 y0=0.6087 饱和蒸汽压f(x) x7= 70 x6= 60 x5= 50 x4= 4