語音信號處理相關(guān)論文

1、摘 要語音信號處理包括語音通信、語音增強(qiáng)、語音合成、語音識別和說話人識別等方面。只有通過語音信號的數(shù)字處理，語音信號的好壞、語音識別率的高低，都取決于語音信號處理的好壞。因此，語音信號處理是一項(xiàng)非常有意義的研究課程。語音，即語言的聲音，是語言符號系統(tǒng)的載體。它由人的發(fā)音器官發(fā)出，負(fù)載著一定的語言意義。語言依靠語音實(shí)現(xiàn)它的社會功能。語言是音義結(jié)合的符號系統(tǒng)，語言的聲音和語言的意義是緊密聯(lián)系著的，因此，語言雖是一種聲音，但又與一般的聲音有著本質(zhì)的區(qū)別語音是人們交流思想和進(jìn)行社會活動的最基本手段,因此我們要對語音信號進(jìn)行處理分析，優(yōu)化人類通信交流。本次實(shí)訓(xùn)用windows的錄音機(jī)錄制一段自己說的話（

2、語音信號），并畫出語音信號的時(shí)域波形和頻譜圖；在MATLAB軟件平臺下，給語音信號加上噪聲，繪出加噪聲后的語音信號時(shí)域波形和頻譜；采用雙線性變換法，設(shè)計(jì)一個(gè)濾波器，并畫出濾波器的頻率響應(yīng)；用自己設(shè)計(jì)的濾波器對語音信號進(jìn)行濾波，畫出濾波后信號的時(shí)域波形和頻譜圖，并對濾波前、后的信號進(jìn)行對比，畫出時(shí)域波形差，分析語言信號濾波前后的變化。本次實(shí)訓(xùn)是利用MATLAB對語音信號進(jìn)行時(shí)域和頻域分析，通過MATLAB對語音進(jìn)行數(shù)字處理。 關(guān)鍵詞：語音；噪聲；MATLAB；濾波器；信號處理AbstractSpeech signal processing including voice communicatio

3、ns, speech enhancement, speech synthesis, speech recognition and speaker recognition and so on. Only through the digital processing of voice signals, voice signals good or bad, high and low speech recognition rate are good or bad depends on the voice signal process .Therefore, the voice signal proce

4、ssing is a very significant research programs.Voice, the voice of the language is the carrier of linguistic sign systems. It is issued by the human vocal organs, load a certain amount of linguistic meaning. Language dependent voice to achieve its social function. Language is the combination of sound

5、 and meaning of the symbol system, language, the meaning of the voice and language are closely linked, and therefore, although the language is a voice, but the sound in general is essentially different. Voice that people exchange ideas and conduct the most basic means of social activities, so we hav

6、e to deal with the speech signal analysis and optimization of human communication exchange. The training record a tape recorder with windows own words (audio signal), and draw the speech signal in time domain waveform and frequency spectrum; In the MATLAB software platform, to the speech signal with

7、 noise, drawn after the speech signal plus noise time domain waveform and frequency spectrum; Bilinear transform method, to design a filter, and draw the frequency response filter; Designed to filter their own filtering of the speech signal, draw the filtered signal in time domain waveform and frequ

8、ency spectrum, and filter before and after the signals are compared to draw difference time-domain waveform analysis language signal changes before and after filtering . The training is the use of MATLAB on the voice signal in time domain and frequency domain analysis, carried out by MATLAB digital

9、processing of speech. Key word: Voice ; Noise; MATLAB; Filter ;Signal Processing引言語言是人類持有的功能聲音是人類常用的工具，是相互傳遞信息的最主要的手段。語音信號是人們構(gòu)成思想疏通和感情交流的最主要的途徑，通過語音傳遞倍息是人類最重要、最有效、最常用和最方便的交換信息的形式。現(xiàn)在隨著時(shí)代的發(fā)展，人類已開始進(jìn)入了信息化時(shí)代，用現(xiàn)代手段研究語音信號，使人們能更加有效地產(chǎn)生、傳輸、存儲、獲取和應(yīng)用語音信息，這對于促進(jìn)社會的發(fā)展具有十分重要的意義。在本次實(shí)訓(xùn)中我設(shè)計(jì)的題目是：語音信號的數(shù)字處理。實(shí)訓(xùn)的設(shè)計(jì)是對語音信號的數(shù)

