数值研究报告Lagrange差值和牛顿插值.docVIP

实验一 一、实验名称 Lagrange插值多项式和牛顿插值多项式 二、实验目的与要求： 实验目的：掌握Lagrange插值多项式和牛顿插值多项式的算法。 实验要求：1. 给出Lagrange插值和牛顿插值算法思路， 2. 用C语言实现算法，运行环境为Microsoft Visual C++， 3. 计算误差（这里只要求给出（-5,5）内101个点的误差）。 三、实验内容： 1. 对Lagrange插值多项式算法作编程练习和上机运算， 2. 对牛顿插值多项式算法作编程练习和上机运算， 3. 比较两种方法。 算法思路： 1. Lagrange算法是把多项式p写成如下形式： ， 其中称为Lagrange基函数。 计算Lagrange基函数的方法： fx=0.0; for(i=0;i=n;i++){ tmp=1.0; for(j=0;ji;j++) tmp=tmp*(x-x[j])/(x[i]-x[j]); for(j=i+1;jn;j++) tmp=tmp*(x-x[j])/(x[i]-x[j]); fx=fx+tmp*y[i]; } return(fx); 2. 牛顿算法是把多项式p写成如下形式： 其中是插值点，是待定系数。可以通过插值点和插值点处的函数值算出待定系数，方法如下： (1) (2) 将(1) - (2)并利用，得