Newton 插值
x1=;y1=;
x=0:0.5:3;
Newton(x1,y1,x);function y=Newton(x1,y1,x);
m=length(x);
n=length(x1);
for i=2:n
for j=n:-1:i
y1(j)=(y1(j)-y1(j-1))/(x1(j)-x1(j-i+1));
end
end
for j=1:m
p(j)=y1(1);
for k=2:n
l=1;
for i=1:k-1
l=l*(x(j)-x1(i));
end
p(j)=p(j)+y1(k)*l;
end
end
re=这个函数接受三个参数：
[*]x1：给定的数据点的 x 坐标。
[*]y1：给定的数据点的 y 坐标。
[*]x：要进行插值的目标 x 值。
该函数首先计算了差商，然后使用 Newton 插值公式计算在目标 x 值处的插值结果。最后，将结果以矩阵的形式返回，其中包含输入的 x 值和对应的插值结果。
页:
[1]