【正文】 語音信號加入單頻噪聲前后的時域和頻域波形圖 濾波器設計 在 Matlab 中，可以利用 矩形窗、三角窗、漢寧窗、漢明窗、布萊克曼窗、凱塞窗等 設計 FIR濾波器， 在本次課程設計中主要應用漢寧 窗設計出 FIR濾波器。 % 將 Hz為單位的模擬頻率換算為 rad 為單位的數(shù)字頻率 wsd=fsd/fs*2*pi。 其調用格式： sound (x， fs， bits)； 可以感覺濾波前后的聲音有變化。掌握了基于 DSP的應用開發(fā)，開闊了視野，增長了知識，是進入現(xiàn)代數(shù)字信號處理領域重要技能，乃至大規(guī)模集成電路的開發(fā)也是會用到的基礎，今后要予以重視并積極努力去學習。 Y=fft(y)。 grid on。加噪聲后的時域波形圖 39。)。 dw=df/fs*2*pi。w/pi39。)。 title(39。 title(39。f39。)。 title(39。f39。 title(39。)。t39。 title(39。Y_fil=Y_fil(1:length(Y_fil)/2)。)。w/pi39。 title(39。fcu=(fpu+fsu)/2。加噪聲后的頻域波形圖 39。 grid on。)。 % 單頻噪 聲頻率設為 1800Hz t=0:1/fs:(length(x)1)/fs。在編寫公式時，應該單擊“插入”，然后選擇“對象”，在選擇“公式 ”就可以編寫公式。 % 用設計好 的濾波器對 y進行濾波 Y_fil=abs(fft(y_fil))。 14 df=min((fsdfpd),(fpufsu))。 % 截取前半部分 deltaf=fs/2/length(X)。這就是說采樣頻率是衡量聲卡采集、記錄和還原聲音文件的質量標準。 常見的窗函數(shù)，可以分為以下主要類型， （ 1） 冪窗 采用時間變量某種冪次的函數(shù)，如矩形、三角形、梯形或其它時間（ t）的高次冪； （ 2） 三角函數(shù)窗 應用三角函數(shù)，即正弦或余弦函數(shù)等組合成復合函數(shù)，例如漢寧窗、海明窗等； （ 3） 指數(shù)窗 采用指數(shù)時間函數(shù) ， 例如高斯窗等 [5]。在新的版本中也加入了對 C， FORTRAN， C++ ， JAVA 的支持。 課程設計目的 數(shù)字信號處理（ Digital Signal Processing，簡稱 DSP）是一門涉及許多學科而又廣泛應用于許多領域的新興學科。 本次的課程設計主要采用的是第一類設計方法，是利用漢寧窗函數(shù)法設計FIR 濾波器對一段語音進行濾波去噪，通過這一過程，對濾波前后波形進行對比分析得到結論。 設 FIRDF 的單位脈沖響應 )(nh 的長度為 N ，則其頻率響應函數(shù)為 ??? ?? 10 )()( Nn njj enheH ?? （ 21） 一般將 )( ?jeH 表示成如下形式： )()()( ??? ? jgj eHeH ? （ 22） 式中， )(?gH 是 ? 的實函 數(shù)（可以去負值）。而 ??nhd 關于 ??n 點偶對稱，所 ? ? 2/1?? N? ，同時要求 ??nw 關于 ? ? 2/1?N 點偶對稱。 sound(x,fs,bits)。利用 Matlab 中的函數(shù) freqz 畫出各濾波器的頻率響應 ，首先利用數(shù)字信號處理里面學過的知識，根據(jù)自己選定的參數(shù)，用漢寧窗函數(shù)法設計 FIR 數(shù)字濾波器，得到數(shù)字濾波器的參數(shù) b， a。 wsu=fsu/fs*2*pi。 聲音變得不會有刺耳的聲音，而且比以前更加地平滑。 在此次課程中主要的要求是用麥克風采集一段語音信號，在 MATLAB 軟件平臺上，利用編程加入噪聲函數(shù)， 繪制 出前后的 波形并觀察其頻譜，給定相應技術指標，用 漢寧窗 設計一個滿足指標的 FIR 濾波器， 自己編寫出符合漢寧窗函數(shù)性質的函數(shù)， 對該語音信號進行濾波去噪處理，比較濾波前后 的波形和頻譜并進行分析，根據(jù)結果和學過的理論得出合理的結論 ，寫出符合要求的課程設計報告。 % 對原始信號和加噪信號進行 fft變換 X=X(1:length(X)/2)。 %加上標題和橫坐標名稱 subplot(2,2,2)。)。ylabel(39。 % 將 Hz 為單位的模擬頻率換算為 rad 為單位的數(shù)字頻率 wsd=fsd/fs*2*pi。)。 grid on。脈沖響應 39。加噪聲前的時域波形 39。)。 grid on。濾波后的時域波形 39。)。plot(f,Y_fil)。Y39。 xlabel(39。plot(f,X)。 % 用設計好的濾波器對 y進行濾波 Y_fil=abs(fft(y_fil))。pha39。 xlabel(39。plot(w/pi,db)。 % 帶阻濾波器設計指標 fcd=(fpd+fsd)/2。 title(39。)。t39。 % 計算信號 x的長度 fn=1800。然后自己查了一些相關書籍得知，繪制流程圖時 一般用橢圓表示 “ 開始 ” 與 “ 結束 ” ， 行動方案普通工作環(huán)節(jié)用矩形表示 ， 問題判斷或判定（審核 /審批 /評審）環(huán)節(jié)用菱形表示 ，箭頭代表工作流方向 ， 流程圖實例圖 ， 輸入輸出為平行四邊形 。 在同一張大圖里分別繪制原始信號 x，加噪信號 x’ ，濾波去噪信號 y的時域波形和頻譜，以便比較和分析； 主程序如下： y_fil=fftfilt(h_bs,y)。fcu=(fpu+fsu)/2。 Y=Y(1:length(Y)/2)。采樣頻率與聲音頻率之間有一 定的關系，根據(jù)奎斯特理論，只有采樣頻率高于聲音信號最高頻率的兩倍時，才能把數(shù)字信號表示的聲音還原成為原來的聲音。窗口設計法就是截取 ??nhd 為有限長 8 的一段因果序列，并用合適的窗口函數(shù)進行加權作為 FIRDF 的單位脈沖響應??nh 。 6 圖 11 MATLAB 的設計平臺 MATLAB 的基本數(shù)據(jù)單位是矩陣，它的指令表達式與數(shù)學、工程中常用的形式十分 相似，故用 MATLAB 來解算問題要比用 C， FORTRAN 等語言完成相同的事情簡捷得多，并且 MathWork也吸收了像 Maple 等軟件的優(yōu)點 ,使 MATLAB 成為一個強大的數(shù)學軟件。 關鍵詞 DSP 應用實例； 濾波去噪； FIR濾波器；漢寧窗； MATLAB 使用