matlab绘制系统函数波特图

已知系统函数：H(s)=30/(s^2+31s+30),先画出幅频特性和相频特性

w=-8*pi:0.01:8*pi; b=[30]; a=[1,31,30]; H=freqs(b,a,w); subplot(211) plot(w,abs(H));grid on xlabel('/omega(rad/s)'),ylabel('|H(/omega)|'); title('H(s)的幅频特性') subplot(212) plot(w,angle(H));grid on xlabel('/omega(rad/s)'),ylabel('/phi(/omega)'); title('H(s)的相频特性')

 

然后再画出波特图

figure sys=tf(b,a); bode(sys);grid on