matlab 插值出错,MATLAB插值问题

一、一元函数插值

已知函数y=f(x)在区间[a,b]上的n+1个不同点

的函数值为

，若存在一个简单函数F(x), 使

，称F(x)为f(x)在区间[a,b]上的插值函数，称(xi, yi)为插值节点。若F(x)为多项式，称为多项式插值(或代数插值) ；常用的代数插值方法有：拉格朗日插值，牛顿插值。

n次代数插值：已知f(x)在n+1个点x0,x1,…,xn处的函数值为 y0,y1,…,yn ， 求一个n次多项式函数Pn(x)，使其满足： Pn(xi)=yi, (i=0,1,…,n). 若Pn(x)按下述方式构造，称为拉格朗日插值

其中Li(x) 为n次多项式：

称为拉格朗日插值基函数．

特别地:

(1)已知两个节点时，得线性插值多项式:

(2)已知三个节点时，得抛物插值多项式:

(3)已知n+1个节点时，可得n次拉格朗日插值多项式。

关于代数插值：

可以看出，当节点较多时，多项式的次数增高，插值函数出现振荡，精度变低。因此，为了保证精度，在节点较多时，一般采用分段插值，但这样在分段点光滑性较差。Matlab采用的多项式插值都是分段插值法。从图形还可以看出，对解析函数，插值精度高；对有奇点的函数，插值精度低。多项式插值对靠近插值区间中点的部分插值精度高，远离中点部分精度低。

规则网点的插值

y=interp1(x0,y0,x,'method')，y是对应x插值后获得的因变量，x0和y0为初始数据的自变量和因变量

Method 的选项有 ‘nearest’, ‘next’, ‘previous’, ‘linear’,‘spline’,‘pchip’, 和 'cubic'. 缺省的机器设置为'linear'.

插值效果：