--> t=to/5; --> prd=dscr(pr,t); --> spec(prd(2)) ans = ! 0.8191692 ! ! 0.8183269 ! ! 0.8189164 + 0.0001594i ! ! 0.8189164 - 0.0001594i ! ! 0.8184613 ! ! 0.8186610 + 0.0002228i ! ! 0.8186610 - 0.0002228i ! ! 0.8187340 ! -->//impulse response --> x0=[0;0;0;0;0;0;0;0]; --> u(1,180)=0;u(1,1)=1; --> y=flts(u,prd,x0); --> draw(0)