% ej4_distorsion.m
% Ejercicio 4 (dia 4) - Documenta tu propia distorsion.  SOLUCION.
%
% Con la onda limpia los numeros son exactos (la ventana tiene 101 muestras):
%   sin filtrar    6.300 uV
%   paso alto 0,1  5.817 uV  (-7.7%)
%   paso alto 1,0  3.043 uV  (-51.7%)
% Con el promedio ruidoso se mueven unos pocos puntos porcentuales.

clear; clc;
rng(42);

srate = 500;
t     = (-100:399) / srate;          % enteros / srate: t(101) es 0 exacto
onda  = 8 * exp(-((t - 0.400).^2) / (2*0.080^2));
vent  = t >= 0.300 & t <= 0.500;

prom = mean(onda + 25 * randn(40, numel(t)), 1);   % el promedio del dia 2

entradas = {onda, prom};
nombres  = {'onda limpia', 'promedio de 40 ensayos'};

for e = 1:numel(entradas)
    x  = entradas{e};
    a0 = mean(x(vent));
    fprintf('\n%s\n', nombres{e});
    fprintf('sin filtrar   : %6.3f uV\n', a0);
    for fc = [0.1 1.0]
        y = pasoalto1(x, fc, srate);
        a = mean(y(vent));
        fprintf('paso alto %3.1f : %6.3f uV  (%.1f%% menos)\n', ...
            fc, a, 100*(1 - a/a0));
    end
end

figure; hold on
plot(t, prom)
plot(t, pasoalto1(prom, 0.1, srate))
plot(t, pasoalto1(prom, 1.0, srate))
plot([t(1) t(end)], [0 0], 'k')
legend('sin filtrar', '0,1 Hz', '1,0 Hz')
xlabel('Tiempo (s)'), ylabel('Amplitud (\muV)')

% Dos frases para la seccion de metodos:
%
% "El EEG continuo se filtro con un paso alto Butterworth no causal de segundo
%  orden, corte a 0,1 Hz (12 dB/octava). Sobre una simulacion del componente
%  medido, ese filtro reduce la amplitud media entre 300 y 500 ms en un 7,6%, y
%  un corte de 1 Hz la habria reducido a la mitad."
%
% Lo segundo casi nunca se reporta, y es lo que permite comparar amplitudes
% entre estudios: un efecto de 2 uV filtrado a 0,1 Hz y uno de 2 uV filtrado a
% 1 Hz no son el mismo efecto.