scipy.signal.cheb1ord¶
- scipy.signal.cheb1ord(wp, ws, gpass, gstop, analog=False, fs=None)[源代码]¶
切比雪夫I型过滤订单选择。
返回最低阶的数字或模拟切比雪夫I类过滤的顺序,该顺序的损失不超过 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 新版功能.
- 退货
参见
示例
设计一种数字低通过滤,使其通带在3dB以内,最高可达0.2dB (fs/2), while rejecting at least -40 dB above 0.3 (FS/2)。绘制其频率响应,以灰色显示通带和阻带约束。
>>> from scipy import signal >>> import matplotlib.pyplot as plt
>>> 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()