10、字處理，通過對所采集的語音信號進(jìn)行時(shí)域波形和頻譜圖分析，采用雙線性變換法設(shè)計(jì)濾波器，語音信號通過對濾波器進(jìn)行濾波，然后回放語音信號，并對其濾波前和濾波后的語音變化進(jìn)行分析?，F(xiàn)在社會衍生了很多現(xiàn)代的語音通信方式，手機(jī)語音、語音聊天軟件以及語音小說等，涉及語音的方面幾乎遍及了我們的大部分生活，和我們生活息息相關(guān)，因此語音信號處理是具有現(xiàn)實(shí)意義的研究課程。1. 語音信號降噪處理分析MATLAB介紹1.1 語音信號 語音信號的采集與分析作為一個(gè)重要的研究領(lǐng)域, 20世紀(jì)60年代中期形成的一系列數(shù)字信號處理方法和技術(shù)，如數(shù)字濾波器、快速博里葉變換(FFT)等成為語音信號數(shù)字處理的理論和技術(shù)基礎(chǔ)。在方法上

11、，隨著電子計(jì)算機(jī)的發(fā)展，以往的以硬件為中心的研究逐漸轉(zhuǎn)化為以軟件為主的處理研究。然而，在語音識別領(lǐng)域內(nèi)，初期有幾種語音打字機(jī)的研究也很活躍，但后來已全部停了下來，這說明了當(dāng)時(shí)人們對話音識別難度的認(rèn)識得到了加深。所以1969年美國貝爾研究所的Pierce感嘆地說“語音識別向何處去?”1971年，以美國ARPA(American Research Projects Agency)為主導(dǎo)的“語音理解系統(tǒng)”的研究計(jì)劃也開始起步。這個(gè)研究計(jì)劃不僅在美國園內(nèi)，而且對世界各國都產(chǎn)生了很大的影響，它促進(jìn)了連續(xù)語音識別研究的興起。 進(jìn)入20世紀(jì)90年代以來，語音信號采集與分析在實(shí)用化方面取得了許多實(shí)質(zhì)性的研究進(jìn)

12、展。其中，語音識別逐漸由實(shí)驗(yàn)室走向?qū)嵱没?。一方面，對聲學(xué)語音學(xué)統(tǒng)計(jì)模型的研究逐漸深入，魯棒的語音識別、基于語音段的建模方法及隱馬爾可夫模型與人工種經(jīng)網(wǎng)絡(luò)的結(jié)合成為研究的熱點(diǎn)。另一方面，為了語音識別實(shí)用化的需要，講者自適應(yīng)、聽覺模型、快速搜索識別算法以及進(jìn)一步的語言模型的研究等課題倍受關(guān)注。 其中MATLAB音頻信號的分析與合成，采用MATLAB分析WAV文件。首先選擇一個(gè)WAV文件作為分析的對象 。繼而進(jìn)行 FFT變換并畫頻域圖形 ，進(jìn)行該聲波主要頻譜的分析，根據(jù)該聲音的頻譜，反演時(shí)域圖形 ，最后進(jìn)行付立葉逆變換IFFT并畫頻域圖形 。1.2 語音信號的特點(diǎn)通過對大量語音信號的觀察和分析發(fā)現(xiàn)，

13、語音信號主要有下面兩個(gè)特點(diǎn)：第一、在頻域內(nèi)，語音信號的頻譜分量主要集中在3003400Hz的范圍內(nèi)。利用這個(gè)特點(diǎn)，可以用一個(gè)防混迭的帶通濾波器將此范圍內(nèi)的語音信號頻率分量取出，然后按8kHz的采樣率對語音信號進(jìn)行采樣，就可以得到離散的語音信號。第二、在時(shí)域內(nèi)，語音信號具有“短時(shí)性”的特點(diǎn)，即在總體上，語音信號的特征是隨著時(shí)間而變化的，但在一段較短的時(shí)間間隔內(nèi)，語音信號保持平穩(wěn)。在濁音段表現(xiàn)出周期信號的特征，在清音段表現(xiàn)出隨機(jī)噪聲的特征。下面是一段語音信號的時(shí)域波形圖(圖2-1)和頻域圖(圖2-2)，由這兩個(gè)圖可以看出語音信號的兩個(gè)特點(diǎn)。圖1-1語音信號時(shí)域波形圖1-2 語音信號頻域波形1.3

