在這里,我試圖創建一個截止頻率為 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
提交
取消