数值分析课件 第四章 数值积分.pptVIP

§2 Composite Quadrature ? 收敛速度与误差估计： 定义 　　 若一个复化积分公式的误差满足 且C ? 0，则称该公式是 p 阶收敛的。 /*中值定理*/ 类似的，可得 2阶收敛 4阶收敛 6阶收敛 例1：计算 解： 其中 = 3.138988494 其中 = 3.141592502 运算量基本相同 Q: 给定精度 ?，如何取 n ? 例如：要求 ，如何判断 n = ？ 上例中若要求 ，则 即：取 n = 409 §2 Composite Quadrature 事后误差估计式，可用来判断迭代 是否停止。 始步长h 0.5?10-2 §4 龙贝格积分 /* Romberg Integration */ 复化梯形公式算法简单，但精度较差，收敛速度（2阶收敛）较慢，如何提高收敛速度？ 注：按上面规律，可以构造线性组合系数为 的新的积分公式，但当m?4时，前一个系数接近于1，后一个系数接近于0，这样构造出的新公式与前一个公式结果差别不大，反而增加计算量，因此实际上常做到Romberg公式为止。 * 第四章 数值积分与数值微分 /* Numerical Integration and differentiation*/ 近似计算 §1 引言 ? 对f(?)采用不同的近似计算方法，从而得到各种不同的求积公式。 ? 以上三种方法都是用被积函数值的线性组合来表示积分值。推广，一般地有 求积节点 求积系数，与被积函数无关 像这样，将积分用若干节点上被积函数值的线性组合来表示的数值积分公式称为机械求积公式。 ? 求积误差 ? 机械型求积公式的构造归结为，确定求积节点xk和求积系数Ak，使在某种意义下精确度较高。总之，要解决三个问题： 精确度的度量标准； 如何构造具体的求积公式； 具体求积公式构造出来后，误差如何估计？ 定义：代数精度 　　　若某个求积公式对次数? m 阶的多项式准确成立，而对m+1 阶的多项式不一定准确成立。即对应的误差满足：R[ Pk ]=0 对任意 k ? m 阶的多项式成立，且 R[ Pm+1 ] ? 0 对某个 m+1 阶多项式成立，则称此求积公式的代数精度为 m 。 ?代数精度与误差的关系：代数精度越高，求积误差越小。 结论： 问题1 由上面代数精度条件确定求积公式可分两种情形： 若事先给定求积节点xk(k=0,…,n),例如被积函数以表的形式给出时xk确定，可令m=n,由上式确定n+1个系数Ak即可----待定系数法和插值法。 若xk和Ak都可选择，令m=2n +1，确定xk和法Ak ---Gauss法 要使求积公式具有m阶代数精度，则它对1,x,…,xm均准确成立，即 m+1个方程，2n+2个未知数 问题2 Case 1---方法1 §1 插值型求积 公式 思路 利用插值多项式 则积分易算。 ? 在[a, b]上取 a ? x0 x1 … xn ? b，做 f 的 n 次插值多项式 ，即得到 Ak 由 决定， 与 无关。 节点 f (x) 插值型积分公式 /*interpolatory quadrature*/ 误差 Case 1---方法2 §1 Newton-Cotes Formulae 例：对于[a, b]上1次插值，有 考察其代数精度。 f(x) a b f(a) f(b) 梯形公式 /* trapezoidal rule*/ 解：逐次检查公式是否精确成立 代入 P0 = 1： = 代入 P1 = x ： = 代入 P2 = x2 ： ? 代数精度 = 1 §1 Newton-Cotes Formulae Th1.形如 的求积公式至少有 n 次代数精度 ? 该公式为插值型（即： ） ? 当节点等距分布时： 令 Cotes系数 注：Cotes 系数仅取决于 n 和 i，可查表得到。与 f (x) 及区间[a, b]均无关。 ? §2 Newton--Cotes 公式 Newton—Cotes formula §1 Newton-Cotes Formulae n = 1: Trapezoidal Rule /* 令 x = a+th, h = b?a, 用中值定理 */ 代数精度 = 1 n = 2: Simpson’s Rule 代数精度 = 3 n = 4: Cotes Rule, 代数精度 =