wr,wc]=margin(G) title(' 系统 10/s(0.5s+1) 幅值相角裕度图') grid 运行结果为: Kg=Inf, r=17.9642, Wr=Inf, Wc=6.1685 6 幅值裕度 h=20log Inf=Inf, 相角裕度 y=17.9642 3.3 绘制系统的根轨迹绘制系统的根轨迹程序如下所示: den=[10]; num=[0.5 10]; G=tf(den,num); rlocus(G) 系统校正前的根轨迹图 7 3.4.2 超前校正后的验证 1.用MATLAB 求校正后系统的幅值裕度和相角裕度程序: den=[5 10]; num=conv([0.29 1],[0.5 10]); G=tf(den,num); margin(G); [kg,r,wr,wc]=margin(G); grid 由图可知 Wc=5.38rad/s, 幅值裕度 Kg=Inf, 相角裕度 y=32.6dB 8 3.4.3 系统校正后和校正装置的 Bode 图系统校正后 Bode 图MATLAB 程序: den=[5 10]; num=conv([0.29 1],[0.5 10]); G=tf(den,num); bode(G); title(' 系统校正后 bode 图'); grid 9 校正装置 Bode 图MATLAB 程序: den=[0.50 1]; num=[0.29 1]; G=tf(den,num); bode(G); title(' 校正装置 bode 图'); grid 10 3.5 系统校正前、后的开环的 Nyquist 图 3.5.1 系统校正前开环 Niqust 图 MATLAB 程序: den=[10]; num=[0.5 10]; G=tf(den,num); nyquist(G) title(' 系统校正前 nyqust 图')