1、數字信號處理課程設計報告數字信號處理課 程 設 計課程名稱數字信號處理題目名稱基于 MATLAB 的語音去噪處理專業班級13級通信工程本一學生姓名余悅婷，牛牧青，陳艷艷，白洋帆學 號 余悅婷 4704130098 牛牧青 4704130102 陳艷艷 4704130093 白洋帆 4704130101指導教師周雅 二一五年十二月二十七日引 言濾波器設計在數字信號處理中占有極其重要的地位，FIR數字濾波器和IIR濾波器是濾波器設計的重要組成部分。利用MATLAB信號處理工具箱可以快速有效地設計各種數字濾波器。課題基于MATLAB有噪音語音信號處理的設計與實現，綜合運用數字信號處理的理論知識對加噪

2、聲語音信號進行時域、頻域分析和濾波。通過理論推導得出相應結論，再利用 MATLAB 作為編程工具進行計算機實現。在設計實現的過程中，使用窗函數法來設計FIR數字濾波器，用巴特沃斯、切比雪夫和雙線性變法設計IIR數字濾波器，并利用MATLAB 作為輔助工具完成設計中的計算與圖形的繪制。通過對對所設計濾波器的仿真和頻率特性分析，可知利用MATLAB信號處理工具箱可以有效快捷地設計FIR和IIR數字濾波器，過程簡單方便，結果的各項性能指標均達到指定要求。關鍵詞 數字濾波器 MATLAB 窗函數法 巴特沃斯 切比雪夫 雙線性變換 目錄1 緒論12 課程設計內容23 課程設計的具體實現23.1 語音信號

3、的采集43.2 語音信號的時頻分析43.3 語音信號加噪與頻譜分析63.4 利用雙線性變換法設計低通濾波器83.5 用濾波器對加噪語音信號進行濾波93.6 分析濾波前后語音信號波形及頻譜的變化103.7回放語音信號103.8小結11結論22致謝23參考文獻24 1 緒論 數字信號處理是利用計算機或專用處理設備，以數值計算的方法對信號進行采集、抽樣、變換、綜合、估值與識別等加工處理，借以達到提取信息和便于應用的目的。它在語音、雷達、圖像、系統控制、通信、航空航天、生物醫學等眾多領域都獲得了極其廣泛的應用。具有靈活、精確、抗干擾強、設備尺寸小、造價低、速度快等優點。數字濾波器, 是數字信號處理中及

4、其重要的一部分。隨著信息時代和數字技術的發展，受到人們越來越多的重視。數字濾波器可以通過數值運算實現濾波，所以數字濾波器處理精度高、穩定、體積小、重量輕、靈活不存在阻抗匹配問題，可以實現模擬濾波器無法實現的特殊功能。數字濾波器種類很多，根據其實現的網絡結構或者其沖激響應函數的時域特性，可分為兩種，即有限沖激響應( FIR，Finite Impulse Response)濾波器和無限沖激響應( IIR，Infinite Impulse Response)濾波器。FIR濾波器結構上主要是非遞歸結構，沒有輸出到輸入的反饋，系統函數H (z)在處收斂，極點全部在z = 0處（因果系統），因而只能用較高

5、的階數達到高的選擇性。FIR數字濾波器的幅頻特性精度較之于IIR數字濾波器低，但是線性相位，就是不同頻率分量的信號經過fir濾波器后他們的時間差不變，這是很好的性質。FIR數字濾波器是有限的單位響應也有利于對數字信號的處理，便于編程，用于計算的時延也小，這對實時的信號處理很重要。 FIR濾波器因具有系統穩定，易實現相位控制，允許設計多通帶（或多阻帶）濾波器等優點收到人們的青睞。IIR濾波器采用遞歸型結構，即結構上帶有反饋環路。IIR濾波器運算結構通常由延時、乘以系數和相加等基本運算組成，可以組合成直接型、正準型、級聯型、并聯型四種結構形式，都具有反饋回路。同時，IIR數字濾波器在設計上可以借助

