计算物理课件 Chapter2.pptVIP

第2章 数值微分和数值积分 2.1 数值微分 2.2 数值积分 2.3 复化数值积分 2.4 重积分计算 2.5 高斯型积分公式 2.6 程序示例 2.1 数值微分 导数与差商的关系 数值微分取成导数的近似值，即差商。 2.1.1 差商型数值微分 数值微分用向前差商表示 截断误差：用Taylor展式， 存在ζ∈[x0,x0+h]， 数值微分用向后差商表示 同理得截断误差 数值微分用中心差商表示 截断误差：由Taylor展式， 数值微分的几何意义 设定最佳步长 计算数值微分时产生的误差有截断误差和舍入误差两部分组成。由数值微分公式产生截断误差，有原始数据产生舍入误差。 一般情况，步长 h 越小，误差也越小，但步长太小，会引起误差的增长。因此，实际计算时，需要选择一个最佳步长。 用中心差商分析数值微分的误差。 ① 截断误差： ② 舍入误差： 设 e 为计算 f(x0-h) 和 f(x0+h) 时的最大舍入误差，可以证明，舍入误差量不超过 e/h。 于是得到计算数值微分的总误差为 实际计算中，采用事后估计方法来选去步长h/2：若 例2.2 对函数 y = ex,选取不同的步长进行计算 f (1.15),观察 误差的变化规律，确定最佳步长。 解 用中心差商表示的数值微分计算公式得到： h f (1.15) error h f (1.15) error 0.1 3.1630 -0.0048 0.05 3.1590 -0.0008 0.09 3.1622 -0.0040 0.04 3.1588 -0.0006 0.08 3.1613 -0.0031 0.03 3.1583 -0.0001 0.07 3.1607 -0.0025 0.02 3.1575 0.0007 0.06 3.1600 -0.0018 0.01 3.1550 0.0032 实际计算时，准确误差难以得到，∵|D(0.05)-D(0.04)|ε, 其中， ε = 2×10-4 最小。所以取h = 0.04。 2.1.2 差值型数值微分 给定 (xi, f(xi)), i = 0,1,2,…,n, 构造Ln(x),则 代入节点xj, j = 0,1,2,…, n, 得到 例2.3 给定三个差值点(xi, f(xi)), i = 0,1,2, 求过三个节点的数 值微分。 解 记 将 x = xi 代入 f (x), 得到三点数值微分公式： 取 n =2，由(2.5)1式，得三点公式的截断误差分别为： 将 f(x1)和 f(x2)在 x0点Taylor展开，也可得截断误差为O(h2)。 2.1.3 样条插值数值微分 给定差值点(xi, f(xi)), 假定 S (xi) = mi, i = 0,1,2,…, n, 由“m关系式”得三次样条插值函数： 2.2 数值积分 数值积分是用来求解定积分中被积函数是以离散点形式(xi, f(xi)) 给出；或被积函数无法求出的情形。它是定积分计算的一种近似方法。 在微积分中，定积分是Riemann和的极限，即 数值积分就是取定积分极限中的有限项的和，即 其中，xi 称为积分节点,αi 称为积分系数。 我们的任务就是确定积分系数αi ，使 I( f )≈In( f )。具体地，就是用插值多项式近似代替被积函数 f(x) 来确定αi。 ★ 代数精度 设[a, b]上以 xi 为积分节点的数值积分公式为 若对任意次数不超过 m 次的多项式成立 In(xk) = I(xk), k = 0,1,2,…, m, 而对 m +1次多项式有 In(xm+1) ≠ I(xm+1), 则称公式