2011计算方法课程的设计汇报.docVIP

《计算方法》课程设计报告 学生姓名: 吴志强 学 号： 0909290124 学 院: 理学院 班 级: 信计091班 题 目: 拉格朗日插值及牛顿插值的通用程序 指导教师： 李鹏松 职称: 教 授 朱秀丽 讲 师 2011年 12 月 30 日 一、摘要 1 二、应用计算方法的基本原理 2 1.拉格朗日插值法 2 1.1算法描述 2 1.2插值余项与误差估计 2 2.牛顿插值法 2 2.1算法描述 2 三、例题的计算结果及误差分析 3 四、总结及心得体会 4 五、参考文献 7 六、附录程序 7 一、摘要 本设计报告主要围绕拉格朗日插值及牛顿插值多项式展开。 首先分别给出了拉格朗日和牛顿插值的基本原理。拉格朗日插值是在给出线性插值基函数的基础上得出插值公式： 利用插值基函数很容易求得拉格朗日插值多项式，公式结构紧凑，理论分析甚为方便，但是当差值节点增减时，全部差值基函数均要随之改变，给实际计算带来不便，而牛顿插值法很好的解决了这个问题。 由基本原理给出程序的算法描述并由此编写出通用程序。给出实例拉格朗日插值节点u=[0 1 4 9 16 25 36 49 64] ，带入程序得出结果：；牛顿插值节点x=[0.2 0.4 0.6 0.8 1.0]，y=[0.98,0.92,0.81,0.64,0.38]；代入程序得出结果： 通过所得出的插值多项式，代入具体数值，分别求出两种插值多项式的误差估计。并估算出误差分别为：则拉格朗日插值法绝对误差为0.0092；牛顿插值绝对误差0.000052。通过误差分析我们可以清楚的看到两种插值多项式的误差是可以接受的。 最后我们对本次课程设计进行了总结，总结了程序的优缺点并对本次试验过程中遇到的问题及困难进行了解答，此外我们还写出了对本次课程设计的心得体会。 二、应用计算方法的基本原理 拉格朗日插值法 1.1算法描述 拉格朗日插值多项式是一种最常见的多项式插值法，也是一种最常用的逼近工具。设是定义在区间[a,b]上的函数，又设，..... 是[a,b]上的n个互不相同的点，可以先构造基函数（i=0,1,2......n）如下： (1) 且有 (2) 显然 (3） 满足（3）的插值多项式可表示为 (4) 由的定义知，我们称为拉格朗日插值多项式。同时引入记号，我们容易求得，于是我们可将公式（4）改 写为。 1.2插值余项与误差估计 若在[a，b]上连续，在（a，b）内存在，节点， 是满足（3）的插值多项式，对任何，插余项，，余项表达式只有在f（x）高阶 导数存在时才能应用，在（a,b）具体位置不能具体给出，如果我们能求出，那么插值多项式的截断误差限为 2.牛顿插值法 2.1算法描述 利用插值基函数很容易求得拉格朗日插值多项式，公式结构紧凑，理论分析甚为方便，但是当差值节点增减时，全部差值基函数均要随之改变，给实际计算带来不便，而牛顿插值法很好的解决了这个问题。在这里我们先引入差商的概念：设有函数，为一系列互不相等的点，我们称为k阶差商这时我们可以求得 其中为牛顿插值的余项。其中为牛顿均差插值多项式，为插值余项。 三、例题的计算结果及误差分析 一．拉格朗日插值法 1.下列数据点的插值 x 1.5 11 28 30 56 101 132 y 6 7.8 2.5 3.2 7.6 6.3 7.5 试用拉格朗日插值法分别求出插值多项式 结果： 图像如下： 图一 例1拉格朗日插值多项式的图像 误差分析：由所求得的插值多项式 例：将x=1.5代入，求得L（4）=6.199968 精确解y（4）=6.2。则绝对误差为0.000012，相对误差为0.00016。 一．牛顿插值法 1.已知函数在下列各点的值为 x 1.0 2 3 4 5 y 0.98 0.92 0.81 0.64 0.38 结果为: 图像： 图二 例2牛顿插值多项式图像 误差为：0.00025 四、总结及心得体会 总结：这次试验使我进一步认识到算法的重要性，为进一步掌握编程能力提供了一次很好的锻炼机会。在学习其他语言的基础上，进一步加强编程能力。 （1）进一步熟悉了matlab编译环境。 （2）完成一个程序的要经过设计，输入 ，