《Matlab在《自动控制原理》中的应用实例》由会员分享，可在线阅读，更多相关《Matlab在《自动控制原理》中的应用实例（9页珍藏版）》请在装配图网上搜索。
1、G1=tf(1,110);G2=tf(1,11);G3=tf(101,144);numg4=11;deng4=16;G4=tf(numg4,deng4);H1=zpk(-1,-2,1);numh2=2;denh2=1;H3=1;nh2=conv(numh2,deng4);dh2=conv(denh2,numg4);H2=tf(nh2,dh2);sys1=series(G3,G4);sys2=feedback(sys1,H1,+1);(结构图化简及闭环传递函数的求解)sys3=series(G2,sys2);sys4=feedback(sys3,H2);sys5=(G1,sys4);sys5=s
2、eries(G1,sys4);sys=feedback(sys5,H3)/zeta=0.707;num=16;den=18*zeta16;sys=tf(num,den);p=roots(den)t=0:0.01:3;(特征根)figure(1)impluse(sys,t);gridimpulse(sys,t);grid(单位脉冲响应)xlabel(t);ylable(c(t);title(stepresponse);xlabel(t);ylabel(c(t);title(stepresponse);figure(2)step(sys,t);grid(单位阶跃响应)xlabel(t);ylabe
3、l(c(t);title(stepresponse);figure(3)u=t;lsim(sys,u,t,0);grid(任意输入响应)xlabel(t);ylabel(c(t);title(stepresponse);steprssponsestepresponseK=20,40;Gc=tf(1,conv(1,0,1,2,10);fori=1:2(开环增益对系统的影响)G1=tf(K(i)*1,1,1,5);G0=series(G1,Gc);G=feedback(G0,1);figure(i);bode(G);gridEnd.rpiiiIIiI11IIi|ihiiifiIaIIliIIII|
4、IIIIIIIIII4IIIIIIIIIIiiIIMilliij|iIaiI14ITI1iiiinii11i|4|i111I1ii|in|i|i|ii|iiiiaana1il11I1llIllI141ii41I4111IIH1111111II1IIIII1II1l1111l1M1Il4111114Iluiihiiiiiiiinciicicicia1i1i1HilHil1ii1i111H11114111H1IIIIIIIIIIIIIII1!*C.inhiaifnVrrnliirIiFrri1111irsjhiiiianiihiiiiiiiiiiirIrit11H1IIII4IIlIII14l11
5、Mi11111111114111114t111iIIInI11Iiiiiiiihi111II11F-LJ11111M1111911111IIIIPIIaIlllI!1iil41Mi1111l11V11la1P411*l111iiibiiqiiii1i:iiiiiijiiAiiiii1I!1I：!1!LI一7-IJ.LLLJl.-.X-i.LLLJ1J.3.|j_|_fl_L11111iiiiaihiqiiiihi11iiiki11iiiiiia)|iiij|ij|iti|ii|iiTjtiiiitititijHi|iil1iIiIIIiiI1riiiii!iI1i11111iirs.Ii111
6、41Iii1iII|I0|II|I111|11HI11111|l111|111TXI1I11IIHIIiI1l14l11l111l41I44III*!iiiini11ii11II111111iiiiiiiiiiiiiiiiiha11iiiirXji111H1!14l11Miiiiiii11if1111ri巾FlV1iiaiiaii1iiiiiia1iiiihiiiiiiiiiihaiiinciiiiIIIIII11IIii11ni11V1V11IlVIi11411I11ilJ1AIIUIIII_11jjiiiji1ii,iiiiliiii.r5010So190sso1(rasp心君一lCLFr
7、equency(rad/secnu-2-80o o -9I-8 岩)芸mild-270-2-10110101010Fr&ciueney (rscJfeae)1DG=tf(1280,640,1,2421604.81,320.24,16);figuremargin(G);(对数幅频曲线与实轴交点横坐标为截止频率，相频曲线与-180交点为幅值1)figure(2)nyquist(G);(奈奎斯特曲线，在开环幅相曲线或对数幅频曲线基础上补画使闭合)Axisequalbo de DiagramGm = 29.5 dS (art 30.9 rad/sec) b Pm = 72.9 deg (at 0.90
8、4 radJsec3o o D I 8 8 7 - 1 2一 (我s翳罡d231010-2a1io- w io mMs u_bee-Frequency Crgd/c)/G1=tf(2000,-4000,conv(1Q0,1,1);伯德图即对数频率分布图，包括对数幅频和对数相频G2=tf(1,1,10,400);G=series(G1,G2);bode(G);grid/BodeDiagram701I101010Frequencvfrad/sec)90Joo8712-b63(6岂asmudG=tf(1,4,1,4,20);(零极点分布图和开环根轨迹图)figurepzmap(G);figure(2)rlocus(G);Pole-ZeroWapRealAxisRootLocus9-8-7-60RealAxis
Matlab在《自动控制原理》中的应用实例
