《拉格朗日插值公式数值分析实验报告》由会员分享，可在线阅读，更多相关《拉格朗日插值公式数值分析实验报告（6页珍藏版）》请在金锄头文库上搜索。
1、学生学号0121414670327实验课成绩学生实验报告书实验课程名称数值分析A开课学院理学院指导教师姓名金升平教授学生姓名陶玮学生专业班级统计140120152016 学年第 2 学期实验项目名称拉格朗日插值公式实验成绩实验者陶玮专业班级统计1401组别同组者实验日期年 月 日部分：实验预习报告（包括实验目的、意义，实验基本原理与方法，主要仪器设备及耗材，实验方案与技术路线等）实验目的：1、学习和掌握拉格朗日插值多项式。2、运用拉格朗日插值多项式进行计算。实验基本原理：拉格朗日插值基函数的一般形式：g) = nn=oXXjxk-xj也即是：l （x） =（xx0）（xxki）（xxk+i）（
2、xxn） k = 0,123k（xkx0）（xkxk-1）（xkxk+J（xkxn）所以可以得出拉格朗日插值公式的一般形式：Pn（x）=k=0lk（x）yk其中，n=1 时，称为线性插值，P1（x） = y0 l0（x） + l/x）n=2时，称为二次插值或抛物插值，精确度相对高些, p2（x）= y0 l0（x）+ y1 l1（x） + y2 l2 （x）主要仪器设备：计算机，MATLAB软件实验内容:输入所求节点X和已知的节点数n,然后输出插值结果s实验方案与步骤:1. 输入 n 对数(xi, yi) (i=0n)2. 令 s=03. 对于 i=0-nT=yi对 j=0- n 但 j!二i
3、 T二T*(X-xj)/(xi-xj)令 s=s+T4、输出结果流程图:第二部分：实验过程记录 (可加页)(包括实验原始数据记录，实验现象记录，实验过 程发现的问题等)#includemain()double X;int n;double lgrr(i nt n,double X);pri ntf( inp ut the aim numberX);printf( input the number of pairs of numbers n:); scanf(%lf%d, &X, &n);lgrr(n,X);return 0;double lgrr(i nt n,double X)double
4、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* (Xxj)/(xixj);s=s+t;pri ntf (s=%fn,s);return 0;第三部分结果与讨论(可加页)一、实验结果分析(包括数据处理、实验现象分析、影响因素讨论、综合分析和结论等)二、小结、建议及体会三、思考题对于现在的许多实际问题来说，我们并不知道f(x)的具体形式，所对 应的函数
5、值可能是由测量仪器或其他设备中直接读出来的，f(x)只是一 个数学概念意义下的函数。(比如：图像的方法处理，天气预报，机床 加工等方面)解答这类问题的方法就是插值方法。泰勒插值要求提供f(x)在点x0处的各阶导数值，这项要求很苛刻， 函数f(x)的表达式必须相当简单才行。如果仅仅给出一系列节点上的函 数值f(xi) = yi(i=0,1,2,n),则插值问题可表述如下：求作 n次多项 式Pn(x)，使满足条件Pn(x)= yi，i = 0，1，n。这就是所谓拉格朗日 (Lagrange)插值。通过本次实验，我不仅学会了如何用程序实现拉格朗日插值的算 法，而且更深刻的理解了拉格朗日插值的原理及方法。教师签字
