2016年北航數(shù)字信號(hào)處理上機(jī)實(shí)驗(yàn)二實(shí)驗(yàn)報(bào)告.doc

離散時(shí)間信號(hào)處理實(shí)驗(yàn)報(bào)告實(shí)驗(yàn)二FIR數(shù)字濾波器的設(shè)計(jì)與實(shí)現(xiàn)班 級(jí) 學(xué) 號(hào) 姓 名 日 期 實(shí)驗(yàn)二FIR數(shù)字濾波器的設(shè)計(jì)與實(shí)現(xiàn)實(shí)驗(yàn)二FIR數(shù)字濾波器的設(shè)計(jì)與實(shí)現(xiàn)濾波器設(shè)計(jì)與實(shí)現(xiàn)是數(shù)字信號(hào)處理中最基本的內(nèi)容之一。我們知道濾波器分IIR濾波器（無(wú)限沖激響應(yīng)濾波器）和FIR濾波器（有限沖激響應(yīng)濾波器）兩種。FIR濾波器的設(shè)計(jì)方法以直接逼近所需離散設(shè)計(jì)系統(tǒng)的頻率響應(yīng)為基礎(chǔ)。FIR濾波器可以很容易地獲得線性相移特性，不存在不穩(wěn)定的問(wèn)題，是實(shí)際系統(tǒng)中廣為采用的一種數(shù)字濾波器。FIR濾波器的設(shè)計(jì)，通常有窗函數(shù)設(shè)計(jì)法、頻率抽樣設(shè)計(jì)法和最佳逼近設(shè)計(jì)法。窗函數(shù)設(shè)計(jì)法比較簡(jiǎn)單，它的頻率特性是理想濾波器頻譜與窗的頻譜的卷積，因而，其頻率特性取決于窗的類型和長(zhǎng)度。頻率抽樣設(shè)計(jì)法比較直觀，但由于頻域的采樣會(huì)造成時(shí)域的混疊，從而濾波器疊性能不可能很高，為提高濾波器的性能，可以在過(guò)渡帶加上01之間的過(guò)渡點(diǎn)。本實(shí)驗(yàn)主要采用窗函數(shù)設(shè)計(jì)法。FIR濾波器在Matlab中的實(shí)現(xiàn)包括時(shí)域卷積和filter濾波等方法，通過(guò)本實(shí)驗(yàn)可以對(duì)兩種方法的實(shí)現(xiàn)過(guò)程和實(shí)現(xiàn)效果進(jìn)行分析比較。相關(guān)Matlab函數(shù)說(shuō)明1、有限沖激響應(yīng)數(shù)字濾波器設(shè)計(jì)函數(shù)fir1和fir2Matlab函數(shù)fir1和fir2可以用來(lái)設(shè)計(jì)加窗的有限沖激響應(yīng)數(shù)字濾波器。兩個(gè)函數(shù)均產(chǎn)生一個(gè)線性相位設(shè)計(jì)。函數(shù)fir1可用于設(shè)計(jì)常規(guī)的低通、高通、帶通和帶阻線性相位有限沖激響應(yīng)濾波器。對(duì)于抽樣頻率為2Hz的情況，命令b = fir1(N,Wn)在向量b中返回以的升冪排列的N階低通或帶通濾波器的沖激響應(yīng)系數(shù)。對(duì)于低通設(shè)計(jì)，歸一化截止頻率由標(biāo)量Wn給定，它是在0和1之間的一個(gè)數(shù)。對(duì)于帶通設(shè)計(jì)，Wn是包含指定通帶邊界的一個(gè)雙元素向量Wn1,Wn2，其中0Wn1Wn21。命令b = fir1(N,Wn,high)其中N為一個(gè)偶數(shù)，用于設(shè)計(jì)高通濾波器。命令b = fir1(N,Wn,stop)其中Wn是一個(gè)雙元素向量，用于設(shè)計(jì)帶阻有限沖激響應(yīng)濾波器。若沒(méi)有指定，用海明窗作為默認(rèn)值。命令b = fir(N,Wn,taper)使用在向量taper中給定的長(zhǎng)度為N+1的窗系數(shù)。然而，窗系數(shù)必須事先用一個(gè)適當(dāng)?shù)腗ATLAB函數(shù)來(lái)產(chǎn)生，如blackman，hamming，hanning，chebwin或kaiser。用到的命令為taper = blackman(N) taper = hamming(N) taper = hanning(N) taper = chebwin(N) taper = kaiser(N,beta)函數(shù)fir2可用于設(shè)計(jì)帶有任意形狀幅度響應(yīng)的線性相位有限沖激響應(yīng)濾波器。在其基本形狀中，命令是b = fir2(N,fpts,mval)它在長(zhǎng)度為N+1的向量b中返回以的升冪排列的沖激響應(yīng)系數(shù)。fpts是給定的頻率點(diǎn)向量，在0到1范圍內(nèi)以遞增的順序排列，其中第一個(gè)頻率點(diǎn)是0，而最后的頻率點(diǎn)是1。和前面一樣，抽樣頻率被假定為2Hz。mval是給定頻率點(diǎn)的給定幅度值的向量，因此也必須和fpts的長(zhǎng)度一樣，默認(rèn)使用海明窗。為使用其他窗口，所用到的命令是b = fir2(N,fpts,mval,taper)其中向量taper包含給定的窗系數(shù)。2、FIR濾波器的實(shí)現(xiàn)(1) FIR濾波器的時(shí)域卷積實(shí)現(xiàn)濾波器是一種能夠讓某些頻率分量通過(guò)而拒絕其他頻率成分或?qū)δ承╊l率成分進(jìn)行修正的系統(tǒng)。對(duì)于一個(gè)線性時(shí)不變系統(tǒng)，其時(shí)域輸入和輸出關(guān)系可用卷積的方法描述，。對(duì)時(shí)域而言，均為數(shù)字信號(hào)，可在數(shù)字系統(tǒng)中實(shí)現(xiàn)。濾波器按照單位脈沖響應(yīng)在某一有限區(qū)間外是否為零分為有限沖擊響應(yīng)（FIR）濾波器和無(wú)限沖擊響應(yīng)（IIR）濾波器。若單位脈沖響應(yīng)在某一有限區(qū)間以外都是零，稱為有限沖擊響應(yīng)FIR系統(tǒng)；若單位脈沖響應(yīng)在某一有限區(qū)間外不是零，則稱為無(wú)限沖擊響應(yīng)IIR系統(tǒng)。在線性時(shí)不變系統(tǒng)的卷積描述中，對(duì)于FIR濾波器，為有限項(xiàng)，范圍從0到M，故。即通過(guò)卷積運(yùn)算，可實(shí)現(xiàn)對(duì)輸入信號(hào)的濾波處理。(2) FIR濾波器的filter濾波實(shí)現(xiàn)MATLAB中有濾波函數(shù)filter，用法為y=filter(b,a,x)。其中b和a為差分方程系數(shù)，x為輸入序列，y為濾波輸出序列。對(duì)于FIR濾波器，a設(shè)置為1即可。MATLAB實(shí)現(xiàn)主要包括兩個(gè)部分，即濾波器設(shè)計(jì)和濾波。代碼如下：%fir濾波器系數(shù)設(shè)計(jì)fl_kaiser = 750 1150;fl_mag = 1 0;fl_dev = 0.01 0.05;fl_n_kaiser,fl_wn,fl_beta,fl_ftype = kaiserord(fl_kaiser,fl_mag,fl_dev,fs);h = fir1(fl_n_kaiser,fl_wn,fl_ftype,kaiser(fl_n_kaiser+1,fl_beta);%濾波函數(shù)y = filter(h,1,x);3、wav格式音頻文件讀取函數(shù) y = wavread(file)y,fs,nBits = wavread(file)輸入?yún)?shù)：file為wav格式的音頻文件。輸出參數(shù)：y為存放音頻數(shù)據(jù)的矩陣；fs為該音頻文件的采樣頻率；nBits表示語(yǔ)音信號(hào)的每個(gè)采樣點(diǎn)所占用的位數(shù)。4、wav格式音頻文件寫入函數(shù)wavwrite(y,fs,nBits,wavfile)wavwrite(y,fs,wavfile)wavwrite(y,wavfile)輸入?yún)?shù)：y是存放要寫入文件的音頻數(shù)據(jù)，fs是寫入文件的音頻信號(hào)采樣頻率，缺省的情況下默認(rèn)fs為8000Hz。nBits是表示語(yǔ)音信號(hào)的每個(gè)取樣點(diǎn)所占用的位數(shù)，它必須是8,16,24或32。nBits缺省的情況下默認(rèn)為16。wavfile是生成的WAV文件的文件名。實(shí)驗(yàn)題目1、已知兩個(gè)信號(hào)序列：用conv函數(shù)求兩個(gè)序列的卷積和，并繪制三個(gè)序列的波形。（1）MATLAB 源代碼：clc;n1=0:20;n2=0:10;y1=0:30;x1n=0.8.n1;x2n=ones(1,11);y=conv(x1n,x2n);subplot(3,1,1)stem(n1,x1n);subplot(3,1,2)stem(n2,x2n);subplot(3,1,3)stem(y1,y);（2）實(shí)驗(yàn)結(jié)果：第一張圖為f1，第二張圖為f2，第三張圖為f1*f22、利用wavread函數(shù)（或者audioread函數(shù)）將加噪聲的音頻文件“第一場(chǎng)雪加噪聲.wav”讀入Matlab的工作空間，利用load函數(shù)從h40.mat文件中讀取40階離散時(shí)間帶通濾波器的脈沖響應(yīng)，用conv函數(shù)求兩個(gè)序列的卷積和，并利用wavwrite函數(shù)（或者audiowrite函數(shù)）將結(jié)果寫入音頻文件“第一場(chǎng)雪濾波40.wav”中。繪制兩個(gè)音頻信號(hào)的波形，并通過(guò)音頻播放器比較兩個(gè)音頻信號(hào)的區(qū)別。（1）MATLAB 源代碼：a,Fs,bits = wavread(第一場(chǎng)雪加噪聲);load(h40.mat);a=a;a1=a(1,1:end);a2=a(2,1:end);h1n=h40;Y1=conv(a1,h1n);subplot(6,1,1)plot(a1);xlabel(n1);ylabel(a1);subplot(6,1,2)plot(h1n);xlabel(n2);ylabel(h1n);subplot(6,1,3)plot(Y1);xlabel(y1);ylabel(a1*h1n);Y2=conv(a2,h1n);subplot(6,1,4)plot(a2);xlabel(n1);ylabel(a2);subplot(6,1,5)plot(h1n);xlabel(n2);ylabel(h1n);subplot(6,1,6)plot(Y2);xlabel(y2);ylabel(a2*h1n);ak=zeros(2,3109334);ak(1,1:end)=Y1;ak(2,1:end)=Y2;ak=ak;wavwrite(ak,Fs,bits,第一場(chǎng)雪濾波40.wav);（2）實(shí)驗(yàn)結(jié)果：說(shuō)明：n1為左聲道信號(hào)，n2為右聲道信號(hào)3、利用wavread函數(shù)（或者audioread函數(shù)）將加噪聲的音頻文件“第一場(chǎng)雪加噪聲.wav”讀入Matlab的工作空間，利用load函數(shù)從h392.mat文件中讀取392階離散時(shí)間帶通濾波器的脈沖響應(yīng)，用conv函數(shù)求兩個(gè)序列的卷積和，并利用wavwrite函數(shù)（或者audiowrite函數(shù)）將結(jié)果寫入音頻文件“第一場(chǎng)雪濾波392.wav”中。繪制兩個(gè)音頻信號(hào)的波形，并通過(guò)音頻播放器比較兩個(gè)音頻信號(hào)的區(qū)別。（1）MATLAB 源代碼：a,Fs,bits = wavread(第一場(chǎng)雪加噪聲);load(h392.mat);a=a;a1=a(1,1:end);a2=a(2,1:end);n1=1:3109295;n2=1:393;y1=1:(3109295+393-1);h1n=h392;Y3=conv(a1,h1n);subplot(6,1,1)plot(n1,a1);xlabel(n1);ylabel(a1);subplot(6,1,2)plot(n2,h1n);xlabel(n2);ylabel(h1n);subplot(6,1,3)plot(y1,Y3);xlabel(y1);ylabel(a1*h1n);Y4=conv(a2,h1n);subplot(6,1,4)plot(n1,a2);xlabel(n1);ylabel(a2);subplot(6,1,5)plot(n2,h1n);xlabel(n2);ylabel(h1n);subplot(6,1,6)plot(y1,Y4);xlabel(y2);ylabel(a2*h1n);ak=zeros(2,3109687);ak(1,1:end)=Y3;ak(2,1:end)=Y4;ak=ak;wavwrite