《龙贝格求积》PPT课件

《《龙贝格求积》PPT课件》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、第7.4-7.5节龙贝格求积公式龙贝格求积公式复合求积公式随着n的增加可以减少积分误差，但高阶N-C公式又会造成数值不稳定，因而采用复合求积公式。复化梯形公式复化辛普森公式上面的求积公式都是定长的，要达到某个精度，则必须选取适当的长度，但是这是一件不容易达到的事情。复合求积公式的截断误差随着n的增大而减少。但每给一种积分方法之后，如何选择n，使精度达到预先选定的精度？1用误差估计式子。但是要求高阶导数，一般实比较困难的。2在实际中一般采用自动选择积分步长。即在求积的过程中，将步长逐步折半，反复利用复合求积公式，直到相邻两次的计算结果之差小于容许的范围。梯形法的逐次分

2、半算法逐次分半算法。。。。。。。。。。。。。解例1例题1例2解例题2例题2在等距节点的情况下，用计算机计算积分值通常都采用把区间逐次分半的方法进行。这样，前一次分割得到的函数值在分半以后仍可被利用,且易于编程。逐次分半算法把区间二等分，每个小区间长度为h/2=(b-a）/2，于是T2=T1/2+[h/2]f(a+h/2）把区间四(22)等分，每个小区间长度为h/22=（b-a）/4，于是T4=T2/2+[h/22][f(a+h/4)+f(a+3h/4)]················把[a,b]2k等分，分点xi=a+(b-a)/2k·i（i=0，1，2···2k

3、）每个小区间长度为(b-a)/2k，由归纳法可得逐次分半算法复合求积方法是用于被积函数变化不太大的积分.如果在求积区间中被积函数变化很大，有的部分函数值变化剧烈，另一部分变化平缓，这时统一将区间等份用复合求积公式计算工作量就会很大.要达到误差要求对变化剧烈部分必须将区间细分，而平缓部分则可用大步长，即针对被积函数在区间上不同情形采用不同的步长，使得在满足精度前提下积分计算的工作量尽可能小.7.4.2自适应simpson公式针对这类问题的算法技巧是在不同区间上预测被积函数变化的剧烈程度确定相应的步长.这种方法称为自适应积分方法.设给定精度要求，计算积分的近似值.先取步

4、长，应用辛普森公式有其中若把区间对分，步长，在每个小区间上用辛普森公式，则得（7.31）（7.32）实际上(7.32)即为与(7.31)比较，若在上变化不大，可假定其中（7.32）’从而可得若不等式(7.33)不成立，则应分别对子区间及进行误差分析，若每个区间误差是否近似在这里.如果有则可期望得到与(7.32)比较，则得此时可取作为的近似，则可达到给定的误差精度.(7.33)对满足要求的区间不再细分，对不满足要求的还要继续上述过程，直到满足要求为止，7.5、龙贝格算法事后估计法利用计算结果来估计误差的方法龙贝格算法当n=1时,我们计算上式右端这恰好是辛普森公式的结果

5、，即有比梯形公式有更好的精确度龙贝格算法类似地可验证：即龙贝格算法龙贝格算法可以验证事实上C1=(42S2-S1)/(42-1)=16S2/15-S1/15=(7y0+32y1+12y2+32y3+7y4)/90恰为柯斯特公式。同理，C2=(42S4-S2)/(42-1)，...龙贝格算法即,R1=(43C2-C1)/(43-1),R2=(43C4-C2)/(43-1),...Rn=(43C2n-Cn)/(43-1);上式即为龙贝格公式，得龙贝格值序列龙贝格算法用若干个积分近似值推算出更为精确的积分近似值的方法，称为外推方法。龙贝格算法计算过程见表：龙贝格算法T1T

6、2S1T4S2C1T8S4C2R1T16S8C4R2T32S16C8R4.................................................上面是Romberg的计算表若则计算停止龙贝格算法或龙贝格算法解：用Romberg方法计算积分近似值例3例题3例题3例题3