scipy.signal.

ellipord

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

橢圓 (Cauer) 濾波器階數選擇。

傳回最低階數位或類比橢圓濾波器的階數，該濾波器在通帶中損失不超過 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，選用

若為 True，則傳回類比濾波器，否則傳回數位濾波器。

fsfloat，選用

數位系統的取樣頻率。

在 1.2.0 版本中新增。

傳回值:

ordint

符合規格的最低階橢圓 (Cauer) 濾波器階數。

wnndarray 或 float

用於 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 [rad/s]')
>>> 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()

返回在新分頁中開啟