Question

lagrange插值的MATLAB实现

Answer 1

function f=Language(x,y,x0)
% 拉格朗日
syms t;
if(length(x)==length(y))
n=length(x);
else
disp('x，和y的维数不相等');
return;
end
f=0.0;
for(i=1:n)
l=y(i);
for(j=1:i-1)
l=l*(t-x(j))/(x(i)-x(j));
end;
for(j=i+1:n)
l=l*(t-x(j))/(x(i)-x(j));
end;
f=f+l;
simplify(f);
if(i==n)
if(nargin==3)
f=subs(f,'t',x0);
else
f=collect(f);
f=vpa(f,6);
end
end
end
看看对不对,很早做的,自己也不太清楚了,呵呵.

Answer 2

用计算器算好像更简单呢，呵呵