8 数值积分.pptVIP

第8章 数值积分 代数精度 例如，有积分公式: Ck(n) §2 复合求积 §3 龙贝格积分 * 软件+桥梁08（12周-星期一） Cotes系数 7/90 16/45 2/15 16/45 7/90 xi 0 1/8 1/4 3/8 1/2 5/8 3/4 7/8 1 f (xi) 1 0.997397867 0.989615837 0.976726744 0.958851077 0.936155637 00.877192574 0.841470985 xi f (xi) 0 1 1/8 0.997397867 1/4 0.989615837 3/8 0.976726744 1/2 0.958851077 5/8 0.936155637 3/4 07/8 0.877192574 1 0.841470985 复化simpson公式： 1 0.9461458822735868 2 0.9460869339517938 4 0.946083310888472 8 0.9460830853849478 16 0.946083071305562 32 0.9460830704258281 64 0.9460830703708484 128 0.9460830703674122 256 0.9460830703671972 512 0.946083070367183 0 0.9207354924039483 1 0.9397932848061772 0.9461458822735868 2 0.9445135216653896 0.9460869339517938 0.9460830040636742 3 0.9456908635827014 0.946083310888472 0.9460830693509172 0.9460830703872227 4 0.9459850299343862 0.9460830853849478 0.9460830703513795 0.9460830703672598 /* Numerical Integration */ 近似计算 §1 Newton-Cotes 公式 思路 利用插值多项式 则积分易算。 ? 在[a, b]上取 a ? x0 x1 … xn ? b，做 f 的 n 次插值多项式 ，即得到 Ak 由 决定， 与 无关。 节点 f (x) 插值型积分公式 /*interpolatory quadrature*/ 误差 §1 Newton-Cotes Formulae 　　　若某个求积公式所对应的误差R[ f ]满足：R[ Pk ]=0 对任意 k ? n 阶的多项式成立，且 R[ Pn+1 ] ? 0 对某个 n+1 阶多项式成立，则称此求积公式的代数精度为 n 。 定义 例：对于[a, b]上1次插值，有 f(x) a b f(a) f(b) 梯形公式 /* trapezoidal rule*/ 解：逐次检查公式是否精确成立 代入 P0 = 1： = 代入 P1 = x ： = 代入 P2 = x2 ： ? 代数精度 = 1 考察其代数精度。 求该积分公式的代数精确度。 对于任意一个一次多项式，求积公式都是精确成立的； 至少存在一个二次多项式使求积公式不精确成立； 故该求积公式的代数精确度为1。 解：取f(x)=1， 取f(x)=x ， 取f(x)=x2 ， ? = = §1 Newton-Cotes Formulae 注：形如 的求积公式至少有 n 次代数精度 ? 该公式为插值型（即： ） ? 当节点等距分布时： 令 Cotes系数 注：Cotes 系数仅取决于 n 和 i，可查表得到。与 f (x) 及区间[a, b]均无关。 ? §1 Newton-Cotes Formulae n = 1: Trapezoidal Rule /* 令 x = a+th, h = b?a, 用中值定理 */ 代数精度 = 1 n = 2: Simpson’s Rule 代数精度 = 3 n = 3: Simpson’s 3/8-Rule, 代数精度 = 3, n = 4: Cotes Rule, 代数精度 =