问题描述
我想在音频文件中实现高通滤波器,以滤除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,order=5):
b,a = butter_highpass(cutoff,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,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开始
解决方法
暂无找到可以解决该程序问题的有效方法,小编努力寻找整理中!
如果你已经找到好的解决方法,欢迎将解决方案带上本链接一起发送给小编。
小编邮箱:dio#foxmail.com (将#修改为@)