关于gammatone的FIR实现

[seacj]

你好 Gammatone滤波器组：性质、实现和应用写的非常好 我想请教关于FIR的实现。 你写道：

g = @(t, f0, erb) t .^ 3 .* exp(- 2 * pi * 1.019 * erb * t) .* cos(2 * pi * f0 * t);
h = g((0:1 / fs:signal_length / fs)', central_freq, band_width);
y = filter(h, 1, x); % y = filter(B, A, x);
y = conv(h, x); % 相当于和输入信号卷积

我想请问这里central_freq，band_width，f0 如何确定