计算方法课件第五章插值法.pptVIP

3.差商的计算第30页，共72页，星期日，2025年，2月5日根据插值多项式的存在唯一性知，如果f(x)充分光滑，则有估计不足：对函数的光滑性要求高；需估计导函数的最值；偏保守。导数型误差估计三、Newton插值余项第31页，共72页，星期日，2025年，2月5日差商型误差估计导数和差商的关系差商型误差估计特点：对被插值函数光滑性要求不高；但不适用于实际计算。第32页，共72页，星期日，2025年，2月5日四、例题解1）建立差商表1.01.52.00.84150.99750.90930.312-0.1764-0.48842）插值第33页，共72页，星期日，2025年，2月5日Newton插值多项式适用于节点任意分布的情形。但当节点等距分布时，可以简化Newton插值公式。§5.3等距节点插值设a=x0x1…xn=b，yi=f(xi)为等距节点xi=x0+h(i=0,1,…,n)上的函数值，其中h=(b-a)/n称为步长。在此基础上我们先定义差分，用差分表示Newton插值多项式，从而得到等距节点的插值公式。第34页，共72页，星期日，2025年，2月5日一、差分的定义与性质定义：称?yi=yi+1-yi (i=0,1,…,n-1)为f(x)在xi处以h为步长的一阶向前差分。?2yi＝?yi＋1-?yi=yi+2-2yi+1+yi(i=0,1,…,n-2)称为f(x)在xi处以h为步长的二阶向前差分。一般地，?myi＝?m-1yi＋1-?m-1yi(i=0,1,…,n-m)称为f(x)在xi处以h为步长的m阶向前差分。第35页，共72页，星期日，2025年，2月5日差分的性质性质1：各阶差分可用函数值线性表示，其计算公式为：其中性质2：差分与差商满足下述关系：证明：利用数学归纳法当k＝1时，有即结论成立。第36页，共72页，星期日，2025年，2月5日设k＝m-1时结论成立，即则当k＝m时，有由数学归纳法知，结论成立。证毕第37页，共72页，星期日，2025年，2月5日Remark：类似地可以定义向后差分与中心差分：性质3：差分与导数满足关系：证明：利用差商与导数、差分的关系，有：证毕第38页，共72页，星期日，2025年，2月5日二、Newton向前插值公式令x=x0+th，由xi=x0+ih(i=0,1,…,n)得：x-xi=(t-i)h，则有:将差商与差分的关系式带入Newton插值多项式，得:第39页，共72页，星期日，2025年，2月5日从而可得Newton向前插值多项式及其余项为：第40页，共72页，星期日，2025年，2月5日三、差分表Newton向前插值公式，又称表初公式，它利用差分表的最上面一个斜行的数值进行计算。第41页，共72页，星期日，2025年，2月5日四、例题解第42页，共72页，星期日，2025年，2月5日#第43页，共72页，星期日，2025年，2月5日五、Newton向后插值公式类似于向前差分，也可以得到差商与向后差分的关系：将插值节点从大到小排列，即类似于向前插值公式，可得到Newton向后插值公式，又称表末公式，它利用差分表的最下面一个斜行的数值进行计算。同样，还可以利用中心差分，构造插值公式，称为贝塞尔（Bessel）插值公式。第44页，共72页，星期日，2025年，2月5日这一类插值问题为埃尔米特(Hermite)插值问题。其几何意义是在插值点上插值曲线与被插值曲线有公共切线。由这2n+2个条件可以唯一确定一个2n+1次的插值多项式。具体我们采用基函数的方法来确定。§5.4埃尔米特插值一、问题第45页，共72页，星期日，2025年，2月5日1.辅助问题及Hermit插值二、一般情形第46页，共72页，星期日，2025年，2月5日2.辅助问题的求解第47页，共72页，星期日，2025年，2月5日第48页，共72页，星期日，2025年，2月5日3.Hermite插值问题解的存在唯一性存在性：唯一性：0第49页，共72页，星期日，2025年，2月5日4.插值余项分析：第50页，共72页，星期日，2025年，2月5日定理证明函数零点（从小到大）…………至少2n+1个零点……至少1个零点ξ证毕第51页，共72页，星期日，2025年，2月5日计算方法课件第五章插值法第1页，共72页，星期日，2025年，2月5日1函数表达式过于复杂不便于计算,而又需要计算许多点处的