当前位置: 首页>>代码示例 >>用法及示例精选 >>正文


Python SciPy signal.cheb2ord用法及代码示例


本文简要介绍 python 语言中 scipy.signal.cheb2ord 的用法。

用法:

scipy.signal.cheb2ord(wp, ws, gpass, gstop, analog=False, fs=None)#

Chebyshev II 型滤波器阶数选择。

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

参数

wp, ws 浮点数

通带和阻带边频率。

对于数字滤波器,它们的单位与 fs 相同。默认情况下,fs 为 2 half-cycles/sample,因此这些从 0 归一化为 1,其中 1 是奈奎斯特频率。 (因此,wp 和 ws 在 half-cycles /样本中。)例如:

  • Lowpass: wp = 0.2, ws = 0.3

  • Highpass: wp = 0.3, ws = 0.2

  • Bandpass: wp = [0.2, 0.5], ws = [0.1, 0.6]

  • Bandstop: wp = [0.1, 0.6], ws = [0.2, 0.5]

对于模拟滤波器,wp 和 ws 是角频率(例如,rad/s)。

gpass 浮点数

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

gstop 浮点数

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

analog 布尔型,可选

如果为 True,则返回模拟滤波器,否则返回数字滤波器。

fs 浮点数,可选

数字系统的采样频率。

返回

ord int

符合规格的切比雪夫 II 型过滤器的最低阶。

wn ndarray 或浮点数

切比雪夫固有频率(“3dB 频率”)用于scipy.signal.cheby2给出过滤结果。如果fs被指定,这是在相同的单位,并且fs也必须传递给scipy.signal.cheby2.

例子

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

相关用法


注:本文由纯净天空筛选整理自scipy.org大神的英文原创作品 scipy.signal.cheb2ord。非经特殊声明,原始代码版权归原作者所有,本译文未经允许或授权,请勿转载或复制。