Python SciPy interpolate.pade用法及代碼示例

本文簡要介紹 python 語言中 scipy.interpolate.pade 的用法。

用法:

scipy.interpolate.pade(an, m, n=None)#

將 Pade 逼近作為兩個多項式的比率返回多項式。

參數 ：：

an： (N,) 數組

泰勒級數係數。

m： int

返回的近似多項式 q 的階數。

n： 整數，可選

返回的近似多項式的階數p。默認情況下，順序是len(an)-1-m.

返回 ：：

p, q： 多項式類

多項式的 Pade 近似由下式定義一個是p(x)/q(x).

例子：

>>> import numpy as np
>>> from scipy.interpolate import pade
>>> e_exp = [1.0, 1.0, 1.0/2.0, 1.0/6.0, 1.0/24.0, 1.0/120.0]
>>> p, q = pade(e_exp, 2)

>>> e_exp.reverse()
>>> e_poly = np.poly1d(e_exp)

比較 e_poly(x) 和 Pade 近似值 p(x)/q(x)

>>> e_poly(1)
2.7166666666666668

>>> p(1)/q(1)
2.7179487179487181