数值分析20定积分计算积分和式.pptVIP

* 定积分计算与积分和式 插值型求积公式 求积公式的代数精度 复合梯形公式 MATLAB求积分命令 《数值分析》 20 椭圆: x =a cos t y =b sin t 0≤ t ≤2? ? 2/18 其中 周长: 问题: 1. 定积分与线积分的计算? 2.重积分的数值计算? 3. 椭球面积的计算？ 4. 由离散数据计算三维体积？ 3/18 4/18 定积分与积分和式 右矩形和 h 1 0.5 0.2 ······ Sn 5.2908 5.1044 4.9835 ······ 4.8999 5/18 左矩形 梯形 右矩形 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/18 A0 = (b – a )/2 A1 =(b – a )/2 插值型求积公式 对 [a，b]做分划: a≤ x0 x1 x2 …… xn≤b 令 7/18 Lagrange插值 插值型求积公式的余项 例2. 梯形公式的误差( 余项 ) 即 8/18 ? 例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/18 定义: 对不高于m次的多项式P(x),求积公式余项 例. 梯形公式 代数精度为1 具有m阶的代数精确度 且有m+1次多项式不具有这样的性质, 则称 10/18 (n+1)点插值型求积公式代数精度至少为n阶. 所以, R[xk] = 0, (k = 0,1,2,···,n) 例4 确定公式 使代数精度尽可能高. 类似有: Simpson公式具有3阶代数精度 对于n次Lagrange插值基函数,有恒等式 11/18 ? 解: 取f(x)= 1, x, x2 若求积公式准确成立,则有 A1 = 0 求积公式 具有至少2阶代数精度 容易验证, 对f (x) = x3 求积公式式不能准确成立. 因此这一公式只具有2次代数精度 12/18 取等距结点xj = a + jh时,插值型求积公式称为Newton-Cotes公式 定理: 当n为偶数时, n阶Newton-Cotes公式至少有(n+1)阶代数精确度。 Newton-Cotes公式代数精度至少为n 13/18 ? *