Question

MATLAB 问题求救~

function [yy,c]=lagrange(x,y,z)
nx=length(x);ny=length(y);
n=min(nx,ny);
yy=0;
for i=1:n
p=1.0;
for j=1:n
if i～=j
if abs(x(i)-x(j))<eps
error('error');
return;
end
ll=[1.0-x(j)]/(x(i)-x(j));
p=conv(p,ll);
end
end
yy=pplus(yy,p*y(i));
end
if nargout=2
c=polyval(yy,2);
end
x=[1 2 3 4 5]
y=1.30 1.80 2.20 2.9 3.5]
p=lagrange(x,y)

已运行它，就会出现??? Strings passed to EVAL cannot contain function declarations

谁帮我修改一下，拜托了。。。。。

Answer 1

07年问的！我现在才挖出来，楼主，不好意思回答晚了

正确的应该如下

function [yy,c]=lagrange(x,y,z)
nx=length(x);ny=length(y);
if nx~=ny
warning('The lengths of X and Y should be equal');
end
n=min(nx,ny);
if n<2
error('the number of the DATA should be greater than 1');
return;
end
yy=0;
for i=1:n
p=1.0;
for j=1:n
if i~=j
if abs(x(i)-x(j))<eps
error('the DATA is error!');
return;
end
ll=[1 0-x(j)]/(x(i)-x(j)); %这边是1 0 不是1.0
p=conv(p,ll);
end
end
yy=plus(yy,p*y(i)); %现在pplus函数matlab里面没有了，变为plus。我以前看到的一些比较老的书里写的也是pplus，实际上就是相加，可以写成yy=yy+p*y(i)
end
if nargout==2
c=polyval(yy,z);
end

楼主，我验算过，应该没问题了！