hermite插值法和拉格朗日插值法的matlab程序实现

2024-11-04 08:51:02
推荐回答(1个)
回答1:

Matlab函数M文件Lagrange程序
function yy=lagrange(x,y,xi)
m=length(x);n=length(y);
if m~=n,error('向量x与y的长度必须一致');end
s=0;
for i=1:n
z=ones(1,length(xi));
for j=1:n
if j~=i
z=z .*(xi-x(j))/(x(i)-x(j));
end
end
s=s+z*y(i);
end
yy=s;
上面是拉格朗日插值法,其中xi为要计算的数值比如
x=[0 3 5 9 31];
y=[2 7 10 12 15];
xi=[1 4 7] ;
yi=Lagrange(x,y,xi)
这是要求在1,4,7处的值。
hermite程序我没有。