function filteredX = Filter(X,f,limits)


% X, f -> [X,f] = Espetro(x,Ta);
% limits -> [400 600] (ex)

N=length(f);

%freqs que queremos do sinal
min = limits(1);
max = limits(2);

H = zeros(N,1); %or (1,N) depende da situacao
H(f > min & f < max) = 1; %o que queremos
H(f > -max & f < -min) = 1; % parte negativa

figure
plot(f, H)
title('Filter');
xlim([-2e3 2e3]);
ylim([0 2]);
  
filteredX = H.*X;

end