14、語音信號的采集在將語音信號進(jìn)行數(shù)字化前，必須先進(jìn)行防混疊預(yù)濾波，預(yù)濾波的目的有兩個(gè)：（1）抑制輸入信導(dǎo)各領(lǐng)域分量中頻率超出的所有分量（為采樣頻率），以防止混疊干擾。（2）抑制50Hz的電源工頻干擾。為了將原始模擬語音信號變?yōu)閿?shù)字信號，必須經(jīng)過采樣和量化兩個(gè)步驟，從而得到時(shí)間和幅度上均為離散的數(shù)字語音信號。采樣也稱抽樣，是信號在時(shí)間上的離散化，即按照一定時(shí)間間隔在模擬信號上逐點(diǎn)采取其瞬時(shí)值。采樣時(shí)必須要注意滿足奈奎斯特定理，即采樣頻率必須以高于受測信號的最高頻率兩倍以上的速度進(jìn)行取樣，才能正確地重建波它是通過采樣脈沖和模擬信號相乘來實(shí)現(xiàn)的。下圖是一段語音信號在采樣頻率44.1KHz情況下的頻譜圖

15、。圖1-3 語音信號頻譜圖 由圖可知，這段語音信號的頻率主要集中在1KHz左右，當(dāng)采樣頻率為44.1KHz時(shí)，由于采樣頻率比較大，所以采樣點(diǎn)數(shù)就越密，所得離散信號就越逼近于原信號，頻譜也沒有發(fā)生混疊。對上述信號進(jìn)行1/80采樣頻率抽取，即采樣頻率變?yōu)閷⒔?00Hz時(shí)，由于采樣頻率比較小，所以采樣點(diǎn)數(shù)就稀疏，所得離散信號就越偏離于原信號，頻譜也發(fā)生了混疊。在采樣的過程中應(yīng)注意采樣間隔的選擇和信號混淆：對模擬信號采樣首先要確定采樣間隔。如何合理選擇涉及到許多需要考慮的技術(shù)因素。一般而言，采樣頻率越高，采樣點(diǎn)數(shù)就越密，所得離散信號就越逼近于原信號。但過高的采樣頻率并不可取，對固定長度（T）的信號，采

16、集到過大的數(shù)據(jù)量，給計(jì)算機(jī)增加不必要的計(jì)算工作量和存儲空間；若數(shù)據(jù)量（N）限定，則采樣時(shí)間過短，會導(dǎo)致一些數(shù)據(jù)信息被排斥在外。采樣頻率過低，采樣點(diǎn)間隔過遠(yuǎn)，則離散信號不足以反映原有信號波形特征，無法使信號復(fù)原，造成信號混淆。根據(jù)采樣定理，當(dāng)采樣頻率大于信號的兩倍帶寬時(shí)，采樣過程不會丟失信息，利用理想濾波器可從采樣信號中不失真地重構(gòu)原始信號波形。量化是對幅值進(jìn)行離散化，即將振動幅值用二進(jìn)制量化電平來表示。量化電平按級數(shù)變化，實(shí)際的振動值是連續(xù)的物理量。具體振值用舍入法歸到靠近的量化電平上。 語音信號經(jīng)過預(yù)濾波和采樣后，由A/D變換器變換為二址制數(shù)字碼。這種防混疊濾波通常與模數(shù)轉(zhuǎn)換器做在一個(gè)集成塊

