% A=[1 -0.5 0.7];B=[0 1 0.5];
% th0=poly2th(A,B);
% u=idinput(300, 'rbs');
% y=idsim([u,randn(300,1)],th0);
% z=[y u];
% ge = etfe(z);
% gs = spa(z);
% bode(ge,gs)

m = idpoly([1 -1.5 0.7],[0 1 0.5]);
u = iddata([],idinput([50,1,10],'sine'));
u.Period = 50;
y = sim(m,u);
me = etfe([y u])
bode(me,'b*',m)