полосовой фильтр с щебечут входного сигнала

[magnetra]

Я пытался подражать выходе из фильтра Gammatone для гиперболического щебечут качестве входных данных.MATLAB сценария показан ниже.
Я вызвала гиперболического щебечут с частотой сметая с 100kHz до 20 кГц в 1 мс.И фильтр Гамма тон центральной частотой 50 кГц и полосой пропускания 20 кГц.Это значит, теоретически выходе фильтра должно пика, когда входной частоты показов щебечут 50 КГц и это происходит время @ T = 250us.
Вместе с тем, запутанные результат (выход из GTF), представляется на пик гораздо раньше, чем 250us которые совершенно неожиданным.

Пожалуйста, помогите мне решить эту проблему.

Код:Fs = 500e3;% Частота Freq

T = 1e-3;% 1ms время развертки

T = 0:1 / Fs: T;% Входной Гиперболические щебечут поколения

% Щебечут Sweeps частот от F1 до F2

% Начало и конец щебечут амплитуды формы с окном Ханнинг

f1 = 100e3;

F2 = 20e3;

B = F2 * T / (F1-F2);

A = B * F1;

Freq = a / (T B);

wHann = Ханн (101) ';

WAMP = [wHann (1:51) Ones (1, длина (T) -101) wHann (52:101)];

X = COS (2 * PI * Freq .* T) .* WAMP;% Гиперболические щебечут с утра по Ханнинг окно% GTF impuse ответ

% Gammatone фильтр с полосой пропускания B и центр ФК Freq

K = 1;% прироста

N = 4;% Выбранная надлежащим прав слуховой модели (Patterson 90)

B = 20e3;% от ширины канала в Гц

FC = 50e3; частота% ЦентрОтвет% impuse

H = K * T. ^ (N-1) .* EXP (-2 * PI * B * T) .* COS (2 * PI * FC * T);% Выход из GTF

у ^ = Conv (H, X) / Ф;

ТГ = (0: длина (У) -1) / Ф;

Рисунок;

участок (ТГ, Ю. Г.);

Название ( 'GTF Output');