17、內(nèi)，因此目前來說，語音信號的數(shù)字化的質(zhì)量還是有保證的。市面上購買到的普通聲卡在這方面做的都很好，語音聲波通過話筒輸入到聲卡后直接獲得的是經(jīng)過防混疊濾波、A/D變換、量化處理的離散的數(shù)字信號。在實(shí)際工作中，我們可以利用windows自帶的錄音機(jī)錄制語音文件，基于PC機(jī)的語音信號采集過程，聲卡可以完成語音波形的A/D轉(zhuǎn)換，獲得WAVE文件，為后續(xù)的處理儲備原材料。調(diào)節(jié)錄音機(jī)保存界面的“更改”選項(xiàng)，可以存儲各種格式的WAV文件。采集到語音信號之后，需要對語音信號進(jìn)行分析，如畫出采樣后語音信號的時(shí)域波形和頻譜圖；雙線性變換設(shè)計(jì)濾波器，并畫出濾波器的頻率響應(yīng)；然后用自己設(shè)計(jì)的濾波器對采集的語音信號進(jìn)行濾

18、波，畫出濾波后信號的時(shí)域波形和頻譜，并對濾波前后的信號進(jìn)行對比，分析信號的變化；回放語音信號。1.4 信號頻譜分析 信號的頻譜分析是將信號源發(fā)出的信號強(qiáng)度按頻率順序展開，使其成為頻率的函數(shù)，并考察變化規(guī)律，稱為頻譜分析。頻譜分析主要分析信號是由哪些頻率的正弦信號疊加得到的，以及這些正弦信號的振幅。信號的頻譜可分為幅值譜、相位譜、功率譜、對數(shù)譜等等。對信號作頻譜分析的設(shè)備主要是頻譜分析儀，它把信號按數(shù)學(xué)關(guān)系作為頻率的函數(shù)顯示出來，其工作方式有模擬式和數(shù)字式二種。模擬式頻譜分析儀以模擬濾波器為基礎(chǔ)，從信號中選出各個(gè)頻率成分的量值；數(shù)字式頻譜分析儀以數(shù)字濾波器或快速傅立葉變換為基礎(chǔ)，實(shí)現(xiàn)信號的時(shí)頻關(guān)

19、系轉(zhuǎn)換分析。傅立葉變換是信號頻譜分析中常用的一個(gè)工具，它把一些復(fù)雜的信號分解為無窮多個(gè)相互之間具有一定關(guān)系的正弦信號之和，并通過對各個(gè)正弦信號的研究來了解復(fù)雜信號的頻率成分和幅值。信號頻譜分析是采用傅立葉變換將時(shí)域信號x(t)變換為頻域信號X(f)，從而幫助人們從另一個(gè)角度來了解信號的特征。時(shí)域信號x(t)的傅氏變換為： 式中X(f)為信號的頻域表示，x(t)為信號的時(shí)域表示，f為頻率。頻譜分析主要用于識別信號中的周期分量，是信號分析中最常用的一種手段。典型信號正弦波、方波、三角波和白噪聲信號是實(shí)際工程測試中常見的典型信號，這些信號時(shí)域、頻域之間的關(guān)系很明確，并且都具有一定的特性，通過對這些典

20、型信號的頻譜進(jìn)行分析，對掌握信號的特性，熟悉信號的分析方法大有益處，并且這些典型信號也可以作為實(shí)際工程信號分析時(shí)的參照資料。本次實(shí)驗(yàn)利用DRVI快速可重組虛擬儀器平臺可以很方便的對上述典型信號作頻譜分析。1.5 信號的時(shí)域分析語音信號的時(shí)域分析就是分析和提取語音信號的時(shí)域參數(shù)。進(jìn)行語音分析時(shí)，最先接觸到并且也是最直觀的是它的時(shí)域波形。語音信號本身就是時(shí)域信號，因而時(shí)域分析是最早使用，也是應(yīng)用最廣泛的一種分析方法，這種方法直接利用語音信號的時(shí)域波形。時(shí)域分析通常用于最基本的參數(shù)分析及應(yīng)用，如語音的分割、預(yù)處理、大分類等。這種分析方法的特點(diǎn)是：（1）表示語音信號比較直觀、物理意義明確；（2）實(shí)現(xiàn)起

