load dryer2
datos=iddata(y2,u2,0.08);
datosd=detrend(datos,'constant');
datose=datosd(1:100);
datosv=datosd(501:1000);
NN=struc(1:10,1:10,1:5);
V=arxstruc(datose,datosv,NN);
nn=selstruc(V,0);
modelo=arx(datose,nn);
%modelo2=arx(datose,[4 3 0]);
nn2=selstruc(V,'aic');
modelo3=arx(datose,nn2);
compare(datosv,modelo,modelo3)