num=[1 35 291 1093 1700];den=[1 9 66 0 1029 2541 4684 5856 4629 1700];G=tf(num,den);bode(G); % 伯德图（对数频率特性曲线）[z,p,k]=zpkdata(G,'v') % 求零极点
matlab里面有画频幅曲线的公式 先得把这个系统处理下