21、來比較簡單、運(yùn)算且少；（3）可以得到語音的一些重要的參數(shù)；（4）只使用示波器等通用設(shè)備，使用較為簡單等。語音信號的時(shí)域參數(shù)有短時(shí)能量、短時(shí)過零率、短時(shí)白相關(guān)函數(shù)和短時(shí)平均幅度差函數(shù)等，這是語音信號的一組最基本的短時(shí)參數(shù)，在各種語音信號數(shù)字處理技術(shù)中都要應(yīng)用。在計(jì)算這些參數(shù)時(shí)使用的一般是凱澤窗或漢明窗。2. MATLAB軟件功能介紹2.1 matlab簡介MATLAB是矩陣實(shí)驗(yàn)室（Matrix Laboratory）的簡稱，是美國MathWorks公司出品的商業(yè)數(shù)學(xué)軟件，用于算法開發(fā)、數(shù)據(jù)可視化、數(shù)據(jù)分析以及數(shù)值計(jì)算的高級技術(shù)計(jì)算語言和交互式環(huán)境，主要包括MATLAB和Simulink兩大部分。

22、MATLAB和Mathematica、Maple并稱為三大數(shù)學(xué)軟件。它在數(shù)學(xué)類科技應(yīng)用軟件中在數(shù)值計(jì)算方面首屈一指。MATLAB可以進(jìn)行矩陣運(yùn)算、繪制函數(shù)和數(shù)據(jù)、實(shí)現(xiàn)算法、創(chuàng)建用戶界面、連接其他編程語言的程序等，主要應(yīng)用于工程計(jì)算、控制設(shè)計(jì)、信號處理與通訊、圖像處理、信號檢測、金融建模設(shè)計(jì)與分析等領(lǐng)域。MATLAB 的應(yīng)用范圍非常廣，包括信號和圖像處理、通訊、控制系統(tǒng)設(shè)計(jì)、測試和測量、財(cái)務(wù)建模和分析以及計(jì)算生物學(xué)等眾多應(yīng)用領(lǐng)域。附加的工具箱（單獨(dú)提供的專用 MATLAB 函數(shù)集）擴(kuò)展了 MATLAB 環(huán)境，以解決這些應(yīng)用領(lǐng)域內(nèi)特定類型的問題。Matlab優(yōu)勢：（1）友好的工作平臺和編程環(huán)境；（

23、2）簡單易用的程序語言；（3）強(qiáng)大的科學(xué)計(jì)算機(jī)數(shù)據(jù)處理能力；（4）出色的圖形處理功能；（5）應(yīng)用廣泛的模塊集合工具箱；（6）實(shí)用的程序接口和發(fā)布平臺；（7）應(yīng)用軟件開發(fā)（包括用戶界面）；2.2基于MATLAB的語音處理Matlab是一個(gè)非常好的語音處理工具箱,對學(xué)習(xí)語音處理的讀者非常有用。基于DSP和MATLAB的語音數(shù)據(jù)采集和處理系統(tǒng)-MATLAB，該工具包主要應(yīng)用于語音處理中聽覺模型的建立和分析，其中也包含了語音處理的各種基本功能如語譜圖，線性預(yù)測，識別和合成等等語音信號的頻域處理,語音雖然是一個(gè)時(shí)變、非平穩(wěn)的隨機(jī)過程。但在短時(shí)間內(nèi)可近似看作是平穩(wěn)的。因此如果能從帶噪語音的短時(shí)譜中估計(jì)出“

