我有一个音频wav文件,其中记录了一些人声。我已经使用scipy.read将wav文件数字化。我已经计算了FFT,然后从中得到了逆FFT。 当我绘制原始信号和逆FFT重叠时。他们完全匹配。 但是当我绘制原始信号和逆FFT的数据集时,它们根本不匹配。 我观察到,这种情况只发生在有人声的文件上,乐器录音也可以
以下是我尝试过的代码:
import scipy.io.wavfile as wavfile
import scipy
import scipy.fftpack as fftpack
from matplotlib import pyplot as plt
import pandas as pd
import csv
import numpy as np
fs_rate, signal = wavfile.read("Sa.wav")
print ("Frequency sampling", fs_rate," Hz")
l_audio = len(signal.shape)
print ("Channels", l_audio,"Audio data shape",signal.shape,"l_audio",l_audio)
if l_audio == 2:
signal = signal.sum(axis=1) / 2
N = signal.shape[0]
length = N / fs_rate
print ("Duration of audio wav file in secs", length,"Samples chosen",signal.shape[0])
time =np.linspace(0, length, signal.shape[0])
sampling_interval=time[1]-time[0]
inverse_fft_amplitude = scipy.ifft(scipy.fft(signal))
plt.xlabel('Time[s]')
plt.ylabel('Amplitude')
plt.plot(time, signal, 'blue', time, inverse_fft_amplitude , 'red')
plt.legend(('Original Signal', 'Reconstructed Signal using Inverse FFT'))
# They match exactly same
#Getting the data points for original signal
df = pd.DataFrame({"Time" : time, "Amplitude" : signal})
df.to_csv("OriginalData.csv", header=None, index=False)
#Getting the data points for inverse FFT signal
df = pd.DataFrame({"Time" : time, "Amplitude" : abs(inverse_fft_amplitude)})
df.to_csv("ReconstructedDataSet.csv", header=None,index=False)
目前没有回答
相关问题 更多 >
编程相关推荐