
close all;
load dataicann;

% Obtencion de coeficientes del filtro fir
fm = 5000;                  % Frecuencia de muestreo
tm = 1/fm;                  % Periodo de muestreo
fc = 50;                    % Frecuencia de corte (Hz)
wc = fc/(fm/2);             % Frecuencia de corte (normalizada 0=0, 1=frec_shannon)
b  = fir1(250,[wc]);        % Coeficientes del fir


% CALIBRADO (PCA)
Registro = 3;
ac = z{Registro}(:,1); 
ax = z{Registro}(:,2); 
ay = z{Registro}(:,3);
ax = ax - mean(ax);
ay = ay - mean(ay);
ac = ac - mean(ac);
axf = filter(b,1,ax);
ayf = filter(b,1,ay);

axf = axf - mean(axf);
ayf = ayf - mean(ayf);
p     = [axf ayf]';
Cpp   = p*p'/length(p);
[v,d] = eig(Cpp);


% ANALISIS DE ORBITAS
figure;
for Registro = 1:9
    ac = z{Registro}(:,1); 
    ax = z{Registro}(:,2); 
    ay = z{Registro}(:,3);
    ax = ax - mean(ax);
    ay = ay - mean(ay);
    ac = ac - mean(ac);
    axf = filter(b,1,ax);
    ayf = filter(b,1,ay);    
    p    = [axf ayf]';
    ppca = sqrt(d^-1)*v*p;    
    subplot(3,3,Registro);    
    plot(ppca(1,:),ppca(2,:)); 
    title(nom(Registro));
    grid on;
    axis equal;
    axis(15*[-1 1 -1 1]);
end


