
clc;

% confronto ingressi (fa una struttura e poi in uscita un txt)

load Dati_Matrice_chiuso_sseh_Mnpe.mat;
Dati_sseh = Dati;

clear Dati
load Dati_Matrice_sseh_chiuso_vinc
Dati_sseh2_vinc = Dati;

clear Dati
load Dati_Matrice_sseh_chiuso_vinc00
Dati_sseh2 = Dati;

clear Dati
load Dati_Matrice_sseh_chiuso_vinc00
Dati_sseh2 = Dati;


% stampiamo per ogni sensore (la norma di) gli ingressi usati (per gli Mnpe)
for k=1:length(Dati)
    
    in_sel_tmp = [];
    for i=1:23
        norm_i = norm(Dati_sseh(1,k).t.sseh_yy_norm.Mnpe(1,i).b);
        in_sel_tmp = [in_sel_tmp norm_i];
        
%         %calola le rms onvece di norme di risposta impulsiva
%         rms_i = sqrt(mean((Dati_sseh(1,k).t.sseh_yy_norm.Mnpe(1,i).b).^2));
%         if isnan(rms_i)
%            rms_i=0; 
%         end
%         in_sel_tmp = [in_sel_tmp rms_i];
    end
    
    Conf_ingressi.in_sseh(1,k).in_selez = in_sel_tmp;

end



% % stampa di tutte le risp impulsive non nulle
% for j=1:23
%     [y t x]=impulse(Dati(1,3).t.sseh_yy_norm.Mnpe(1,j));
%     hold on
%     figure();
%  
%     if sum(y)==0
%        close(); 
%     else
%         plot(y)
%         legend(num2str(j))
%         pause;
%         close();
%     end
%    
% end


% per sseh2 con ingressi non vincolati
for k=1:length(Dati)
    in_sel_tmp = [];
    for i=1:23
        norm_i = norm(Dati_sseh2(1,k).t.sseh_yy_norm.Mnpe(1,i).b);
        in_sel_tmp = [in_sel_tmp norm_i];
    end
    
    Conf_ingressi.in_sseh2(1,k).in_selez = in_sel_tmp;

end


% per sseh2 con ingressi vincolati
for k=1:length(Dati)
    in_sel_tmp = [];
    for i=1:23
        norm_i = norm(Dati_sseh2_vinc(1,k).t.sseh_yy_norm.Mnpe(1,i).b);
        in_sel_tmp = [in_sel_tmp norm_i];
    end
    
    Conf_ingressi.in_sseh2_vinc(1,k).in_selez = in_sel_tmp;

end


%% 
% stampa in txt

fid = fopen('ValidazioneChiuso\conf_norme_ingressi.txt', 'w+');

caso = ['Mnpe sseh' '\t'];
fprintf(fid, [ caso, '\n']);    
for i=1:length(Dati)
    v = ['Sens \t'];
    for k=1:23
       v = [v num2str(Dati_sseh(1,i).t.Colonne_U(k)) '\t' ]; 
    end
    fprintf(fid, [v  '\n']);
    v=[];
    for j=1:23
       v = [v num2str(Conf_ingressi.in_sseh(1,i).in_selez(j),3) '\t' ]; 
    end
    fprintf(fid, [num2str(i) '\t' v  '\n']);
end 

caso = ['\n' 'Mnpe sseh2' '\t'];
fprintf(fid, [ caso, '\n']);
for i=1:length(Dati)
    v = ['Sens \t'];
    for k=1:23
       v = [v num2str(Dati_sseh2(1,i).t.Colonne_U(k)) '\t' ]; 
    end
    fprintf(fid, [v  '\n']);
    v=[];
    for j=1:23
       v = [v num2str(Conf_ingressi.in_sseh2(1,i).in_selez(j),3) '\t' ]; 
    end
    fprintf(fid, [num2str(i) '\t' v   '\n']);
end

caso = ['\n' 'Mnpe sseh2 vincolati' '\t'];
fprintf(fid, [ caso, '\n']);
for i=1:length(Dati)
    v = ['Sens \t'];
    for k=1:23
       v = [v num2str(Dati_sseh2_vinc(1,i).t.Colonne_U(k)) '\t' ]; 
    end
    fprintf(fid, [v  '\n']);
    v=[];
    for j=1:23
       v = [v num2str(Conf_ingressi.in_sseh2_vinc(1,i).in_selez(j),3) '\t' ]; 
    end
    fprintf(fid, [num2str(i) '\t' v '\n']);
end

fclose(fid);

disp('fine');








    