Вам нужно несколько дополнительных шагов в дополнение к БПФ. Об этом уже много раз говорилось в предыдущих подобных вопросах на SO, и вы можете найти дополнительный материал, выполнив поиск по «dsp», «fft», «spectrum», «спектрограмма» и т. Д., Но по сути вам нужно сделать следующее :
- применить оконную функцию к входным данным (например, Hann (ing) )
- применить БПФ к оконным входным данным (для комплексного БПФ мнимые входные данные должны быть равны нулю)
- рассчитать квадратную величину первых
N / 2
выходных лотков FFT (re * re + im * im
)
- преобразование квадрата в шкалу дБ (
10 * log10(squared_magnitude)
)