ch3数值积分与数值微分.pptVIP

Ch4 数值积分与数值微分 Numerical Integration 解：逐次检查公式是否精确成立 代入 P0 = 1： = 代入 P1 = x ： = 代入 P2 = x2 ：  代数精度 = 1 例 求梯形公式的代数精确度 插值型求积公式 节点 f (x) 插值型积分公式 /*interpolatory quadrature*/ 注：Cotes 系数仅取决于 n 和 i，可查表得到。与 f (x) 及区间[a, b]均无关。 Newton-Cotes 公式 Simpson’s formula 可以验证代数精度 = 3 定理3 当阶n为偶数时，netow-cotes公式至少 具有n+1次代数精确度. 求积公式的余项 复化求积公式 高次插值不稳定，故采用分段低次插值  分段低次合成的 Newton-Cotes 复合求积公式。  复合梯形公式：  复化 Simpson 公式： = Sn 龙贝格求积公式 Romberg Integration 考察复合梯形公式，区间二分前后： 理查森（Richardson）外推法 理查森外推加速方法 高斯求积公式 代入 f (x) = 1, x, x2, x3 不是线性方程组，不易求解。 §4 Gaussian Quadrature 证明： “” 对任意次数不大于n 的多项式 Pm(x)， Pm(x) w(x)的次数不大于2n+1，则代入公式应精确成立： = 0   求 Gauss 点  求w(x) Step 1：构造正交多项式2   §4 Gaussian Quadrature Step 2：求2 = 0 的 2 个根，即为 Gauss 点 x0 ，x1 Step 3：代入 f (x) = 1, x 以求解 A0 ，A1 解线性方程组，简单。 §4 Gaussian Quadrature  特殊正交多项式族： ① Legendre 多项式族： 以 Pn+1 的根为节点的求积公式称为Gauss-Legendre 公式。 ② Chebyshev 多项式族： 注意到积分端点 1 可能是积分的奇点，用普通Newton-Cotes公式在端点会出问题。而Gauss公式可能避免此问题的发生。 数值微分 中点公式，精度高 注：h不宜太小，也不宜太大 插值型求导公式： 给定函数表，计算函数f(x)的导数 x x0 x1 x2 … xn y y0 y1 y2 … yn 利用数值微分 3次样条插值 外推法