为什么逆fft与原始信号不匹配

2024-09-10 08:09:41 发布

您现在位置:Python中文网/ 问答频道 /正文

我有一个音频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)

原始信号与逆FFT的重叠 enter image description here

重叠原始信号和逆FFT数据集 enter image description here


Tags: csvimportfftdfsignaltime信号as