在MATLAB编程实现拉格朗日插值多项式
请懂得朋友帮忙一下,用MATLAB语言写出拉格朗日插值多项式。
function Ys=langr(X,Y,Xs)
%拉格朗日插值
n=length(X);
nl=length(Xs);
for i=1:nl
x=Xs(i);
s=0;
for j=1:n
p=1;
for k=1:n
if k~=j
p=p*(x-X(k))/(X(j)-X(k));
end
end
s=s+p*Y(j);
end
Ys(i)=s;
end
这个应该行吧,对给定的一组数据,求插值点的值