matlab課程設(shè)計報告--基于matlab有噪聲語音信號處理_第1頁
已閱讀1頁,還剩23頁未讀, 繼續(xù)免費閱讀

下載本文檔

版權(quán)說明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請進行舉報或認領(lǐng)

文檔簡介

1、<p>  Matlab課程設(shè)計報告</p><p>  題目:基于MATLAB有噪聲語音信號處理</p><p>  系 (院): 計算機與信息工程學院 </p><p>  專 業(yè): 通信工程 </p><p>  班 級:

2、 </p><p>  指導教師: </p><p>  學年學期: 2011 ~ 2012 學年 第 2 學期</p><p><b>  簡介:</b></p><p>  我們通信工程專業(yè)在實踐中經(jīng)常碰到需要對已接收信號進行處理的情況

3、,而濾波器設(shè)計在數(shù)字信號處理中占有極其重要的地位。本課題基于MATLAB有噪音語音信號處理的設(shè)計與實現(xiàn),綜合運用數(shù)字信號處理的理論知識對加噪語音信號進行時域、頻域分析和濾波。通過理論推導得出相應結(jié)論,再利用MATLAB作為編程工具進行計算機實現(xiàn)。在設(shè)計實現(xiàn)的過程中,我們使用雙線性變換法設(shè)計IIR數(shù)字濾波器,對模擬加噪語音信號進行低通濾波、高通濾波及帶通濾波,并利用MATLAB作為輔助工具完成設(shè)計中的計算與圖形的繪制。</p>

4、<p><b>  1 緒論:</b></p><p>  數(shù)字信號處理是利用計算機或?qū)S锰幚碓O(shè)備,以數(shù)值計算的方法對信號進行采集、抽樣、變換、綜合、估值與識別等加工處理,借以達到提取信息和便于應用的目的。數(shù)字濾波器, 是數(shù)字信號處理中及其重要的一部分。本課題采用IIR濾波器對加噪聲音信號進行處理。</p><p>  IIR濾波器采用遞歸型結(jié)構(gòu),即結(jié)構(gòu)

5、上帶有反饋環(huán)路。IIR濾波器運算結(jié)構(gòu)通常由延時、乘以系數(shù)和相加等基本運算組成,可以組合成直接型、正準型、級聯(lián)型、并聯(lián)型四種結(jié)構(gòu)形式,都具有反饋回路。同時,IIR數(shù)字濾波器在設(shè)計上可以借助成熟的模擬濾波器的成果,如巴特沃斯、契比雪夫和橢圓濾波器等,有現(xiàn)成的設(shè)計數(shù)據(jù)或圖表可查,在設(shè)計一個IIR數(shù)字濾波器時,我們根據(jù)指標先寫出模擬濾波器的公式,然后通過一定的變換,將模擬濾波器的公式轉(zhuǎn)換成數(shù)字濾波器的公式。</p><p&g

6、t;  2.原始語音信號采集與處理</p><p>  2.1語音信號的采集</p><p>  由于MATLAB只識別格式為.wav的聲音文件,我們利用PC機上的聲卡和WINDOWS操作系統(tǒng)進行數(shù)字信號的采集。啟動錄音機進行錄音,以文件名“Orisound”保存入原程序所屬的文件夾中??梢钥吹?文件存儲器的后綴默認為.wav ,這是WINDOWS操作系統(tǒng)規(guī)定的聲音文件存的標準。</

7、p><p><b>  程序流程圖:</b></p><p>  2.2語音信號的時頻分析</p><p>  利用MATLAB中的“wavread”命令來讀入(采集)語音信號,將它賦值給某一向量。再對其進行采樣,記住采樣頻率和采樣點數(shù)。</p><p>  對語音信號Orisound.wav進行采樣其程序如下:</p

