%scelta degli orizzonti futuri e passati per PBSID_opt
%procedura da effettuare dopo almeno una simulazione del main
   %numero di esperimenti 
   l = 2;
    dat=iddata(YidC',Ucp');
    %uscite
    p=54;
    %ingressi
    m=3;
    for i=1:l,
        order =[i*ones(p,p) i*ones(p,m) ones(p,m)];
        M = arx(dat,order);
        Aic_crit(i)=aic(M);
        fpe_crit(i)=fpe(M);
    end ;
     figure
     plot(Aic_crit)
     hold on
     plot(fpe_crit,':')
     hold off
     drawnow
 
%questi sono gli ordini da inserire in PBSID_opt
     past = floor(find(Aic_crit==min(Aic_crit)));
     futuro = past;