【正文】 信號(hào)濾波處理對(duì)語音信號(hào)信號(hào)進(jìn)行濾波處理主要是濾掉加入的噪聲。wsu=fsu/fs*2*pi。As=70。當(dāng)所有的準(zhǔn)備工作完成后就可以調(diào)用自編的函數(shù)計(jì)算理想帶阻濾波器的脈沖響應(yīng)和用窗函數(shù)法計(jì)算實(shí)際的濾波器的脈沖響應(yīng)。而FIR濾波器的設(shè)計(jì)也有很多方法。Y=abs(fft(y))。 % 單頻噪聲頻率，此參數(shù)可改t=0:1/fs:(N1)/fs。采樣完后再語音信號(hào)中加入一個(gè)單頻噪聲，單頻的噪聲的頻率可以自己設(shè)置。布萊克曼窗函數(shù)的最大旁瓣之比主瓣值低57db，但是主瓣寬度是矩形窗函數(shù)的主瓣寬度的三倍。窗口函數(shù)對(duì)理想特性的影響：改變了理想頻響的邊沿特性，形成過渡帶，寬為 , 等于WR(ω)的主瓣寬度；過渡帶兩旁產(chǎn)生肩峰和余振（帶內(nèi)、帶外起伏），取決于 WR(ω)的旁瓣，旁瓣多，余振多；旁瓣相對(duì)值大，肩峰強(qiáng)，與 N無關(guān)；N增加,過渡帶寬減小,肩峰值不變。由h（n）求H (ejω) 是否在誤差容限之內(nèi)。下面介紹窗函數(shù)設(shè)計(jì)法的基本設(shè)計(jì)過程。做法是從信號(hào)中截取一個(gè)時(shí)間片段，然后用觀察的信號(hào)時(shí)間片段進(jìn)行周期延拓處理，得到虛擬的無限長(zhǎng)的信號(hào)，然后就可以對(duì)信號(hào)進(jìn)行傅里葉變換、相關(guān)分析等數(shù)學(xué)處理。有限長(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)）也包含有反饋的遞歸部分。在保持內(nèi)核不變的情況下，MATLAB可以針對(duì)不同的應(yīng)用學(xué)科推出相應(yīng)的工具箱，目前已經(jīng)推出了圖像處理工具箱、信號(hào)處理工具箱、小波工具箱、神經(jīng)網(wǎng)絡(luò)工具箱以及通信工具箱等多個(gè)學(xué)科的專用工具箱，極大的方便了不同學(xué)科的研究工作。MATLAB名稱是有兩個(gè)英文單詞Matrix和Laboratory的前三個(gè)字母組成。綜合運(yùn)用本課程的理論知識(shí)進(jìn)行頻譜分析以及濾波器設(shè)計(jì)，通過理論推導(dǎo)得出相應(yīng)結(jié)論，并利用MATLAB作為編程工具進(jìn)行計(jì)算機(jī)實(shí)現(xiàn)，從而復(fù)習(xí)鞏固了課堂所學(xué)的理論知識(shí)，提高了對(duì)所學(xué)知識(shí)的綜合應(yīng)用能力，并從實(shí)踐上初步實(shí)現(xiàn)了對(duì)數(shù)字信號(hào)的處理。在本課程設(shè)計(jì)中，是用MATLAB的集成環(huán)境完成一系列的設(shè)計(jì)。然后在該語言信號(hào)中加一個(gè)噪音，利用布萊克曼窗設(shè)計(jì)一個(gè)FIR濾波器，對(duì)該語音信號(hào)進(jìn)行慮噪處理，然后比較濾波前后的波形與頻譜。 課程設(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)中占有重要地位。（4）獨(dú)立完成課程設(shè)計(jì)并按要求編寫課程設(shè)計(jì)報(bào)告書。MATLAB軟件有很強(qiáng)的開放性和適應(yīng)性。FIR濾波器系數(shù)計(jì)算較為繁瑣，在設(shè)計(jì)時(shí)借助Matlab工具箱，選擇合適的窗函數(shù)，可以方便地計(jì)算濾波器系數(shù)，并分析其幅頻、相頻特性。不過，當(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)行分析。窗函數(shù)設(shè)計(jì)法就是截取hd(n)為有限長(zhǎng)的一段因果序列，并用合適的窗函數(shù)進(jìn)行加權(quán)做為FIRDF的單位脈沖響應(yīng)h(n)。 （4）、檢驗(yàn)濾波器的性能。邊界頻率因?yàn)榧哟暗挠绊懚灰卓刂?。加N可減少過渡帶。再調(diào)用函數(shù)sound此時(shí)可以聽見錄制的語音。 % 計(jì)算信號(hào)x的長(zhǎng)度fn=2100。下面是對(duì)未加和加入噪聲的頻譜分析的程序：X=abs(fft(x))。而本課程設(shè)計(jì)采用的是窗函數(shù)法設(shè)計(jì)FIR濾波器。當(dāng)所有的指標(biāo)都設(shè)置完后，可以用這些數(shù)字來計(jì)算上下邊帶的中心頻率和頻率間隔，并計(jì)算布萊克曼窗設(shè)計(jì)該濾波器所需要的階數(shù)和產(chǎn)生幾階的布萊克曼窗。Rp=1。 % 將Hz為單位的模擬頻率換算為rad為單位的數(shù)字頻率wsd=fsd/fs*2*pi。 % 用窗口法計(jì)算實(shí)際濾波器脈沖響應(yīng)[db,mag,pha,grd,w]=freqz_m(h_bs,1)。下面是對(duì)語音信號(hào)進(jìn)行濾波的程序：y_fil=filter(h_bs,1,y)。而第二圖與第四圖相比較在頻率f=2100時(shí)多了一個(gè)尖銳脈沖。有幾次因?yàn)榍懊娴脑肼曨l率設(shè)置為f=2100，而后面的帶阻濾波器設(shè)計(jì)指標(biāo)沒有與前面的頻率相匹配導(dǎo)致畫出來的圖怎么也不能濾掉噪音。5 結(jié)束語為期兩周的數(shù)字信號(hào)處理課程設(shè)計(jì)已經(jīng)結(jié)束了，但在這次設(shè)計(jì)中我學(xué)到了許多的東西。雖然在這次課程設(shè)計(jì)中遇到很多的困難，但通過自己查找有關(guān)資料以及老師和同學(xué)的幫助下都一一解決了，而且在與同學(xué)交流的過程中使同學(xué)之間的感情更進(jìn)一步。 % 按指定的采樣率和每樣本編碼位數(shù)回放N=length(x)。 % 對(duì)原始信號(hào)和加噪信號(hào)進(jìn)行fft變換X=X(1:N/2)。)。plot(f,X)。語音信號(hào)幅度譜39。幅度39。頻率（單位：Hz）39。fpd=1800。df=min((fsdfpd),(fpufsu))。 % 定義時(shí)間范圍w_black=blackman(M)。濾波器的db39。subplot(2,2,2)。ylabel(39。xlabel(39。title(39。)。)。plot(f,X)。X39。t39。加噪聲后的頻域圖39。subplot(3,2,5)。ylabel(39。)。)