scipy.signal.
cheb2ord#
- scipy.signal.cheb2ord(wp, ws, gpass, gstop, analog=False, fs=None)[源代码]#
切比雪夫 II 型滤波器阶数选择。
返回最低阶数字或模拟切比雪夫 II 型滤波器的阶数,该滤波器在通带内的衰减不超过 gpass dB,并且在阻带内至少有 gstop dB 的衰减。
- 参数:
- wp, wsfloat
通带和阻带边缘频率。
对于数字滤波器,这些频率的单位与 fs 相同。默认情况下,fs 为 2 个半周期/采样,因此这些频率被归一化为 0 到 1,其中 1 是奈奎斯特频率。(因此,wp 和 ws 的单位为半周期/采样。)例如
低通:wp = 0.2, ws = 0.3
高通:wp = 0.3, ws = 0.2
带通:wp = [0.2, 0.5], ws = [0.1, 0.6]
带阻:wp = [0.1, 0.6], ws = [0.2, 0.5]
对于模拟滤波器,wp 和 ws 是角频率(例如,弧度/秒)。
- gpassfloat
通带中的最大损耗 (dB)。
- gstopfloat
阻带中的最小衰减 (dB)。
- analogbool, 可选
当为 True 时,返回一个模拟滤波器,否则返回一个数字滤波器。
- fsfloat, 可选
数字系统的采样频率。
在 1.2.0 版本中新增。
- 返回:
另请参阅
示例
设计一个数字带阻滤波器,该滤波器在 0.2*(fs/2) 到 0.5*(fs/2) 之间抑制 -60 dB,同时在 0.1*(fs/2) 以下或 0.6*(fs/2) 以上保持在 3 dB 以内。绘制其频率响应,以灰色显示通带和阻带约束。
>>> from scipy import signal >>> import matplotlib.pyplot as plt >>> import numpy as np
>>> N, Wn = signal.cheb2ord([0.1, 0.6], [0.2, 0.5], 3, 60) >>> b, a = signal.cheby2(N, 60, Wn, 'stop') >>> w, h = signal.freqz(b, a) >>> plt.semilogx(w / np.pi, 20 * np.log10(abs(h))) >>> plt.title('Chebyshev II bandstop filter fit to constraints') >>> plt.xlabel('Normalized frequency') >>> plt.ylabel('Amplitude [dB]') >>> plt.grid(which='both', axis='both') >>> plt.fill([.01, .1, .1, .01], [-3, -3, -99, -99], '0.9', lw=0) # stop >>> plt.fill([.2, .2, .5, .5], [ 9, -60, -60, 9], '0.9', lw=0) # pass >>> plt.fill([.6, .6, 2, 2], [-99, -3, -3, -99], '0.9', lw=0) # stop >>> plt.axis([0.06, 1, -80, 3]) >>> plt.show()