拉格朗日插值法_matlab.doc

文档介绍：function f = Language(x,y,x0)
%求已知数据点的拉格朗日插值多项式
%已知数据点的x坐标向量: x
%已知数据点的y坐标向量: y
%插值点的x坐标: x0
%求得的拉格朗日插值多项式或在x0处的插值: f
x=[ ];%input x data(可替换为自己的数据)
y=[0 ];%input y data(可替换为自己的数据)
x0=[ ];%input x0 data(可替换为自己的数据)
syms t l;
if(length(x) == length(y))
n = length(x);
else
disp('x和y的维数不相等!');
return; %检错
end
p=sym(0);
for (i=1:n)
l=sym(y(i));
for(k=1:i-1)
l=l*(t-x(k))/(x(i)-x(k));
end;
for(k=i+1:n)
l=l*(t-x(k))/(x(i)-x(k));
end;
p=p+l;
end
simplify(p);%简化多项式
f = subs (p,'t',x0); %计算插值点的函数值
f = vpa(f,6); %将插值多项式的值化成6位精度的小数
end