matlab绘制系统函数波特图
- 格式:doc
- 大小:68.50 KB
- 文档页数:3
已知系统函数:H(s)=30/(s^2+31s+30),先画出幅频特性和相频特性view plaincopy to clipboardprint?
1. w=-8*pi:0.01:8*pi;
2. b=[30];
3. a=[1,31,30];
4. H=freqs(b,a,w);
5. s ubplot(211)
6. p lot(w,abs(H));grid on
7. x label('/omega(rad/s)'),ylabel('|H(/omega)|');
8. t itle('H(s)的幅频特性')
9. s ubplot(212)
10. plot(w,angle(H));grid on
11. xlabel('/omega(rad/s)'),ylabel('/phi(/omega)');
12. title('H(s)的相频特性')
然后再画出波特图
view plaincopy to clipboardprint?
1. f igure
2. s ys=tf(b,a);
3. b ode(sys);grid on