matlab控制系统仿真课程设计.docVIP

上机作业： 请将下面的传递函数模型输入到matlab环境。 s=tf(s); G=(s^3+4*s+2)/((s^3)*(s^2+2)*((s^2+1)^3+2*s+5)); G Transfer function: s^3 + 4 s + 2 ------------------------------------------------------ s^11 + 5 s^9 + 9 s^7 + 2 s^6 + 12 s^5 + 4 s^4 + 12 s^3 ，T=0.1s num=[1 0 0.568]; den=conv([1,-1],[1 -0.2 0.99]); H=tf(num,den,Ts,0.1) Transfer function: z^2 + 0.568 ----------------------------- z^3 - 1.2 z^2 + 1.19 z - 0.99 Sampling time: 0.1 请将下面的零极点模型输入到matlab环境。请求出上述模型的零极点，并绘制其位置。 P=[0;0;-5;-6;-j;j];Z=[-1-j;-1+j]; G=zpk(Z,P,8) Zero/pole/gain: 8 (s^2 + 2s + 2) -------------------------- s^2 (s+5) (s+6) (s^2 + 1) pzmap(G) ，T=0.05s pzmap(H) 设描述系统时域行为的高阶微分方程为，试建立其状态空间模型并将其输入到matlab环境。 y=dsolve(D3y+13*D2y+4*Dy+5*y=2*u,t) y=2/5*u+C1*exp(1/591576*(-98596*(16244+12*112461^(1/2))^(1/3)-4061*(16244+12*112461^(1/2))^(2/3)+3*(16244+12*112461^(1/2))^(2/3)*112461^(1/2)-2563496)*t)+C2*exp(-1/1183152*(-98596*(16244+12*37487^(1/2)*3^(1/2))^(1/3)+3*(16244+12*37487^(1/2)*3^(1/2))^(2/3)*37487^(1/2)*3^(1/2)-4061*(16244+12*37487^(1/2)*3^(1/2))^(2/3)+5126992-9*i*(16244+12*37487^(1/2)*3^(1/2))^(2/3)*37487^(1/2)+4061*i*3^(1/2)*(16244+12*37487^(1/2)*3^(1/2))^(2/3)-98596*i*3^(1/2)*(16244+12*37487^(1/2)*3^(1/2))^(1/3))*t)+C3*exp(-1/1183152*(-98596*(16244+12*37487^(1/2)*3^(1/2))^(1/3)+3*(16244+12*37487^(1/2)*3^(1/2))^(2/3)*37487^(1/2)*3^(1/2)-4061*(16244+12*37487^(1/2)*3^(1/2))^(2/3)+5126992+9*i*(16244+12*37487^(1/2)*3^(1/2))^(2/3)*37487^(1/2)-4061*i*3^(1/2)*(16244+12*37487^(1/2)*3^(1/2))^(2/3)+98596*i*3^(1/2)*(16244+12*37487^(1/2)*3^(1/2))^(1/3))*t) 上机作业： 请分析下面传递函数模型的稳定性。 s=tf(s); G=1/(s^3+2*s^2+s+2); G Transfer function: 1 ------------------- s^3 + 2 s^2 + s + 2 pzmap(G) 由于系统的两个极点位于S平面虚轴上，所以该系统是临界稳定状态。 s=tf(s); G=(3*s+1)/((s^2)*(300*s^2+600*s+50)+3*s+1) Transfer function: 3 s + 1 ------------------------------------ 300 s^4 + 600 s^3 + 50 s^2 + 3 s + 1 pzmap(G) 由于系统有两个极点位于S平面虚轴左侧，所以该系统是不