



for i=1:23
   figure(i);
   
   [y, t, d] =  impulse(tf(Dati(1,4).t.sseh_yy_norm.Mnpf.b(i,:),Dati(1,4).t.sseh_yy_norm.Mnpf.a,1));
   [y1, t, d] =  impulse(tf(Dati_PEM(1,4).t.pem_norm.mod(1,6).Mpem_BIC.b(i,:),Dati_PEM(1,4).t.pem_norm.mod(1,6).Mpem_BIC.a,1));

   plot(y,'b');
   hold on
   plot(y1,':r');

   
   
%   if sum(y)==0
%       close(i)
%    else
%        S = tf(Dati(1,3).t.sseh_yy_norm(3,1).Mnpf.b(i,:),Dati(1,3).t.sseh_yy_norm(3,1).Mnpf.a,1);
%        [S G] = balreal(ss(S));
%        figure
%        bar(G)
%        pause
%        close
%        
%   end
   
    
end
