《实验一 拉格朗日插值Matlab实验报告》由会员分享，可在线阅读，更多相关《实验一 拉格朗日插值Matlab实验报告（3页珍藏版）》请在金锄头文库上搜索。
1、北京理工大学珠海学院实验报告ZHUHAI CAMPAUS OF BEIJING INSTITUTE OF TECHNOLOGY班级2012电气2班 学号xxxxxxxxx姓名 陈冲 指导教师 张凯 成绩 实验题目(实验一)拉格朗日插值 实验地点及时间JD501 2013/12/26(6-7节)一、实验目的1. 掌握用程序语言来编辑函数。2. 学会用MATLAB编写Lagrange.m函数。二、实验环境Matlab软件三、实验内容1、以书中第55页题目13为例编辑程序来实现计算结果。2、使用MATLAB进行编写：第一步：编写Lagrange.m函数，代码如下function s=Lagrange
2、(xJy1x|O)%1招既切驴m函数为求已知数据点的拉格朗日插值多项式 为数据点的x坐标向量 为数据点的y坐标向量为插值的x坐标为求得的拉格朗日插值多项式或在湖处的插值symsp;nlength(x);%读取x向量维数s=;for(k=l:n)la=y(k);end;fbiG=k+l ;门)5a=la*(p-xQ)/(x(k)-x )；end;s=&+Ia;simplify(s);end.葺对输入参数个数作判断如果只有两个参数.则直接给出插值多项式%如果3个参数j则给出插值点的插值结臬 第三个参数可以为向量 ；f(nargin=2)s=subs(矿 p?K)；s=w lleGt(s);%多项式
3、展开尸vp&(M)；%把系数取到6位精度表达elsem=length(xO);%读取 xO 度%分别对湘的每个分量插值for i=l :mtemp(i)=subs(s,pxO(i);sd%得到的是系列插僵点的插值结果stemp;end第二步：利用这个函数来编辑命令：(可见实验结果中的截图)x=0.32,0.34,0.36;y=sin(0.32),sin(0.34),sin(0.36);x0=0.3367;yt=Lagrange(x,y,x0)得出抛物线插值为：0.3304以及x=0.32,0.34;y=sin(0.32),sin(0.34);x0=0.3367;yt=Lagrange(x,y,
4、x0)得出线性插值为：0.3304四、实验题目尤i0.320.340.36sin工0.3145670.3334870.352274试用线性插值和抛物线插值分别计算sin0.3367的近似值并估计误差。五、实验结果 x= 0. 32, 0. 54.0. 36: y= siritO. 32)j sim(0. S4)j sin(0. 3S); z0=0. 3367; yt=Lagrangey, k0)yt =0.3304 disp C电气二班陈冲电气二班陈冲0. 32, a. 34:y= s in (0. 32) s sin (. 34);x0=0.336?；Yt=LagrarB(Mj yjsOl dispC电气二班陈冲yt -0. 3304电气二班陈冲O六、总结通过这次实验我学会用MATLAB软件编辑口令进行计算，实验结果是正确的，我相 信在以后的实验中，我可以做好每一步，练习好每一次的上机。实验难度不是很大，主要要注意标点符号的正确性。
