清华第五版数值分析第4章ppt课件.pptVIP

Simpson公式的余项 定理　若f(x)? C4[a,b] ,则Simpson公式余项为 证明：1）确定代数精度为3. 2）令 3） 柯特思公式余项 例 用梯形公式和Simpson公式计算积分 ,并估计误差. 解 记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 复合求积法 当积分区间较大时，直接使用牛顿-柯特斯公式所得积分近似值的精度是很难得到保证的。因此在实际应用中，为了既能提高结果的精度，又使算法简便且易在计算机上实现，往往采用复合求积的方法。 将区间[a,b]适当分割成若干个子区间，对每个子区间使用低阶求积公式，构成所谓的复合求积公式，这是提高积分精度的一个常用的方法。 复合梯形公式 将区间[a, b]分成n等份，其分点为xi=a+ih (i=0,1,2,…n), 步长h=(b-a)/n. 在每个小区间[xk-1, xk] (k=1,2,…n)上利 用梯形公式 则 称 为复合梯形公式 复合梯形公式的误差 定理：若 ,则 证 ( ) 由于f(x)?C2[a,b] ，利用闭区间上连续函数的介值定理知存在一点??[a,b],使 因为求积系数为正，所以复合梯形公式稳定。 复合辛甫生公式 复合Simpson公式 复合Simpson公式的截断误差 定理 若f(x)? C4[a,b] ,则 可以看出，复合Simpson公式是收敛的。 证明 ( ) ( ) 实际上，只要f(x)可积，就可得到收敛性。 复合辛普森公式计算稳定。 例：用n=8的复合梯形公式和复合Simpson公式计算 解： = 3.138988494 = 3.141592502 结论：相同节点个数时，辛甫生求积公式的精度更高 第三节 龙贝格(Romberg)算法 一. 梯形法的步长逐次分半算法 将[a,b]分成n等分,有复合梯形公式 将[a,b]分成2n等分,有复合梯形公式 即 由以上递推公式可以看出，在已经算出Tn的基础上再计算T2n时，只要计算n个新分点上的函数值就行了。与直接利用复合梯形公式求T2n相比较，计算工作量几乎节省了一半。 为编程方便,常取 ,则有 本题考查学生综合知识解决问题的能力。最早研制成功的火车、飞机分别是在第一、二次工业革命中。火车动力来源于蒸汽机，飞机动力来源于内燃机 第四章 数值积分和数值微分 为什么要数值积分？ 要求被积函数f(x) ? 有解析表达式； ? f(x)的原函数F(x)为初等函数． 问题 1) f(x)没有解析表达式，只有数表形式 e.g. x 1 2 3 4 5 f(x) 4 4.5 6 8 8.5 2) f(x)有表达式，但原函数不是初等函数 e.g.，　　　　　　　　 它们的原函数都不是初等函数。 2 x e dx - ò 解决办法 1) 我们用不同的办法近似 可得到不同的积分公式。 2）用简单曲线的积分代替复杂曲线的积分。 求积公式举例 1）梯形公式 2）中矩形公式 3）一般公式 求积节点， 求积系数，也称为节点 的权，权仅与节点的选取有关，不依赖于被积函数的具体形式。 这种数值积分方法称为机械求积，特点是积分问题转变为被积函数值的计算。 求积公式的代数精度 定义　若求积公式对所有次数不超过 m 的代数多项式都精确成立，而对于某个 m+1次多项式不能精确成立，则称此求积公式具有m次代数精度. 上述定义等价于：若求积公式对 f(x)=1,x,x2,…,xm 均精确成立，而对 f(x)=xm+1 不精确成立,则称此求积公式具有 m 次代数精度. 求积公式的代数精度概念是衡量公式逼近好坏的标准之一. 例1 判别下列求积公式的代数精度 例2 试确定一个具有3次代数精度