scipy.signal.
cheb1ord#
- scipy.signal.cheb1ord(wp, ws, gpass, gstop, analog=False, fs=None)[source]#
切比雪夫 I 型滤波器阶次选择。
返回最低阶的数字或模拟切比雪夫 I 型滤波器阶次,该滤波器在通带中的衰减不超过 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 是角频率(例如,rad/s)。
- gpassfloat
通带中的最大衰减(dB)。
- gstopfloat
阻带中的最小衰减(dB)。
- analog布尔值,可选
如果为 True,则返回模拟滤波器,否则返回数字滤波器。
- fs浮点数,可选
数字系统的采样频率。
在 1.2.0 版中添加。
- 返回值:
另请参见
示例
设计一个数字低通滤波器,其通带在高达 0.2*(fs/2) 的范围内不超过 3 dB,同时抑制至少 -40 dB 的 0.3*(fs/2) 以上频率。绘制其频率响应,以灰色显示通带和阻带约束。
>>> from scipy import signal >>> import matplotlib.pyplot as plt >>> import numpy as np
>>> N, Wn = signal.cheb1ord(0.2, 0.3, 3, 40) >>> b, a = signal.cheby1(N, 3, Wn, 'low') >>> w, h = signal.freqz(b, a) >>> plt.semilogx(w / np.pi, 20 * np.log10(abs(h))) >>> plt.title('Chebyshev I lowpass filter fit to constraints') >>> plt.xlabel('Normalized frequency') >>> plt.ylabel('Amplitude [dB]') >>> plt.grid(which='both', axis='both') >>> plt.fill([.01, 0.2, 0.2, .01], [-3, -3, -99, -99], '0.9', lw=0) # stop >>> plt.fill([0.3, 0.3, 2, 2], [ 9, -40, -40, 9], '0.9', lw=0) # pass >>> plt.axis([0.08, 1, -60, 3]) >>> plt.show()