化工计算方法化工计算方法数值微分与积分.pptVIP

5 数值微分与积分 5 数值微分与积分 数值微分2 数值微分2 数值微分2 例5-1 数值微分结果显示（二阶精度） 5.2 数值积分基础 选择什么样的简单函数去近似被积函数呢？ 5.2.1 梯形公式（n=1） 5.2.2 辛普森（Simpson）公式（n=2） 5.2.3 牛顿-科特斯（Newton-Cotes）公式（n等分） 求Ak 5.3 复化求积公式 5.3.2 复化辛普生（Simpson）公式 5.3.3 变步长辛普生求积 5.4 龙贝格（Romberg)积分法 5.4 龙贝格（Romberg)积分法 5.4 龙贝格（Romberg)积分法 例题5-2：分别用定步长辛普生积分、变步长辛普生积分、龙贝格求积法求定积分的 值，定步长求积时取N = 6。 例题 5-3 (本例选自《化工原理 》) 用水吸收混合气中的丙酮。进气丙酮含量0.05(mol%)，塔内吸收98%，吸收系数Kv=1.1111×10-4，常压和20℃时气体流量Vt=0.5556 m3/s。相平衡关系式为Y=MX，其中M=1.68。求： (1）出塔吸收液的最大浓度和吸收剂最小用量； （2）当吸收液的浓度为最大浓度的60％时的吸收剂用量； （3）若填料比表面积 a0=190 m2/m3，空塔气速 u=0.8m/s，计算塔径和填料层高度。 例题 5-3 (本例选自《化工原理》) (b) 计算X1m，Lmin，X1，L 惰性气体流量(kmol/s) V=[273×(1-0.05)/293×22.4]×Vt 吸收速率(kmol/s) GA=V ×(Y1-Y2) 吸收液出口最大浓度 Xlm=Y1/M (与Y1平衡的液体浓度） 吸收剂最小用量(kmol/s) Lmin = GA/(X1m-X2) 吸收液出口实际浓度 X1= 0.6 Xlm 吸收剂实际用量(kmol/s) L = GA×(X1-X2) 例题 5-3 (本例选自《化工原理》) （2）计算程序及结果 例题：已知函数值如下表 计算结果比较 例题 5-4 用脉冲法测量某反应器的停留时间分布，由实验数据计算寿命分布密度与时间t的关系如下表 本章要求 清楚数值积分的原理 掌握数值积分的梯形公式和辛普生公式、复化梯形公式和复化辛普生公式 了解龙贝格积分的主要思想 能用数值积分求解实际问题 * * 5.1 数值微分 ——根据离散点上的函数值求取某点导数近似值的方法。 1）差商代替导数 2)插值型数值求导 向前差商 用差商近似代替导数，f（x）在xi处的差商形式为 例:已知自变量和函数值如下， 求各点的导数值。 0.35 0.8 0.29 0.9 0.41 0.5 0.51 0.6 0.43 0.7 0.31 0.2 0.33 0.3 0.36 0.4 0.28 0.38 0.48 f (x) 1.0 0.1 0 x 向后差商 中心差商 yi yi-1 yi+1 y x xi-1 xi xi+1 向前 向后 整理后得向前差商 右端除第一项外其后各项省略 ?x=h (?x)2=h2 (O) ?x # f(xi)在其相邻点 xi+1处可展开成泰勒级数 考查差商代替导数的误差 O(h)——意为截断误差是步长 h 的一次方量级（一阶精度） 整理可得向后差商 在相邻点 xi-1处可展开成泰勒级数 O(h)——一阶精度 将 xi+1和xi-1处式泰勒级数展开式相减得中心差商 O(h2)——二阶精度 对 i=1和 i=n（端点），二阶精度的数值微分式分别为 2)插值型数值求导 用插值多项式Ln(x)近似f(x)，对Ln(x)求导作为f(x)的导数的近似值，称为插值型求导公式。 线性插值型公式对x求导 两点公式 两个点的求导公式形式相同，两点公式的截断误差可用插值余项公式(3-18)计算，为O(h) （一阶精度） 同理可由三点插值公式求导得三点公式 ＃ 上式求导可得三点求导公式 二阶导数的三点公式 三点求导公式的截断误差为O(h2) （二阶精度） 从截断误差的角度看，减小步长h会提高计算的精度，但另一方面，h的减小会使舍入误差增大。由于数值微分计算对舍入误差非常敏感，因此步长h不宜太小。 例题5-1 实验测得一组数据如下，求x=0.8处及各点的导数 解 先求 x=0.8处 导数 df/dx 0.35 0.8 0.29 0.9 0.41 0.5 0.51 0.6 0.43 0.7 0.31 0.2 0.33 0.3 0.36 0.4 0.28 0.3