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