matlab求拉格朗日插值多项式函数
时间: 2023-12-12 16:34:46 浏览: 126
拉格朗日插值多项式
以下是使用Matlab实现拉格朗日插值算法的代码,可以通过给定的平面上的n个数据点,计算拉格朗日多项式Pn(x)的值,并将其作为实际函数f(x)的估计值。用户输入了数据点(xi,f(xi))、插值点之后,程序能够输出插值点对应的函数估值。
```matlab
function = Lagrange(x,f,x0)
%得出数据点的个数
n = length(x) ;
%得出插值点的个数
m = length(x0);
for i = 1:m
D = x0(i);
y = 0.0;
for k = 1:n
%置l(x)的初值
l = 1.0;
for j = 1:n
%迭代表示出l(xj)
%记住,第k个是取不到的!!!!
if j~=k
l = l*(D-x(j))/(x(k)-x(j));
end
end
%Pn(x)
y = y + l*f(k);
end
xx = num2str(D,'%.4f');
y =num2str(y,'%.4f');
disp('f(x)的近似值点坐标为:');
disp(['(',xx,',',y,')']);
end
end
```
阅读全文