matlab課程設(shè)計(jì)報(bào)告--基于matlab有噪聲語音信號處理

上傳人：奔*** IP屬地：河北更新時(shí)間：2024-03-01 格式：doc 頁數(shù)：24 大?。?.27MB 人氣指數(shù)：12 舉報(bào) 版權(quán)申訴

matlab課程設(shè)計(jì)報(bào)告--基于matlab有噪聲語音信號處理_第1頁

已閱讀1頁，還剩23頁未讀，繼續(xù)免費(fèi)閱讀

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

文檔簡介

1、　　Matlab課程設(shè)計(jì)報(bào)告　　題目:基于MATLAB有噪聲語音信號處理　　系（院）：計(jì)算機(jī)與信息工程學(xué)院 　　專業(yè)：通信工程 　　班級：

2、 　　指導(dǎo)教師： 　　學(xué)年學(xué)期: 2011 ~ 2012 學(xué)年第 2 學(xué)期　　簡介：　　我們通信工程專業(yè)在實(shí)踐中經(jīng)常碰到需要對已接收信號進(jìn)行處理的情況

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

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

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

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

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

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

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

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

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

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

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

14、　　s=y+Noise;　　sound(s);　　figure;　　subplot(2,1,1);　　plot(s);

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

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

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

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

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

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

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

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

23、帶截止頻率、阻帶最小衰減。　　（2）將數(shù)字低通濾波器的技術(shù)指標(biāo)轉(zhuǎn)換成相應(yīng)的模擬低通濾波器的技術(shù)指標(biāo)。　?。?）按照模擬低通濾波器的技術(shù)指標(biāo)設(shè)計(jì)及過渡模擬低通濾波器。　　（4）用雙線性變換法，模擬濾波器系統(tǒng)函數(shù)轉(zhuǎn)換成數(shù)字低通濾波器系統(tǒng)函數(shù)。　　程序如下：&

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

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

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

27、gend(‘IIR低通濾波器’,’Location’,’NorthWest’);　　grid;　　程序結(jié)果如下圖：　　3.2驗(yàn)證所設(shè)計(jì)的濾波器：　　為了驗(yàn)證濾波器的可使用性，我們用常用的sin函數(shù)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

47、的加噪處理后，可在Matlab中用函數(shù)sound對聲音進(jìn)行回放。其調(diào)用格式：sound(y,Fs)，sound(y)和sound(y,Fs,bits)?？梢圆煊X濾波前后的聲音有明顯的變化。　　5.創(chuàng)建GUI界面　　總結(jié)　　本次設(shè)計(jì)通過選擇

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

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

50、310208　　課程報(bào)告：賀靜文　　參考文獻(xiàn)　　1、張志涌楊祖櫻 MATLAB教程R2011a 北京航空航天大學(xué)出版社　　2、高西全，丁玉美.數(shù)字信號處理（第三版）[M].西安：西安電子科技

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

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

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

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

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

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

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

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

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

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

61、0 0 1.5]);　　grid;　　程序結(jié)果如下圖：　　%======================IIR帶通濾波器==========================　　F

62、p1=1200;　　Fp2=3000;　　Fs1=1000;　　Fs2=3200;　　Ft=8000;

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

眾賞文庫> 全部分類> 畢業(yè)設(shè)計(jì)

溫馨提示

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

matlab課程設(shè)計(jì)報(bào)告--基于matlab有噪聲語音信號處理

文檔簡介

溫馨提示

最新文檔

評論

matlab課程設(shè)計(jì)報(bào)告--基于matlab有噪聲語音信號處理

文檔簡介

溫馨提示

最新文檔

評論

免費(fèi)下載