6、成熟的模擬濾波器的成果，如巴特沃斯、契比雪夫和橢圓濾波器等，有現成的設計數據或圖表可查，在設計一個IIR數字濾波器時，我們根據指標先寫出模擬濾波器的公式，然后通過一定的變換，將模擬濾波器的公式轉換成數字濾波器的公式。濾波器的設計可以通過軟件或設計專用的硬件兩種方式來實現。隨著MATLAB軟件及信號處理工具箱的不斷完善，MATLAB很快成為應用學科等領域不可或缺的基礎軟件。它可以快速有效地實現數字濾波器的設計、分析和仿真，極大地減輕了工作量,有利于濾波器設計的最優化。2 課程設計內容錄制一段自己的語音信號（可以錄制含有噪音的信號，或者錄制語音后再加進噪音信號），要求：（1）對語音信號進行采樣；（

7、2）畫出采樣后語音信號的時域波形和頻譜圖；（3）設計一個合適的濾波器，并畫出濾波器的頻率響應曲線；（4）用設計的濾波器對語音信號進行濾波，畫出濾波后信號的時域波形和頻譜圖；（5）對濾波前后的信號進行對比，分析信號的變化；（6）回放語音信號。3 課程設計的具體實現3.1 語音信號的采集利用PC 機上的聲卡和WINDOWS 操作系統可以進行數字信號的采集。將話筒輸入計算機的語音輸入插口上,啟動錄音機。按下錄音按鈕,接著對話筒說一段話,說完后停止錄音,屏幕左側將顯示所錄聲音的長度。點擊放音按鈕,可以實現所錄音的重現。以文件名“speech”保存入D: MATLAB work 中。可以看到,文件存儲器

8、的后綴默認為. wav ,這是WINDOWS 操作系統規定的聲音文件存的標準。3.2 語音信號的時頻分析利用MATLAB中的“wavread”命令來讀入（采集）語音信號，將它賦值給某一向量。再對其進行采樣，記住采樣頻率和采樣點數。下面介紹Wavread 函數幾種調用格式。（1）y=wavread（file）功能說明：讀取file所規定的wav文件，返回采樣值放在向量y中。（2）y,fs,nbits=wavread(file) 功能說明：采樣值放在向量y中，fs表示采樣頻率（hz），nbits表示采樣位數。（3）y=wavread（file，N）功能說明：讀取錢N點的采樣值放在向量y中。（4）y

9、=wavread（file，N1,N2）功能說明：讀取從N1到N2點的采樣值放在向量y中。接下來，對語音信號speech.wav進行采樣。其程序如下：>> y,fs,nbits=wavered (speech.wav); %把語音信號加載入Matlab 仿真軟件平臺中。然后，畫出語音信號的時域波形，再對語音信號進行頻譜分析。MATLAB提供了快速傅里葉變換算法FFT計算DFT的函數fft,其調用格式如下：Xk=fft(xn,N)參數xn為被變換的時域序列向量，N是DFT變換區間長度，當N大于xn的長度時，fft函數自動在xn后面補零。，當N小于xn的長度時，fft函數計算xn的前N