8、><p>  [y,fs,nbits]=wavread (‘OriSound’); %把語音信號加載入MATLAB仿真軟件平臺中</p><p>  畫出語音信號的時域波形,再對語音信號進行頻譜分析。在本次設(shè)計中,我們利用fft函數(shù)對語音信號進行快速傅里葉變換,就可以得到信號的頻譜特性。程序如下:</p><p>  [y,fs,nbits]=wavread(‘Oris

9、ound’); %語音信號的采集</p><p>  sound(y,fs,nbits); %語音信號的播放</p><p>  n=length(y) ; </p><p>  Y=fft(y,n); %快速傅里葉變換</p><p>&l

10、t;b>  figure;</b></p><p>  subplot(2,1,1);</p><p><b>  plot(y);</b></p><p>  title(‘原始信號波形’,’fontweight’,’bold’);</p><p>  axis([ 78000 80000 -1 1]

11、);</p><p><b>  grid;</b></p><p>  subplot(2,1,2);</p><p>  plot(abs(Y));</p><p>  title(‘原始信號頻譜’,’fontweight’,’bold’);</p><p>  axis([ 0 150000

12、 0 4000]);</p><p><b>  grid;</b></p><p><b>  程序結(jié)果如下圖:</b></p><p>  2.3語音信號加噪與頻譜分析</p><p>  利用MATLAB中的隨機函數(shù)Randn(m,n)產(chǎn)生噪聲加入到語音信號中,模仿語音信號被污染,并對其頻譜分

13、析。主要程序如下:</p><p>  [y,fs,nbits]=wavread('OriSound');</p><p>  sound(y,fs,nbits); </p><p>  n = length (y) ; </p><p>  Noise=0.2*randn(n,2);</p><p>

14、  s=y+Noise;</p><p><b>  sound(s);</b></p><p><b>  figure;</b></p><p>  subplot(2,1,1);</p><p><b>  plot(s);</b></p><p>

15、;  title('加噪語音信號的時域波形','fontweight','bold');</p><p>  axis([ 78000 80000 -1 1]);</p><p><b>  grid;</b></p><p>  S=fft(s); </p>

16、<p>  subplot(2,1,2);</p><p>  plot(abs(S));</p><p>  title('加噪語音信號的頻域波形','fontweight','bold');</p><p>  axis([ 0 150000 0 4000]);</p><p>&

17、lt;b>  grid;</b></p><p><b>  程序結(jié)果如下圖:</b></p><p><b>  3設(shè)計數(shù)字濾波器</b></p><p>  3.1 數(shù)字濾波器設(shè)計的基本思路</p><p>  數(shù)字濾波器的實現(xiàn)有兩個關(guān)鍵步驟:一個從數(shù)字域到模擬域間的變換,這個

18、變換實現(xiàn)了數(shù)字濾波器技術(shù)指標到模擬濾波器技術(shù)指標的轉(zhuǎn)換,同樣也實現(xiàn)了模擬濾波器系統(tǒng)函數(shù)到數(shù)字濾波器系統(tǒng)函數(shù)的轉(zhuǎn)換;另一個是從模擬濾波器技術(shù)指標到滿足該指標的模擬濾波器的設(shè)計。</p><p>  3.2 模擬濾波器概述</p><p>  用模擬—數(shù)字變換法設(shè)計IIR數(shù)字濾波器,首先必須設(shè)計一個模擬濾波器,它有許多不同的類型,主要有以下兩種類型:</p><p> 

19、?。?)、巴特沃思(Botterworth簡寫B(tài)W)濾波器。BW濾波器是根據(jù)幅頻特性在通帶內(nèi)具有最拼圖特性而定義的濾波器,對一個N階低通濾波器來說,所謂最平坦特性就是模擬函數(shù)的前(2N-1)階導數(shù)在處都為零。BW濾波器的另一個特性是在通帶和阻帶內(nèi)的幅頻特性始終是頻率的單調(diào)下降函數(shù),且其模擬函數(shù)隨階次N 的增大而更接近于理想低通濾波器。</p><p> ?。?)、切比雪夫(Chbyshev簡寫為CB)濾波器。CB

