数值计算8.docVIP

数值计算实验报告 班级： 实验题目 (1)用单项式基函数对上述数据进行插值 (2)对上述数据进行Lagrange插值 (3) 对上述数据进行Newton插值 并画出对应的插值曲线。 实验目的 多项式插值 编程环境 采用C++或JAVA实现 实验内容（算法、程序、步骤和方法） 算法简述： 单项式基函数插值法： 对n个数据点，应取k=n-1，以使空间的维数与数据点的个数相匹配，次数不超过n-1次的多项式形成的向量空间的最直接的基底由前n个单项式组成 j=1,...n多项式为： 拉格朗日插值法： 对给定的数据点，i=1,...n,中的拉格朗日基函数为： j=1,...,n 数据点的多项式插值为： 牛顿插值法： 对给定的数据集，i=1，2，3，...，n 中的牛顿基函数为： j=1，2，...，n 当下标范围无意义时，取乘积的值为1，在牛顿基底下，多项式为： 程序代码如下： public class Interpolation { /** * 插值法拟合离散数据，包括： * 单项式基函数插值法 * 拉格朗日（Lagrange）插值法 * 牛顿（Newton）插值法 * */ /**单项式基函数插值法*/ //计算单项式基数函数的系数 public float[] monomial(float t[],float y[]){ float[] x = null; float[][]a=new float[t.length][t.length]; for(int i=0;ia.length;i++) for(int j=0;ja[i].length;j++) a[i][j]=(float)Math.pow(t[i], j); GaussElimination ge=new GaussElimination(); x=ge.gauss(a, y); return x; } //通过秦九韶算法计算单项式的值 protected float calculateMonomial(float[] a,float x){ float y=0; for(int i=a.length-1;i=0;i--) y=y*x+a[i]; return y; } /**Lagrange插值法*/ //计算Lagrange插值法的多项式的系数 public float[] lagrange(float t[],float y[]){ float[] x=y; return x; } //计算Lagrange插值法的多项式的值 public float calculateLagrange(float t[],float a[],float x){ float y=0; float z=1; for(int i=0;it.length;i++){ for(int j=0;jt.length;j++){ if(i!=j) z=z*((x-t[j])/(t[i]-t[j])); } y=y+a[i]*z; z=1; } return y; } /**Newton插值法*/ //计算Newton插值法求得的多项式的系数 public float[] newton(float t[],float y[]){ float[][] a=new float[t.length][t.length]; float[] x=null; for(int i=0;it.length;i++){ for(int j=0;jt.length;j++){ if(ji) a[i][j]=0; else if(j==0) a[i][j]=1; else{ a[i][j]=1; for(int n=0;nj;n++) a[i][j]=a[i][j]*(t[i]-t[n]); } } } GaussElimination gs=new GaussElimination(); x=gs.fs(a, y); return x; } //计算Newton插值法求得的多项式的值 public float calculateNewton(float t[],float a[],float x){ float y=0; for(int i=0;ia.length;i++){