拉格朗日插值法_matlab

《拉格朗日插值法_matlab》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、MATLAB结课作业姓名：郭海阳班级：机械093学号：2009111006成绩：时间：2012/6/8一．任务。用matlab编写拉格朗日插值算法的程序并且以（x=-2.00,f(x)=17.00x=0.00,f(x)=1.00x=1.00,f(x)=2.00x=2.00,f（x)=17.00)为数据基础，在整个插值区间上采用拉格朗日插值算法计算f(x=0.6)，写出程序源代码，输出计算结果二．算法。x0=-2.00;x1=0.00;x2=1.00;x3=2.00;y0=17.00;y1=1.00;y2=2.0

2、0;y3=17.00;x=0.6y=(x-x1).*(x-x2).*(x-x3)/((x0-x1).*(x0-x2).*(x0-x3))*y0+(x-x0).*(x-x2).*(x-x3)/((x1-x0).*(x1-x2).*(x1-x3))*y1+(x-x0).*(x-x1).*(x-x3)/((x2-x0).*(x2-x1).*(x2-x3))*y2+(x-x0).*(x-x1).*(x-x2)/((x3-x0).*(x3-x1).*(x3-x2))*y3;disp('y=');disp(y);结果为：x

3、=0.6000y=0.2560三.程序。functions=Lagrange(x,y,x0)%lagrange插值，x，y为已知的插值点及其函数值%x0为要求的插值点的x值nx=length(x);ny=length(y);ifnx~=nywarning('矢量x与y的长度应该相等')returnendm=length(x0);%按照公式，对要求的插值点矢量x0的每个元素进行计算fori=1:mt=0.0;forj=1:nxu=1.0;fork=1:nxifk~=ju=u*(x0(i)-x(k))/(x(j)-

4、x(k));endendt=t+u*y(j);ends(m)=t;endreturn四．算例及分析。给定x,y,x0↓nx←length(x)ny←length(y)n←length(x0)t←0U=1↓nx=ny↓j=1…nx↓k=1:nx