python 多项式求系数_在Python中用于计算“多项式系数”的numpy / scipy函数

是否有任何

python函数(可能来自numpy或scipy)计算扩展中x ** r的系数(1 xx ** 2 x ** 3 … x **(k-1))** n ,其中k> = 1,n> = 0且0 <= r <= n(k-1)？这有时被称为多项式系数(PC)(参见,例如,

here).

如果没有,你能想到一种有效的计算方法吗？ (我对天真/贪婪的方式不感兴趣).

你实际上正在进行[1,1,1,…,1,1,1]的n次卷积.

因此,您是否考虑在足够零填充的阵列上使用FFT,将其元素提升到幂n并使用逆FFT来恢复所有系数

(1+x+x**2+x**3+...+x**(k-1))**n

然后只是阅读你感兴趣的那些？

更新：

由于FFT是循环的,因此您需要一个不小于其中的项数的数组

(1+x+x**2+x**3+...+x**(k-1))**n

或者换句话说,(k-1)* n 1使结果不会在末端环绕(或者至少在它们这样做时它们只向受影响的元素添加零).通常它的长度也应该是2的幂,因为这是FFT算法所要求的(实现不需要它将用零填充输入直到它).

在类似C的伪代码中：

unsigned int m = 1;

while(m

complex c[m];

for(unsigned int i=0;i!=k;++i) c[i] = complex(1.0, 0.0);

for(unsigned int i=k;i!=m;++i) c[i] = complex(0.0, 0.0);

c = fft(c);

for(unsigned int i=0;i!=m;++i) c[i] = pow(c[i], double(n));

c = inv_fft(c);

在此结束时,复数组c的第r个元素的实部等于x ** r的系数和零的虚部.

现在,由于这都是以浮点形式完成的,因此您应该知道这些元素将累积舍入错误.您可以通过将它们四舍五入到最接近的整数来部分纠正这一点,但要注意,对于足够大的k和n,这些误差将超过0.5,因此这会产生因某些小的相对误差而关闭的结果.

在网上快速搜索显示,numpy分别具有FFT及其逆,numpy.fft.rfft和numpy.fft.irfft的实现,您可以在输入数据为真时使用.

python 多项式求系数_在Python中用于计算“多项式系数”的numpy / scipy函数相关推荐

python多项式运算_python – 用于计算“多项式系数”的numpy / scipy函数
是否有任何python函数(可能来自numpy或scipy)计算扩展中x ** r的系数(1 xx ** 2 x ** 3 - x **(k-1))** n ,其中k> = 1,n> = ...
python 多项式求系数_Python 从一元多项式中提取系数和次数，并进行简单的运算...
import re #正则 from collections import defaultdict #defaultdict #提取一元多项式(type: str)中的次数和系数并转化为字典 -> ...
python 多项式求系数_多项式系数的计算
def evalPoly(lst, x): total = 0 for power, coeff in enumerate(lst): # starts at 0 by default total + ...
python编程求导数_在python中如何计算导数？
这个问题的答案很简单.当然,在另一个答案中给出的subs选项可用于计算某个数的导数,但如果要绘制导数,则不起作用.有一种方法可以解决这个问题:lambdify,如下所述. 使用lambdify将所有s ...
python字典求平均值_查找字典中值的平均值
9 个答案: 答案 0 :(得分:5) 如果您使用numpy: import numpy as np np.array(list(dict.values())).mean() 答案 1 :(得分:4) ...
python多项式求导_链表一元多项式求导
[c++]代码库#include using namespace std; typedef struct PNode { int data1; int data2; struct PNode *nex ...
python 多项式求系数_python – sklearn：如何获得多项式特征的系数
顺便说一句,现在有更合适的功能: PolynomialFeatures.get_feature_names. from sklearn.preprocessing import PolynomialF ...
python惰性求值_让Python中类的属性具有惰性求值的能力
起步我们希望将一个只读的属性定义为 property 属性方法,只有在访问它时才进行计算,但是,又希望把计算出的值缓存起来,不要每次访问它时都重新计算. 解决方案定义一个惰性属性最有效的方法就是利 ...
python迭代法求极值_用Python实现最速下降法求极值的方法
对于一个多元函数 ,用最速下降法(又称梯度下降法)求其极小值的迭代格式为其中为负梯度方向,即最速下降方向,αkαk为搜索步长. 一般情况下,最优步长αkαk的确定要用到线性搜索技术,比如精确线性搜 ...

python 多项式求系数_在Python中用于计算“多项式系数”的numpy / scipy函数

python 多项式求系数_在Python中用于计算“多项式系数”的numpy / scipy函数相关推荐

最新文章

热门文章