scipy.signal.

cheb2ord#

scipy.signal.cheb2ord(wp, ws, gpass, gstop, аналог=False, fs=None)[источник]#

Выбор порядка фильтра Чебышева типа II.

Возвращает порядок цифрового или аналогового фильтра Чебышёва типа II самого низкого порядка, который теряет не более gpass дБ в полосе пропускания и имеет не менее gstop Затухание в дБ в полосе задерживания.

Параметры:

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 являются угловыми частотами (например, рад/с).

gpassfloat

Максимальные потери в полосе пропускания (дБ).

gstopfloat

Минимальное затухание в полосе задерживания (дБ).

аналогbool, необязательно

Если True, возвращает аналоговый фильтр, в противном случае возвращается цифровой фильтр.

fsfloat, опционально

Частота дискретизации цифровой системы.

Добавлено в версии 1.2.0.

Возвращает:

ordint: Наименьший порядок для фильтра Чебышева типа II, удовлетворяющего спецификациям.
wnndarray или float: Натуральная частота Чебышёва (частота «3 дБ») для использования с cheby2 для получения результатов фильтрации. Если fs указан, это в тех же единицах, и fs также должен быть передан в cheby2.

Смотрите также

cheby2: Проектирование фильтра с использованием порядка и критических точек
buttord: Найти порядок и критические точки по спецификациям полосы пропускания и полосы задерживания
cheb1ord, ellipord
iirfilter: Общее проектирование фильтров с использованием порядка и критических частот
iirdesign: Общее проектирование фильтра с использованием спецификаций полосы пропускания и задерживания

Примеры

Спроектируйте цифровой полосно-заграждающий фильтр, который ослабляет сигнал на -60 дБ в диапазоне от 0.2*(fs/2) до 0.5*(fs/2), оставаясь в пределах 3 дБ ниже 0.1*(fs/2) или выше 0.6*(fs/2). Постройте его частотную характеристику, отобразив ограничения полосы пропускания и задерживания серым цветом.

>>> 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