20、低通濾波器的模擬函數(shù)由切比雪夫多項式定義,且在通帶內(nèi)的幅頻響應是波動的,在阻帶則單調(diào)變化。</p><p>  3.3設(shè)計IIR濾波器</p><p>  目前IIR數(shù)字濾波器設(shè)計的最通用的方法是借助于模擬濾波器的設(shè)計方法。模擬濾波器設(shè)計已經(jīng)有了一套相當成熟的方法,它不但有完整的設(shè)計公式,而且還有較為完整的圖表供查詢,因此,充分利用這些已有的資源將會給數(shù)字濾波器的設(shè)計帶來很大方便。IIR數(shù)

21、字濾波器的設(shè)計步驟是:</p><p> ?。?)、按一定規(guī)則將給出的數(shù)字濾波器的技術(shù)指標轉(zhuǎn)換為模擬低通濾波器的技術(shù)指標;</p><p> ?。?)、根據(jù)轉(zhuǎn)換后的技術(shù)指標設(shè)計模擬低通濾波器;</p><p> ?。?)、再按一定規(guī)則將G(s)轉(zhuǎn)換成H(z)。</p><p>  若設(shè)計的數(shù)字濾波器是低通的,那么上述設(shè)計工作可以結(jié)束,若所設(shè)計

22、的是高通,帶通或帶阻濾波器,那么還有步驟(4):</p><p>  (4)、將高通、帶通、或帶阻數(shù)字l不去的技術(shù)指標先轉(zhuǎn)化為低通模擬濾波器的技術(shù)指標,然后按照上述步驟(2)設(shè)計出低通,再將轉(zhuǎn)換為所需的H(z)。</p><p>  利用模擬濾波器設(shè)計IIR數(shù)字低通濾波器的步驟:</p><p> ?。?)確定數(shù)字低通濾波器的技術(shù)指標:通帶邊界頻率、通帶最大衰減,阻

23、帶截止頻率、阻帶最小衰減。</p><p>  (2)將數(shù)字低通濾波器的技術(shù)指標轉(zhuǎn)換成相應的模擬低通濾波器的技術(shù)指標。</p><p>  (3)按照模擬低通濾波器的技術(shù)指標設(shè)計及過渡模擬低通濾波器。</p><p>  (4)用雙線性變換法,模擬濾波器系統(tǒng)函數(shù)轉(zhuǎn)換成數(shù)字低通濾波器系統(tǒng)函數(shù)。</p><p><b>  程序如下:&

24、lt;/b></p><p><b>  Ft=8000;</b></p><p><b>  Fp=1000;</b></p><p><b>  Fs=1200;</b></p><p>  wp=2*pi*Fp/Ft;</p><p>  w

25、s=2*pi*Fs/Ft;</p><p>  fp=2*Ft*tan(wp/2);</p><p>  fs=2*Fs*tan(wp/2);</p><p>  [n11,wn11]=buttord(wp,ws,1,50,’s’); </p><p>  [b11,a11]=butter(n11,wn11,’s’); </p>

26、<p>  [num11,den11]=bilinear(b11,a11,0.5);</p><p>  [h,w]=freqz(num11,den11);</p><p><b>  figure;</b></p><p>  plot(w*8000*0.5/pi,abs(h));</p><p>  le

27、gend(‘IIR低通濾波器’,’Location’,’NorthWest’);</p><p><b>  grid;</b></p><p><b>  程序結(jié)果如下圖:</b></p><p>  3.2驗證所設(shè)計的濾波器:</p><p>  為了驗證濾波器的可使用性,我們用常用的sin函數(shù)

28、來進行驗證。其具體程序及運行結(jié)果如下:</p><p>  t=[0:1/1023:1];</p><p>  s=sin(2*pi*t);</p><p>  N=length(s);</p><p>  y=s+0.5*rand(1,N);</p><p>  subplot(2,1,1);</p>

