我想在我的音频文件中实现高通滤波器,以滤除低于6kHz的任何信号。代码如下:
def butter_highpass(cutoff, fs, order=5):
nyq = 0.5 * fs
normal_cutoff = cutoff / nyq
b, a = signal.butter(order, normal_cutoff, btype='high', analog=False)
return b, a
def butter_highpass_filter(data, cutoff, fs, order=5):
b, a = butter_highpass(cutoff, fs, order=order)
y = signal.filtfilt(b, a, data)
return y
(fs,rate)=wav.read('2020_09_02-21_21_16_51.wav')
y=butter_highpass_filter(rate,6000, fs, order=5)
plt.figure(3)
frequencies, times, spectrogram =
signal.spectrogram(rate,fs)
plt.pcolormesh(times, frequencies, spectrogram)
plt.ylabel('Frequency [Hz]')
plt.xlabel('Time [sec]')
plt.show()
然而,结果输出是不正确的,即信号分量应在6kHz以下移除,相反,它们仍然处于活动状态,并更改为从kHz开始
目前没有回答
相关问题 更多 >
编程相关推荐