資源描述:
《Matlab數(shù)字信號處理實驗報告材料》由會員上傳分享,免費在線閱讀,更多相關內(nèi)容在工程資料-天天文庫。
1、實用數(shù)字信號處理實驗報告基礎實驗篇實驗一離散時間系統(tǒng)及離散卷積一、實驗原理利用Matlab軟件計算出系統(tǒng)函數(shù)的零極點分布、單位脈沖響應和系統(tǒng)頻率響應等的圖像并于筆算結果進行比較,找出異同。編譯合適程序能計算取值范圍不同的離散卷積。二、實驗目的(1)熟悉MATLAB軟件的使用方法。(2)熟悉系統(tǒng)函數(shù)的零極點分布、單位脈沖響應和系統(tǒng)頻率響應等概念。(3)利用MATLAB繪制系統(tǒng)函數(shù)的零極點分布圖、系統(tǒng)頻率響應和單位脈沖響應。三、實驗步驟(1)自編并調(diào)試實驗程序,并且,給實驗程序加注釋;(2)按照實驗內(nèi)容完成筆算結
2、果;(3)驗證計算程序的正確性,記錄實驗結果。(4)至少要求一個除參考實例以外的實驗結果,在實驗報告中,要描述清楚實驗結果對應的系統(tǒng),并對實驗結果進行解釋說明。四、實驗源程序及實驗結果實驗1-1運行結果xlabel('n');文檔實用ylabel('h(n)');figure(2)[z,p,g]=tf2zp(b,a);zplane(z,p)title('零極點');function[x,n]=chongji(n1,n2)n=[n1:n2];x=[n==0];functionshiyan1()a=[1,-1,0
3、.9];b=1;x=chongji(-20,120);n=-20:120;h=filter(b,a,x);figure(1)stem(n,h);title('沖擊響應');實驗1-2運行結果b=[0.0181,0.0543,0.0543,0.0181];a=[1.000,-1.76,文檔實用1.1829,-0.2781];w=pi*freqspace(500);H=freqz(b,a,w);MH=abs(H);AH=angle(H);subplot(2,1,1);plot(w/pi,MH);grid;axis
4、([0,1,0,1]);xlabel('w(pi)');ylabel('
5、H
6、');title('幅度、相位響應');subplot(2,1,2);plot(w/pi,AH);grid;xlabel('w(pi)');ylabel('angle(H)');實驗1-3運行結果文檔實用n=0:30;%輸入x(n)和沖激響應h(n)x=zeros(1,length(n));h=zeros(1,length(n));x([find((n>=0)&(n<=4))])=1;h([find((n>=0)&(n<=8))])
7、=0.5;figure(1)subplot(3,1,1);stem(n,x);axis([0,30,0,2]);title('輸入序列');xlabel('n');ylabel('x(n)');subplot(3,1,2);stem(n,h);axis([0,30,0,2]);title('沖激響應序列');xlabel('n');ylabel('h(n)');%輸出響應文檔實用y=conv(x,h);subplot(3,1,3);n=0:length(y)-1;stem(n,y);title('輸出響應')
8、;xlabel('n');ylabel('y(n)');實驗二離散傅立葉變換與快速傅立葉變換一、實驗原理對有限長序列使用離散Fouier變換(DFT)可以很好的反映序列的頻譜特性,而且易于用快速算法在計算機上實現(xiàn),當序列x(n)的長度為N時,它的DFT定義為反變換為??有限長序列的DFT是其Z變換在單位圓上的等距采樣,或者說是序列Fourier變換的等距采樣,因此可以用于序列的譜分析。??FFT是為了減少DFT運算次數(shù)的一種快速算法。它是對變換式進行一次次分解,使其成為若干較短序列的組合,從而減少運算量。常用
9、的FFT是以2文檔實用為基數(shù)的,其長度。它的效率高,程序簡單,使用非常方便,當要變換的序列長度不等于2的整數(shù)次方時,為了使用以2為基數(shù)的FFT,可以用末位補零的方法,使其長度延長至2的整數(shù)次方。??用FFT可以實現(xiàn)兩個序列的圓周卷積。在一定的條件下,可以使圓周卷積等于線性卷積。一般情況,設兩個序列的長度分別為N1和N2,要使圓周卷積等于線性卷積的充要條件是FFT的長度N≥N1+N2對于長度不足N的兩個序列,分別將他們補零延長到N。二、實驗目的1、加深理解離散傅立葉變換及快速傅立葉變換概念;2、學會應用FFT對
10、典型信號進行頻譜分析的方法;3、研究如何利用FFT程序分析確定性時間連續(xù)信號;4、熟悉應用FFT實現(xiàn)兩個序列的線性卷積的方法。三、實驗步驟1、調(diào)試實驗程序,并且,給參考程序加注釋;2、利用編制的計算卷積的計算程序,分別給出一下三組函數(shù)的卷積結果一、實驗源程序及結果實驗2-1運行結果文檔實用b=[0.0181,0.0543,0.0543,0.0181];a=[1.000,-1.76,1.1829,-0