function Lagrange(X,Y,x0)
Len_y=length(Y);
Len_x=length(X);
if Len_y~=Len_x
errorlog('维数不等');
else
syms x;
p=0;
for i=1:Len_y
q=1;
for j=1:Len_x
if j~=i
q=q*(x-X(j))/(X(i)-X(j));
end
end
p=p+q*Y(i);
end
simplify(p);
if nargin==3
sprintf('给定点处的函数值得近似值为: %f',subs(p,'x',x0))
else
sprintf('所求的拉格朗日的插值函数为: ')
p=collect(p);
p=vpa(p,6),
end
end
end
clc;
X=[2 2.5 4];
Y=[0.5 0.4 0.25];
Lagrange(X,Y);
matlab 拉格朗日插值法
最新推荐文章于 2024-06-24 16:37:34 发布