Code666 (代码贴、代码片段)

[Python] python计算对角线有理函数插值 →→→→→进入此内容的聊天室

来自 , 2020-01-02, 写在 Python, 查看 114 次.
URL http://www.code666.cn/view/c862b029
  1. ''' p = rational(xData,yData,x)
  2.    Evaluates the diagonal rational function interpolant p(x)
  3.    that passes through he data points
  4. '''    
  5. from numpy import zeros
  6.  
  7. def rational(xData,yData,x):
  8.     m = len(xData)
  9.     r = yData.copy()
  10.     rOld = zeros(m)
  11.     for k in range(m-1):
  12.         for i in range(m-k-1):
  13.             if abs(x - xData[i+k+1]) < 1.0e-9:
  14.                 return yData[i+k+1]
  15.             else:
  16.                 c1 = r[i+1] - r[i]
  17.                 c2 = r[i+1] - rOld[i+1]
  18.                 c3 = (x - xData[i])/(x - xData[i+k+1])
  19.                 r[i] = r[i+1] + c1/(c3*(1.0 - c1/c2) - 1.0)
  20.                 rOld[i+1] = r[i+1]
  21.     return r[0]
  22. #//python/7402

回复 "python计算对角线有理函数插值"

这儿你可以回复上面这条便签

captcha

贴代码片段 - Code666