% load dryer2
% dat = iddata(y2,u2,0.08);
% dat1 = merge(dat(1:498),dat(502:719),dat(719:1000))
% m = pem(dat1);
% m1 = pem(getexp(dat1,[1,2]));
% compare(getexp(dat1,3),m1)

load dryer2
z = iddata(y2,u2,0.08);
m = armax(z,[2 2 2 1]);
g = spa(z)
% g = spafdr(z,[],{0,10})
bode(g,m)

g = idfrd(m('m'))
phi = idfrd(m('n'))



% A = [1 -1.5 0.7];
% B = [0 1 0.5];
% C = [1 -1 0.2];
% m0 = idpoly(A,B,C);
% e = iddata([],randn(300,1)+sin(randn(300,1)));
% u1 = iddata([],idinput(300,'rbs',[0 1],[2 4]));
% y1 = sim(m0,[u1 e]);
% z1=[y1,u1];
% A = [1 -1.2 0.5];
% B = [0 0.5 0.6];
% C = [1 -0.8 0.7];
% m0 = idpoly(A,B,C);
% e = iddata([],randn(300,1)+sin(randn(300,1)));
% u2 = iddata([],idinput(300,'rbs',[0 1],[2 4]));
% y2 = sim(m0,[u2 e]);
% z2=[y2,u2];
% m1 = arx(z1,[2 3 4]);
% m2 = arx(z2,[2 3 4]);
% ma = merge(m1,m2)