function y = pasoalto1(x, fc, srate)
% PASOALTO1  Filtro paso alto de un polo, causal, escrito a mano.
%   y = pasoalto1(x, fc, srate) filtra el vector x con frecuencia de corte fc.
%
%   Sirve para VER lo que hace un filtro, no para analizar datos: recorre la
%   senal solo hacia adelante, asi que ademas de bajar la amplitud corre el
%   componente en el tiempo. Los filtros de ERPLAB pasan hacia adelante y hacia
%   atras (no causales), lo que cancela ese corrimiento pero no la perdida de
%   amplitud ni el valle que el filtro fabrica despues del pico.

    a = 1 / (1 + 2*pi*fc/srate);
    y = zeros(size(x));
    for n = 2:numel(x)
        y(n) = a * (y(n-1) + x(n) - x(n-1));
    end
end