资源描述：
《matlab与控制系统仿真部分习题答案》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。
1、【4.2】程序：num=[5,0];den=conv([1,1],conv([1,2],[1,3]));[numc,denc]=cloop(num,den);[z,p,k]=tf2zp(numc,denc);[A,B,C,D]=tf2ss(numc,denc);g_zp=zpk(z,p,k)g_tf=tf(numc,denc)g_ss=ss(A,B,C,D)运行结果：Zero/pole/gain:5s----------------------------------(s+0.4432)(s^2+5.557s+13.54)Transferfunction
2、:5s----------------------s^3+6s^2+16s+6a=x1x2x3x1-6-16-6x2100x3010b=u1x11x20x30c=x1x2x3y1050d=u1y10【4.3】程序：A=[000-1;100-2;010-3;001-4];B=[0;0;0;1];C=[1000];g_ss=ss(A,B,C,D)[num,den]=ss2tf(A,B,C,D);g_tf=tf(num,den)[z,p,k]=ss2zp(A,B,C,D);g_zpk=zpk(z,p,k)运行结果：a=x1x2x3x4x1000-1x2100-
3、2x3010-3x4001-4b=u1x10x20x30x41c=x1x2x3x4y11000d=u1y10Continuous-timemodel.Transferfunction:-3.109e-015s^3-s^2-3.331e-015s-4.441e-016-------------------------------------------------s^4+4s^3+3s^2+2s+1Zero/pole/gain:-s^2----------------------------------------------(s+0.6724)(s+3.
4、234)(s^2+0.0936s+0.4599)【5.1】（1）程序num=[0,10];den=conv([1,0],[1,7,17]);[numc,denc]=cloop(num,den,-1);G=tf(numc,denc)[y,t]=step(G);plot(t,y,'b-')C=dcgain(G);n=1;whiley(n)<0.1*Cn=n+1;endm=1;whiley(m)<0.9*Cm=m+1;endrisetime=t(m)-t(n)[Y,k]=max(y);percentovershoot=100*(Y-C)/Ci=length(t
5、);while(y(i)>0.98*C)&(y(i)<1.02*C)i=i-1;endsettlingtime=t(i)运行结果：Transferfunction:10-----------------------s^3+7s^2+17s+10risetime=2.7312percentovershoot=-0.4399settlingtime=5.1372图：（1）程序k=[10,100,1000];t=linspace(1,20,200);num=1;den=conv([1,0],[1,7,17]);forj=1:3;s1=tf(num*k(j),d
6、en);sys=feedback(s1,1)y(:,j)=step(sys,t);endplot(t,y(:,1),'r',t,y(:,2),'b',t,y(:,3),'g')gtext('k=10');gtext('k=100');gtext('k=1000')运行结果：Transferfunction:10-----------------------s^3+7s^2+17s+10Transferfunction:100------------------------s^3+7s^2+17s+100Transferfunction:1000-----
7、--------------------s^3+7s^2+17s+1000图：图：【6.1】程序：（1）num1=[1,1];den1=conv([1,0,0],conv([1,2],[1,4]));sys1=tf(num1,den1)rlocus(sys1)运行结果：（2）num2=[1,1];den2=conv([1,0],conv([1,-1],[1,4,16]));sys2=tf(num2,den2)rlocus(sys2)运行结果：（3）num3=[1,8];den3=conv([1,0,0],conv([1,3],conv([1,5],con
8、v([1,7],[1,15]))));sys3=tf(num3,den3)rlo
