有一条长长的石头,外形像一条龙,石头上长满斑斑点点的东西,像龙身上的鳞。人生就是一场旅行,不在乎目的地,在乎的应该是沿途的风景以及看风景的心情。春色满园关不住,一枝红杏出墙来。
本文实例讲述了python计算对角线有理函数插值的方法。分享给大家供大家参考。具体实现方法如下:
''' p = rational(xData,yData,x)
Evaluates the diagonal rational function interpolant p(x)
that passes through he data points
'''
from numpy import zeros
def rational(xData,yData,x):
m = len(xData)
r = yData.copy()
rOld = zeros(m)
for k in range(m-1):
for i in range(m-k-1):
if abs(x - xData[i+k+1]) < 1.0e-9:
return yData[i+k+1]
else:
c1 = r[i+1] - r[i]
c2 = r[i+1] - rOld[i+1]
c3 = (x - xData[i])/(x - xData[i+k+1])
r[i] = r[i+1] + c1/(c3*(1.0 - c1/c2) - 1.0)
rOld[i+1] = r[i+1]
return r[0]
希望本文所述对大家的Python程序设计有所帮助。
到此这篇关于python计算对角线有理函数插值的方法就介绍到这了。学的时候就要想着现在和将来怎么去用,否则是浪费时间。更多相关python计算对角线有理函数插值的方法内容请查看相关栏目,小编编辑不易,再次感谢大家的支持!