scipy.signal.cheb2ord¶
- scipy.signal.cheb2ord(wp, ws, gpass, gstop, analog=False, fs=None)[源代码]¶
切比雪夫第二类过滤订单选择。
返回最低阶的数字或模拟切比雪夫II型过滤的顺序,其损失不超过 gpass DB在通带内,并且至少具有 gstop 阻带中的dB衰减。
- 参数
- wp、ws浮动
通带和阻带边缘频率。
对于数字滤波器,它们的单位与 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)。
- gpass浮动
通带内的最大损耗(DB)。
- gstop浮动
阻带内的最小衰减(DB)。
- analog布尔值,可选
如果为True,则返回模拟过滤,否则返回数字过滤。
- fs浮动,可选
数字系统的采样频率。
1.2.0 新版功能.
- 退货
参见
示例
设计一种抗0.2dB-60dB的数字带阻过滤 (fs/2) to 0.5 (fS/2),同时保持在低于0.1的3 dB内 (fs/2) or above 0.6 (FS/2)。绘制其频率响应,以灰色显示通带和阻带约束。
>>> from scipy import signal >>> import matplotlib.pyplot as plt
>>> 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()