拉格朗日插值和牛顿插值多项式的C程序算法.docVIP

本科生毕业论文 题 目: 专业代码: 作者姓名: 学 号: 单 位: 指导教师: 2012年 月 20 日 拉格朗日插值多项式的C程序算法 1 1引言 1 1.1插值问题的提出 1 1.2插值法 2 1.3插值法思想 2 2拉格朗日插值法 3 2.1拉格朗日插值法的由来 3 2.2n次插值基函数 4 2.3拉格朗日插值多项式 4 3牛顿插值法 5 3.1均差： 5 3.2牛顿插值多项式： 5 4C程序设计 6 4.1算法设计： 7 4.2程序源码编写 7 5程序检测 12 5.1对拉格朗日插值的检测 12 5.2对牛顿插值的检测 13 总结 15 参考文献 16 致谢 17 摘 Abstract This article discuss the method to calculate Lagrange interpolation and Newton interpolation with C program. Base on the results of predecessors research, firstly, this article introduces the thoughts and methods of Lagrange interpolation and Newton interpolation. Improving the old method by adding functions which can repeatedly computing interpolation and correct illegal data. Then spreading it and making it more reasonable and perfect, checking it with some examples. Finally, summing up the main results of this article and application prospect. Key words：Lagrange interpolation; Newton interpolation ; C program;拉格朗日插值多项式的C程序算法 1引言 插值法是一种古老的数学研究方法，他的产生来自与社会的生产实践活动。在我国，早在一千多年前的隋唐时期，制定历法时，就应用了二次插值的方法。隋朝刘焯将等距节点二次插值应用于天文计算。但是，终究没有形成系统的理论。插值理论都是10世纪微积分产生以后渐渐发展起来的。拉格朗日插值和牛顿插值都是优秀的重要研究成果。《数值分析》[1]对此作了详细介绍，最近50多年来计算机技术的飞速发展和广泛应用，以及轻重工业等各方面实际问题的需要，促使插值法得到了更进一步的发展。 之前也有不少关于拉格朗日插值和牛顿插值的C程序算法，但是，经过实际运用发现都有各种各样的缺点，主要分为以下两种： 1、每次只能执行一次，算完一次之后，就会出现“press anykey to continue”，从而没法在进行下一次的计算； 2、没有纠错功能，通常情况下，为了计算的精确，我们这一个程序一般只用于计算20组以内的（即不超过20个节点的），当超过之后，会产生较大误差，甚至用户输入负组数之后，程序崩溃，即程序的健壮性没有设计好； 本算法在尽量弥补这两个不足的同时，也注意尽量优化程序，使占用的资源和运算的时间不会明显增加。 1.1插值问题的提出 在实际生活中，我们常用来表示某种内在的数量关系，其中很多数据可以通过实验或观测得到。这样虽然在给定的区间上是存在的，但是也仅仅能够得到上的一系列点的函数值。但这也只能刻画有限的情形。 为了研究函数整体的变化规律，以及实际的需要，我们往往要求出不再上的情形。因此我们常常会试着找出一个既能方便运算，又能和比较接近的函数。为了计算的方便，我们一般选一类比较简单的函数作为，使得满足=。这样确定的函数就是我们想要得到的插值函数。 1.2插值法 《拉格朗日平均插值法》[2]介绍，当一些实际问题用数学函数关系来描述时,往往没有明显的解析表达式,只能根据实验观测或其他途径提供一些离散点处的函数值和导数,有时尽管有表达式,却比较复杂,不便于研究和使用对此,人们希望构造一个简单的连续函数p(x)来近似替代所考察的函数f(x),使问题得到简化+1个点的函数表以后，我们要构造一个函数，这样的满足两个条件：第一，是一个不超过次的多项式；第二，在给定的点（）上的值必须与给定的值相同。 设函数在区间上有定义，且已知在点上