        %prova FPE o AIC

        d = iddata(y',u');
        %m = n4sid(d,size(d.y,2),'N4Weight','CVA'); é instabile!!
        %m_n4 = n4sid(d,size(d.y,2),'focus','simulation');
        
        %d5 = iddata(y_id',u_id');
        %m = n4sid(d,size(d.y,2),'N4Weight','CVA'); é instabile!!
        %m1_n4 = n4sid(d5,size(d.y,2),'focus','simulation');
        %risposta in frequenza e spettro
        %ffplot(m_n4)
        %mappa zero/poli per LTI
        %pzmap(m_n4)
        
        m_5 = n4sid(d,5,'focus','simulation');
        m_7 = n4sid(d,7,'focus','simulation');
        m_8 = n4sid(d,8,'focus','simulation');
        m_9 = n4sid(d,9,'focus','simulation');
        m_13 = n4sid(d,13,'focus','simulation');
        m_20 = n4sid(d,20,'focus','simulation');
        
        aic(m_5, m_6, m_7, m_8, m_9)
        fpe(m_5, m_6, m_7, m_8, m_9)
        