基于matlab有噪聲語音信號(hào)的處理-預(yù)覽頁
【正文】 功能說明：讀取file所規(guī)定的wav文件，返回采樣值放在向量y中。接下來。當(dāng)N小于xn的長度時(shí)，fft函數(shù)計(jì)算xn的前N個(gè)元素，忽略其后面的元素。)。 %傅里葉變換subplot(2,1,1)。)。原始信號(hào)頻譜39。也可直接用randn函數(shù)產(chǎn)生高斯分布序列。語音信號(hào)添加噪聲及其頻譜分析的主要程序如下：[y,fs,nbits]=wavread (39。 %隨機(jī)函數(shù)產(chǎn)生噪聲Si=y+Noise。title(39。 %傅里葉變換subplot(2,1,2)。)。FIR濾波器比鞥采用間接法，常用的方法有窗函數(shù)法、頻率采樣發(fā)和切比雪夫等波紋逼近法。具體設(shè)計(jì)步驟如下：（1）確定所需類型數(shù)字濾波器的技術(shù)指標(biāo)。（6）采用雙線性變換法將相應(yīng)類型的過渡模擬濾波器轉(zhuǎn)換成所需類型的數(shù)字濾波器。（2）將數(shù)字低通濾波器的技術(shù)指標(biāo)轉(zhuǎn)換成相應(yīng)的模擬低通濾波器的技術(shù)指標(biāo)。下面我們著重介紹用窗函數(shù)法設(shè)計(jì)FIR濾波器的步驟。（2）構(gòu)造希望逼近的頻率響應(yīng)函數(shù)。hn=fir1(M，wc，window)，可以指定窗函數(shù)向量window。MATLAB信號(hào)處理工具箱函數(shù)cheblap,cheblord和cheeby1是切比雪夫I型濾波器設(shè)計(jì)函數(shù)。ws=2*pi*Fs/Ft。s39。)。legend(39。 %帶通到低通濾波器的轉(zhuǎn)換wp2=tan(pi*Fp2/Ft)。bw=wp2wp1。s39。)。 %根據(jù)參數(shù)求出頻率響應(yīng)plot(w*8000*,abs(h))。)。wp1=tan(pi*Fp/Ft)。[n13,wn13]=cheb1ord(wp,ws,1,50,39。s39。 %利用雙線性變換實(shí)現(xiàn)S域到Z域轉(zhuǎn)換[h,w]=freqz(num13,den13)。)。圖5 IIR 高通濾波器%**************************FIR低通濾波器*****************************% 用窗函數(shù)設(shè)計(jì)低通濾波器的程序如下rp=1。fpts=[wp ws]。%kaiserord求階數(shù)截止頻率b21=fir1(n21,wn21,Kaiser(n21+1,beta))。FIR低通濾波器39。ws1=tan(pi*Fs1/Ft)。wp=1。)。 %求S域的頻率響應(yīng)的參數(shù)[num2,den2]=lp2bp(b22,a22,sqrt(wp1*wp2),bw)。axis([0 4000 0 ])。圖7 FIR帶通濾波器%***************************%FIR高通濾波器*******************p=110.^(rp/20)。dev=[p s]。,Kaiser(n23+1,beta))。FIR高通濾波器39。函數(shù)fftfilt用的是重疊相加法實(shí)現(xiàn)線性卷積的計(jì)算。其中xn是輸入信號(hào)向量，yn輸出信號(hào)向量。n = length (y) 。 %傅里葉變換z11=filter(num11,den11,s)。plot(abs(S),39。濾波前信號(hào)的頻譜39。plot(abs(m11),39。濾波后信號(hào)的頻譜39。plot(s)。grid。濾波后的信號(hào)波形39。)。 %語音信號(hào)加入噪聲S=fft(s)。 %求濾波后的信號(hào)subplot(2,2,1)。title(39。plot(abs(m12),39。濾波后信號(hào)的頻譜39。title(39。plot(z12)。圖10 雙線性法帶通濾波%*****************************高通濾波器*****************************[y,fs,nbits]=wavread (39。 %求出語音信號(hào)的長度noise=*randn(n,2)。sound(z13)。g39。)。)。subplot(2,2,3)。)。濾波后的信號(hào)波形39。)。 %語音信號(hào)加入噪聲 S=fft(s)。 %求濾波后的信號(hào)subplot(2,2,1)。title(39。plot(abs(m21),39。濾波后信號(hào)的頻譜39。 title(39。plot(z21)。 圖12 窗函數(shù)法低通濾波%****************************帶通濾波器***************************[y,fs,nbits]=wavread (39。 %求出語音信號(hào)的長度noise=*randn(n,2)。sound(z22)。g39。)。)。subplot(2,2,3)。)。濾波后的信號(hào)波形39。)。 %語音信號(hào)加入噪聲 S=fft(s)。 %求濾波后的信號(hào)subplot(2,2,1)。subplot(2,2,2)。title(39。plot(s)。subplot(2,2,4)。)。 系統(tǒng)界面設(shè)計(jì)工具—GUI概述圖形用戶界面（graphical user interfaces ，GUI）則是由窗口、光標(biāo)、按鍵、菜單、文字說明等對(duì)象（objects）構(gòu)成的一個(gè)用戶界面。GUIDE將設(shè)計(jì)好的GUI保存在一個(gè)FIG文件中，同時(shí)生成M文件框架。 GUI創(chuàng)建包括界面設(shè)計(jì)和控件編程兩部分，主要步驟如下。最后，使用說明。也就是說，課題更多的還是體現(xiàn)了數(shù)字信號(hào)處理技術(shù)。通過MATLAB里幾個(gè)命令函數(shù)的調(diào)用，很輕易的在實(shí)際語音與數(shù)字信號(hào)的理論之間搭了一座橋。通過比較加噪前后，語音的頻譜和語音回放，能明顯的感覺到加入噪聲后回放的聲音與原始的語音信號(hào)有很大的不同,前者隨較尖銳的干擾嘯叫聲。由此可見，語音信號(hào)主要分布在低頻段，而噪聲主要分布在高頻段。對(duì)報(bào)告的書寫格式及內(nèi)容，老師多次幫助我分析思路，開拓視角。同時(shí)還要學(xué)感謝我的同，尤其是我們同一課題的幾個(gè)同學(xué)，我們花費(fèi)課很多的時(shí)間和精力。在一周的課程設(shè)計(jì)過程中， 學(xué)院的機(jī)房工作人員給我們提供的便利的條件，天氣寒冷，實(shí)驗(yàn)室空調(diào)一直開放，我們覺得很溫暖，在此，表達(dá)對(duì)工作人員的謝意。 最后，再次感謝所有幫助過我的老師和同學(xué)！參考文獻(xiàn)1高西全，丁玉美．?dāng)?shù)字信號(hào)處理．第3版．北京：西安電子科技大學(xué)出版社，20082 劉泉，闕大順．?dāng)?shù)字信號(hào)處理原理與實(shí)現(xiàn)．北京：電子工業(yè)出版社，20053 張磊，畢靖，郭蓮英．MATLAB實(shí)用教程．北京：人民郵電出版社，20084 張威．MATLAB基礎(chǔ)與編程入門．西安：西安電子科技大學(xué)出版社，20065 周利清，蘇菲．?dāng)?shù)字信號(hào)處理基礎(chǔ)．北京：北京郵電大學(xué)出版社，200525 指導(dǎo)教師評(píng)語學(xué)號(hào)1071302126姓名班級(jí)通信1071選題名稱基于MATLAB有噪聲語音信號(hào)的處理序號(hào)評(píng)價(jià)內(nèi)容權(quán)重（%）得分1考勤記錄、學(xué)習(xí)態(tài)度、工作作風(fēng)與表現(xiàn)。205報(bào)告的格式規(guī)范程度、是否圖文并茂、語言規(guī)范及流暢程度；主題是否鮮明、重心是否突出、論述是否充分、結(jié)論是否正確；是否提出了自己的獨(dú)到見
環(huán)評(píng)公示相關(guān)推薦
鄂ICP備17016276號(hào)-1