一、一元函数插值
已知函数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'.
插值效果: