信號與系統(tǒng)課程實驗報告實驗四實驗名稱教師姓名學生姓名一、實驗內容

信號的時域分析及Matlab實現(xiàn)實驗地點

系別5309學號

實驗日期 2011-06-20預習實驗原理。對實驗內容編寫程序(M文件)，上機運行。記錄并整理實驗數(shù)據(jù)。二、實驗目的Matlab分析系統(tǒng)時間響應的方法。Matlab分析連續(xù)、離散系統(tǒng)的沖激響應的方法。理解系統(tǒng)零、極點分布與系統(tǒng)穩(wěn)定性關系。三、涉及實驗的相關情況介紹（包含使用軟件或實驗設備等情況）計算機一臺(安裝MATLAB6.5版本或以上版本)四、實驗試做記錄（含程序、數(shù)據(jù)記錄及分析）設H(s) s(sp1)(sp2)設①p1=-2，p2=-30；②p1=-2，p2=3針對極點參數(shù)①②，畫出系統(tǒng)零、極點分布圖，判斷該系統(tǒng)穩(wěn)定性。程序：clearnum=[1];den=[1,32,60];[r,poles,k]=residue(num,den);p=roots(den);z=roots(num);subplot(2,2,1);plot(real(p),imag(p),'*');holdon;plot(real(z),imag(z')o,');gridonT=0:0.1:10;y1=impulse(num,den,T);subplot(2,2,2);plot(T,y1);grid;title('脈沖響應')結果：位于s左半平面，穩(wěn)定針對極點參數(shù)①②，繪出系統(tǒng)的脈沖響應曲線，并觀察t趨勢。程序：clearnum=[1];den=[1,-1,-6];[r,poles,k]=residue(num,den);p=roots(den);z=roots(num);subplot(2,2,1);plot(real(p),imag(p),'*');holdon;plot(real(z),imag(z')o,');gridonT=0:0.1:10;y1=impulse(num,den,T);subplot(2,2,2);plot(T,y1);grid;title('脈沖響應')結果：由系統(tǒng)零、極點分布圖可知，該系統(tǒng)有一極點位于s由系統(tǒng)零、極點分布圖可知，該系統(tǒng)有一極點位于s右半平面，故系統(tǒng)是不穩(wěn)定的。（3）（必做）程序：w=linspace(0,5,200);b=[10];a=[13260];H=freqs(b,a,w);subplot(2,1,1);plot(w,abs(H));set(gca,'xtick',[012345]);set(gca,'ytick',[00.40.7071]);xlabel('\omega');ylabel('|H|(j\omega)|');subplot(2,1,2);plot(w,angle(H));set(gca,'xtick',[012345]);xlabel('\omega');ylabel('\phi(\omega)');結果：利用filter命令求下面系統(tǒng)的沖激響應：y(n)-0.7y(n-1)-0.6y(n-2)+y(n-3)=x(n)+0.5x(n-1)（必做）程序：A=[1,-0.7,-0.6,1];B=[1,0.5];n=[-20:100];X=[n==0];y=filter(B,A,X);stem(n,y);title('沖擊響應');grid結果：五、實驗總結通過本次實驗，我學習了如何用五、實驗總結通過本次實驗，我學習了如何用Matlab軟件產生基本信號（連續(xù)/離散的正弦、