【正文】 % 播放加噪語音信號sound (y_fil,fs,bits)。 %（7）回放原始語音信號，加噪前后語音信號sound(x,fs,bits)。濾波后語音信號幅度譜39。)。ylabel(39。頻率f(單位: Hz)39。plot(abs(f),abs(Y_fil))。)。title(39。幅度39。)。grid onxlabel(39。 subplot(3,2,5)。加干擾語音信號幅度譜39。)。ylabel(39。頻率f(單位: Hz)39。plot(f,Y)。)。title(39。幅度39。)。 grid onxlabel(39。 subplot(3,2,3)。原始語音信號幅度譜39。)。ylabel(39。頻率f(單位: Hz)39。plot( f,X)。)。title(39。幅度39。)。 grid onxlabel(39。 % 計算頻譜取前一半%（6）畫出原始信號，加噪信號，濾波去噪信號的時域波形和頻譜subplot(3,2,1)。% 用設計好的濾波器對y進行濾波Y_fil=fft(y_fil)。)%（5）用濾波器濾波器對信號進行濾波去噪處理。title(39。h(n)39。)。xlabel(39。plot(n,hd)。濾波器相位響應圖39。)。ylabel(39。w/pi39。grid on。)subplot(2,2,3)。title(39。幅度mag39。)。xlabel(39。plot(w/pi,mag)。濾波器幅度響應圖39。)。ylabel(39。w/pi39。grid on。Rp=(min(db(1:1:wp/delta_w+1))) % 真正的通帶波紋As=round(max(db(ws/delta_w+1:1:501))) % 真正的阻帶衰減%（4）畫出濾波器的頻譜特性圖subplot(2,2,1)。 % 計算實際濾波器脈沖響應[db,mag,pha,grd,w]=freqz_m(h,[1])。 % 調(diào)用自編函數(shù)計算理想低通濾波器的脈沖響應w_bohmanwin=(bohmanwin(M))39。 % 定義時間范圍wc=(ws+wp)/2。tr_width=wswp。 % 低通濾波器設計指標wp=fp/fs*2*pi。 fp=4000。As=45。加入干擾后的語音信號幅度譜圖39。)。ylabel(39。頻率（單位： Hz）39。plot(f,magY)。加入單頻干擾后的語音信號39。)。ylabel(39。時間（單位： s）39。plot(t,y)。語音信號幅度譜圖39。)。ylabel(39。頻率（單位： Hz）39。plot(f,magX)。原始語音信號39。)。ylabel(39。時間（單位： s）39。plot(t,x)。magY=abs(Y)。 % 計算頻譜的譜線間隔f=0:deltaf:fs/2deltaf。 Y=Y(1:N/2)。 Y=abs(fft(y))。 % 加噪聲sound(y,fs,bits)。 % 計算時間范圍，樣本數(shù)除以采樣頻率x=x(:,1)39。 % 計算信號x的長度fn=9000。sound(x,fs,bits)。)。% 程序作者：王羚% 最后修改日期：201335%（1）對原始語音信號加噪[x,fs,bits]=wavread(39。普羅克斯．劉樹棠．數(shù)字信號處理(MATLAB版)．西安：西安交通大學出版社，2008[6] 趙紅怡．數(shù)字信號處理及其MATLAB實現(xiàn)．英格爾，約翰參考文獻[1] 高西全，丁玉美，闊永紅．數(shù)字信號處理——原理、實現(xiàn)及應用．北京：電子工業(yè)出版社，2009[2] 張圣勤．．北京：機械工程出版社，2008[3] 綠波杰能．FIR濾波器．百度文庫：，2013227[4] 張立材, 王民．數(shù)字信號處理．北京：人民郵電出版社，2008[5] 維納俗話說柴多火焰高，如果只憑自己的知識，我會浪費許多的時間和精力，可是滴水成海，在同學和老師的幫助下，我不僅把學到的知識運用的課程設計中，還體會到設計成功后帶來的成就感。讓我意識到將理論與實際結(jié)合起來的重要性，只有理論知識是遠遠不夠的，只有把所學的理論知識與實踐相結(jié)合起來，從理論中得出結(jié)論，才能將知識為我所用，從而運用到實際生活中，提高自己的實際動手能力和獨立思考的能力。當然，在設計過程中或多或少的遇到了一些問題，但在老師和同學的幫助下找出了原因所在。數(shù)字信號處理課程設計也讓我受益匪淺，通過此次課程設計，讓我更深入的了解了語音信號濾波去噪，窗口設計法設置FIR濾波器的相關知識，注意到很多以前沒有注意到的細節(jié)，學到了很多在書本上所沒有學到過的知識。數(shù)字信號處理的核心算法是離散傅里葉變換(DFT)，是DFT使信號在數(shù)字域和頻域都實現(xiàn)了離散化，從而可以用通用計算機處理離散信號。數(shù)字信號處理的算法需要利用計算機或?qū)Ｓ锰幚碓O備如數(shù)字信號處理器（DSP）和專用集成電路（ASIC）等。因此在進行數(shù)字信號處理之前需要將信號從模擬域轉(zhuǎn)換到數(shù)字域，這通常通過模數(shù)轉(zhuǎn)換器實現(xiàn)。數(shù)字信號處理是通信、電子、信息領域中最重要的專業(yè)課之一，數(shù)字信號處理對于我們通信工程的學生來說是一門很重要的科目，不僅僅因為它涵蓋了通信工程的許多理論知識，是考研的科目，而且其實踐性很強。5 結(jié)束語兩周的課程設計說長不長說短不短，如果好好利用，我們不僅可以鞏固學過的知識，也能拓寬自己的視野，提升自己的能力。程序中存在的問題：程序中有一些警告存在，多次修改調(diào)試程序還是沒有解決，最后尋求了老師的幫助，老師解釋到如果不影響整個程序的運行和濾波器的設計效果可以不多做考慮。為了盡量減少BOHMAM窗的階數(shù)，減小M的值，低通濾波器的過渡帶寬參數(shù)設置的過大，加噪語音信號通過低通濾波器后不能將噪聲完全濾除，經(jīng)過多次調(diào)試后將參數(shù)設置為了合適的值，能完全濾除噪聲。老師建議增大濾波器的過渡帶寬，或者對錄制的語音信號的屬性進行適當?shù)恼{(diào)整。但是運行程序后得出的真正的通帶波紋Rp和真正的阻帶衰減As是總是不達標，濾波器不符合設計要求，老師建議將過渡帶寬增大，經(jīng)過多次調(diào)試后，最終將過渡帶寬定為，設計的低通濾波器也達到了指標。由于坐標設置不適當，最初顯示的加噪前后語音信號的時域和頻域波形對比不是很明顯，對坐標進行調(diào)整后，通過對比加噪前后的頻域圖，可以明顯觀察到在頻率為9000Hz處，有一個尖脈沖，也就是噪聲。由于數(shù)字信號處理是上個學期學習的課程，所以理論知識對我來說已經(jīng)有點陌生了，最開始拿到課程設計題目時不知道從哪下手，在復習并掌握了FIR濾波器，窗口設計法相關的知識點后，我對本課程設計有了一定了理解。 繪制濾波器結(jié)構(gòu)圖本課程設計為頻率采樣型FIR濾波器，先求出濾波器的系數(shù)，對應語句如下： [C,B,A]=dir2fs(h) % 計算濾波器系數(shù)得出濾波器的系數(shù)為：C =