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