拉格朗日插值公式数值分析实验报告

1、学生学号0实验课成绩或修俎7次浮学生实验报告书数值分析A理学院金升平教授陶玮统计1401实验课程名称 开课学院 指导教师姓名 学生姓名 学生专业班级20152016 学年 第 2 学期实验课程名称：数值分析实验项目名称拉格朗日插值公式实验成绩实验者陶玮专业班级统计1401组 别同组者实验日期年 月 日一部分：实验预习报:告( 包括实验目的、意义，实验基本原理与方法，主要仪器设备及耗材，实验方案与技术路线等 )实验目的：1、学习和掌握拉格朗日插值多项式。?2、运用拉格朗日插值多项式进行计算实验基本原理：拉格朗日插值基函数的一般形式：ugo = n彘j*心Xj也即是：/ _ (X - Xo)?(x

2、 - Kj, - i)(x - Xk 4 1)?(X - X)_(xk - Xo)?(Kk -叫-仇办-4 l)?(Xk - Xn)所以可以得出拉格朗日插值公式的一般形式：Pn(x)= 2： . olk(X)Vk其中，n=1时,称为线性插值，Pi(x)=y例式K)+ Vi?li(x)n=2时,称为二次插值或抛物插值，精确度相对高些,冉(x) = vo?loG) + yi?h(x) + ya? 12 W主要仪器设备：计算机，MATLAB实验内容：输入所求节点x和已知的节点数n,然后输出插值结果s实验方案与步骤：1. 输入 n 对数(xi,yi)(i=0-n )2. 令 s=03. 对于 i=0-

3、nT=yi对 j=0- n 但 j!=i T=T*(X-xj)/(xi-xj)令 s=s+T4. 、输出结果流程图：开始第二部分：实验过程t己录(可加页)(包括实验原始数据记录，实验现象记录，实验过程发现的问题等)#includemain()double X;int n;double lgrr(int n,double X);printf(input the aim numberX);printf(input the number of pairs of numbers n:);scanf(%lf%d,&X,&n);lgrr(n,X);return 0;double lgrr(int n,do

4、uble X)double s=0,t;double x81,y81;int i,j;printf(input n pairs of numbers:n);for(i=0;in;i+)scanf(%lf%lf,&xi,&yi);for(i=0;in;i+)t=yi;for(j=0;jn;j+)if(j!=i)t=t*(X-xj)/(xi-xj);s=s+t;printf(s=%fn,s);return 0;第三部分结果与讨论(可加页)一、实验结果分析(包括数据处理、实验现象分析、影响因素讨论、综合分析和结论等)二、小结、建议及体会三、思考题对于现在的许多实际问题来说，我们并不知道f(x)的具体形式，所对应的函数值可能是由测量仪器或其他设备中直接读出来的，f(x)只是一个数学概念意义下的函数。(比如：图像的方法处理，天气预报，机 床加工等方面)解答这类问题的方法就是插值方法。？泰勒插值要求提供f(x)在点X0处的各阶导数值，这项要求很苛刻， 函数f(x)的表达式必须相当简单才行。如果仅仅给出一系列节点上的函 数值f(x i)?=?yi(i=0,1,2,n),则插值问题可表述如下：求作？n?次多项 式？Pn(x),使满足条件Pn(x尸?yi, i?=?0, 1,，n?。这就是所谓拉