scipy.signal.

ellipord#

scipy.signal.ellipord(wp, ws, gpass, gstop, analog=False, fs=None)[source]#

椭圆函数（科尔）滤波器阶数选择。

返回最低阶数数字或模拟椭圆滤波器的阶数，该阶数在通带内损失不超过 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）。

模拟量布尔值，可选

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

fs浮点数，可选

数字系统的采样频率。

在版本 1.2.0 中添加。

返回:
ord整数

满足规范的椭圆（Cauer）滤波器的最低阶。

wnndarray 或浮点数

与 ellip 配合使用以实现滤波结果的切比雪夫自然频率（“3dB 频率”）。如果指定了 fs，则采用相同单位，并且也必须将 fs 传递给 ellip。

另请参见

ellip

使用阶次和临界点进行滤波器设计

buttord

从通带和阻带规格中查找阶次和临界点

cheb1ord、cheb2ord
iirfilter

使用阶次和临界频率进行通用滤波器设计

iirdesign

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

示例

设计一个模拟高通滤波器，使通带在 30 rad/s 以上为 3 dB，而在 10 rad/s 时衰减 -60 dB。绘制其频率响应，以灰色显示通带和阻带约束。

>>> from scipy import signal
>>> import matplotlib.pyplot as plt
>>> import numpy as np

>>> N, Wn = signal.ellipord(30, 10, 3, 60, True)
>>> b, a = signal.ellip(N, 3, 60, Wn, 'high', True)
>>> w, h = signal.freqs(b, a, np.logspace(0, 3, 500))
>>> plt.semilogx(w, 20 * np.log10(abs(h)))
>>> plt.title('Elliptical highpass filter fit to constraints')
>>> plt.xlabel('Frequency [radians / second]')
>>> plt.ylabel('Amplitude [dB]')
>>> plt.grid(which='both', axis='both')
>>> plt.fill([.1, 10,  10,  .1], [1e4, 1e4, -60, -60], '0.9', lw=0) # stop
>>> plt.fill([30, 30, 1e9, 1e9], [-99,  -3,  -3, -99], '0.9', lw=0) # pass
>>> plt.axis([1, 300, -80, 3])
>>> plt.show()
../../_images/scipy-signal-ellipord-1.png
本页面