load pilotdata.mat
datos1=iddata(Y,U,5);
datos=datos1(300:1200);
datosd=detrend(datos,'constant');
datose=datosd(1:500);
datosv=datosd(501:length(datosd.u));
[model, fit, yhat]=orthobasesid(datose,datosv);
NN=struc(1:10,1:10,1:5);
V=arxstruc(datose,datosv,NN);
nn=selstruc(V,0);
modelo=arx(datose,nn);
figure
compare(datosv,model,modelo)
[model2, fit, yhat]=orthobasesid(datose,datosv,zeros(1,200));
figure
compare(datosv,model,modelo,model2)
[model3, fit, yhat]=orthobasesid(datose,datosv,0.985*ones(1,7));
figure
compare(datosv,model,modelo,model2,model3)