24、純凈”語音的短時(shí)譜，即可達(dá)到語音增強(qiáng)的目的。由于噪聲也是隨機(jī)過程，因此這種估計(jì)只能建立在統(tǒng)計(jì)模型基礎(chǔ)上。利用人耳感知對語音頻譜分量的相位不敏感的特性，這類語音增強(qiáng)算法主要針對短時(shí)譜的幅度估計(jì)。3.MATLAB濾波器設(shè)計(jì)3.1 matlab濾波器設(shè)計(jì)法介紹3.3 模擬濾波器設(shè)計(jì)3.3.1設(shè)計(jì)指標(biāo)模擬低通濾波器的設(shè)計(jì)指標(biāo)有，和。其中和分別稱為通帶截止頻率和阻帶截止頻率，是通帶)中的最大衰減系數(shù)，是阻帶的最小衰減系數(shù)，和一般用dB數(shù)表示。在本次實(shí)訓(xùn)中%用切比雪夫低通濾波器 t=1;fs=1/t; fp=0.2; fr=0.3;Ap=1;Ar=25;wp=2*pi*fp/fs;wr=2*pi*fr/f

25、s;wpp=2/t*tan(wp/2);wrr=2/t*tan(wr/2);n,wn=cheb1ord(wpp,wrr,Ap,Ar,s);b,a=cheby1(n,Ap,wn,low,s);bz,az=bilinear(b,a,fs);freqz(bz,az,512,fs)3.3.2巴特沃斯濾低通濾波器的設(shè)計(jì)%用巴特沃斯低通濾波器 t=1;fs=1/t;fp=0.2;fr=0.3;Ap=1;Ar=20;wp=2*pi*fp/fs;wr=2*pi*fr/fs;wpp=2/t*tan(wp/2);wrr=2/t*tan(wr/2);n,wn=buttord(wpp,wrr,Ap,Ar,s);b,a

26、=butter(n,wn,low,s);bz,az=bilinear(b,a,fs);freqz(bz,az,512,fs)典型的模擬濾波器有多種，常用的有巴特沃斯濾波器和切比雪夫?yàn)V波器。其設(shè)計(jì)步驟如下：根據(jù)技術(shù)指標(biāo)，和，求出濾波器的階數(shù)N。求出歸一化極點(diǎn),將代入傳輸函數(shù)中，得到歸一化傳輸函數(shù)Ha(p)。 將Ha(p)去歸一化。將代入Ha(p)，得到實(shí)際的濾波器傳輸函數(shù)Ha(s)。以上是巴特沃斯低通濾波器，通過頻率變換法我們還可以靈活地設(shè)計(jì)高通，帶通等模擬濾波器。3.4 數(shù)字濾波器3.4.1IIR濾波器的設(shè)計(jì)IIR濾波器常用的一類設(shè)計(jì)方法是借助于模擬濾波器的設(shè)計(jì)方法進(jìn)行的。通過模擬濾波器得到傳

