scipy.signal.

cheb2ord#

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

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

返回最低阶的数字或模拟切比雪夫 II 类滤波器的阶数，使其在通带上损失不超过 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）。

analogbool，可选

为真时返回模拟滤波器，否则返回数字滤波器。

fsfloat，可选

数字系统的采样频率。

1.2.0 版本中新增。

返回:
ordint

满足规格的 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
本页内容