本文简要介绍 python 语言中 scipy.signal.gammatone
的用法。
用法:
scipy.signal.gammatone(freq, ftype, order=None, numtaps=None, fs=None)#
Gammatone 滤波器设计。
此函数计算 FIR 或 IIR 伽马调数字滤波器 [1] 的系数。
- freq: 浮点数
滤波器的中心频率(以与 fs 相同的单位表示)。
- ftype: {‘fir’, ‘iir’}
函数生成的过滤器类型。如果‘fir’,该函数将生成一个 N 阶 FIR gammatone 滤波器。如果‘iir’,该函数将生成一个 8 阶数字 IIR 滤波器,建模为 4 阶 gammatone 滤波器。
- order: 整数,可选
过滤器的顺序。仅在
ftype='fir'
时使用。默认值为 4 以模拟人类听觉系统。必须介于 0 和 24 之间。- numtaps: 整数,可选
过滤器的长度。仅在使用时
ftype='fir'
.默认为fs*0.015
如果fs大于 1000, 15 如果fs小于或等于 1000。- fs: 浮点数,可选
信号的采样频率。频率必须介于 0 和
fs/2
.默认值为 2。
- b, a: 数组,数组
滤波器的分子 (
b
) 和分母 (a
) 多项式。
- ValueError
如果频率小于或等于 0 或大于或等于
fs/2
, 如果类型不是 ‘fir’ 或 ‘iir’,如果次序小于或等于 0 或大于 24 时ftype='fir'
参数 ::
返回 ::
抛出 ::
参考:
[1]Slaney, Malcolm,“Patterson-Holdsworth 听觉滤波器组的有效实现”,Apple Computer Technical Report 35, 1993, pp.3-8, 34-39。
例子:
以 440 Hz 为中心的 16 个样本四阶 FIR Gammatone 滤波器
>>> from scipy import signal >>> signal.gammatone(440, 'fir', numtaps=16, fs=16000) (array([ 0.00000000e+00, 2.22196719e-07, 1.64942101e-06, 4.99298227e-06, 1.01993969e-05, 1.63125770e-05, 2.14648940e-05, 2.29947263e-05, 1.76776931e-05, 2.04980537e-06, -2.72062858e-05, -7.28455299e-05, -1.36651076e-04, -2.19066855e-04, -3.18905076e-04, -4.33156712e-04]), [1.0])
IIR Gammatone 滤波器以 440 Hz 为中心
>>> import matplotlib.pyplot as plt >>> import numpy as np
>>> b, a = signal.gammatone(440, 'iir', fs=16000) >>> w, h = signal.freqz(b, a) >>> plt.plot(w / ((2 * np.pi) / 16000), 20 * np.log10(abs(h))) >>> plt.xscale('log') >>> plt.title('Gammatone filter frequency response') >>> plt.xlabel('Frequency') >>> plt.ylabel('Amplitude [dB]') >>> plt.margins(0, 0.1) >>> plt.grid(which='both', axis='both') >>> plt.axvline(440, color='green') # cutoff frequency >>> plt.show()
相关用法
- Python SciPy signal.gausspulse用法及代码示例
- Python SciPy signal.gauss_spline用法及代码示例
- Python SciPy signal.group_delay用法及代码示例
- Python SciPy signal.get_window用法及代码示例
- Python SciPy signal.czt_points用法及代码示例
- Python SciPy signal.chirp用法及代码示例
- Python SciPy signal.residue用法及代码示例
- Python SciPy signal.iirdesign用法及代码示例
- Python SciPy signal.max_len_seq用法及代码示例
- Python SciPy signal.kaiser_atten用法及代码示例
- Python SciPy signal.oaconvolve用法及代码示例
- Python SciPy signal.hilbert用法及代码示例
- Python SciPy signal.ricker用法及代码示例
- Python SciPy signal.cheb2ord用法及代码示例
- Python SciPy signal.lfilter用法及代码示例
- Python SciPy signal.morlet用法及代码示例
- Python SciPy signal.coherence用法及代码示例
- Python SciPy signal.dfreqresp用法及代码示例
- Python SciPy signal.TransferFunction用法及代码示例
- Python SciPy signal.dbode用法及代码示例
- Python SciPy signal.residuez用法及代码示例
- Python SciPy signal.bilinear_zpk用法及代码示例
- Python SciPy signal.firls用法及代码示例
- Python SciPy signal.impulse用法及代码示例
- Python SciPy signal.buttord用法及代码示例
注:本文由纯净天空筛选整理自scipy.org大神的英文原创作品 scipy.signal.gammatone。非经特殊声明,原始代码版权归原作者所有,本译文未经允许或授权,请勿转载或复制。