27、輸函數(shù)Ha(s)后，再將Ha(s)按某種方法轉(zhuǎn)換成數(shù)字濾波器的系統(tǒng)函數(shù)H(z)。轉(zhuǎn)換方法有脈沖響應(yīng)不變法和雙線性變換法。脈沖響應(yīng)不變法所謂脈沖響應(yīng)不變法就是使數(shù)字濾波器的單位脈沖響應(yīng)序列h(n)等于模擬濾波器的單位沖激響應(yīng)和(t)的采樣值，即：，其中，T為采樣周期。在MATLAB中，可用函數(shù)impinvar實(shí)現(xiàn)從模擬濾波器到數(shù)字濾波器的脈沖響應(yīng)不變映射，調(diào)用格式為：Bz,Az=impinvar(B,A,Fs(3.6)Bz,Az=impinvar(B,A)(3.7)其中，B、A分別為模擬濾波器的分子和分母多項(xiàng)式系數(shù)向量；Fs為采樣頻率（Hz），缺省值Fs=1Hz；Bz、Az分別為數(shù)字濾波器分子和

28、分母多項(xiàng)式系數(shù)向量。雙線性變換法圖3.3 s平面和z平面的單值雙線性映射關(guān)系示意圖平面和z平面的單值雙線性映射關(guān)系如圖3.3所示由于s平面和z平面的單值雙線性映射關(guān)系為s=，其中T為采樣周期。因此，若已知模擬濾波器的傳遞函數(shù)，將上式代入即可得到數(shù)字濾波器的系統(tǒng)函數(shù)H(z)。在雙線性變換中，模擬角頻率和數(shù)字角頻率的變換關(guān)系為：(3.8)可見，和w之間的變換關(guān)系為非線性的。設(shè)計(jì)步驟3.2.1 IIR數(shù)字濾波器的設(shè)計(jì)步驟由于采用了級聯(lián)結(jié)構(gòu)，因此如何將濾波器的每一個(gè)極點(diǎn)和零點(diǎn)相組合，從而使得數(shù)字濾波器輸出所含的噪聲最小是個(gè)十分關(guān)鍵的問題。為了產(chǎn)生最優(yōu)的量化后的IIR數(shù)字濾波器，采用如下步驟進(jìn)行設(shè)計(jì)。

29、(1)、首先計(jì)算整體傳遞函數(shù)的零極點(diǎn)； (2)、選取具有最大幅度的極點(diǎn)以及距離它最近的零點(diǎn)，使用它們組成一個(gè)二階基本節(jié)的傳遞函數(shù)； (3)、對于剩下的極點(diǎn)和零點(diǎn)采用與3.2相類似的步驟，直至形成所有的二階基本節(jié)。 通過上面三步法進(jìn)行的設(shè)計(jì)可以保證IIR數(shù)字濾波器中N位乘法器產(chǎn)生的量化舍入誤差最小。利用模擬濾波器設(shè)計(jì)IIR數(shù)字低通濾波器的步驟如下：確定數(shù)字低通濾波器的技術(shù)指標(biāo)：通帶截止頻率p、通帶衰減p、阻帶截止頻率s、阻帶衰減s。將數(shù)字低通濾波器的技術(shù)指標(biāo)轉(zhuǎn)換成模擬低通濾波器的技術(shù)指標(biāo)。(3.9)如果采用雙線性變換法，邊界頻率的轉(zhuǎn)換關(guān)系為（3.10)按照模擬低通濾波器的技術(shù)指標(biāo)設(shè)計(jì)模擬低通濾波

30、器。 將模擬濾波器Ha(s)，從s平面轉(zhuǎn)換到z平面，得到數(shù)字低通濾波器系統(tǒng)函數(shù)H(z)?；贐utterworth模擬濾波器原型，使用雙線性狀換設(shè)計(jì)數(shù)字濾波器：各參數(shù)值為：通帶截止頻率Ap=0.2*pi,阻帶截止頻率Ar=0.3*pi,通帶波動值Rp=1dB,阻帶波動值Rs=15dB,設(shè)Fs=3400Hz。代碼：%用巴特沃斯數(shù)字低通濾波器 wp=0.2*pi;ws=0.3*pi;Fs=3400;T=1/Fs; Ap=(2/T)*tan(wp/2);Ar=(2/T)*tan(ws/2);rp=0.6;rs=8;as=50;ripple=10(-rp/20);attn=10(-rs/20);n,w

31、n=buttord(Ap,Ar,rp,rs,s);z,p,k=Buttap(n);b,a=zp2tf(z,p,k);bt,at=lp2lp(b,a,wn);b,a=bilinear(bt,at,Fs);3.4.3 濾波器的MATLAB實(shí)現(xiàn)MATLAB設(shè)計(jì)FIR濾波器的相關(guān)文件buttord.m用來確定數(shù)字低通或模擬低通濾波器的階次，其調(diào)用格式分別是 N,Wn=buttord(Wp,Ws,Rp,Rs)(3.14) N,Wn=buttord(Wp,Ws,Rp,Rs, s)(3.15)格式3.14中的Wn，Wp，Ws為歸一化頻率，其值在01之間，1對應(yīng)抽樣頻率的一半；格式3.15中Wn，Wp，Ws單

32、位為rad/s，因此，它們實(shí)際上是頻率。對于不同類型的濾波器，參數(shù)Wp和Ws有一些限制：對于低通濾波器，WpWs；對于帶通濾波器，wp和ws分別為具有兩個(gè)元素的矢量，Wp=Wp1,Wp2和Ws=Ws1,Ws2，并且Ws1Wp1Wp2Ws2；對于帶阻濾波器Wp1Ws1Ws2 wp=0.2*pi;ws=0.3*pi;Fs=3400;T=1/Fs; Ap=(2/T)*tan(wp/2);Ar=(2/T)*tan(ws/2);rp=0.6;rs=8;as=50;ripple=10(-rp/20);attn=10(-rs/20);n,wn=buttord(Ap,Ar,rp,rs,s);z,p,k=But

