【正文】 。然后在Matlab軟件平臺下，利用函數(shù)wavread對語音信號進行采樣，記住采樣頻率和采樣點數(shù)。q = abs(linspace(1,1,n))。窗函數(shù)頻譜中的主瓣寬度較窄，獲得的過渡帶較陡；旁瓣衰減越大，阻帶的衰減越快。為了從hd(n)得到一個FIR濾波器必須在hd(n)兩邊將它截斷。(3)實現(xiàn)：上面一步的結(jié)果是一個濾波器的表述，它可能是一個差分方程的形式，或者是某一系統(tǒng)函數(shù)H(z)，或者是某一脈沖響應(yīng)h(n)。[3] FIR濾波器的設(shè)計方法有多種，如窗口設(shè)計法，頻率采樣法，最優(yōu)等波紋設(shè)計法等。MATLAB的基本數(shù)據(jù)單位是矩陣，它的指令表達式與數(shù)學(xué)、工程中常用的形式十分相似，故用MATLAB來解算問題要比用其他語言完成相同的事情簡捷得多，并且MATLAB也吸收了像Maple等軟件的優(yōu)點，使MATLAB成為一個強大的數(shù)學(xué)軟件。(2)設(shè)計完后應(yīng)檢查其頻率響應(yīng)曲線是否滿足指標(biāo)。隨著現(xiàn)代通信的數(shù)字化，數(shù)字濾波器變得更加重要。開發(fā)工具是M語言編程。數(shù)字濾波器是一種用來過濾時間離散信號的數(shù)字系統(tǒng)，它是通過對抽樣數(shù)據(jù)進行數(shù)字處理來達到頻域濾波的目的。[1] 課程設(shè)計目的 數(shù)字信號處理課程設(shè)計是數(shù)字信號處理課程的重要實踐性環(huán)節(jié)，是通信工程一次較全面的工程師能力訓(xùn)練。它將數(shù)值分析、矩陣計算、科學(xué)數(shù)據(jù)可視化以及非線性動態(tài)系統(tǒng)的建模和仿真等諸多強大功能集成在一個易于使用的視窗環(huán)境中，為科學(xué)研究、工程設(shè)計以及必須進行有效數(shù)值計算的眾多科學(xué)領(lǐng)域提供了一種全面的解決方案，并在很大程度上擺脫了傳統(tǒng)非交互式程序設(shè)計語言的編輯模式，代表了當(dāng)今國際科學(xué)計算軟件的先進水平。FIR濾波器相對于IIR濾波器具有易實現(xiàn)線性相位，系統(tǒng)總是穩(wěn)定的，允許設(shè)計多通帶（或多阻帶）濾波器的優(yōu)點，F(xiàn)IR濾波器主要有4種結(jié)構(gòu)，直接型，級聯(lián)型，線性相位型，頻率采樣型，與IIR結(jié)構(gòu)相比，結(jié)構(gòu)是相對簡單的。這些技術(shù)要求是由用途決定的。一理想帶寬為wcπ的LPF由下式給出為 1ejw , |w|≤wc Hd(ejw)= (22) 0 , wc|w|≤ π 其中wc 也稱為截止頻率，α稱為樣本延遲（注意，根據(jù)DTFT性質(zhì)，ejαw意味者在正n方向的位移或延遲）。 hd(n)= (25) (2)由性能指標(biāo)確定窗函數(shù)w(n)和窗口長度N。[6]Bohman窗對應(yīng)的Bohmanwin函數(shù)定義源程序如下：function w = bohmanwin(n)error(nargchk(1,1,nargin))。 圖21 Bohwin窗3 設(shè)計步驟 設(shè)計流程圖語音信號濾波去噪——使用BOHMAN窗設(shè)計的頻率采樣型FIR濾波器的設(shè)計流程圖如圖31所示。圖33 聲音屬性圖34 聲音選定 對應(yīng)程序如下： [x,fs,bits]=wavread(39。 % 計算時間范圍，樣本數(shù)除以采樣頻率 x=x(:,1)39。 Y=Y(1:N/2)。由于BOHMAN窗的過渡帶寬定為。最后畫出該低通濾波器的頻譜特性圖。ws=fo/fs*2*pi。 % 調(diào)用自編函數(shù)計算濾波器的頻率特性 delta_w=2*pi/1000。對應(yīng)程序如下： sound(x,fs,bits)。由于坐標(biāo)設(shè)置不適當(dāng)，最初顯示的加噪前后語音信號的時域和頻域波形對比不是很明顯，對坐標(biāo)進行調(diào)整后，通過對比加噪前后的頻域圖，可以明顯觀察到在頻率為9000Hz處，有一個尖脈沖，也就是噪聲。程序中存在的問題：程序中有一些警告存在，多次修改調(diào)試程序還是沒有解決，最后尋求了老師的幫助，老師解釋到如果不影響整個程序的運行和濾波器的設(shè)計效果可以不多做考慮。數(shù)字信號處理的算法需要利用計算機或?qū)Ｓ锰幚碓O(shè)備如數(shù)字信號處理器（DSP）和專用集成電路（ASIC）等。讓我意識到將理論與實際結(jié)合起來的重要性，只有理論知識是遠遠不夠的，只有把所學(xué)的理論知識與實踐相結(jié)合起來，從理論中得出結(jié)論，才能將知識為我所用，從而運用到實際生活中，提高自己的實際動手能力和獨立思考的能力。普羅克斯．劉樹棠．?dāng)?shù)字信號處理(MATLAB版)．西安：西安交通大學(xué)出版社，2008[6] 趙紅怡．?dāng)?shù)字信號處理及其MATLAB實現(xiàn)． % 計算信號x的長度fn=9000。 Y=Y(1:N/2)。時間（單位： s）39。plot(f,magX)。語音信號幅度譜圖39。)。ylabel(39。 fp=4000。 % 調(diào)用自編函數(shù)計算理想低通濾波器的脈沖響應(yīng)w_bohmanwin=(bohmanwin(M))39。w/pi39。plot(w/pi,mag)。title(39。ylabel(39。xlabel(39。)%（5）用濾波器濾波器對信號進行濾波去噪處理。)。plot( f,X)。原始語音信號幅度譜39。幅度39。頻率f(單位: Hz)39。 subplot(3,2,5)。title(39。ylabel(39。 % 播放加噪語音信號sound (y_fil,fs,bits)。)。)。grid onxlabel(39。ylabel(39。title(39。 subplot(3,2,3)。頻率f(單位: Hz)39。幅度39。% 用設(shè)計好的濾波器對y進行濾波Y_fil=fft(y_fil)。)。)。)subplot(2,2,3)。xlabel(39。ylabel(39。 % 計算實際濾波器脈沖響應(yīng)[db,mag,pha,grd,w]=freqz_m(h,[1])。 % 低通濾波器設(shè)計指標(biāo)wp=fp/fs*2*pi。)。加入單頻干擾后的語音信號39。plot(t,y)。頻率（單位： Hz）39。ylabel(39。 % 計算頻譜的譜線間隔f=0:deltaf:fs/2deltaf。 % 計算時間范圍，樣本數(shù)除以采樣頻率x=x(:,1)39。% 程序作者：王羚% 最后修改日期：201335%（1）對原始語音信號加噪[x,fs,bits]=wavread(39。俗話說柴多火焰高，如果只憑自己的知識，我會浪費許多的時間和精力，可是滴水成海，在同學(xué)和老師的幫助下，我不僅把學(xué)到的知識運用的課程設(shè)計中，還體會到設(shè)計成功后帶來的成就感。數(shù)字信號處理的核心算法是離散傅里葉變換(DFT)，是DFT使信號在數(shù)字域