可能的“最佳”解决方案确实是 DFT,丢弃前 3/4 的频率,并执行逆 DFT,域限制在底部的 1/4。在这种情况下,丢弃前 3/4 是低通滤波器。填充到 2 个样本的幂可能会给您带来速度优势。请注意您的 FFT 包如何存储样本。如果它是一个复杂的 FFT(更容易分析,并且通常具有更好的特性),频率将从 -22 到 22,或 0 到 44。在第一种情况下,您需要中间的 1/4。在后者中,最外层的 1/4。
您可以通过将样本值平均在一起来完成足够的工作。以四乘四的方式抓取样本并进行相等加权平均的天真方法是可行的,但并不太好。相反,您需要使用“内核”函数以非直观的方式将它们平均在一起。
在数学上,丢弃低频带之外的所有内容是乘以频率空间中的盒函数。(逆)傅立叶变换将逐点乘法转换为函数的(逆)傅立叶变换的卷积,反之亦然。因此,如果我们想在时域工作,我们需要使用盒函数的(逆)傅里叶变换进行卷积。结果证明这与“sinc”函数(sin at)/at 成正比,其中 a 是频率空间中框的宽度。因此,在每第 4 个位置(因为您要进行 4 倍的下采样),您可以将其附近的点相加,乘以 sin (a dt) / a dt,其中 dt 是到该位置的时间距离。有多近?好吧,这取决于您希望它听起来有多好。它'
最后,还有一种很糟糕(但很快)的方法,就是丢弃大部分样本,只保留第零个、第四个等等。
老实说,如果它适合记忆,我建议只走 DFT 路线。如果它不使用其他人推荐的软件过滤器包之一为您构建过滤器。