在这里,我试图创建一个截止频率为 0.1 Hz 的高通巴特沃斯数字滤波器。我已经实现了以下代码,但我不确定它是否正确#%% creating the filter # filter parameters order=6btype='highpass'cutoff_frequency=0.1*2*np.pianalog=Falseb, a= signal.butter(order,cutoff_frequency,btype, analog)w, h = signal.freqs(b, a)plt.figure() plt.plot(w, 20 * np.log10(abs(h)))plt.xscale('log')plt.title('Butterworth filter frequency response')plt.xlabel('Frequency [radians / second]')plt.ylabel('Amplitude [dB]')plt.margins(0, 0.1)plt.grid(which='both', axis='both')plt.axvline(0.1*2*np.pi, color='green') # cutoff frequencyplt.show()我的困惑是关于截止频率,我将其乘以 2*pi,因为据我了解,scipy.signal.butter 的 cutoff_frequency 对应于以 rad/s 为单位的角频率。
添加回答
举报
0/150
提交
取消