(转载作者请注明出处)
废话不多说,直接上货....
在算法分析过程中,我们经常会遇到数据需要处理插值的过程,为了方便理解,我们这里给出相关概念和源程序,希望能帮助到您!
案例一:
已知坐标 (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)])) / floa