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