计算方法课件11月7日.pptVIP

Today’s Plan2011-11-7 上节课重点回顾： 数值积分， 梯形求积公式，抛物线求积公式， Newton-Cotes求积公式 今日新内容： 梯形求积公式和抛物线求积公式 的误差估计 上节课(2011-10-31)重点回顾 谢谢陆鋆同学 带领大家做回顾总结！ 梯形求积公式 以x0=a , x1=b为插值节点,作线性插值多项式， 对插值多项式求积分。 几何意义:用梯形面积代替曲边梯形面积. 抛物线求积公式 将区间[a,b]二等分，以分点x0=a , x1=(a+b)/2, x2=b为插 值节点,作二次插值多项式，对插值多项式求积分。 几何意义:用抛物线围成的曲边梯形的面积代替y=f(x)围成的曲边梯形面积 当n=1时， Newton-Cotes公式为梯形公式 当n=２时， Newton-Cotes公式为抛物线公式 今日新内容 梯形求积公式/抛物线求积公式 的误差估计 代数精确度的定义 若用区间左端点高度为曲边梯形平均高度 则得到左矩形公式 ? 若用区间右端点高度为曲边梯形平均高度, 则得到右矩形公式 ? 若用梯形的面积近似地代替曲边梯形的面积,则得到计算定积分的梯形求积公式 当n=1时， Newton-Cotes公式为梯形公式 当n=２时， Newton-Cotes公式为抛物线公式 求积公式对f(x) =1， x， x2 ， x3准确成立： 但对f(x) =x4不成立 所以该求积公式具有3次代数精度。 n=1: 梯形公式 具有 1次 代数精度 n=2: 抛物线公式 具有3 次代数精度 对于n次的Newton-Cotes求积公式 当n为奇数时，具有n次代数精度； 当n为偶数时，具有n+1次代数精度. 梯形公式的余项估计 定理　若f(x)? C2[a,b] ,则梯形公式余项为 证明: 由插值余项定理知 等式两边积分得 由于f(x)?C2[a,b] ，且(x-a)(x-b)在[a,b]内非正 (不变号),故根据积分中值定理知,至少存在一点??[a,b],使 抛物线公式的余项 定理　若f(x)? C4[a,b] , 则抛物线公式余项为 证明: 抛物线公式的代数精度为3,为此构造三次多项式φ3(x)满足插值条件 φ3(a)=f(a), 则余项为 等式两边从a 到b 积分得 由于φ3(x)是三次多项式，故抛物线公式对它准确成立，即 这样 由于f(x)?C4[a,b] ，且 在[a,b] 非正(不变号),故根据积分中值定理知,至少存在一 点??[a,b],使 这样 由于f(x)?C4[a,b] ，且 在[a,b] 非正(不变号),故根据积分中值定理知,至少存在一 点??[a,b],使 （代数精度2）抛物线公式余项为 例1. 用梯形公式和抛物线公式计算积分 ,并估计误差. 解 记a=0, b=1, f(x)=e-x ,则f (x)=-e-x f(x)=e-x , f(x)=-e-x , f(4)(x)=e-x 例2 设有求积公式 求A0，A1，A2，使其代数精度尽量高，并问此时求积公式的代数精度 解：（3个未知系数需三个方程） 令求积公式分别对f(x) = 1、x、x2准确成立。即 ? 解之得A0 = A2 = 1/3，A1 = 4/3，即有 当n为偶数时，具有n+1次代数精度. 当n为奇数时，具有n次代数精度， 但不具有n+1次代数精度？ 因为被积函数 n=2k为偶数，则 为整数， 为奇函数，所以 再令 n为奇数则以下不再成立： 结论： 一个求积公式的代数精度越高， 它就越能对更多的被积函数f(x) 准确地成立， 从而具有更好的实际计算意义。 求积公式的代数精度概念是衡量 公式逼近好坏的标准之一. 误差/余项估计 真实值与近似计算值的差---截断误差 b a b a x g x f b a x g x f ] , [ ] , [ ) ( ) ( ] , [ ) ( ), ( ò ò = b a b a dx x g f dx x g x f ) ( ) ( ) ( ) ( x x ，使得 一个数 上至少存在 上不变号,则在区间 在区间 连续， 上有界可积， 在区间 若函数 积分的第一中值定理： 回忆： a b x (x-a)(x-b) ≤0 梯形求积公式余项定理证明 抛物线插值余项： 是否足够？ 抛物线公式的代数精度为3。 注意 在这里，f(x)可以为任意函数。 （Hermite插值） 表