MATLAB校正系统程序.docxVIP

clear; Gl=tf (20, [1ZO.4,0]) %系统的开环传递函数 G=feedback (Gl, lf-l)殆原系统的闭环传递函数 figure (1) step(G) %原系统的单位阶跃响应 figure (2) rlocus (G1) $原系统的根轨迹 figure(3) margin (G1) $原系统的bode图 [GmO,PmO,WgO,WcO]=margin(G1) deta=5 %补偿角 Phim=4 5-Pm0+deta匕校正装置最大超前角 a=sin(Phim*pi/180) alpha=(1+a)/(1-a) Wm=WcO*alpha/s (1/4) T=1/(Wm*sqrt(alpha)); alpha_T=alpha*T Gc=tf ([alpha_Tzl]z [T,l]) %校正装置 Ga=Gl*Gc %校正后的传递函数 Gll=feedback (Glf 1) ; $原系统闭环传递函数 Gaa=feedback (Ga, 1) ; $校正后系统闭环传递函数 figure(4) rlocus (Ga) %校正后系统的根轨迹 legend(Ga) figure(5) margin (Ga)衛校正后系统的频域特性 [Gm,Pm,Wcm,Wpm]=margin(Ga) figure(6) step (Gil, - . k z Gaa) %原系统与校正后系统的时域特性的比较 legend(Gl1,Ga) figure(7) bode(Glz -.kzGc, 1 : \ Ga)%原系统与校正后系统的频域特性的比较 legend ( G1 \ Gc\ ) G1 = Continuous-time transfer function? sA2 + 0.4 s + 20 Continuous-time transfer function? GmO = Inf PmO = 5.1214 WgO = Inf WcO = 4.4631 deta = 5 Phim = 44.8786 0.7056 alpha = 5.7936 Wm = 6.9243 alpha_T = 0.3476 Gc = 0.3476 s + 1 0.06 s + 1 function.0 ? 4 sfunction? function. 0 ? 4 s function? Ga = 6.952 s + 20 0.06 sA3 + 1.024 s入2 + Continuous-time transfer Gm = Inf E6?9 HUIdM JUI H UIOM SZ (spuooes) eluH 5Z 0 Aspuooes) 8sfi LoinleA leuiz o melsxs L6L (spuogs)oE-?-￡as o WSSAS TT 7h~ co o Bsuodsaad￡s 0 945 0 89 -0 976 Root Locus 0 68 03 0.05谕-0.05 (「spuooas)曲 xv AJeuoew- l 0 994 -0 994 r 0 976 0 945 0 35 03 0 89 .45 114 -0.4 -0 35 0 15 0 05 OJ 1 0 68 05 03 -0.3 -025 -02 -0.15 -0.1 Real Axis (seconds 1) 005 0 0.05 (8P) epnuu6es Bode Diagram Gm = Inf dB (at Inf rad⑸，Pm = 5.12 deg (at 4.46 rad/s) 35 1 ? (Bap) aseqd 101 10°Frequency (racVs) 101 校正后 T o S o Rea-Axis (seconds. Rea-Axis (seconds. s_ep Response System Ga seeing 5? (seconds) 0 908 06 Time (seconds) System Ga 58 Ume (seconds^ 0」62 System Ga F-E vaErcx1 100500510090 (8P) 0-90 Bode Diagram Gm = Inf dB (at Inf rad/s), Pm = 48.2 deg (at 6.94 rad/s) 101 101 100500-50100-90 (8P) 9pn_c6es 35 (Bop) heqd 103 Frequency (rad/s)