源代码：1nump=[4];2denp=[1684];3sysGp=tf(nump,denp);4den=[10];5t=0:0.01:8;6forK=3:0.2:57fora=0.1:0.1:38num=[K2*K*aK*a^2];9sysG=tf(num,den);10sysGz=feedback(series(sysG,sysGp),1);11sys

代码如下，非常简单，不多说：1k=1;2s=tf('s');3sysG=1/(s+k);45w=logspace(-2,2);67[mag,phase]=bode(sysG,w);8subplot(311)9bode(sysG,{0.01,100})10gridon11subplot(312)12loglog(w,squeeze(mag));13gridon14subplot(313)15semilog