华科材院大二上工控matlab实验代码2华科材院大二上工控matlab实验代码2.doc

1-8(a)1 t=0:0.01:10; for sigma=[0.1,0.2,0.3,0.7,1,2,3] num=1; den=[1 sigma 1]; G=tf(num,den); step(G,t) hold on; end grid; hold off; a= ξ=0.1;text(3,1.7,a), a= ξ=0.3;text(3,1.4,a), a= ξ=0.5;text(3,1.15,a), a= ξ=0.7;text(3,1,a), a= ξ=1;text(3,0.8,a), a= ξ=2;text(4,0.6,a), a= ξ=3;text(3,0.3,a), title(阶跃响应) 1-8(a)2 num=[1]; t=0:0.1:20; den=[1 0.6 1]; G=tf(num,den); y=step(G,t);plot(t,y); grid; [ym,km]=max(y);line(t(km),ym,marker,.,markeredgecolor,r,markersize,20); ystr=[ymax=,sprintf(%1.6g,ym)];tstr=[tp=,sprintf(%1.4g,t(km))]; text(t(km)+1,ym-0.05,{ystr,tstr}); ttt=t(find(abs(y-1)0.05));ts=max(ttt); hold on;plot(ts,1.05,bo,markersize,10); hold off; text(ts+1.2,1.05,[ts=,num2str(ts)]); tt=t(find(y-10));tr=min(tt); hold on;plot(tr,1,bo,markersize,10); hold off; text(tr,0.9,[tr=,num2str(tr)]); 1-8(b)1 t=0:0.01:20; for x=0.1:0.1:1 num=[x*x] den=[1 x x*x]; G=tf(num,den); step(G,t) hold on; end grid; hold off; 1-8(b)2 num=[0.49]; t=0:0.1:20; den=[1 0.7 0.49]; G=tf(num,den); y=step(G,t);plot(t,y); grid; [ym,km]=max(y);line(t(km),ym,marker,.,markeredgecolor,r,markersize,20); ystr=[ymax=,sprintf(%1.6g,ym)];tstr=[tp=,sprintf(%1.4g,t(km))]; text(t(km)+1,ym-0.05,{ystr,tstr}); ttt=t(find(abs(y-1)0.05));ts=max(ttt); hold on;plot(ts,1.05,bo,markersize,10); hold off; text(ts+1.2,1.05,[ts=,num2str(ts)]); tt=t(find(y-10));tr=min(tt); hold on;plot(tr,1,bo,markersize,10); hold off; text(tr,0.9,[tr=,num2str(tr)]); 2-2(1) for i=0.1:0.2:0.9; num=[1]; den=[1 2*i 1]; bode(num,den); hold on; end; hold off; legend(tao=0.1,tao=0.3,tao=0.5,’tao=0.7’,’tao=0.9’); 2-2(2) for i=0.1:0.2:0.9; y=tf(i*i,[1 2*0.707*i i*i]); hold on; bode(g); hold off;end; grid; legend(w=0.1,w=0.3,w=0.5,w=0.7,w=0.9); 2-5 num=100*conv([1 5],[1 5]); den=conv([1 1],[1 1 9]); G=tf(num,den);[Gm,Pm,Wcg,Wcp]=margin(G); G_close=feedback(G,1); step(G_close); 运行结果：Gm =Inf Pm =85.4365 Wcg =NaN Wcp =100.3285 1