数值积分基本概插值型求积公式求积公式的代数精度复合梯.pptVIP

* * 数值积分基本概念 插值型求积公式 求积公式的代数精度 复合梯形公式 格林公式中曲线积分处理 数值积分 定积分与积分和式 右矩形和 h 1 0.5 0.2 ······ Sn 5.2908 5.1044 4.9835 ······ 4.8999 数值求积公式的一般形式 R[f ] —— 数值求积公式余项 x0, x1, ···, xn —— 求积结点 A0, A1, ···, An —— 求积系数 插值型求积公式 对 [a，b]做分划: a≤ x0 x1 x2 …… xn≤b 令 Lagrange插值 插值型求积公式的余项 梯形公式的误差( 余项 ) 即 ? 例1. 梯形公式 线型插值 a b 左矩形 梯形 右矩形 4.4429 4.8669 5.2908 4.6804 4.8924 5.1044 4.8139 4.8987 4.9835 4.8572 4.8996 4.9420 取 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 公式 例2. Simpson 公式 L(x)=l0(x)y0+l1(x)y1+l2(x)y2 定义: 对不高于m次的多项式P(x),求积公式余项 例. 梯形公式 代数精度为1 具有m阶的代数精确度 且有m+1次多项式不具有这样的性质, 则称 求积公式的代数精度 (n+1)点插值型求积公式代数精度至少为n阶. 所以, R[xk] = 0, (k = 0,1,2,···,n) 例3 确定公式 使代数精度尽可能高. 类似有: Simpson公式具有3阶代数精度 对于n次Lagrange插值基函数,有恒等式 ? 解: 取f(x)= 1, x, x2 若求积公式准确成立,则有 A1 = 0 求积公式 具有至少2阶代数精度 容易验证, 对f (x) = x3 求积公式式不能准确成立. 因此这一公式只具有2阶代数精度。 取等距结点xj = a + jh时,插值型求积公式称为Newton-Cotes公式 定理: 当n为偶数时, n阶Newton-Cotes公式至少有(n+1)阶代数精确度。 Newton-Cotes公式代数精度至少为n ? 复合梯形求积公式 将积分区间[a,b] n 等分.取 h=(b-a)/n . xj=a+jh 取 递推，得 给定允许误差界ε0,当 时，结束计算并以T2n作为定积分的近似值. T1 ? T2 ? T4 ? ········ ? Tn ? T2n 格林公式的离散化——多边形面积计算 令 Q=x, P= –y 得区域 D 的面积计算公式 将? 分划, ?j : (xj，yj)?(xj+1，yj+1) ( j = 1,2,···,n) t∈( 0，1 ) 参数方程: D 格林公式中曲线积分处理 ? 面积公式 顶点按逆时针排列,且 (xn+1，yn+1)=(x1，y1) 蒙特卡罗法求积分 N =2000: q= 4.8975, 4.9256, 4.7550, 4.9800 ····· 在D中投入N个点,落入曲边梯形内的点数为n 复合梯形公式 误差估计 外推计算 龙贝格外推计算公式 记 则 欧拉-麦克劳林公式 有:T2(h) = I + O(h6) 记 记 所以 ↓T0(0) ↓T0(1) → T1(1) ↓T0(2) → T1(2) → T2(2) ↓T0(3) → T1(3) → T2(3) → T3(3) 记 龙贝格外推计算公式 高斯型数值求积公式 插值型求积公式 代数精度为3,取 f(x)=1, x, x2, x3 (1) (2) (3) (4) (4)-(2)×x02 x12= x02 (3)-(1) ×x02 x02=1/3 ? ? ? 代数精度为 3 的数值求积公式 对于[a, b]区间上的定积分,构造变换 t∈[-1, 1] 定义 如果求积结点x0, x1,······,xn,使插值型求积公式 的代数精度为2n+1,则称该求积公式为Gauss型求积公式. 称这些求积结点为Gauss点. 定理7.2 如果多项式wn+1(x)=(x – x0) (x – x1)···(x – xn) 与任意的不超过n次的多项式P(x