if ~exist('data')
    importa_tutto()
    % TOGLI LA MEDIA? (xcov lo fa da solo, solo x i grafici)
    % data.T = data.T - repmat(mean(data.T),length(data.names),1);
end

durata = 120; % un ora
R = zeros(length(data.names),length(data.names));
ins = 1:length(data.names);
for ix = ins
for k = ins
    [r,l] = xcov(data.T(ix,1:20:end),data.T(k,1:20:end),durata,'coeff');
    R(ix,k) = r(durata);
end
end
%%
imagesc(exp(R))
colormap

%%
[R2,ix_R2] = sort(R);

%%
global startTime dt
startTime = datenum('Jun 15, 2007  15:45:00.0000'); % anno mese giorno ora minuti secondi
dt = datenum('Jan 0, 0000  0:00:30.0000'); % tempo di campionamento


t = 0:length(data.T)-1;
figure(2)
for k = 1:length(data.names)
    indici = flipud(ix_R2(end-3:end,k));
plot(t,data.T(indici,:))
xlabel('n')
ylabel('°C')
title(sprintf('Temperature, andamenti più simili a %s',data.names{k}))
legend(data.names(indici))
axis tight
dcm_obj = datacursormode(gcf);
set(dcm_obj,'UpdateFcn',@myfunction)
pause
end