拉格朗日插值 MATLAB

Question

这个程序为什么运行不了
function f=lagrange(x,y)
syms x
n=length(x);
for i=1:n
for j=2:n
a=x-x(j);
ji1=1;
ji1=ji1*a;
if i~=j
b=x(i)-x(j);
else
b=1;
end
ji2=1;
ji2=b*ji2;
end
f=sum(y(i)*ji1/ji2);
end

Answer 1

额，你这是lagrange插值？
这里有个lagrange插值的例子，你参照一下吧。
function y=lagrange_interp(xdata,ydata,x)
% Lagrange插值
% 输入参数：
% ---xdata：给定的节点横坐标
% ---ydata：给定的节点纵坐标
% ---x：需要进行插值的节点横坐标
% 输出参数：
% ---y：Lagrange插值函数在x处的函数值
n=length(xdata);m=length(ydata);
if n~=m
error('插值数据长度不等！');
end
ii=1:n;y=zeros(size(x));
for i=ii
ij=find(ii~=i);V=1;
for j=1:length(ij)
if abs(xdata(i)-xdata(ij(j)))<eps
error('输入的n+1个节点不是互异的。');
end
V=V.*(x-xdata(ij(j)));
end
y=y+V*ydata(i)/prod(xdata(i)-xdata(ij));
end