function g = filtersinc(PR); % filtersinc.m % % Written by Waqas Akram % % "a": This parameter varies the filter magnitude response. % When "a" is very small (a<<1), the response approximates |w| % As "a" is increased, the filter response starts to % roll off at high frequencies. a = 1; [Length, Count] = size(PR); w = [-pi:(2*pi)/Length:pi-(2*pi)/Length]; rn1 = abs(2/a*sin(a.*w./2)); rn2 = sin(a.*w./2); rd = (a*w)./2; r = rn1*(rn2/rd)^2; f = fftshift(r); for i = 1:Count IMG = fft(PR(:,i)); fimg = IMG.*f'; g(:,i) = ifft(fimg); end g = real(g);