数值第2章插值程序.ppt

例 给定f(x)在等距节点上的函数值表如下: xi 0.4 0.6 0.8 1.0 f(xi) 1.5 1.8 2.2 2.8 分别用Newton向前和向后差分公式求f(0.5)及f(0.9)的近似值. 解 先构造向前差分表如下: xi fi △fi △2fi △3fi 0.4 1.5 0.6 1.8 0.3 0.8 2.2 0.4 0.1 1.0 2.8 0.6 0.2 0.1 x0=0.4, h=0.2, x3=1.0. 分别用差分表中对角线上的值和最后一行的值,得Newton向前和向后插值公式如下: 由此可以提出一种新的法则为：当速率介于为 英里/小时之间时遵循“2秒法则”， 当速率介于为 英里/小时之间时遵循“3秒法则”， 当速率介于为 英里/小时之间时遵循“4秒法则”。 历史与注记 1756年，拉格朗日被任命为普鲁士科学院通讯院士。1766年任普鲁士科学院数学部主任。1786年出任法国米制委员会主任。1795年拉格朗日被选为法兰西研究院科学院数理委员会主席。1813年4月3日，拿破仑授予他帝国大十字勋章。 拉格朗日在数学上最突出的贡献是使数学分析与几何与力学脱离开来，使数学的独立性更为清楚。他在数值计算上的主要贡献是发展了欧拉所开创的变分法，为变分法奠定了理论基础。近百余年来，数学领域的许多新成就都可以直接或间接地溯源于拉格朗日的工作，在数学史上被认为是对数学的发展产生全面影响的数学家之一。 拉格朗日(Joseph-Louis Lagrange，1736～1813) 拉格朗日1736年生于意大利都灵，1813年卒于巴黎。 在近代，插值法是观测数据处理和函数制表所常用的工具，又是导出其它许多 数值方法（如数值积分、非线性方程求根、微分方程数值解等）的依据。插值法的一般参考资料见文献[1,2]，关于样条的一篇有重要影响的论文参见文献[3]。 插值一词最早是由Wallis提出的，公元6世纪，中国刘焯已将等距二次插值用于天文计算。17世纪，牛顿和格雷果里建立了等距节点上的插值公式。18世纪，拉格朗日给出了更一般的非等距节点上的插值公式。1946年, Schoenberg 首先提出了样条插值函数。 [1] P. J. Davies. The Finite Element Method: A First Approach. Oxford University Press, New York, 1980. [2] M. J. D. Powell. Approximation Theory and Methods.Cambridge University Press, New York, 1981. [3] I. J. Schoenberg. Contributions to the problem of approximation of equidistant data by analytic functions. Quarterly of Applied Mathematics 4, 1946. 参考文献 ? 问 题 ? 分段低次插值虽然具有简单、收敛性、整体连续性及数值计算的稳定性等优点，但在节点处常有“尖点”出现，光滑性较差。特别是需要给出节点处的导数值，这在多数问题中是不实际的。如何在没有节点导数数据时也能达到上述目的？为此引入样条插值函数。 七、三次样条插值 1.引入 定义：设对y = f (x)在区间[a, b]上给定一组节点 a = x0 x1 x2 … xn = b和相应的函数值y0, y1,…, yn， 如果s(x)具有如下性质： （1）在每个子区间[xi-1, xi] (i = 1, 2,…, n)上s (x)是不高 于三次的多项式; （2）s (x)，s ’ (x)，s? (x)在 [a, b]上连续；则称s (x)为 三次样条函数.如再有 （3） s (xi ) = f (xi) (i = 0, 1, 2,…, n)， 则称