Lagrange插值及Newton插值.docxVIP

PAGE PAGE 1 实验报告 实验项目 插值法 实验日期 2016/9/30 理论内容 Lagrange插值与Newton插值 授课日期 实验室名称 文理管203 微机编号 E1 实验目的及要求： 1、了解多项式差值公式的存在唯一性条件及其余项表达式的推导。 2、了解拉格朗日插值多项式的构造、计算及其基函数的特点，牛顿插值多项式的构造与应用，差商、差分的计算及基本性质。 实验内容： 编写Lagrange插值法及Newton插值法通用子程序，依据数据表 0.32 0.34 0.36 0.314567 0.333487 0.352274 构造一个抛物插值多项式及，计算的近似值并估计误差。 实验步骤及程序： Lagrange插值公式算法流程图 开始 开始 输入x 0?y 0?k 1?t x- y+t? k=0? k+1?k 输出y 结束 ≠ Newton插值公式算法流程图 调用 调用函数ChaShang（） f=0;temp=0; f=f+temp i 返回i阶差商f 调用函数Newton（） return=0 I阶差商f=ChaShang（i,X,Y） temp=1 temp=tremp* result=result+f*temp iX.Sie() 返回结果result i=i+1 i=i+1 ≥ ≥ 2、Lagrange插值源程序 import java.util.Scanner; /*拉格朗日插值*/ public class Lagrange_interpolation { /*拉格朗日插值法*/ private static double[] Lag_method(double X[], double Y[], double X0[]) { int m = X.length; int n = X0.length; double Y0[] = new double[n]; for (int i1 = 0; i1 n; i1++) {//遍历X0 double t = 0; for (int i2 = 0; i2 m; i2++) {//遍历Y double u = 1; for (int i3 = 0; i3 m; i3++) {//遍历X if (i2 != i3) { u = u * (X0[i1] - X[i3]) / (X[i2] - X[i3]); } } u = u * Y[i2]; t = t + u; } Y0[i1] = t; } return Y0; } public static void main(String[] args) { /*输入插值点横纵坐标*/ System.out.println(Input number of interpolation point:); Scanner scan = new Scanner(System.in); int m = scan.nextInt(); System.out.println(Input number of test point:); int n = scan.nextInt(); double X[] = new double[m]; double Y[] = new double[m]; double X0[] = new double[n]; System.out.println(Input the elements of X:);//已知插值点 for (int i = 0; i m; i++) { X[i] = scan.nextDouble(); } System.out.println(Input the elements of Y:);//已知插值点的函数值 for (int i = 0; i m; i++) { Y[i] = scan.nextDouble(); } System.out.println(Input