均差定义及性质(python实现)

简述

k阶均差

f[x0,x1,...,xk]=f[x0,x1,...,xk−2,xk]−f[x0,x1,...,xk−1]xk−xk−1f[x0,x1,...,xk]=f[x0,x1,...,xk−2,xk]−f[x0,x1,...,xk−1]xk−xk−1

一阶均差

f[x0,x1]=f(x1)−f(x0)x1−x0f[x0,x1]=f(x1)−f(x0)x1−x0

上述式子可以看到递推性质。
最终都会收敛到一阶均差。

性质

• f[x0,x1,...,xk]=∑j=0kf(xj)∏t!=j(xj−xt)f[x0,x1,...,xk]=∑j=0kf(xj)∏t!=j(xj−xt) 这个性质，带来的就是简化均差的简化运算（其实也很好记的，这里，分母不就是之前w′n(xj)wn′(xj)）
• f[x0,x1,...,xk]=f[x1,x2,...,xk]−f[x0,x1,...,xk−1]xk−x0f[x0,x1,...,xk]=f[x1,x2,...,xk]−f[x0,x1,...,xk−1]xk−x0
• f[x0,x1,...,xk]=f(n)(ξ)n!f[x0,x1,...,xk]=f(n)(ξ)n!

python实现

下面是用python实现的均差，ff()就是f[,,]
其中f(x)=1x2+1f(x)=1x2+1

def f(X):return 1 / (x ** 2 + 1)def ff(X=list()):if len(X) < 2:raise ValueError('X\'s length must be bigger than 2')ans = 0for i in range(len(X)):temp = 1.0for j in range(len(X)):if j == i:continuetemp *= (X[i] - X[j])ans += (f(X[i]) / temp)return ans

总结

以上是生活随笔为你收集整理的均差定义及性质(python实现)的全部内容，希望文章能够帮你解决所遇到的问题。

如果觉得生活随笔网站内容还不错，欢迎将生活随笔推荐给好友。