function yy=lagrange(x1,y1,xx) %本程序为Lagrange1插值,其中x1,y1 %为插值节点和节点上的函数值,输出为插值点xx的函数值, %xx可以是向量。 syms x n=length(x1); for i=1:n t=x1;t(i)=[];L(i)=prod((x-t)./(x1(i)-t));% L向量用来存放插值基函数 end u=sum(L.*y1); p=simplify(u) % p是简化后的Lagrange插值函数(字符串) yy=subs(p,x,xx); %p是以x为自变量的函数,并求xx处的函数值 end
Matlab实现拉格朗日插值函数
猜你喜欢
转载自blog.csdn.net/zhanshen112/article/details/79778545
今日推荐
周排行