scipy.signal.

cheb2ord#

scipy.signal.cheb2ord(wp, ws, gpass, gstop, analog=False, fs=None)[源代码][源代码]#

切比雪夫II型滤波器阶数选择。

返回在通带中损失不超过 gpass dB 且在阻带中至少有 gstop dB 衰减的最低阶数字或模拟切比雪夫 II 型滤波器的阶数。

参数:
wp, ws浮动

通带和阻带边缘频率。

对于数字滤波器,这些单位与 fs 相同。默认情况下,fs 是 2 半周期/样本,因此这些值从 0 到 1 归一化,其中 1 是奈奎斯特频率。(wpws 因此以半周期/样本为单位。)例如:

  • 低通滤波器: 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]

对于模拟滤波器,wpws 是角频率(例如,弧度/秒)。

gpass浮动

通带中的最大损耗(dB)。

gstop浮动

阻带中的最小衰减(dB)。

模拟bool, 可选

当为 True 时,返回一个模拟滤波器,否则返回一个数字滤波器。

fsfloat, 可选

数字系统的采样频率。

Added in version 1.2.0.

返回:
ord整数

满足规格的切比雪夫II型滤波器的最低阶数。

wnndarray 或 float

用于 cheby2 的切比雪夫自然频率(“3dB频率”)以获得滤波器结果。如果指定了 fs,则此频率单位相同,并且 fs 也必须传递给 cheby2

参见

cheby2

使用阶数和临界点的滤波器设计

buttord

从通带和阻带规格中找到阶数和临界点

cheb1ord, ellipord
iirfilter

使用阶数和临界频率进行一般滤波器设计

iirdesign

使用通带和阻带规格进行一般滤波器设计

示例

设计一个数字带阻滤波器,它在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()
../../_images/scipy-signal-cheb2ord-1.png