33、tap(n);b,a=zp2tf(z,p,k);bt,at=lp2lp(b,a,wn);b,a=bilinear(bt,at,Fs);濾波器的頻率響應(yīng)圖如下圖1-8 濾波器的頻率響應(yīng)圖 4.2.4 語音信號進(jìn)行濾波對比分析用自己設(shè)計(jì)的濾波器對語音信號進(jìn)行濾波，畫出濾波后信號的時(shí)域波形和頻譜圖，并對濾波前、后的信號進(jìn)行對比，畫出時(shí)域波形差，分析信號的變化。在MATLAB中，用巴特沃斯數(shù)字低通濾波器，其程序設(shè)計(jì)如下：%用巴特沃斯數(shù)字低通濾波器 wp=0.2*pi;ws=0.3*pi;Fs=3400;T=1/Fs; Ap=(2/T)*tan(wp/2);Ar=(2/T)*tan(ws/2);rp=0

34、.6;rs=8;as=50;ripple=10(-rp/20);attn=10(-rs/20);n,wn=buttord(Ap,Ar,rp,rs,s);z,p,k=Buttap(n);b,a=zp2tf(z,p,k);bt,at=lp2lp(b,a,wn);b,a=bilinear(bt,at,Fs);y,fs,bits=wavread(C:Documents and SettingsAdministrator.COMPUTER桌面1.wav);d=0.02*randn(length(y),2);x1=y+d;y1=fft(x1,2048);d=filter(b,1,x1); D=fft(d)

35、; subplot(221);plot(x1);grid on;axis tight;title(濾波前的加噪語音信號);xlabel(time(s);ylabel(幅度);subplot(222);plot(abs(y1);grid on;axis tight;title(濾波前的加噪語音信號頻譜);xlabel(Hz);ylabel(幅度);subplot(223);plot(d);grid on;axis tight;title(濾波后的語音信號); xlabel(time(s);ylabel(幅度);subplot(224);plot(abs(D);grid on;axis tight

36、;title(濾波后的語音信號頻譜);xlabel(Hz);ylabel(幅度);sound(y,fs,bits);濾波前后信號的時(shí)域波形和頻譜圖對比圖如下：5. 總結(jié)通過本次數(shù)字信號處理實(shí)訓(xùn)，我對MATLAB的相關(guān)知識有了更深了學(xué)習(xí)了解，MATLA是一種以數(shù)值計(jì)算和數(shù)據(jù)圖示為主的計(jì)算機(jī)軟件，并包含適應(yīng)多個(gè)學(xué)科的專業(yè)軟件包，以及完善程序開發(fā)功能。本實(shí)訓(xùn)設(shè)計(jì)旨在要求學(xué)生熟悉Matlab的功能，掌握Matlab程序設(shè)計(jì)，培養(yǎng)學(xué)生的應(yīng)用開發(fā)能力，為以后的畢業(yè)設(shè)計(jì)奠定一定的基礎(chǔ)。首先通過自己對各類MATLAB知識的了解和興趣，選擇了合適自己這次的實(shí)訓(xùn)題目，我之所以選擇語言信號的降噪，是覺得我在Matlab語音方面了解的很少，借助這次實(shí)訓(xùn)的機(jī)會，讓自己學(xué)習(xí)更多關(guān)于Matlab語音方面的知識。 在這次的實(shí)訓(xùn)中，我對Matlab很多的知識進(jìn)行了再一次或是更深入的學(xué)習(xí)，對一些不是很懂，課本介紹又少的相關(guān)知識，我借助了網(wǎng)絡(luò)進(jìn)行了搜索學(xué)習(xí)，達(dá)到了更全面的學(xué)習(xí)和掌握，本次實(shí)訓(xùn)我涉及到Matlab一些基礎(chǔ)的知識，主要的就涉及了語音信號