使用布萊克曼窗設(shè)計(jì)的fir濾波器正文(文件)
【正文】 利用函數(shù)filter對(duì)y進(jìn)行濾波，然后對(duì)y進(jìn)行傅里葉變換。當(dāng)將設(shè)計(jì)好的濾波器濾掉噪聲后我們也可以再一次調(diào)用函數(shù)sound，聽此時(shí)的聲音是否與原始語音信號(hào)基本一樣，若沒有單頻嘯叫聲則說明此次設(shè)計(jì)是成功的，否則應(yīng)重新設(shè)置指標(biāo)。 % 計(jì)算頻譜取前一半sound (x，fs，bits)。從圖中可以看出，第一圖和第三圖相比因?yàn)榧尤朐肼暤木壒仕缘谌龍Dy軸的幅度要比第一圖要大，但其形狀還是基本沒有改變。從圖中可以清楚的了解濾波器的幅頻和相頻特性?？偨Y(jié)起來大概有以下幾個(gè)方面：首先，最主要的是要把設(shè)計(jì)濾波器的參數(shù)設(shè)置正確才能濾除語音信號(hào)中的噪聲。其次，在對(duì)加入噪聲后的信號(hào)進(jìn)行濾波時(shí)沒有使用正確的形式也不能濾去噪聲。因?yàn)樵诰幊痰倪^程過有可能出現(xiàn)錯(cuò)誤，如果建一個(gè)文檔有助于程序出現(xiàn)錯(cuò)誤時(shí)可以在文檔中直接修改，這樣可以省很多的時(shí)間，又這個(gè)課程設(shè)計(jì)的程序多而繁雜，一不小心就有可能寫錯(cuò)，如果在工作環(huán)境下修改這樣利于將所有的程序復(fù)制在課程設(shè)計(jì)中，而且還要對(duì)復(fù)制后的程序進(jìn)行刪除。雖然這次課程設(shè)計(jì)已經(jīng)完成了，但是遇到的困難也是很多的。其次就是一些函數(shù)的細(xì)節(jié)問題。設(shè)計(jì)的過程中，我通過查閱大量有關(guān)資料，與同學(xué)交流經(jīng)驗(yàn)和自學(xué)，并向老師請(qǐng)教等方式，使我學(xué)到了不少的東西，雖然有許多的辛酸，但是看到自己課程設(shè)計(jì)完成后心中的那份激動(dòng)是無法用言語來形容的。 % 輸入?yún)?shù)為文件的全路徑和文件名，輸出的第一個(gè)參數(shù)是每個(gè)樣本的值，fs是生成該波形文件時(shí)的采樣率，bits是波形文件每樣本的編碼位數(shù)sound(x,fs,bits)。 % 計(jì)算時(shí)間范圍，樣本數(shù)除以采樣頻率x=x(:,1)39。Y=abs(fft(y))。 % 計(jì)算頻譜的譜線間隔f=0:deltaf:fs/2deltaf。時(shí)間（單位：s）39。)。subplot(2,2,2)。)。title(39。plot(t,y)。ylabel(39。 加入單頻干擾后的語音信號(hào)39。xlabel(39。幅度譜39。)。fpu=2000。fcu=(fpu+fsu)/2。dw=df/fs*2*pi。 % 計(jì)算布萊克曼窗設(shè)計(jì)該濾波器時(shí)需要的階數(shù)n=0:M1。.*hd_bs。title(39。w/pi39。)。濾波器的幅頻特性39。)。subplot(2,2,3)。)。ylabel(39。plot(h_bs)。xlabel(39。h_bs39。 %對(duì)y進(jìn)行傅里葉變化Y_fil=Y_fil(1:length(Y_fil)/2)。未加噪聲的時(shí)域圖39。)。subplot(3,2,2)。)。ylabel(39。plot(t,y)。xlabel(39。y39。title(39。f39。)。加噪聲后的頻域圖39。)。title(39。t39。)。濾去噪聲后的頻域圖39。)。Y_fil39。xlabel(39。plot(f,Y_fil)。ylabel(39。)。Y39。xlabel(39。plot(t,y_fil)。ylabel(39。)。subplot(3,2,4)。)。加噪聲后的時(shí)域圖39。)。f39。title(39。x39。xlabel(39。plot(t,x)。y_fil=filter(h_bs,1,y)。)。濾波器的脈沖響應(yīng) 39。)。w/pi39。title(39。mag39。xlabel(39。plot(w/pi,mag)。ylabel(39。)。 % 調(diào)用自編函數(shù)計(jì)算濾波器的頻率特性subplot(2,2,1)。 % 產(chǎn)生M階的布萊克曼窗hd_bs=ideal_lp(wcd,M)+ideal_lp(pi,M)ideal_lp(wcu,M)。wsu=fsu/fs*2*pi。 % 計(jì)算上下邊帶中心頻率，和頻率間隔wcd=fcd/fs*2*pi。As=70。fsd=2050。title(39。)。subplot(2,2,4)。)。時(shí)間（單位：s）39。)。幅度譜39。xlabel(39。原始語音信號(hào)39。ylabel(39。plot(t,x)。Y=Y(1:N/2)。sound (y,fs,bits)。 % 計(jì)算信號(hào)x的長度fn=2100。e:39。這次設(shè)計(jì)不僅讓我學(xué)會(huì)如何獨(dú)立完成一項(xiàng)工作，而且提高了獨(dú)立解決問題的能力，為以后的課程設(shè)計(jì)打下良好的基礎(chǔ)。在設(shè)置過程中有很多次因?yàn)樵O(shè)置的參數(shù)不合適而導(dǎo)致設(shè)計(jì)的濾波器不能慮出單頻噪聲信號(hào)。通過這次的設(shè)計(jì)，不僅加深了我對(duì)課本基礎(chǔ)理論知識(shí)的理解，而且增強(qiáng)了我的實(shí)踐能力，同時(shí)更加認(rèn)識(shí)到理論知識(shí)和實(shí)踐結(jié)合的重要性。再次，在利用自編函數(shù)計(jì)算理想帶阻濾波器的脈沖響應(yīng)，如果在MATLAB軟件中的Work下沒有定義，那么不能調(diào)用自編的函數(shù)，否則將會(huì)報(bào)錯(cuò)。后來慢慢的改變指標(biāo)使其在噪聲頻率的左右，這時(shí)才能濾掉噪聲得到原始的語音信號(hào)。所以用布萊克曼窗設(shè)計(jì)的濾波器是符合要求的，也就是說該課程設(shè)計(jì)是成功的。說明原始語音信號(hào)加入噪聲是成功的。 結(jié)果分析 要確定本課程設(shè)計(jì)是否成功就得看原始信號(hào)的頻域圖和時(shí)域圖與經(jīng)過濾波器后的語音信號(hào)的頻域圖和時(shí)域圖是否一樣，若一樣則表示該設(shè)計(jì)是成功的，否則是不成功的。 %用設(shè)計(jì)好的濾波器對(duì)y進(jìn)行濾波Y_fil=fft(y_fil)。最后在同一個(gè)圖中畫出原始信號(hào)的
電大資料相關(guān)推薦
鄂ICP備17016276號(hào)-1