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

Question

我有一个离散的粒子轨迹表示q(t), t=0..T，即我有q_n, n=0..N。我希望准确地计算衍生物的光谱表示，但信号不是周期性的。鉴于应用程序，甚至扩展似乎是最好的方法，但我如何有效地实现这一点（在Python / scipy中）？在我看来，必须存在一些方法来执行此操作，而首先使用2N条目构建偶数扩展数组。此外，这样做会在2N之后产生dct余弦，但其中一半将为零。