matlab中,已知原函数和插值点,怎么求三次拉格朗日插值多项式
发布网友
发布时间：2022-09-06 19:57
共1个回答
热心网友
时间：2023-11-10 10:17
%本程序为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);
clf
plot(x1,y1,'ro',xx,yy,'*')
