设系统为sys = 40*(s^2+s+1)/(s*(2s+1)(0.2s+1)(0.05*s+1))，代码如下

s = tf('s');
sys = 40*(s^2+s+1)/(s*(2*s+1)*(0.2*s+1)*(0.05*s+1))
res = freqresp(sys, w);
amp_true = abs(res)
amp_db = 20*log10(abs(res))
angle_true = angle(res)*180/pi

其中res为复数形式的传函结果，amp_true和angle_true分别为幅值和相角（单位为度deg），amp_db为以dB为单位的幅值