freson - peak frequencies
returns the vector of peak frequencies in Hz for the SISO plant h
h=syslin('c',-1+%s,(3+2*%s+%s^2)*(50+0.1*%s+%s^2)) fr=freson(h) bode(h) g=20*log(abs(repfreq(h,fr)))/log(10)