29、<p><b>  plot(y);</b></p><p>  title('加噪語音信號的時域波形','fontweight','bold');</p><p>  S=fft(y); </p><p>  subplot(2,1,2);</p>

30、<p>  plot(abs(S));</p><p>  title('加噪語音信號的頻域波形','fontweight','bold');</p><p><b>  Ft=8000;</b></p><p><b>  Fp=1000;</b></p&g

31、t;<p><b>  Fs=1200;</b></p><p>  wp=2*pi*Fp/Ft;</p><p>  ws=2*pi*Fs/Ft;</p><p>  [n11,wn11]=buttord(wp,ws,1,50,'s'); %求低通濾波器的階數(shù)和截止頻率</p><p> 

32、 [b11,a11]=butter(n11,wn11,'s'); %求S域的頻率響應的參數(shù) </p><p>  [num11,den11]=bilinear(b11,a11,0.5); %利用雙線性變換實現(xiàn)頻率響應S域到Z域的變換 </p><p>  z11=filter(num11,den11,s);</p><p>  sound(z1

33、1);</p><p>  m11=fft(z11); %求濾波后的信號</p><p><b>  figure;</b></p><p>  subplot(2,2,1);</p><p>  plot(abs(S),'g');</p><p>  title('濾波

34、前信號的頻譜','fontweight','bold');</p><p><b>  grid;</b></p><p>  subplot(2,2,2);</p><p>  plot(abs(m11),'r');</p><p>  title('濾波

35、后信號的頻譜','fontweight','bold');</p><p><b>  grid;</b></p><p>  subplot(2,2,3);</p><p><b>  plot(y);</b></p><p>  title('濾波

36、前信號的波形','fontweight','bold');</p><p><b>  grid;</b></p><p>  subplot(2,2,4);</p><p>  plot(z11);</p><p>  title('濾波后的信號波形','

37、fontweight','bold');</p><p><b>  grid;</b></p><p><b>  程序結(jié)果如下圖:</b></p><p>  由所得結(jié)果可知,所設(shè)計的濾波器符合要求。</p><p><b>  4 濾波</b>&l

38、t;/p><p>  用設(shè)計好的IIR低通濾波器對加噪的語音信號進行濾波,程序如下:</p><p>  [y,fs,nbits]=wavread (‘OriSound’); %IIR低通</p><p>  n = length (y) ; %求出語音信號的長度</p><p>  Noise=0.2*randn(n,2);

39、 %隨機函數(shù)產(chǎn)生噪聲</p><p>  s=y+Noise; %語音信號加入噪聲</p><p>  S=fft(s); </p><p><b>  Ft=8000;</b></p><p><b>  Fp=1000;</b></p>

40、;<p><b>  Fs=1200;</b></p><p>  wp=2*pi*Fp/Ft;</p><p>  ws=2*pi*Fs/Ft;</p><p>  [n11,wn11]=buttord(wp,ws,1,50,’s’);%求低通濾波器的階數(shù)和截止頻率</p><p>  [b11,a11]

41、=butter(n11,wn11,’s’); %求S域的頻率響應的參數(shù) </p><p>  [num11,den11]=bilinear(b11,a11,0.5); %利用雙線性變換實現(xiàn)頻率響應S域到Z域的變換 </p><p>  z11=filter(num11,den11,s);</p><p>  sound(z11);</p><

42、;p>  m11=fft(z11); %求濾波后的信號</p><p><b>  figure;</b></p><p>  subplot(2,2,1);</p><p>  plot(abs(S),’g’);</p><p>  title(‘濾波前信號的頻譜’,’fontweight’,’bold’);&

43、lt;/p><p>  axis([ 0 150000 0 4000]);</p><p><b>  grid;</b></p><p>  subplot(2,2,2);</p><p>  plot(abs(m11),’r’);</p><p>  title(‘濾波后信號的頻譜’,’fontw

44、eight’,’bold’);</p><p>  axis([ 0 150000 0 4000]);</p><p><b>  grid;</b></p><p>  subplot(2,2,3);</p><p><b>  plot(s);</b></p><p> 

45、 title(‘濾波前信號的波形’,’fontweight’,’bold’);</p><p>  axis([95000 100000 -1 1]);</p><p><b>  grid;</b></p><p>  subplot(2,2,4);</p><p>  plot(z11);</p>&l

46、t;p>  title(‘濾波后的信號波形’,’fontweight’,’bold’);</p><p>  axis([95000 100000 -1 1]);</p><p><b>  grid;</b></p><p><b>  程序結(jié)果如下圖:</b></p><p>  經(jīng)過以上

47、的加噪處理后,可在Matlab中用函數(shù)sound對聲音進行回放。其調(diào)用格式:sound(y,Fs),sound(y)和sound(y,Fs,bits)。可以察覺濾波前后的聲音有明顯的變化。</p><p><b>  5.創(chuàng)建GUI界面</b></p><p><b>  總結(jié)</b></p><p>  本次設(shè)計通過選擇

48、一個語音信號作為分析的對象,對其進行頻譜分析;利用MATLAB中的隨機函數(shù)產(chǎn)生噪聲加入到語音信號中,模仿被噪聲干擾的語音信號,并對其進行頻譜分析;運用數(shù)字信號處理理論設(shè)計IIR數(shù)字濾波器, 并對被噪聲污染的語音信號進行濾波,分析濾波后信號的時域和頻域特征,回放語音信號。最后,設(shè)計一個信號處理系統(tǒng)界面返回。</p><p>  我們四個人共同完成這個課程設(shè)計,分工是很重要的,首先我們確定了大致的方向,我們要做低通,

49、帶通以及高通三種濾波器,此外還有界面的設(shè)計,程序的整合,以及課程報告的編寫等方面的內(nèi)容,我們的具體分工如下:</p><p>  聲音信號的讀入以及ppt的制作:翁淦泉1062310209</p><p>  低通,帶通,高通濾波器的設(shè)計及相關(guān)程序的編寫:馮竹青1062310202,賀靜文1062301203</p><p>  最后整合及界面的設(shè)計:王雅青1062

50、310208</p><p><b>  課程報告:賀靜文</b></p><p><b>  參考文獻</b></p><p>  1、張志涌 楊祖櫻 MATLAB教程R2011a 北京航空航天大學出版社</p><p>  2、高西全,丁玉美.數(shù)字信號處理(第三版)[M].西安:西安電子科技

51、大學出版社,2008.8</p><p>  3、蔣瑞艷 振動,噪聲處理系統(tǒng)[J]. 大連理工大學學院學報(自然科學版) 2002</p><p>  4、馬超  高世倫 基于MATLAB的噪聲信號采集與分析系統(tǒng)研究[J]. 華中科技大學能源與動力工程學院院報 2004.5</p><p>  5、李正周.MATLAB數(shù)字信號處理與應用[M].北京:清華

52、大學出版社,2008</p><p>  6、高萍,祖靜.基于MATLAB小波去噪技術(shù)淺析[J].科技信息(學術(shù)版),2006(11)1-3.</p><p>  7、劉智 基于MATLAB的機械噪聲分析與處理[J]. 華中江師范學院學報(自然科學版) 2010(4)</p><p><b>  附錄</b></p><p&

53、gt;  附錄(I) 設(shè)計IIR數(shù)字濾波器</p><p>  %=========================IIR低通濾波器=======================</p><p><b>  Ft=8000;</b></p><p><b>  Fp=1000;</b></p><p&g

54、t;<b>  Fs=1200;</b></p><p>  wp=2*pi*Fp/Ft;</p><p>  ws=2*pi*Fs/Ft;</p><p>  fp=2*Ft*tan(wp/2);</p><p>  fs=2*Fs*tan(wp/2);</p><p>  [n11,wn11]

55、=buttord(wp,ws,1,50,’s’); </p><p>  [b11,a11]=butter(n11,wn11,’s’); </p><p>  [num11,den11]=bilinear(b11,a11,0.5);</p><p>  [h,w]=freqz(num11,den11);</p><p><b>  

56、figure;</b></p><p>  plot(w*8000*0.5/pi,abs(h));</p><p>  legend(‘IIR低通濾波器’,’Location’,’NorthWest’);</p><p><b>  grid;</b></p><p><b>  程序結(jié)果如下圖:&

57、lt;/b></p><p>  %======================IIR高通濾波器========================</p><p><b>  Ft=8000;</b></p><p><b>  Fp=4000;</b></p><p><b>  F

58、s=3500;</b></p><p>  wp1=tan(pi*Fp/Ft);</p><p>  ws1=tan(pi*Fs/Ft);</p><p><b>  wp=1;</b></p><p>  ws=wp1*wp/ws1;</p><p>  [n13,wn13]=che

59、b1ord(wp,ws,1,50,’s’);</p><p>  [b13,a13]=cheby1(n13,1,wn13,’s’); </p><p>  [num,den]=lp2hp(b13,a13,wn13);</p><p>  [num13,den13]=bilinear(num,den,0.5); </p><p>  [h,w]

60、=freqz(num13,den13);</p><p><b>  figure;</b></p><p>  plot(w*21000*0.5/pi,abs(h));</p><p>  legend(‘IIR高通濾波器’,’Location’,’NorthWest’);</p><p>  axis([0 1100

61、0 0 1.5]);</p><p><b>  grid;</b></p><p><b>  程序結(jié)果如下圖:</b></p><p>  %======================IIR帶通濾波器==========================</p><p><b>  F

62、p1=1200;</b></p><p><b>  Fp2=3000;</b></p><p><b>  Fs1=1000;</b></p><p><b>  Fs2=3200;</b></p><p><b>  Ft=8000;</b>

63、;</p><p>  wp1=tan(pi*Fp1/Ft);</p><p>  wp2=tan(pi*Fp2/Ft);</p><p>  ws1=tan(pi*Fs1/Ft);</p><p>  ws2=tan(pi*Fs2/Ft);</p><p>  w=wp1*wp2/ws2;</p>&l

64、t;p>  bw=wp2-wp1; %有效通帶頻率</p><p><b>  wp=1;</b></p><p>  ws=(wp1*wp2-w.^2)/(bw*w);</p><p>  [n12,wn12]=buttord(wp,ws,1,50,’s’); </p><p>  [

65、b12,a12]=butter(n12,wn12,’s’);</p><p>  [num2,den2]=lp2bp(b12,a12,sqrt(wp1*wp2),bw);</p><p>  [num12,den12]=bilinear(num2,den2,0.5);</p><p>  [h,w]=freqz(num12,den12);</p>&l

66、t;p><b>  figure;</b></p><p>  plot(w*8000*0.5/pi,abs(h));</p><p>  axis([0 4500 0 1.5]);</p><p>  legend(‘IIR帶通濾波器’,’Location’,’NorthWest’);</p><p><b

67、>  grid;</b></p><p><b>  程序結(jié)果如下圖:</b></p><p>  附錄(II)比較濾波前后語音信號的波形及頻譜</p><p>  % ======================雙線性變換法=======================</p><p>  %****

68、*********************低通濾波器************************</p><p>  [y,fs,nbits]=wavread (‘OriSound’); %IIR低通</p><p>  n = length (y) ; %求出語音信號的長度</p><p>  Noise=0.2*randn(n,2);

69、%隨機函數(shù)產(chǎn)生噪聲</p><p>  s=y+Noise; %語音信號加入噪聲</p><p>  S=fft(s); </p><p><b>  Ft=8000;</b></p><p><b>  Fp=1000;</b></p>

70、<p><b>  Fs=1200;</b></p><p>  wp=2*pi*Fp/Ft;</p><p>  ws=2*pi*Fs/Ft;</p><p>  [n11,wn11]=buttord(wp,ws,1,50,’s’);%求低通濾波器的階數(shù)和截止頻率</p><p>  [b11,a11]=b

71、utter(n11,wn11,’s’); %求S域的頻率響應的參數(shù) </p><p>  [num11,den11]=bilinear(b11,a11,0.5); %利用雙線性變換實現(xiàn)頻率響應S域到Z域的變換 </p><p>  z11=filter(num11,den11,s);</p><p>  sound(z11);</p><p

72、>  m11=fft(z11); %求濾波后的信號</p><p><b>  figure;</b></p><p>  subplot(2,2,1);</p><p>  plot(abs(S),’g’);</p><p>  title(‘濾波前信號的頻譜’,’fontweight’,’bold’);<

73、;/p><p>  axis([ 0 150000 0 4000]);</p><p><b>  grid;</b></p><p>  subplot(2,2,2);</p><p>  plot(abs(m11),’r’);</p><p>  title(‘濾波后信號的頻譜’,’fontwei

74、ght’,’bold’);</p><p>  axis([ 0 150000 0 4000]);</p><p><b>  grid;</b></p><p>  subplot(2,2,3);</p><p><b>  plot(s);</b></p><p>  t

75、itle(‘濾波前信號的波形’,’fontweight’,’bold’);</p><p>  axis([95000 100000 -1 1]);</p><p><b>  grid;</b></p><p>  subplot(2,2,4);</p><p>  plot(z11);</p><

76、p>  title(‘濾波后的信號波形’,’fontweight’,’bold’);</p><p>  axis([95000 100000 -1 1]);</p><p><b>  grid;</b></p><p><b>  程序結(jié)果如下圖:</b></p><p>  附II-1

77、雙線性低通濾波器比較</p><p>  %**********************高通濾波器*****************************</p><p>  [y,fs,nbits]=wavread (‘OriSound’); %IIR高通</p><p>  n = length (y) ; %求出語音信號的長度</p>&l

78、t;p>  Noise=0.2*randn(n,2); %隨機函數(shù)產(chǎn)生噪聲</p><p>  s=y+Noise; %語音信號加入噪聲</p><p>  S=fft(s); %傅里葉變換</p><p><b>  Fp1=1200;</b></p><p><b> 

79、 Fs1=1000;</b></p><p><b>  Ft=8000;</b></p><p>  wp1=tan(pi*Fp1/Ft);</p><p>  ws1=tan(pi*Fs1/Ft);</p><p><b>  wp=1;</b></p><p&g

80、t;  ws=wp1*wp/ws1;</p><p>  [n13,wn13]=cheb1ord(wp,ws,1,50,’s’); %求模擬的低通濾波器階數(shù)和截止頻率</p><p>  [b13,a13]=cheby1(n13,1,wn13,’s’); %求S域的頻率響應的參數(shù)</p><p>  [num,den]=lp2hp(b13,a13,wn13);

81、 %將S域低通參數(shù)轉(zhuǎn)為高通的</p><p>  [num13,den13]=bilinear(num,den,0.5); %利用雙線性變換實現(xiàn)頻率響應S域到Z域轉(zhuǎn)換</p><p>  z13=filter(num13,den13,s);</p><p>  sound(z13);</p><p>  m13=fft(z13); %

82、求濾波后的信號</p><p><b>  figure;</b></p><p>  subplot(2,2,1);</p><p>  plot(abs(S),’g’);</p><p>  title(‘濾波前信號的頻譜’,’fontweight’,’bold’);</p><p>  ax

83、is([0 150000 0 4000]);</p><p><b>  grid;</b></p><p>  subplot(2,2,2);</p><p>  plot(abs(m13),’r’);</p><p>  title(‘濾波后信號的頻譜’,’fontweight’,’bold’);</p>

84、;<p>  axis([0 150000 0 4000]);</p><p><b>  grid;</b></p><p>  subplot(2,2,3);</p><p><b>  plot(s);</b></p><p>  title(‘濾波前信號的波形’,’fontw

85、eight’,’bold’);</p><p>  axis([95000 100000 -1 1]);</p><p><b>  grid;</b></p><p>  subplot(2,2,4);</p><p>  plot(z13);</p><p>  title(‘濾波后的信號波形

86、’,’fontweight’,’bold’);</p><p>  axis([95000 100000 -1 1]);</p><p><b>  grid;</b></p><p><b>  程序結(jié)果如下圖:</b></p><p>  附II-2 雙線性高通濾波器比較</p>

87、<p>  %**********************帶通濾波器*****************************</p><p>  [y,fs,nbits]=wavread (‘OriSound’); %IIR帶通</p><p>  n = length (y) ; %求出語音信號的長度</p><p>  Noise=0.

88、2*randn(n,2); %隨機函數(shù)產(chǎn)生噪聲</p><p>  s=y+Noise; %語音信號加入噪聲 </p><p>  S=fft(s); %傅里葉變換</p><p><b>  Ft=8000;</b></p><p><b>  Fp=1000;

89、</b></p><p><b>  Fs=1200;</b></p><p>  wp=2*Fp/Ft;</p><p>  ws=2*Fs/Ft;</p><p><b>  rp=1;</b></p><p><b>  rs=50;</b&

90、gt;</p><p>  p=1-10.^(-rp/20); %通帶阻帶波紋</p><p>  q=10.^(-rs/20);</p><p>  fpts=[wp ws];</p><p>  mag=[1 0];</p><p>  dev=[p q];</p><p>  [n

91、21,wn21,beta,ftype]=kaiserord(fpts,mag,dev);%由kaiserord求濾波器的階數(shù)和截止頻率</p><p>  b21=fir1(n21,wn21,kaiser(n21+1,beta)); %由fir1設(shè)計濾波器</p><p>  z21=fftfilt(b21,s);</p><p>  sound(z21);&l

92、t;/p><p>  m21=fft(z21); %求濾波后的信號</p><p>  figure(4);</p><p>  subplot(2,2,1);</p><p>  plot(abs(S),’g’);</p><p>  title(‘濾波前信號的頻譜’,’fontweight’,’bold’);<

93、;/p><p>  axis([0 150000 0 4000]);</p><p><b>  grid;</b></p><p>  subplot(2,2,2);</p><p>  plot(abs(m21),’r’);</p><p>  title(‘濾波后信號的頻譜’,’fontwei

94、ght’,’bold’);</p><p>  axis([0 150000 0 4000]);</p><p><b>  grid;</b></p><p>  subplot(2,2,3);</p><p><b>  plot(s);</b></p><p>  t

95、itle(‘濾波前信號的波形’,’fontweight’,’bold’);</p><p>  axis([95000 100000 -1 1]);</p><p><b>  grid;</b></p><p>  subplot(2,2,4);</p><p>  plot(z21);</p><

96、p>  title(‘濾波后的信號波形’,’fontweight’,’bold’);</p><p>  axis([95000 100000 -1 1]);</p><p><b>  grid;</b></p><p><b>  程序結(jié)果如下圖:</b></p><p>  附II-3

溫馨提示

  • 1. 本站所有資源如無特殊說明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請下載最新的WinRAR軟件解壓。
  • 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
  • 3. 本站RAR壓縮包中若帶圖紙,網(wǎng)頁內(nèi)容里面會有圖紙預覽,若沒有圖紙預覽就沒有圖紙。
  • 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
  • 5. 眾賞文庫僅提供信息存儲空間,僅對用戶上傳內(nèi)容的表現(xiàn)方式做保護處理,對用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對任何下載內(nèi)容負責。
  • 6. 下載文件中如有侵權(quán)或不適當內(nèi)容,請與我們聯(lián)系,我們立即糾正。
  • 7. 本站不保證下載資源的準確性、安全性和完整性, 同時也不承擔用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。

評論

0/150

提交評論