无法使用 FFT 重新创建相同的声音答案

【问题标题】：Not able to recreate same sound using FFT无法使用 FFT 重新创建相同的声音
【发布时间】：2018-12-01 03:35:42
【问题描述】：

我正在尝试使用傅立叶变换 (FFT) 返回的前 10 个频率重新创建音符。产生的声音与原始声音不匹配。不确定我是否没有正确找到频率或没有正确地从中产生声音。此代码的目标是匹配原始声音。这是我的代码：

import numpy as np
from scipy.io import wavfile
from scipy.fftpack import fft
import matplotlib.pyplot as plt

i_framerate = 44100
fs, data = wavfile.read('./Flute.nonvib.ff.A4.stereo.wav') # load the data

def findFrequencies(arr_data, i_framerate = 44100, i_top_n =5):
    a = arr_data.T[0] # this is a two channel soundtrack, I get the first track
#        b=[(ele/2**8.)*2-1 for ele in a] # this is 8-bit track, b is now normalized on [-1,1)
    y = fft(a) # calculate fourier transform (complex numbers list)

    xf = np.linspace(0,int(i_framerate/2.0),int((i_framerate/2.0))+1) /2 # Need to find out this last /2 part
    yf = np.abs(y[:int((i_framerate//2.0))+1])

    plt.plot(xf,yf)

    yf_top_n = np.argsort(yf)[-i_top_n:][::-1]
    amp_top_n =  yf[yf_top_n] / np.max(yf[yf_top_n])
    freq_top_n = xf[yf_top_n]

    return freq_top_n, amp_top_n

def createSoundData(a_freq, a_amp, i_framerate=44100, i_time = 1, f_amp = 1000.0):
    n_samples = i_time * i_framerate

    x = np.linspace(0,i_time, n_samples)
    y = np.zeros(n_samples)
    for i in range(len(a_freq)):
        y += np.sin(2 * np.pi * a_freq[i] * x)* f_amp * a_amp[i]
    data2 = np.c_[y,y] # 2 Channel sound
    return data2

top_freq , top_freq_amp = findFrequencies(data, i_framerate = 44100 , i_top_n = 200)

print('Frequencies: ',top_freq)
print('Amplitudes : ',top_freq_amp) 

soundData = createSoundData(top_freq, top_freq_amp,i_time = 2, f_amp = 50 / len(top_freq))
wavfile.write('createsound_A4_v6.wav',i_framerate,soundData)

【问题讨论】：

您是否尝试过绘制波形和频谱图？
“不匹配”非常模糊 - 而且非常主观。
看起来你正在抛弃灵态。不要那样做。使用相位，卢克。
谢谢@CrisLuengo。我将更改代码以添加阶段并在评论中更新结果。

标签： python scipy fft wav wave

【解决方案1】：

音符中的前 10 个频谱频率与前 10 个 FFT 结果 bin 幅度的中心频率不同。实际频率峰值可以在 FFT 区间之间。

不仅频率峰值信息可以在 FFT bin 之间，重现任何音符瞬态（起音、衰减等）所需的相位信息也可以在 bin 之间。 FFT bin 之间的频谱信息由复数 FFT 结果的跨度（最大宽度）承载。

【讨论】：