8数值积分与数值微分要点.pptVIP

* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * 用归纳法容易证明，由(4)得到的Fm(h)逼近F*的误差为 上面的这种方法，称为Richardson外推法。 Romberg求积公式 注： 外推法不只用在积分的计算，也可用于计算其他量，只要能把该量看成f(h)当 h? 0时的极限为f(0),而f(0)不易直接计算，要通过{hi}序列对应的函数值序列{f(hi)}来外推。 如果f(h)有类似的展式，即 f(h)= 则一切的算法和讨论可类似进行，更一般地设为 f(h)= 其中 0r1r2···rm，ri不必是整数，?I是不依赖于h的常数，?0= limf(h)，可以类似地推出外推计算公式。 例：单位圆内接正n边形的面积为 S= sin 设h=2/n,当n ? ?有h ?0可以类似地推出外推计算公式。 解:将S与h的关系写成 S(h)= 。 当h ? 0时,S(h)的极限是?，即单位圆的面积。将S(h)作Taylor展开： ··· 即S(h)有如上展式的形式，注意?是无理数，所以可按Romberg算法计算积分的过程来计算?到任意精度。 取n=6,12,24,即h=1/3,1/6,1/12,这时S(h)值很容易计算。经过计算S(h)及两次外推法，可计算出S(0)= ?的近似值。 取更多的h值外推过程还可继续下去。我国古代数学家刘徽在公元263年就是用这种“割圆术”加上特殊的外推技巧，计算得??≈3.1416。 例1：取e=0.00001,用龙贝格方法计算积分 I=? dx 4 1+x2 0 1 例2:分别用不同方法计算如下积分,并做比较 令I= 各种做法比较如下： （1）、Newton-Cotes公式 当n=1时：即用梯形公式，I=0.9270354 当n=2时：即用Simpson公式，I=0.9461359 当n=3时：I=0.9461090 当n=4时：I=0.9460830 当n=5时：I=0.9460831 （2）用复化梯形公式： 令h=1/8=0.125 （3）用复化抛物线 令h=1/8=0.125 （4） Romberg公式 K T1 T2 T3 T4 0 0.9207355 1 0.9397933 0.9461459 2 0.9445135 0.9460869 0.9400830 3 0.9456906 0.9460833 0.9460831 0.9460831 比较 此例题的精确值为0.9460831... 由例题的各种算法可知： 对Newton-cotes公式，当n=1时只有1位有效数字，当n=2时有3位有效数字，当n=5时有7位有效数字。 对复化梯形公式有2位有效数字，对复化抛物线公式有6位有效数字。 用复合梯形公式，对积分区间[0，1]二分了11次用2049个函数值，才可得到7位准确数字。 用Romberg公式对区间二分3次，用了9个函数值，得到同样的结果。 外推法在数值微分中的应用 总结 1：梯形求积公式和抛物线求积公式是低精度的方法，但对于光滑性较差的函数有