控制系统仿真-答案-薛定宇-第二版

第一局部：MATLAB必备根底知识、控制系统模型与转换、线性控制系统的计算机辅助分析(4学时)1.2.3.答：〔a〕>>s=tf('s');G=(s^3+4*s+2)/(s^3*(s^2+2)*((s^2+1)^3+2*s+5))Transferfunction:s^3+4s+2------------------------------------------------------s^11+5s^9+9s^7+2s^6+12s^5+4s^4+12s^3〔b〕>>num=[100.568];den=[1-1.21.19-0.99];H=tf(num,den,'Ts',0.1)Transferfunction:-----------------------------4.设y〔t〕=x1；A=[0,1,0;0,0,1;-13,-4,-5];B=[0;0;2];C=[1,0,0];D=0;G=SS(A,B,C,D)a=x1x2x3x1010x2001x3-13-4-5b=u1x10x20x32c=x1x2x3y1100d=u1y10Continuous-timemodel.传递函数模型：对式子两联进行拉氏变换，G〔s〕=2/(s^3+13*s^2+4*s+5)对分子分母进行分解因式处理5.进行Z变换6.s=tf('s');symsJKpsKi;G=(s+1)/(J*s^2+2*s+5);Gc=(Kps+Ki)/s;H=1;GG=feedback(G*Gc,H)7.〔a〕s=tf('s');G=(211.87*s+317.64)/((s+20)*(s+94.34)*(s+0.1684));Gc=(169.6*s+400)/s/(s+4);H=1/(0.01*s+1);GG=feedback(G*Gc,H)Transferfunction:359.3s^3+3.732e004s^2+1.399e005s+127056----------------------------------------------------------------------------------0.01s^6+2.185s^5+142.1s^4+2444s^3+4.389e004s^2+1.399e005s+1270568．Ga=feedback(1/s^2,50);G1=feedback((1/(s+1))*(s/(s^2+2)),(4*s+2)/((s+1)^2));G2=feedback(Ga*G1,(s^2+2)/(s^3+14));G=3*G2Transferfunction:3s^6+6s^5+3s^4+42s^3+84s^2+42s-----------------------------------------------------------------------------------------------------s^10+3s^9+55s^8+175s^7+300s^6+1323s^5+2656s^4+3715s^3+7732s^2+5602s+14009.>>s=tf('s');G=(s+1)^2*(s^2+2*s+400)/(s+5)^2/(s^2+3*s+100)/(s^2+3*s+2500)Transferfunction:s^4+4s^3+405s^2+802s+400-------------------------------------------------------s^6+16s^5+2694s^4+34040s^3+393225s^2>>step(G)s=tf('s');G=(s+1)^2*(s^2+2*s+400)/(s+5)^2/(s^2+3*s+100)/(s^2+3*s+2500);G1=c2d(G,0.01)Transferfunction:------------------------------------------------------------------------------------------>>step(G1)>>s=tf('s');G=(s+1)^2*(s^2+2*s+400)/(s+5)^2/(s^2+3*s+100)/(s^2+3*s+2500);G1=c2d(G,0.1)Transferfunction:0.0003982z^5-0.0003919z^4-0.000336z^3-------------------------------------------------------z^6-2.644z^5+4.044z^4-3.94z^3+2.549z^2>>step(G1)>>s=tf('s');G=(s+1)^2*(s^2+2*s+400)/(s+5)^2/(s^2+3*s+100)/(s^2+3*s+2500);G1=c2d(G,1)Transferfunction:8.625e-005z^5-4.48e-005z^4+6.545e-006z^3-------------------------------------------------------z^6-0.0419z^5-0.07092z^4-0.0004549z^3Samplingtime:1>>step(G1)10.〔a〕>>s=tf('s');G=1/(s^3+2*s^2+s+2);eig(G)ans=-2.0000临界稳定〔b〕>>s=tf('s');G=1/(6*s^4+3*s^3+2*s^2+s+1);eig(G)ans=不稳定〔c〕>>s=tf('s');G=1/(s^4+s^3-3*s^2-s+2);eig(G)ans=不稳定11.〔a〕>>z=tf('z');G=(-3*z+2)/(z^3-0.2*z^2-0.25*z+0.05);abs(eig(G))ans=稳定〔b〕z=tf('z');G=(3*z^2-0.39