資源描述:
《mb機電系統(tǒng)動態(tài)仿真上機備用》由會員上傳分享,免費在線閱讀,更多相關(guān)內(nèi)容在行業(yè)資料-天天文庫。
1、MB機電系統(tǒng)動態(tài)仿真上機備用3-3p41業(yè)t=(0:pi/100:pi)';%有邊界的限定,p41業(yè)。機電系統(tǒng)動態(tài)仿真。y1=sin(t)*[1,-1];y2=sin(t).*sin(9*t);plot(t,y1,'b:',t,y2,'r')axis([0,pi,-1,1]);title('fontsize{14}fontname{斜體}調(diào)制波形y=sintsin9t及其包絡(luò)線');gridontext(2.5,0.6,'leftarrow包絡(luò)線sin(t)*[1,-1]')text(0.96,0.6,'leftarrow波形曲線:sintsin9t')x
2、label('fontsize{28}bft');ylabel('fontsize{28}y')p44;%圖像標示演示。t=0:0.01:pi;B=[1,3]';w=[2,5]';y=sin(w*t).*exp(-B*t);plot(t,y(1,:),'-.',t,y(2,:))legend('w=2,B=1','w=5,B=3')xlabel('fontsize{28}bft'),ylabel('fontsize{28}y')title('fontsize{14}fontname{斜體}y=sin(wt)*exp(-Bt)')text(1,sin
3、(2)*exp(-1),'fontsize{19}fontname{bf}y=sin(2t)*exp(-t)')text(0.5,sin(2.5)*exp(-1.5),'fontsize{14}fontname{it}y=sin(5t)*exp(-3t)')text(2.5,-0.15,'fontsize{14}fontname{黑體}備注:此圖網(wǎng)格線特別密集')gridminor;x=linspace(-2,2,25);y=linspace(-2,2,25);[xx,yy]=meshgrid(x,y);zz=xx.*exp(-xx.^2-yy.^2)
4、;subplot(121),mesh(xx,yy,zz);title('立體網(wǎng)狀圖');xlabel('X');ylabel('Y');zlabel('Z');subplot(122),surf(xx,yy,zz);title('立體曲面圖');xlabel('X');ylabel('Y');zlabel('Z');subplot(241),peaks;[x,y,z]=peaks;title('圖1:對peaks函數(shù)取點');subplot(242),meshc(x,y,z);axis([-inf,inf,-inf,inf,-inf,inf]);title('圖2
5、:網(wǎng)狀圖和等高線');subplot(243),contour3(peaks,25);axis([-inf,inf,-inf,inf,-inf,inf]);title('圖3:曲面在三維空間中的等高線');subplot(244),contour3(peaks,25);axis([-inf,inf,-inf,inf]);title('圖4:等高線在X-Y平面的投影');subplot(245),meshz(x,y,z);axis([-inf,inf,-inf,inf,-inf,inf]);title('圖5:曲面加上圍裙');subplot(246),waterf
6、all(x,y,z);axis([-inf,inf,-inf,inf,-inf,inf]);title('圖6:在X方向加上水流效果');subplot(247),waterfall(x',y',z');axis([-inf,inf,-inf,inf,-inf,inf]);title('圖7:在Y方向加上水流效果');set(gca,'ytick',[0.1.3]);表示在哪些地方顯示set(gca,'yticklabel',{'a°','b°','c°'});表示在這些地方顯示什么內(nèi)容第二次上機1.1w=1;y=zeros(200,1);i=1;num=w*w;
7、forzeta=0.1:0.1:1den=[1,2*zeta*w,w*w];t=[0:0.1:19.9];y(:,1)=step(num,den,t);i=i+1;endmesh(fliplr(y),[12030])zlabel('Y')ylabel('t')xlabel('zeta')grid課本118sys=tf([1,-1],[1,1.5]);step(sys,30)自創(chuàng)forzeta=0.1:0.1:1;den=[1,2*zeta,1]sys=tf([1],den)第二次上機1.2G1=tf([1,0.5],[1,0.1]);x=[1,0];y=[1,2
8、];z=[