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

[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');