牛顿插值多项式例题 matlab中,已知原函数和插值点,怎么求三次拉格朗日插值多项式?

matlab中,已知原函数和插值点,怎么求三次拉格朗日插值多项式?函数YY=Lagrange(x1,Y1,XX)%,本程序为Lagrange 1插值,其中x1,Y1%为插值节点及节点上的函数值,输出为

matlab中,已知原函数和插值点,怎么求三次拉格朗日插值多项式?

函数YY=Lagrange(x1,Y1,XX)%,本程序为Lagrange 1插值,其中x1,Y1%为插值节点及节点上的函数值,输出为插值点XX的函数值,%XX可为矢量。Syms xn=length(x1)for I=1:NT=x1t(I)=[]l(I)=prod((x-t)/(x1(I)-t))%l向量用于存储插值基函数endu=sum(l.*Y1)P=simplify(U)%P是简化的拉格朗日插值函数(字符串)YY=subs(P,x,XX)clfplot(x1,Y1,“ro”,XX,YY,“*”)