【学习笔记】第七章 插值与拟合

7.1 插值

7.1.1 插值方法

1. Lagrange插值

例7.1 编写函数Lag_intp(x,y,x0) ，实现Lagrange插值，其中x和y是两个具有相同长度的Numpy数组。

def w(x, y, i, x0):
    p=1.0
    for j in range(len(x)):
        if j==i: continue
        p *= (x0-x[j]); p /= (x[i]-x[j])
    return p

def Lag_intp(x, y, x0):
    s=0
    for i in range(len(i)): s += w(x, y, i, x0)
    return s

2.分段线性插值

原理：用折现代替曲线y=f(x)，也称折现插值或者分段线性插值。分段线性插值多项式为：

3. 分段二次插值

原理：用分段抛物线代替曲线y=f(x)，也称分段抛物线插值，此时要求有2n+1个节点。其插值公式为：