
for f=1:1
% UMIDITA' MEDIA 0 VAR 1
for k=1:length(Dati)
  
    y = Dati(k).h.y_org; 
    U = Dati(k).h.Udati_org;   
    
    %%%%%%%%%%%%%%%%%%%%%%%%%%%      
    varU = var(U);
    meanU = mean(U);

    for i=1:size(U,2)
       u(:,i) = (U(:,i)-meanU(i))/sqrt(varU(i));    
    end
    Dati(k).h.Udati_norm = u(:,:); % <- aggiunta del campo

    varY = var(y);
    meanY = mean(y);
    
    Dati(k).h.y_norm=(y(:,1)-meanY(1))/sqrt(varY(1));  % <- aggiunta del campo
end

save(str,'Dati')

% TEMPERATURA MEDIA 0 VAR 1
for k=1:length(Dati)
  
    y = Dati(k).t.y_org; 
    U = Dati(k).t.Udati_org;   
    
    %%%%%%%%%%%%%%%%%%%%%%%%%   
    varU = var(U);
    meanU = mean(U);

    for i=1:size(U,2)
       u(:,i) = (U(:,i)-meanU(i))/sqrt(varU(i));    
    end
    Dati(k).t.Udati_norm = u(:,:); % <- aggiunta del campo

    varY = var(y);
    meanY = mean(y);
    
    Dati(k).t.y_norm=(y(:,1)-meanY(1))/sqrt(varY(1));  % <- aggiunta del campo    
end

% ora salvo la Dati_Matrice a cui ho aggiunto i due campi.
save(str,'Dati')

% TEMPERATURA MEDIA 0 VAR !=1
% si aggiunge l'ident. sseh sulla temperatura con modello 'yy' e 'ny'
for k=1:length(Dati)
  
    y = Dati(k).t.y_org; 
    U = Dati(k).t.Udati_org;   
    
    %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%   
    meanU = mean(U);

    for i=1:size(U,2)
       u(:,i) = (U(:,i)-meanU(i));    
    end
    Dati(k).t.Udati= u(:,:); % <- aggiunta del campo

    meanY = mean(y);
    
    Dati(k).t.y = (y(:,1)-meanY(1));  % <- aggiunta del campo    
end

% ora salvo la Dati_Matrice a cui ho aggiunto i due campi.
save(str,'Dati')


% UMIDITA' MEDIA 0 VAR !=1
% si aggiunge l'ident. sseh sulla temperatura con modello 'yy' e 'ny'
for k=1:length(Dati)
  
    y = Dati(k).h.y_org; 
    U = Dati(k).h.Udati_org;   
    
    %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
    meanU = mean(U);

    for i=1:size(U,2)
       u(:,i) = (U(:,i)-meanU(i));    
    end
    Dati(k).h.Udati= u(:,:); % <- aggiunta del campo

    meanY = mean(y);
    
    Dati(k).h.y = (y(:,1)-meanY(1));  % <- aggiunta del campo    
    %%%%%%%%%%%%%%%%%%%%%%%%%%% 
end


save(str,'Dati')

end % for inutile. solo per chiudere questa parte di codice
% fine normalizzazione dei dati.
 

%% Inizio identificazione dei dati:


% contiene una struttura Dati che viene via via aggiornata con i risultati
% dei vari metodi di indentificazione. Poi però viene salvata in un fil
% chiamato Dati_Matrice_id

len = length(Dati);

i_end = 2406;    % fine dei dati utilizzati

Ts=60*5; % tempo di campionamento

TEM = true % mettere true se si vuole effettuare una id sulla temper 
HUM = false % mettere true se si vuole effettuare una id sulla umidità

%Per scegliere algoritmo di identificazione
SSEH = true
PEM = true



% NOR 1: uso dati con media 0 e var 1 (con nome_norm)

NOR=1;
switch NOR
    case 1
        for k=1:len 
            Dati(k).t.y_id = Dati(k).t.y_norm;
            Dati(k).t.Udati_id = Dati(k).t.Udati_norm;   

            Dati(k).h.y_id = Dati(k).h.y_norm;
            Dati(k).h.Udati_id = Dati(k).h.Udati_norm;   
        end

end
% da qui si usano solo i due nuovi campi che contengono i dati con media  e
% varianza a scelta


