数值逼近实习报告一(Lagrange).docxVIP

数 值 逼 近实验报告班级姓名指导教师实习题目给定（，），,构造插值多项式,使得=，目的意义掌握Lagrange插值多项式的构造方法数学模型 已知+1组离散数据 按照二次Lagrange插值多项式的构造方法，令：将插值条件=，=，=代入，容易得到： 可得到公式算法描述Step1: 输入n值及（，），,Step2: For ,对于给定x，计算Lagrange基函数Step3:计算（流程图）程序#includestdio.h#define MAX_N 20int main(){ printf(******lagrange插值计算公式******\n);double x[MAX_N]={0},y[MAX_N]={0};int i,j;int n; //定义插入节点的个数为nprintf(请输入插入节点的个数:\n);printf(n=);scanf(%d,n);if(nMAX_N){printf(输入的结点个数太大，请重新设置MAX_N的值:\n);return 1;}if(n=0){printf(输入的值应从1到一个数%d\n,MAX_N);return 1;}printf(请输入已知点的值(Xi Yi):\n); //在Xi、Yi中存放变量 double xx,lagrange=0; //这里xx表示要插入的点的x坐标printf(请输入已知点的x坐标:\n);for(int k=0;k=n;k++)scanf(%lf,x[k]);printf(请输入已知点的y坐标;\n);for(int m=0;m=n;m++)scanf(%lf,y[m]);printf(请输入插值点的x值:\n); //在m中存放要去求的未知变量printf(x=);scanf(%lf,xx);double t=1; //此处定义一个连乘计算的标记值for(i=0;i=n;i++)//利用Lagrange插值公式求其结果{for(j=0,t=1;j=n;j++){if(j==i)continue;t=t*(xx-x[j])/(x[i]-x[j]);}lagrange=lagrange+t*y[i];}printf(lagrange=%lf\n,lagrange);}数值算例给定，构造Lagrange函数，计算. 对计算结果进行分析评价误差估计式，插值多项式及其误差估计这里.当时误差为当时如果存在，则可以估计误差限：参考文献[1] ***，数值逼近，***：****大学***学院，20**