10、個元素，忽略其后面的元素。在本次設計中，我們利用fft對語音信號進行快速傅里葉變換，就可以得到信號的頻譜特性。其程序如下：y,fs,bits=wavread('speech.wav'); % sound(y,fs,bits); Y=fft(y,40000);%采樣點數subplot(211);plot(y);title('原始信號波形'); subplot(212);plot(abs(Y);title('原始信號頻譜'); 程序運行如下：3.3 語音信號加噪與頻譜分析在原始信號中加入高頻噪音，其具體的程序如下：fc=20000;%載波頻率y1=m

11、odulate(y,fc,fs,'fm');%對原語音信號調制Y1=fft(y1,40000);%采樣點數40000 figure(2); subplot(211);plot(y1);title('調制后信號波形'); subplot(212);plot(abs(Y1);title('調制后信號頻譜'); % sound(y1); 對原信號進行調制得到的波形和頻譜圖如下： 接著對調制后的信號加入原始信號，其程序為：% 調制信號+原始信號波形y2=y+y1; Y2=fft(y2,40000);%采樣點數40000 figure(3)subplot(

12、211);plot(y2);title('調制信號+原始信號波形'); subplot(212);plot(abs(Y2);title('調制信號+原始信號頻譜');% sound(y2); 得到的圖形為：; 3.4 利用雙線性變換法設計低通濾波器 用雙線性變換法設計低通濾波器程序如下所示：figure(4); % 用雙線性變換法設計的低通濾波器fp=1000;fc=1200;As=100;Ap=1;fs1=22050; wc=2*fc/fs1;wp=2*fp/fs1; n,wn=ellipord(wp,wc,Ap,As); b,a=ellip(n,Ap,As,

13、wn); freqz(b,a,512,fs1); x=filter(b,a,y2);%對要進行低通濾波X=fft(x,40000); 運行結果如圖：3.5 用濾波器對加噪語音信號進行濾波利用設計的低通濾波器對信號進行濾波，程序如下：figure(5); subplot(211);plot(x);title('低通濾波后信號波形'); subplot(212);plot(abs(X);title('低通濾波后信號頻譜'); sound(x,40000) 運行結果如下圖：3.6 分析濾波前后語音信號波形及頻譜的變化通過比較加噪前后，語音的頻譜和語音回放，能明顯的感覺

14、到加入噪聲后回放的聲音與原始的語音信號有很大的不同,前者隨較尖銳的干擾嘯叫聲。從含噪語音信號的頻譜圖中可以看出含噪聲的語音信號頻譜,在整個頻域范圍內分是布均勻。其實，這正是干擾所造成的。通過濾波前后的對比，低通濾波后效果最好，高通濾波后的效果最差。由此可見，語音信號主要分布在低頻段，而噪聲主要分布在高頻段。3.7 回放語音信號 經過以上的加噪聲處理后，可在Matlab中用函數sound對聲音進行回放。其調用格式：sound(y,Fs)，sound(y)和sound(y,Fs,bits)。可以察覺濾波前后的聲音有明顯的變化。3.8小結 本設計圓滿的完成了對加噪聲語音信號的讀取與打開，與課題的要求

15、十分相符；也較好的完成了對原始語音信號和加噪聲后的語音信號的頻譜分析，通過fft變換，得出了語音信號的頻譜圖；在濾波這一部分，課題主要是從雙線性變換法入手來設計濾波器等入手，實現了預期的濾波效果；完成了界面的設計，但也存在許多的不足，只是很勉強的達到了打開語音文件、顯示已定濾波前后的波形圖，界面沒有用到其它的一些控件，如列表框、下拉菜單等等。結 論語音信號處理是語音學與數字信號處理技術相結合的交叉學科，課題在這里不討論語音學，而是將語音當做一種特殊的信號，即一種“復雜向量”來看待。也就是說，課題更多的還是體現了數字信號處理技術。從課題的中心來看，課題“基于MATLAB的有噪聲語音信號處理”是希

16、望將數字信號處理技術應用于某一實際領域，這里就是指對語音及加噪處理。作為存儲于計算機中的語音信號，其本身就是離散化了的向量，我們只需將這些離散的量提取出來，就可以對其進行處理了。這一過程的實現，用到了處理數字信號的強有力工具MATLAB。通過MATLAB里幾個命令函數的調用，很輕易的在實際語音與數字信號的理論之間搭了一座橋。課題的特色在于它將語音信號看作一個向量，于是就把語音數字化了。那么，就可以完全利用數字信號處理的知識來解決語音及加噪處理問題。我們可以像給一般信號做頻譜分析一樣，來對語音信號做頻譜分析，也可以較容易的用數字濾波器來對語音進行濾波處理。通過比較加噪前后，語音的頻譜和語音回放，

17、能明顯的感覺到加入噪聲后回放的聲音與原始的語音信號有很大的不同,前者隨較尖銳的干擾嘯叫聲。從含噪語音信號的頻譜圖中可以看出含噪聲的語音信號頻譜,在整個頻域范圍內分是布均勻。其實，這正是干擾所造成的。通過濾波前后的對比，低通濾波后效果最好，高通濾波后的效果最差。由此可見，語音信號主要分布在低頻段，而噪聲主要分布在高頻段。致 謝本次課程設計選題及進行過程中得到同學們的幫助，尤其是我們同一課題的幾個同學，我們花費課很多的時間和精力，相互之間幫忙協作，上網搜索相關資料，到圖書館查閱相關文獻，遇到難題，共同商討。解決不了的問題，我們就向老師和其他同學虛心請教。最終，我們一起解決了一個又一個難題，雖然，我們有過爭吵，但是在真理面前，我們的行動是一致的，在課程設計過程中， 學院的機房工作人員給我們提供的便利的條件，天氣寒