h11：数值计算方法之数值微分与外推方法.pptVIP

2.涅维尔插值的递推格式 定理：对于前面给出的问题及记号，我们有 证明：不难验证，对m=0及任意k值，(*)式均成立。 对k-mjk+1，由于 所以有pk+1,m+1(xj)= pk+1,m(xj)= pk,m(xj)=yj。 显然， pk+1,m+1(xk+1)= pk+1,m(xk+1)=yk+1; pk+1,m+1(xk-m)= pk,m(xk+m)=yk-m; 所以， pk+1,m+1(x)满足m+2个插值条件。 3.几点具体说明 实际数值计算时，应采用教材第313页(4.4)式计算，可节省一点计算时间，程序稍简单一点。 无论是手工列表计算还是编程由计算机计算，均可仿照上一节中的里查逊外推法的实际计算方法。事实上，如果我们把涅维尔插值方法应用于外推，里查逊外推方法便成为这里的特例。 在工程应用中，可以先用涅维尔插值方法计算，并根数值结果得出一个比较好的插值多项式的阶数，在实际应用中再采用具体的方法编程。这样一来，程序显得更简洁一些。实际使用里查逊外推方法也可采用这种模式。 在后面的数值微分领域，寻找比较好的外推方法是很有意义的，涅维尔插值方法为此提供了基础方法。 4.利用涅维尔插值方法构造外推格式 在涅维尔递推插值公式中，如果诸xk是一个单调，不妨假设是单调减，收敛于x的序列，这时相应的计算问题便是一个外推问题。 不失一般性，可在涅维尔递推插值公式(1)中令xk-x=hk,此时，{hk}是单调收敛于零的序列。再 把公式中的诸pk,m,(x)改写为pk,m,则有 在进一步讨论上面的公式的用法之前，我们进一步寻找偶函数在原点处的外推格式。 偶函数在原点处的外推格式 与数值微分方法类似，大量数值计算问题均可等价地变换为下面的“步长”收敛于零的极限问题： 而且，在一般情况下还可以假定F(x0,h)是关于h的偶函数。 偶函数在原点处的泰勒展式只含偶次幂的项，简单的处理方法是把前面的(2)式中所有hk改写为(hk)2,从而有 其中，pk,0=F(x0,hk)。 6. 里查逊外推格式的推导 为了说明如何利用前面给出的外推格式，首先我们利用上面的(3)式来推导里查逊外推格式。 仍假定h0为初始步长，令hk=h0/2k,k=0,1,…,n，即可得到： 11.5 实用外推加速方法讨论 对于寻找某个二元函数F(x0,h)当h趋近于+0市的极限问题，上面一节中已经给出了3个外推格式，首先熟悉一下他们各自的特点是很有好处的： (2)式是最一般性的公式，适用于所有场合； 如果F(x0,h)是关于h的偶函数，虽然也能使用(2)式，但使用(3)式的效果更好一些，包括计算量更少一些，数值性能更好一些，所能得到的最好结果会更精确一些。 (4)式是(3)式的特例，按每次步长减半的规则确定步长序列，所以程序最为简单。 大家也可以在(2)式中按每次步长减半的规则确定步长序列，也可以得到类似于(4)式得简单结果。 (2)式和(3)式的意义在于，我们对外推方法的研究，可以专注于寻找更好的步长序列，或者构造更好的步长序列的策略。 利用等比序列外推加速 我们最容易想到的构造步长序列的方法是产生一个单调收敛于零的等比数列作为步长序列，这样处理的确有很多好处： 方法简单，取0q1，对适当给定的h0，取hk=h0qk,k=0,1,…，即可搞定。 由于几何序列具有无记忆性，与初始步长h0的选取没有关系，在数值计算领域，这是难得的优点。 具体做法： 取q=1/r，则有hk=h0/rk，一般情况下r取1.2~2之间的数为宜，权衡各种细节因素，取r=1.6是比较理想的。 把诸hk分别代入到(2)式和(3)式中，化简，即可得到具体的外推格式。 2.利用单调增的整数序列外推加速 设{an}是一个单调增的整数序列，对于给定的h00，记hk=h0/ak,则{hk}是单调减收敛于零的序列。把他们分别代入到(2)式和(3)式中，也可以得到相应的外推格式。 利用整数序列外推加速的好处是方法简单，可选择的范围更大，所以更容易寻找具体的外推格式。 最简单的整数序列可采用等差数列，由于首项、公差、以及系列的长度均能对算法性能产生影响，再加之与初始步长有密切关系，不宜作为一般性的方法使用。 利用斐波拉奇数列外推也许是个好主意：构造简单，从1、2、3开始，以后每一项都是前面相邻两项之和，相邻两项之比大体为1.6，与初始步长的选取基本没有关系。 11.6 二阶导数的计算方法 对于求f(x)在x0处的二阶导数来说，也可以采用基本的处理方法：构造一个与步长有关的近似公式F(x0,h)，最好是关于步长的偶函数；确定适当的初始步长；构造变步长系列；利用变步长系列外推得到加速收敛序列。 对于求f(x)在x0处的二阶导数