稀疏傅里叶变换（频率降采样）

设复信号 x(n) 长度为 N，X(k) 为其 DFT 结果。参数 B 整除 N，想要以等间隔 N/B对信号频域进行降采样，即求Y(K)=X(K*N/B) K∈[0,B-1].
可以断言，Y(k) 是以下构造信号的 B 点 DFT 结果：

参考文献的公式推导：

参考文献：Simple and Practical Algorithm for Sparse Fourier Transform