二阶导数的五点数值微分公式.docxVIP

1二阶导数的五点数值微分公式卜面给出问题的条件：设f⑴为定义在区间［。,b］上的函数，给定f(%)在等距节点a%%%%%b处的函数值01234f(%k),(k=0,1,2,3,4),且%k+1-%k=h。 在区间［a,b］上作f(%)的4次Lagrange插值函数，将%=%0+th，te［0,4］，%*=%0+kh代入，并将方程两端对t求二次导数，再分别把t=0，1,2,3,4代入，即可得到%(k=0,1,2,3,4)节点二阶导数的5点数值微分公式 f“(%。)=总[35f(%0)-104f(气)+114f(%2)-56f(%3)+11f(%4)](1) f(%)=-1-I11f(%)-20f(%)+6f(%)+4f(%)-f(%)](2)112h2012341f”(%2)=诫[-f(%0)+16f(气)-30f(%2)+16f(%「f(%4)](3)1 f”(%3)=床Lf(%0)+4f(气)+6f(%2)+20f(%3)+11f(%4)](4)f”(%4)=总hlf(%0)-56f(%1)+114f(%2)-104f(%3)+35f(%4)](5)2二阶导数五点微分公式的外推算法2.1中点节点五点微分公式的外推算法分别将f(%)，f(%)，f(%)，f(%)，f(%)在％点作Tay-lor展开并代入(3)012342式的右端整理得到YYYYf”(%)一S(h)=ah4+ah6+ah+....(6)式中：S.k1表示%.节点处第k次外推公式-f(x)+16f(x)-30f(x)+16f(x)-f(x)S(h)=01234 12h2f(5)(x)f(7)(x)f9(x)以1902，以210082，以321600， 对于固定的x2，a.(i=1,2,...)是与h无关的常数，所以上面的误差估计式符合Richardson外推算法，将h缩小一倍得到(7)式： f(x)-S(h)=a(h)4+a(h)6+a(h+...(7)22.12i2i2i2由16x(7)-6式整理得到一hhh f(x)-S(h)=a(=)4+a侦)6+a侦￡+... TOC\o1-5\h\z22.212121216S(h)-S(h)式中：以=-以，以=-以，S(h)=——2^-2*2220231632215由此可见，x2点的五点微分公式外推一次后，精度由Oh4提高到Oh6依次类推，可得到外推算法的递推序列: vS-f(x)+16f(x)-30f(x)+16f(x)-f(x)〕 S(h)=01234-12h28)一h,22(k+1)s(3)-S(h) 8) S(h)=,(k=1,2,...)2.k+12k+2一1式中k为外推次数，S2k+1(h)的截断误差为O(h2(k+1))，既每外推一次收敛阶增长两阶。 2.2其他节点五点微分公式的外推算法对于其他节点x(i=0,1,2,3)，将f(x)(i=0,1,2,3,4)在x点处作Taylor展开并分别代入(1)(2)(3)(4)(5)式整理得到TOC\o1-5\h\z〃，5,〃119,〃43,〃797…八f、(x)-S(h)七h3f⑸(x)+h4f⑹(x)+—h5f()(x)+h6f⑻(x)+...(9)00.16090036010080一一1一一19一一1一一23一一一f(x)-S(h)=h3f⑸(x)+h4f⑹(x)+--h5f⑺(x}+——h6f⑻(x)+...(10)11.1121360136120161一一1一一19一一1一一23一一f”(x3)-S3.1(h)=12时⑸(x3)一疝⑹(x3)+话时⑺b一而h6产8)(x3)+...(11)“，5,“119,“43,“797,“一f(x4)-S41(h)=甘h3f(5)(x4)--90h4f(6)(x)+—h5f⑺(x)-仍彼h6/⑻也)+...(12)x(i=0,1,2,3,4)Oh3hx(i=0,1,3,4)由上式可知，节点x(i=0,1,2,3,4)的5点数值微分公式精度均为Oh3。分别将(9)至(12)式h缩小一倍乘以8减去(9)至(12)】〕[35f(x)-104f(x)+114f(x)-56f(x)+11f(x)]式即得到x(i=0,1,3,4)节点外推公式的递推算法TOC\o1-5\h\zS10.112h2~0-1~2-3S=-^-I11f(x)-20f(x)+6f(x)+4f(x)-f(x)]12h201234S=t—-[-f(x)+4f(x)+6f(x)-20f(x)+11f(x)]12h201234S=-^-\11f(x)-56f(x)+114f(x)-104f(x)+35f(x)]12h201234一h.2k+2s(-)-S(h)S(h)=——2,:2业_,(k=1,2,...)、2.k+12k+2-1(13)经过一次外推，点x,x