python 线性插值处理_详解Python实现线性插值法

在算法分析过程中，我们经常会遇到数据需要处理插值的过程，为了方便理解，我们这里给出相关概念和源程序，希望能帮助到您!

案例一：

已知坐标 (x0, y0) 与 (x1, y1)，要求得区间 [x0, x1] 内某一点位置 x 在直线上的y值。两点间直线方程，我们有

那么，如何实现它呢？

依据数值分析，我们可以发现存在递归情况

import matplotlib.pyplot as plt

"""

详细描述;

1)计算n阶差商 f[x0, x1, x2 ... xn]

2) xi 所有插值节点的横坐标集合 o

3) fi 所有插值节点的纵坐标集合 / \

4) 返回xi的i阶差商(i为xi长度减1) o o

5) a. 必须确保xi与fi长度相等 / \ / \

b. 由于用到了递归，所以留意不要爆栈了. o o o o

c. 递归减递归(每层递归包含两个递归函数), 每层递归次数呈二次幂增长，总次数是一个满二叉树的所有节点数量(所以极易栈溢出)

"""

def get_order_diff_quot(xi = [], fi = []):

if len(xi) > 2 and len(fi) > 2:

return (get_order_diff_quot(xi[:len(xi) - 1], fi[:len(fi) - 1]) - get_order_diff_quot(xi[1:len(xi)], fi[1:len(fi)])) / float(xi[0] - xi[-1])

return (fi[0] - fi[1]) / float(xi[0] - xi[1])

"""

6) 获得Wi(x)函数;

Wi的含义举例 W1 = (x - x0); W2 = (x - x0)(x - x1); W3 = (x - x0)(x - x1)(x - x2)

7) i i阶(i次多项式)

8) xi 所有插值节点的横坐标集合

9) 返回Wi(x)函数

"""

def get_Wi(i = 0, xi