1、實(shí)驗(yàn)一 MATLAB及仿真實(shí)驗(yàn)（控制系統(tǒng)的時(shí)域分析）一、實(shí)驗(yàn)?zāi)康膶W(xué)習(xí)利用MATLAB進(jìn)行控制系統(tǒng)時(shí)域分析，包括典型響應(yīng)、判斷系統(tǒng)穩(wěn)定性和分析系統(tǒng)的動(dòng)態(tài)特性；二、實(shí)驗(yàn)內(nèi)容(一) 穩(wěn)定性1 系統(tǒng)傳函為，試判斷其穩(wěn)定性num1=0 3 2 5 4 6;den1=1 3 4 2 7 2;sys1=tf(num1,den1);figure(1);hold ongm,pm,wcp,wcg=margin(sys1);margin(sys1);title(對(duì)數(shù)頻率特性圖);xlabel(頻率rad/sec);ylabel(Gain dB); 2 用Matlab求出的極點(diǎn)。a=0 0 1 2 2;b=1 7 3

2、 5 2;z,p,k=tf2zpk(a,b) ;（二）階躍響應(yīng)1. 二階系統(tǒng)1）鍵入程序，觀察并記錄單位階躍響應(yīng)曲線num1=10;den1=1 2 10;step(num1,den1);grid on;2）計(jì)算系統(tǒng)的閉環(huán)根、阻尼比、無阻尼振蕩頻率，并記錄wn=sqrt(10);%自然振蕩頻率zunibi=2/wn;%阻尼比syms s;S=solve(s2+2*s+10);%求閉環(huán)根3）修改參數(shù)，分別實(shí)現(xiàn)和的響應(yīng)曲線，并記錄n0=10;d0=1 2 10;step(n0,d0);%原響應(yīng)曲線hold on;n1=10;d1=1 6.32 10;step(n1,d1);n2=10;d2=1 1

3、2.64 10;step(n2,d2);4）修改參數(shù)，分別寫出程序?qū)崿F(xiàn)和的響應(yīng)曲線，并記錄n0=10;d0=1 2 10;step(n0,d0);%原響應(yīng)曲線hold on;n1=2.5;d1=1 1 2.5;step(n1,d1);n2=40;d2=1 4 40;step(n2,d2);2. 作出以下系統(tǒng)的階躍響應(yīng)，并分析結(jié)果（1）（2）（3） （4）n0=2 10;d0=1 2 10;step(n0,d0);hold on;n1=1 0.5 10;d1=1 2 10;step(n1,d1);hold on;n2=1 0.5 0;d2=1 2 10;step(n2,d2);hold on;n3=1 0;d3=1 2 10;step(n3,d3);3. 求該系統(tǒng)單位階躍響應(yīng)曲線，并在所得圖形上加網(wǎng)格線和標(biāo)題num0=25;den0=1 4 25;step(num0,den0);grid on;xlabel(X);ylabel(Y );title(單位階躍曲線);（三）系統(tǒng)動(dòng)態(tài)特性分析用Matlab求二階系統(tǒng)和的峰值時(shí)間，上升時(shí)間，調(diào)整時(shí)間，超調(diào)量。n0=120;d0=1 12 120;step(n0,d0);damp(d0);y,x,t=step(n0,d0);y,t; n1=0.01