Newton插值法的程序设计与应用.doc

引言 在我们生活中，许多实际问题都都可以利用函数来表示某种内在规律的数量关系，但是在很多应用领域，函数有时不能直接写出表达式或者过于复杂不易计算，而只能给出函数在若干个点上的函数值或导数值，通常也是造一张函数表，当遇到要求表中未列出的变量的函数值时，就必须做数值逼近．例如给定了函数在几个特定点上的函数值，为了研究函数的变化规律，就要根据这个表，将其公式化寻求某一函数去逼近，在给定点上等于函数值，在其他点上约等于函数值，这样既能反映的函数特性，又便于计算．称的插值函数，为被插值函数. 我们可以求一个待定函数来近似反映函数的特性，使得待定函数在给定点上等于函数值，在其它点上的函数的值作为函数的近似值，这种方法称为插值法.利用差值函数很容易得到Lagrange插值多项式，公式结构紧凑，在理论分析中甚为方便.由于公式中的每一项与所有的插值结点有关.因此，如果需要增加一个插值结点，则Lagrange插值公式中的每一项都需要改变计算量大，为了克服这一缺点，于是产生了Newton插值法.有的要求更高需要在某些点处的导数也相同于是前两种方法达不到精度于是产生了Hermite插值法. 本题主要讨论牛顿插值多项式.主要是对牛顿插值法的理论总结，利用C++程序和Matlab语言程序编写，以及对其应用讨论. 第一章 插值方法 求得近似函数的方法通常有两种：一种是插值法，另一种是曲线拟合法.本章主要讨论插值法用多项式来逼近列表函数问题，即具有唯一插值函数的多项式插值,对其中的多项式插值主要讨论n次多项式插值的方法，即给定n+1各点处的函数值后，怎样构造一个n次插值多项式的方法. §1.1插值问题的由来 在许多实际问题及科学研究中，因素之间往往存在函数关系.然而这种关系经常很难有明显的解析表达式，通常只是由观察与测试得到一些离散数值.有时由于表达式过于复杂，不仅使用不便，而且不易于进行计算与理论分析.解决此类问题的方法有两种：一种是插值法，另一种是拟合法. 插值法是一种古老的数学方法，他来自生产实践，早在一千多年前，我国科学家在研究历法上就应用了线性插值与二次插值，但它的基本理论却是在微积分产生之后才逐渐完善的，其应用也逐步增多.特别是在计算机软件中，许三多库函数.等的计算实际上归结于它的逼近函数的计算.在工程的实际问题中，由于航空、造船、精密机械加工等实际问题的需要，插值方法在实践上和理论上显得更为重要，并得到了进一步的发展和广泛的应用.因此我们希望能够得到一个“简单函数”逼近被计算函数的函数值.于是就得到了一种方法.叫做插值逼近或者插值法. §1.2 几个基本概念 §1.2.1 插值函数的概念及相关性质 定义1.1 设连续函数在区间上有定义，已知在个互异的点上取值分别为（设）.若在函数类中存在以简单函数，使得（）,则称为的插值函数. 称为插值节点，称为插值区间. 定义1.2 设如下表 求一个，使其满足： ， 则称为的n次代数插值多项式，该问题称为n次代数插值问题. 定理1.1 n次代数插值问题的解存在且唯一 . §1.2.2 Lagrange插值函数的概念及相关性质 定义1.3 若n次多项式（）在个互异节点上满足 则称（）为节点上的n次插值基函数. n次差值多项式为 称之为拉格朗日多项式. 定理1.2 设，则余项为： 其中. §1.2.3 差商的概念及相关性质 定义1.4 称函数于点（）上的平均变化率为在处的一阶差商，记作，即 ，. 一阶差商的差商为在点处的二阶差商.记作，即 一般地，在求出的阶差商后就可以构造的阶差商.称 为在处的阶差商. 特别地，规定零阶差商. 性质1.1 任意调换插值节点的次序其值不变 即 性质1.2 差商是一个数，可表示为某个函数的线性组合. 性质1.3 当的阶导数在包含节点的某个区间上存在时，在 之间至少存在一点使 注1.1 由性质2可以推出，次多项式的一阶差商是多项式. 表1-1差商表 一阶差商 二阶差商 三阶差商 四阶差商 §1.2.3 Newton插值函数与差分的概念及相关性质 定义1.5 我们定义牛顿插值公式为 （1.1） 定义1.6 设函数在等距节点上的函数值 ，其中为常数，称为步长.称函数在每个小区间上的