Gammatone人耳滤波器

人耳滤波器 Gammatone带通滤波器

function [ dataOut ] = gammatone（ dataInfsfi ）
%GAMMATONE Summary of this function goes here
%   Detailed explanation goes here
% fi ： 中心频率
% fs : 采样频率
% damorn 2014/5
[NsMs] = size（dataIn）;
bi = 1.019*24.7*（4.37*fi/1000+1）;
N = 1024;
t = [0:N-1]/fs;
h = t.^3.*exp（-2*pi*bi*t）.*cos（2*pi*fi*t+0）;
dataOut = zeros（Ns+N-1Ms）;
for i = 1:Ms
dataOut（:i） = conv（hdataIn（:i））;
end
dataOut = dataOut/max（max（dataOut））;
%sound（datafs）;
end

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件         467  2014-05-12 21:55  readme.txt
     文件         467  2014-05-12 21:52  gammatone.m