紅河學(xué)院工學(xué)院實(shí)驗(yàn)報(bào)告單《機(jī)械工程控制基礎(chǔ)》Matlab仿真實(shí)驗(yàn)報(bào)告單課程名稱：《機(jī)械工程控制基礎(chǔ)》實(shí)驗(yàn)姓名：日期：2013-12-6成績：年級(jí)專業(yè)：2010級(jí)機(jī)械工程學(xué)號(hào)：201101030415實(shí)驗(yàn)場地：任美福樓222實(shí)驗(yàn)二：一階、二階線性系統(tǒng)時(shí)間響應(yīng)的Matlab仿真實(shí)驗(yàn)一、實(shí)驗(yàn)內(nèi)容。1、一階線性系統(tǒng)，當(dāng)T分別取0.1，2以及5時(shí)，系統(tǒng)單位脈沖響應(yīng)、單位階躍響應(yīng)的Matlab仿真。2、二階線性系統(tǒng)單位脈沖響應(yīng)、單位階躍響應(yīng)、單位正弦輸入響應(yīng)的Matlab仿真。二、實(shí)驗(yàn)?zāi)康摹?、熟悉Matlab中，輸出圖形的操作2、掌握Matlab中一階、二階線性系統(tǒng)的表達(dá)3、常握Matlab中一階、二階線性系統(tǒng)的典型輸入信號(hào)的時(shí)間響應(yīng)圖形的求法。三、相關(guān)Matlab仿真程序與仿真圖形1、程序G1=tf(1,[0.1,1])G2=tf(1,[2,1])G3=tf(1,[5,1])t=0:0.1:20;[y1,t]=step(G1,t);[y2,t]=step(G2,t);[y3,t]=step(G3,t);plot(t,y1,'r')holdonplot(t,y2,'b')holdonplot(t,y3,'g')xlabel('時(shí)間')ylabel('輸出')title('三個(gè)一階系統(tǒng)的單位階躍響應(yīng)')legend('y1T=0.1','y12T=2','y3T=5')axis([0,15,0,10])仿真圖形程序G1=tf(1,[0.1,1])G2=tf(1,[2,1])G3=tf(1,[5,1])t=0:0.1:15;[y1,t]=impulse(G1,t);[y2,t]=impulse(G2,t);[y3,t]=impulse(G3,t);plot(t,y1,'r')holdonplot(t,y2,'b')holdonplot(t,y3,'g')xlabel('時(shí)間')ylabel('輸出')title('三個(gè)一階系統(tǒng)的單位脈沖響應(yīng)')legend('y1T=0.1','y12T=2','y3T=5')axis([0,10,0,5])仿真圖形2、程序clearG=tf([16],[4,4,16]);t=0:0.1:20;figure(1);impulse(G,t)xlabel('時(shí)間')ylabel('輸出')title('二階系統(tǒng)單位脈沖響應(yīng)')figure(2);step(G,t)xlabel('時(shí)間')ylabel('輸出')title('二階系統(tǒng)單位階躍響應(yīng)')figure(3);r=sin(t);lsim(G,r,t)xlabel('時(shí)間')ylabel('輸出')tit