ifft(fft(audio)) 只是噪音

Question

无论我只是嵌套它们（iff(fft(audio))）还是尝试逐个窗口（@98​​7654322@ 音频，执行fft，执行ifft，然后反转窗口，用 eps 替换零，然后将样本合并回来（在管道中到处尝试abs））我只得到噪音。

我知道 ifft 仅与 fft 反向，具有无限精度算术、无限多样本等（对吗？）我正在使用 64 位浮点和 44kHZ 采样率。但是，我希望至少能够听到原始音频。

我的错误是实际的还是理论上的？如果是错误，我可以提供代码。

Answer 1

根据 JoeKington 的评论，我已经下载了this file，并尝试了以下操作

>>> import scipy.io.wavfile
>>> rate, data = scipy.io.wavfile.read('wahoo.wav')
>>> data
array([134, 134, 134, ..., 124, 124, 124], dtype=uint8)
>>> data_bis = np.fft.ifft(np.fft.fft(data))
>>> data_bis
array([ 134. +6.68519934e-14j,  134. -4.57982480e-14j,
        134. -1.78967708e-14j, ...,  124. -2.09835513e-14j,
        124. -1.61750469e-14j,  124. -2.14867343e-14j])
>>> data_bis = data_bis.astype('uint8')
C:\Users\Jaime y Eva\Desktop\stack_exchange.py:1: ComplexWarning: Casting complex values to real discards the imaginary part
  # -*- coding: utf-8 -*-
>>> data_bis
array([134, 133, 133, ..., 123, 123, 123], dtype=uint8)
>>> scipy.io.wavfile.write('wahoo_bis.wav', rate, data_bis)

生成的文件与原始文件完全相同。

因此，将返回的复数值转换为实数只是问题的一半（您可能希望使用 np.abs 而不是 data.real，就像上面的代码隐含的那样），然后您还需要重铸浮动将数字指向适当位深度的uints。