Newton-Cotes公式代数精度至少为n1316.pptVIP

* 定积分计算与积分和式 插值型求积公式 求积公式的代数精度 复合梯形公式 《数值分析》 20 ? ? ? ? 椭圆: x =a cos t y =b sin t 0≤ t ≤2? ? 2/16 其中 周长: 数值积分要解决的问题: 1. 定积分与线积分的计算? 2. 重积分的数值计算? 3. 由离散数据计算面积和体积？ 4. 数值积分的误差和精度？ 3/16 4/16 定积分与积分和式 右矩形和 h 1 0.5 0.2 ······ Sn 5.2908 5.1044 4.9835 ······ 4.8999 5/16 左矩形 梯形 右矩形 4.4429 4.8669 5.2908 4.6804 4.8924 5.1044 4.8139 4.8987 4.9835 4.8572 4.8996 4.9420 数值求积公式的一般形式 R[f ] —— 数值求积公式余项 x0, x1, ···, xn —— 求积结点 A0, A1, ···, An —— 求积系数 例1. 梯形公式: a b 6/16 A0 = (b – a )/2 A1 =(b – a )/2 插值型求积公式 对 [a，b]做分划: a≤ x0 x1 x2 …… xn≤b 令 7/16 Lagrange插值 插值型求积公式的余项 例2. 梯形公式的误差( 余项 ) 即 8/16 ? 例3. 取 x0 =a, x1 =0.5(a+b), x2 = b ,则 h=0.5(b – a ) A0= (b-a)/6 A1=2(b-a)/3 A2= (b-a)/6 即著名的 Simpson 公式 9/16 定义: 对不高于m次的多项式P(x),求积公式余项 例. 梯形公式 代数精度为1 具有m阶的代数精确度 且有m+1次多项式不具有这样的性质, 则称 10/16 (n+1)点插值型求积公式代数精度至少为n阶. 所以, R[xk] = 0, (k = 0,1,2,···,n) 例4 确定公式 使代数精度尽可能高. 类似有: Simpson公式具有3阶代数精度 对于n次Lagrange插值基函数,有恒等式 11/16 ? 解: 取f(x)= 1, x, x2 若求积公式准确成立,则有 A1 = 0 求积公式 具有至少2阶代数精度 容易验证, 对f (x) = x3 求积公式式不能准确成立. 因此这一公式只具有2阶代数精度 12/16 *