Фильтр Хеннинга в пространственной области
Итак, фильтр Хеннинга это:
H(k)=(1+cos(2*pi*k/W))/2 , W=2*kmax
и на основании этого предполагается, что обратный фильтр Фурье Хеннинга:
h(x) = 1/2delt(x) + 1/4(delta(x-1)+delta(x+1))
, где delta(.) является дельта-функцией Дирака. Так что если Hk а также hx обозначает H(k) а также h(x)соответственно:
n = 128;
Hk = hann(n,'symmetric')';
hx = ifft(ifftshift(Hk));
Я проверял hx а также Hk на примере сигнала 1D здесь, и результат, кажется, работает, как я ожидал. Тем не мение, hx не выглядит как то, что я ожидал, что является суммой трех дельта-функций Дирака. Чего не хватает в моем коде?
