Python 执行偶数扩展的离散傅里叶变换的最快方法是什么？

我有一个粒子轨迹的离散表示

q（t），t=0..t

，也就是说，我有

q\n，n=0..n

。我想要一个精确计算导数的谱表示，但信号不是周期性的。对于应用程序来说，均衡扩展似乎是最好的方法，但是如何高效地实现它（在Python/scipy中）？在我看来，一定有某种方法可以做到这一点，而不必首先构造带有

2N

项的偶数扩展数组。此外，这样做将导致在

dct

之后出现

2N

余弦，但其中一半将为零