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 新版功能.

退货
ord集成

符合规格的契比雪夫I型过滤的最低订单。

wnndarray或浮动

与一起使用的切比雪夫固有频率(“3dB频率”) cheby1 给出过滤结果。如果 fs 是指定的,则使用相同的单位,并且 fs 还必须传递给 cheby1 。

参见

cheby1

过滤设计的使用顺序和关键点

buttord

从通带和阻带规格中查找顺序和临界点

cheb2ord, ellipord
iirfilter

使用阶数和临界频率的通用过滤设计

iirdesign

通带和阻带规范的通用过滤设计

示例

设计一种数字低通过滤,使其通带在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()
../../_images/scipy-signal-cheb1ord-1.png