系统仿真上机作业.docxVIP

系统仿真上机作业 姓名：屈佳龙 学号班级：041032 系统仿真上机作业 r :单位阶跃 Gn :非线性器件 Go Go(s)二 s(10s 1)(0.625s 1)(0.025s 1) GN =1 , K=40时，用MATLAB 画出开环Bode图，求出仁c、0po由其估 计出 tr、ts、6% 答： den=conv(conv(conv([1 0],[10 1]),[0.625 1]),[0.025 1]); num=[40 40]; sys=tf(num,den); sys=feedback(sys,1); [y,t]=step(sys); ytr=find(y=1); [ym,tp]=max(y); 0.7772 （上升时间） 0.7772 （上升时间） rise_time = 1.2953 （峰值时间 1.2953 （峰值时间） peak_time = 0.8168 (超调量) 0.8168 (超调量) max_overshoot = s=length(t); while y(s)0.98y(s)1.02 s=s-1; end settling_time=t(s+1) settling_time = 17.2280 (调整时间) bode(num,den) [mag,phase,w]=bode(num,den); [gm,pm,wcg,wcp]=margin(mag,phase,w) gm = 4.3168 (增益裕量) pm = 10.0158 (相位裕量) wcg = 5.1598 (交界频率) wcp = 2.3975 (增益=0dB 时) Bode图如下： I Il I V V I| IIi1 I Il I V V I| II i 1 t VR t ? 1 B ■ II 1 1 1 II j 4 4 t *1 It 9 1 1 4 i| |i| a 1 1 I f 111 4 ? I P IP r 7 I 1 ?ll 4 ■ Illi! 4 * 1 4 , 1 1 1 ll 9 111 p p || II I 1 P VI 1 1 I 11 V i| |i |i i| Pfi ■ l ? ■ 1 ■ 14 IP g ■ivan 4 * | 4 p I i|i| t V I V 1 1 ?l I* .二 11 I 4 N ■ m 1 ? ! 4 4 1 1 11i ? 4111 4 ? 1 ■ t ■ tl V V 9 I i 4 1 VI ?I* 1 p ■I |i | 4 “ ■ l i ■ 1 I 1 1 1 P i 4 1 F1V 4 ? I . R 1 Il 申 ri ■ ■ 4 ■ 4 ■■断 l?M *■ ■ ■? * T皂营I* ■ ■ ■■ irr ? ■ ■ ar I H a i * 4 1 II R ?? * 1 1 1 1 I 4 | | 4 1 1 ■ I 4 4 1111 1 ■ 4 1 tIE 4 ? 1 it t ■ 1l P f 4 F i 1 P fl f ■ ? 4 1 1 11 4 iim l i 1 1 i I f IF i 1 1 MR 4 ? ! ? P | |i| i it a i ? 1 I t i 1 1 1 II V J里 1 1 1 14 11 ■ ■ 41 tit 4 * 1 ■ 1 1 ll I ? I i i 1 ■ II ■ 1 ? 1 1 1 II i i 4 4 rrr— ■1 1 1 1 14 11 1 4 1 MG 1 * 1 ■ 1 1 ll I l 4 I i 1 1 ?l h 1 l 1 1 1 II i 4 ? 1 4 M 4 1 1 V 1 H , ■ ■ 41 III 4 ? 1 4 1 1 ll 卜―■ * 4 , i*. ■ ■. 4 -*■ ? * ■- ■ d A ■ —. 触 * M *■- ■■*■? fi fi i I ? 4 d H 4 1 i 1 1 1 II i 4 a i oa a l ■ 1 1 a 4 1 MG 4 i 1 4 I 1 19 I ii a ii ? 4 ■ , l * 1 g ft 1 e l! 4 4 4 Od d 1 ? ?4 4 ■ ■ @甘、 i , 41+11 4 * 1 i b 1 IV h b I Ii I 1 I fel i F l t * Illi h 4 t 1 rilit d 1 ! I* ri l If Ml tltlt i t i , 1 1 +1 Ii l! 9 l! i 4?” k 1 I* 1 1 1 11 fi 4 a?m 1 b 4 ii I 111 |iV 4 1 IE 4 ? B i I ■ II