x=[1 2 3]
y=[1 4 9]
xi=5
yi=lagrange(x,y,xi)
function yi=lagrange(x,y,xi)
n=length(x);m=length(y);
if n~=m
error('The length of X must be equal!');
return;
end
p=zeros(1,n);
for k=1:n
t=ones(1,n);
for j=1:n
if j~=k
if abs(x(k)-x(j))<eps
error('the DATA is error! ');
return;
end
end
p(k)=prod(t);
end
end
yi=sum(y.*p)
end
MATLAB的拉格朗日算法
最新推荐文章于 2024-09-14 23:37:48 发布