数值微分与数值积分研究.ppt

* n xk Ak 0 0.0000000 2.0000000 1 ±0.5773503 1.0000000 2 ±0.7745967 0.0000000 0.5555556 0.8888889 3 ±0.8611363 ±0.3399810 0.3478548 0.6521452 4 ±0.9061798 ±0.5384693 0.0000000 0.2369269 0.4786287 0.5688889 高斯-勒让德求积公式的节点和系数 * 例题 例17 应用Gauss-Legendre公式计算积分 解 先作变量代换化积分区间[0,?]为[-1,1]，令 则 精确值为 -12.0703463 * 例题 应用两点Gauss－Legendre公式，得 应用五点Gauss－Legendre公式的计算结果为 * 例18 利用三点高斯-勒让德求积公式计算 的近似值。精确值为 2.399529 解:由表可知，得到三点高斯型求积公式为 由所求公式得 例题 * §4.8 数值微分 根据函数在若干个点处的函数值去求该函数的导数的近似值称为数值微分。 * 最简单的数值微分是用差商（向前差商）近似代替导数，即 同样，也可用向后差商近似代替导数，即 差商 * 或中心差商的方法，即 差商 可以看出中心差商是向前差商和向后差商的算术平均值。 上述三种方法的截断误差分别为 、 和 * 可见弦BC的斜率更接近于切线AT的斜率 ，因此从精度方面看，用中心差商近似代替导数值更可取。 如右图所示，前述三种导数的近似值分别表示弦线AB，AC和BC的斜率，将这三条通过A点的弦的斜率与切线AT的斜率进行比较后, 差商 * 例题 例22 已知函数表如下： 1.3 1.5 1.7 1.9 2.1 2.3 2.5 3.669 4.482 5.474 6.686 8.166 9.974 12.182 试求 在x＝1.7处的导数值。 * 例题 解 取h＝0.2,得 应用向前差商公式,得 应用向后差商公式,得 * 例题 真值5.473 应用中心差商公式,得 * 本章介绍了积分和微分的数值计算方法，其基本原理主要是逼近论，即设法构造某个简单函数P(x)近似表示f(x)，然后对P(x)求积得到f(x)的积分的近似值。基于插值原理，推导了数值积分的基本公式。 插值型求积公式介绍了牛顿-柯特斯公式和高斯公式两类。前者取等距节点，算法简单而容易编制程序。但是，由于在n≥8 时出现了负系数，从而影响稳定性和收敛性。因此实用的只是低阶公式。 本章小结 * 解决长区间与低阶公式的矛盾是使用复合求积公式，因此，常用的数值积分法都是复合求积公式。高斯公式不但具有最高代数精度，而且收敛性和稳定性都有保证，因此是高精度的求积公式。高斯公式的主要缺点是节点与系数无规律。所以高阶高斯公式不便于上机使用。 本章小结 * * * 例13 给定求积公式 试确定求积系数A-1, A0 ,A1, 使其有尽可能高的代数精度，并指出其代数精度。 例题 * 解：令求积公式对f(x)=1, x, x2准确成立，则有 例题 解之得 * 其代数精度至少为2,将f(x)=x3代入求积公式两端相等,而将将f(x)=x4代入求积公式两端不相等,所以其代数精度为3次。 例题 * 例 14 确定求积公式 使其具有尽可能高的代数精度。 例题 * 解：不妨设a=0, b=h, b-a=h, 设所求公式的代数 精度为2,则当f(x)=1,x,x2时公式变成等式,即 例题 解之得： 其中h=b-a, 令f(x)=x3代入上式, 两端不等, 说明求积公式只有2次代数精度。 * 例15 求证当节点为n+1个时, 插值求积系数之和为 例题 * 例题 证： 当节点为n+1个时，插值求积公式有n次代数精度，对于 ，上式严格相等，所以取f(x)=1时，上式也严格相等，因此有 即 * 构造插值求积公式有如下特点： 复杂函数f(x)的积分转化为计算多项式的积分 求积系数Ak只与积分区间及节点xk有关，而与被积函数f(x)无关，可以不管f(x)如何，预先算出Ak的值 n+1个节点的插值求积公式至少具有n次代数精度 求积系数之和 可用此检验计算求积系数的正确性 插值求积公式的特点 * 2 高斯求积公式的构造与应用 在构造形如 的两点公式时,如果限定求积节点, 那么所得插值求积公式 高斯求积公式的构造与应用