【正文】 )。ylabel(39。加噪聲后的頻域圖39。X39。)。title(39。ylabel(39。濾波器的db39。df=min((fsdfpd),(fpufsu))。頻率（單位：Hz）39。語音信號(hào)幅度譜39。)。 % 按指定的采樣率和每樣本編碼位數(shù)回放N=length(x)。5 結(jié)束語為期兩周的數(shù)字信號(hào)處理課程設(shè)計(jì)已經(jīng)結(jié)束了，但在這次設(shè)計(jì)中我學(xué)到了許多的東西。而第二圖與第四圖相比較在頻率f=2100時(shí)多了一個(gè)尖銳脈沖。 % 用窗口法計(jì)算實(shí)際濾波器脈沖響應(yīng)[db,mag,pha,grd,w]=freqz_m(h_bs,1)。Rp=1。而本課程設(shè)計(jì)采用的是窗函數(shù)法設(shè)計(jì)FIR濾波器。 % 計(jì)算信號(hào)x的長(zhǎng)度fn=2100。加N可減少過渡帶。 （4）、檢驗(yàn)濾波器的性能。不過，當(dāng)運(yùn)用計(jì)算機(jī)實(shí)現(xiàn)工程測(cè)試信號(hào)處理時(shí)，不可能對(duì)無限長(zhǎng)的信號(hào)進(jìn)行測(cè)量和運(yùn)算，而是取其有限的時(shí)間片段進(jìn)行分析。MATLAB軟件有很強(qiáng)的開放性和適應(yīng)性。 課程設(shè)計(jì)目的數(shù)字信號(hào)處理課程設(shè)計(jì)是數(shù)字信號(hào)處理課程的重要實(shí)踐性環(huán)節(jié)，是學(xué)生在校期間一次較全面的工程師能力訓(xùn)練，在實(shí)現(xiàn)學(xué)生總體培養(yǎng)目標(biāo)中占有重要地位。在本課程設(shè)計(jì)中，是用MATLAB的集成環(huán)境完成一系列的設(shè)計(jì)。MATLAB名稱是有兩個(gè)英文單詞Matrix和Laboratory的前三個(gè)字母組成。有限長(zhǎng)單位沖激響應(yīng)（FIR）濾波器有以下特點(diǎn)：(1)、系統(tǒng)的單位沖激響應(yīng)h (n)在有限個(gè)n值處不為零；(2)、系統(tǒng)函數(shù)H (z)在|z|0處收斂，極點(diǎn)全部在z = 0處（因果系統(tǒng)）；(3)、 結(jié)構(gòu)上主要是非遞歸結(jié)構(gòu)，沒有輸出到輸入的反饋，但有些結(jié)構(gòu)中（例如頻率抽樣結(jié)構(gòu)）也包含有反饋的遞歸部分。下面介紹窗函數(shù)設(shè)計(jì)法的基本設(shè)計(jì)過程。窗口函數(shù)對(duì)理想特性的影響：改變了理想頻響的邊沿特性，形成過渡帶，寬為 , 等于WR(ω)的主瓣寬度；過渡帶兩旁產(chǎn)生肩峰和余振（帶內(nèi)、帶外起伏），取決于 WR(ω)的旁瓣，旁瓣多，余振多；旁瓣相對(duì)值大，肩峰強(qiáng)，與 N無關(guān)；N增加,過渡帶寬減小,肩峰值不變。采樣完后再語音信號(hào)中加入一個(gè)單頻噪聲，單頻的噪聲的頻率可以自己設(shè)置。Y=abs(fft(y))。當(dāng)所有的準(zhǔn)備工作完成后就可以調(diào)用自編的函數(shù)計(jì)算理想帶阻濾波器的脈沖響應(yīng)和用窗函數(shù)法計(jì)算實(shí)際的濾波器的脈沖響應(yīng)。wsu=fsu/fs*2*pi。 %用設(shè)計(jì)好的濾波器對(duì)y進(jìn)行濾波Y_fil=fft(y_fil)。后來慢慢的改變指標(biāo)使其在噪聲頻率的左右，這時(shí)才能濾掉噪聲得到原始的語音信號(hào)。這次設(shè)計(jì)不僅讓我學(xué)會(huì)如何獨(dú)立完成一項(xiàng)工作，而且提高了獨(dú)立解決問題的能力，為以后的課程設(shè)計(jì)打下良好的基礎(chǔ)。Y=Y(1:N/2)。xlabel(39。)。fsd=2050。 % 產(chǎn)生M階的布萊克曼窗hd_bs=ideal_lp(wcd,M)+ideal_lp(pi,M)ideal_lp(wcu,M)。plot(w/pi,mag)。w/pi39。y_fil=filter(h_bs,1,y)。title(39。)。plot(t,y_fil)。ylabel(39。t39。)。xlabel(39。subplot(3,2,2)。h_bs39。)。)。 % 計(jì)算布萊克曼窗設(shè)計(jì)該濾波器時(shí)需要的階數(shù)n=0:M1。)。ylabel(39。subplot(2,2,2)。Y=abs(fft(y))。其次就是一些函數(shù)的細(xì)節(jié)問題。總結(jié)起來大概有以下幾個(gè)方面：首先，最主要的是要把設(shè)計(jì)濾波器的參數(shù)設(shè)置正確才能濾除語音信號(hào)中的噪聲。當(dāng)將設(shè)計(jì)好的濾波器濾掉噪聲后我們也可以再一次調(diào)用函數(shù)sound，聽此時(shí)的聲音是否與原始語音信號(hào)基本一樣，若沒有單頻嘯叫聲則說明此次設(shè)計(jì)是成功的，否則應(yīng)重新設(shè)置指標(biāo)。dw=df/fs*2*pi。若不能按照這些來設(shè)計(jì)濾波器則不可能慮掉噪聲。將所有未加和加入噪聲的時(shí)域圖和頻域圖畫在同一個(gè)圖中，便于比較和分析。錄制好語音信號(hào)后，打開MATLAB軟件平臺(tái)，利用函數(shù)waveread對(duì)語音信號(hào)進(jìn)行采樣，記住采樣頻率的和采樣點(diǎn)數(shù)。缺點(diǎn)是當(dāng)Hd(ejω)較為復(fù)雜時(shí)，hd(n)就不容易由反付里葉變換求得。因此hd(n)是無限長(zhǎng)非因果序列，不能直接作為FIRDF的單位脈沖響應(yīng)。使用具有并行處理特性的FPGA來實(shí)現(xiàn)FIR濾波器，既有很強(qiáng)的實(shí)時(shí)性，又兼顧了靈活性，為數(shù)字信號(hào)處理提供了一種很好的解決方案。（3）處理結(jié)果和分析結(jié)論應(yīng)該一致，而且應(yīng)符合理論。第１４頁共14 頁音樂信號(hào)濾波去噪——使用布萊克曼窗設(shè)計(jì)的FIR濾波器摘 要 本課程設(shè)計(jì)主要是用麥克風(fēng)采集一段語音信號(hào)，繪制其波形并觀察其頻譜。XX本課程設(shè)計(jì)能使學(xué)生對(duì)通信工程領(lǐng)域各種技術(shù)的DSP實(shí)現(xiàn)的設(shè)計(jì)有較熟練的掌握。國(guó)內(nèi)已有越來越多的科研和技術(shù)人員認(rèn)識(shí)到MABLAB的強(qiáng)大作用，并在不同領(lǐng)域內(nèi)使用MATLAB來快速實(shí)現(xiàn)科研構(gòu)想和提高工作效率。無線長(zhǎng)的信號(hào)被截?cái)嘁院?，其頻譜發(fā)生了畸變，原來集中在f(0)處的能量被分散到兩個(gè)較寬的頻帶中去了(這種現(xiàn)象稱之為頻譜能量泄漏)。如果不滿足，則返回第（2）步。布萊克曼窗主瓣寬，旁瓣小，頻率識(shí)別精度最低，但幅值識(shí)別精度最高。 % 計(jì)算時(shí)間范圍，樣本數(shù)除以采樣頻率x=x