用MATLAB进行数据插值.ppt

在MATLAB中编程实现向后牛顿插值法函数为：Newtonback。功能：求已知数据点的向前牛顿插值法多项式；调用格式：f=Newtonback(x,y)或f=Newtonback(x,y,x0)。其中，x为已知数据点的x坐标向量；y为已知数据点的y坐标向量；x0为插值点的x坐标；f为求得的向前牛顿插值法多项式或x0处的插值。第96页,共107页，2024年2月25日，星期天functionf=Newtonback(x,y,x0)%求已知数据点的向后差分牛顿插值多项式%已知数据点的x坐标向量:x%已知数据点的y坐标向量:y%为插值点的x坐标:x0%求得的向前差分牛顿插值多项式或x0处的插值:fsymst;if(length(x)==length(y))n=length(x);c(1:n)=0.0;elsedisp('x和y的维数不相等！');return;end第97页,共107页，2024年2月25日，星期天f=y(n);y1=0;xx=linspace(x(1),x(n),(x(2)-x(1)));if(xx~=x)disp('节点之间不是等距的！');return;endfor(i=1:n-1)for(j=i+1:n)y1(j)=y(j)-y(j-1);endc(i)=y1(n);l=t;for(k=1:i-1)l=l*(t+k);end;f=f+c(i)*l/factorial(i);simplify(f);y=y1;if(i==n-1)if(nargin==3)f=subs(f,'t',(x0-x(n))/(x(2)-x(1)));elsef=collect(f);f=vpa(f,6);endendend第98页,共107页，2024年2月25日，星期天例5根据下表的数据点求出其差分形式的牛顿插值多项式，并计算当x=1.55时y的值。解：>>x=[11.21.41.61.8];>>y=[0.84150.9320 0.98540.99960.9738];>>f=Newtonforward(x,y)f=.841500+.108025*t-.169042e-1*t^2-.675000e-3*t^3+.541667e-4*t^4>>f=Newtonforward(x,y,1.55)f=0.9998f=Newtonback(x,y)f=.973800-.457417e-1*t-.198042e-1*t^2+.191667e-3*t^3+.541667e-4*t^4>>f=Newtonback(x,y,1.55)f=0.9998x11.21.41.61.8y0.84150.93200.98540.99960.9738第99页,共107页，2024年2月25日，星期天Hermite插值Hermite插值满足在节点上等于给定函数值，而且在节点上的导数值也等于给定的导数值。对于高阶导数的情况，Hermite插值多项式比较复杂，在实际中，常常遇到的是函数值与一阶导数给定的情况。在此情况下，n个节点x1，x2，…，xn的Hermite插值多项式的表达式如下：其中，，，第100页,共107页，2024年2月25日，星期天在MATLAB中编程实现Hermite插值法函数为：Hermite。功能：求已知数据点的Hermite插值法多项式；调用格式：f=Hermite(x,y,y_1)或f=Hermite(x,y,y_1,x0)。其中，x为已知数据点的x坐标向量；y为已知数据点的y坐标向量；