Python numpy polynomial.polyint用法及代码示例

本文简要介绍 python 语言中 numpy.polynomial.polynomial.polyint 的用法。

用法:

polynomial.polynomial.polyint(c, m=1, k=[], lbnd=0, scl=1, axis=0)

积分多项式。

返回多项式系数c融合的m次从lbnd沿着轴.在每次迭代中，结果序列是倍增经过scl和一个积分常数，k， 被添加。比例因子用于变量的线性变化。 (“Buyer beware”：请注意，根据一个人在做什么，一个人可能想要scl成为人们可能期望的倒数；有关详细信息，请参阅下面的注释部分。)参数c是一个系数数组，沿每个轴从低到高，例如，[1,2,3] 表示多项式1 + 2*x + 3*x**2而 [[1,2],[1,2]] 表示1 + 1*x + 2*y + 2*x*y如果轴 = 0 是x和轴= 1 是y.

参数：

c： array_like

多项式系数的一维数组，从低到高排序。

m： 整数，可选

积分顺序，必须是正数。 (默认值：1)

k： {[]，列表，标量}，可选

积分常数。第一个零积分的值是列表中的第一个值，第二个零积分的值是第二个值，依此类推。如果k == [](默认)，所有常量都设置为零。如果 m == 1 ，可以给出单个标量而不是列表。

lbnd： 标量，可选

积分的下界。 (默认值：0)

scl： 标量，可选

在每次积分之后，在添加积分常数之前，结果会乘以 scl。 (默认值：1)

axis： 整数，可选

进行积分的轴。 (默认值：0)。

返回：

S： ndarray

积分的系数数组。

抛出：

ValueError

如果 m < 1 、 len(k) > m 、 np.ndim(lbnd) != 0 或 np.ndim(scl) != 0 。

注意：

请注意，每次积分的结果是倍增经过scl.为什么要注意这一点？假设正在对变量进行线性变化\(u = ax + b\)在一个积分相对于x.然后\(dx = du/a\)，所以需要设置scl等于\(1/a\)- 也许不是人们首先想到的。

例子：

>>> from numpy.polynomial import polynomial as P
>>> c = (1,2,3)
>>> P.polyint(c) # should return array([0, 1, 1, 1])
array([0.,  1.,  1.,  1.])
>>> P.polyint(c,3) # should return array([0, 0, 0, 1/6, 1/12, 1/20])
 array([ 0.        ,  0.        ,  0.        ,  0.16666667,  0.08333333, # may vary
         0.05      ])
>>> P.polyint(c,k=3) # should return array([3, 1, 1, 1])
array([3.,  1.,  1.,  1.])
>>> P.polyint(c,lbnd=-2) # should return array([6, 1, 1, 1])
array([6.,  1.,  1.,  1.])
>>> P.polyint(c,scl=-2) # should return array([0, -2, -2, -2])
array([ 0., -2., -2., -2.])