MATLAB程序与数值实验
一、 插值法(课本第二章)
1、拉格朗日插值方法
function s=lag(x,y,t)
%拉格朗日插值方法，可以同时对多点插值，t可以为向量
syms p;
n=length(x);%读取x向量的维数
s=0;
for(k=1:n)
la=y(k);
%构造基函数
for(j=1:k-1)
la=la*(p-x(j))/(x(k)-x(j));
end;
for(j=k+1:n)
la=la*(p-x(j))/(x(k)-x(j));
end;
s=s+la;
simplify(s);
end
if(nargin==2)
s=subs(s,'p','x');
s=collect(s);%展开多项式
s=vpa(s,6);
else
m=length(t);
for i=1:m
temp(i)=subs(s,'p',t(i));
end
s=temp;
End
