bet188金宝搏

进口matplotlib。Pyplot作为PLT
进口scipy。信号作为信号
将numpy导入为np
从数学导入圆周率
plt.close(所有的)
根据巴特沃斯规格设计IIR过滤器
Fs = 1000 #采样频率
n = 7 #过滤器的阶数
fc = 100 #滤波器截止频率(LPF/HPF)
Fc1 = np。array([200,400]) # BPF/BSF的启动和停止频率
w_c = 2 * fc/Fs
w_c1 = 2 * fc1/Fs
#不同类型的IIR过滤器配置
[b1, a1] = sig.butter(n, w_c, btype='lowpass')
[b2, a2] = sig.butter(n, w_c, btype='highpass')
[b3, a3] = sig.butter(n, w_c1, btype='带通')
[b4, a4] = sig.butter(n, w_c1, btype='bandstop')
IIR滤波器的频率响应(LPF型)
[w1, h1] = sig.freqz(b1, a1, worN=2000)
w_lpf = Fs*w1/(2*pi)
H_db_lpf = 20 * np.log10(abs(h1))
# IIR滤波器的频率响应(HPF型)
[w2, h2] = sig.freqz(b2, a2, worN=2000)
w_hpf = Fs*w2/(2*pi)
H_db_hpf = 20 * np.log10(abs(h2))
IIR滤波器的频率响应(BPF型)
[w3, h3] = sig.freqz(b3, a3, worN=2000)
w_bpf = Fs*w3/(2*pi)
H_db_bpf = 20 * np.log10(abs(h3))
# IIR滤波器的频率响应(BSF型)
[w4, h4] = sig.freqz(b4, a4, worN=2000)
w_bsf = Fs*w4/(2*pi)
H_db_bsf = 20 * np.log10(abs(h4))
#地块
图，轴= plt。次要情节(2,2)
轴(0,0)。情节(w_lpf h_db_lpf)
轴(0,0)。set_title("IIR过滤器响应(LPF类型)")
轴(0,0)。set_xlabel(“频率(赫兹)”)
轴(0,0)。set_ylabel(“(dB)级”)
轴[0,1]。情节(w_hpf h_db_hpf)
轴[0,1]。set_title("IIR过滤器响应(HPF类型)")
轴[0,1]。set_xlabel(“频率(赫兹)”)
轴[0,1]。set_ylabel(“(dB)级”)
轴(1,0)。情节(w_bpf h_db_bpf)
轴(1,0)。set_title("IIR过滤器响应(BPF类型)")
轴(1,0)。set_xlabel(“频率(赫兹)”)
轴(1,0)。set_ylabel(“(dB)级”)
轴[1]。情节(w_bsf h_db_bsf)
轴[1]。set_title("IIR过滤器响应(BSF类型)")
轴[1]。set_xlabel(“频率(赫兹)”)
轴[1]。set_ylabel(“(dB)级”)
plt.tight_layout ()
plt.show ()