[理学]数值分析第七章.pptVIP

* 插值型求积公式 复合梯形公式及外推计算 Gauss型数值求积公式 数值微分方法 《数值分析》 Ch7 ? ? ? ? 椭圆: x =a cos t y =b sin t 0≤ t ≤2? ? 其中 周长: 数值积分要解决的问题: 1. 定积分与线积分的计算? 2. 重积分的数值计算? 3. 由离散数据计算面积和体积？ 4. 数值积分的误差和精度？ 定积分与积分和式 右矩形和 h 1 0.5 0.2 ······ Sn 5.2908 5.1044 4.9835 ······ 4.8999 左矩形 梯形 右矩形 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 A0 = (b – a )/2 A1 =(b – a )/2 插值型求积公式 对 [a，b]做分划: a≤ x0 x1 x2 …… xn≤b 令 Lagrange插值 插值型求积公式的余项 例2. 梯形公式的误差( 余项 ) 即 ? 例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 公式 定义: 对不高于m次的多项式P(x),求积公式余项 例. 梯形公式 代数精度为1 具有m阶的代数精确度 且有m+1次多项式不具有这样的性质, 则称 (n+1)点插值型求积公式代数精度至少为n阶. 所以, R[xk] = 0, (k = 0,1,2,···,n) 例4 确定公式 使代数精度尽可能高. 类似有: 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 例1. 证明辛卜生公式代数精度为3. Simpsion公式 证 设 f(x) 是任意的三次多项式,记 x3 的 系数为C. 取插值结点x0=a, x1=(a+b)/2, x2=b,有 f(x) =L2(x) + C(x – x0)(x – x1)(x – x2) 所以 x0=a, x1=(a+b)/2, x2= b. 三次多项式H(x)满足 H(x0)=f(x0), H(x1)=f(x1), H(x2)=f(x2), H’(x1)=f’(x1) 插值误差 令 h =(b – a)/2, x = x0+ t h ,则 令 u = t – 1, 则 Simpson公式的求积误差: 令h=(b – a)/(2m) , 取 xk=a+kh ( k =0,1,···,2m) 复合Simpson公式 例 圆周率计算的外推算法 考虑半径为R=0.5的圆内接正多边形周长。设正多边形每条边的边长为ln，对应三角形顶角为 ? 设 ? ? 其中 误差的级数表达式 ? ? 一次外推 其中 记 二次外推 实验数据 边数 周长误差 一次外推误差 二次外推误差 5.4352e-001 1.4159e-001 7.6181e-003 3.5764e-002 4.8793e-004 1.2590e-005 8.9640e-003 3.0683e-005 1.9969e-007