




版權(quán)說明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請(qǐng)進(jìn)行舉報(bào)或認(rèn)領(lǐng)
文檔簡介
1、語音信號(hào)的數(shù)字濾波FIR數(shù)字濾波器的(漢寧)窗函數(shù)法設(shè)計(jì)設(shè)計(jì)題目:語音信號(hào)的數(shù)字濾波FIR數(shù)字濾波器的(漢寧)窗函數(shù)法設(shè)計(jì)一、課程設(shè)計(jì)的目的通過對(duì)常用數(shù)字濾波器的設(shè)計(jì)和實(shí)現(xiàn),掌握數(shù)字信號(hào)處理的工作原理及設(shè)計(jì)方法;掌握利用數(shù)字濾波器對(duì)信號(hào)進(jìn)行濾波的方法。并能夠?qū)υO(shè)計(jì)結(jié)果加以分析。二、設(shè)計(jì)步驟2.1窗函數(shù)設(shè)計(jì)法的原理窗函數(shù)的基本思想:先選取一個(gè)理想濾波器(它的單位抽樣響應(yīng)是非因果、無限長的),再截?。ɑ蚣哟埃┧膯挝怀闃禹憫?yīng)得到線性相位因果FIR濾波器。這種方法的重點(diǎn)是選擇一個(gè)合適的窗函數(shù)和理想濾波器。設(shè)x(n)是一個(gè)長序列,是長度為N的窗函數(shù),用截?cái)啵玫絅點(diǎn)序列,即 在頻域上則有由此可見,窗函
2、數(shù)不僅僅會(huì)影響原信號(hào)在時(shí)域上的波形,而且也會(huì)影響到頻域內(nèi)的形狀。2.2漢寧窗函數(shù)簡介漢寧窗、海明窗和布萊克曼窗,都可以用一種通用的形式表示,這就是廣義余弦窗。這些窗都是廣義余弦窗的特例,漢寧(Hanning)窗又稱升余弦窗,漢寧窗可以看作是3個(gè)矩形時(shí)間窗的頻譜之和,或者說是 3個(gè) sinc(t)型函數(shù)之和,而括號(hào)中的兩項(xiàng)相對(duì)于第一個(gè)譜窗向左、右各移動(dòng)了/T,從而使旁瓣互相抵消,消去高頻干擾和漏能。適用于非周期性的連續(xù)信號(hào)。公式如下:2.3進(jìn)行語音信號(hào)的采集(1)按“開始”“程序”“附件”“娛樂”“錄音機(jī)”的順序操作打開Windows系統(tǒng)中的錄音機(jī)軟件。如圖1所示。圖1 windows錄音機(jī)(2
3、)用麥克風(fēng)錄入自己的聲音信號(hào)并保存成wav文件。如圖2所示。圖2 保存文件保存的文件按照要求如下: 音信號(hào)文件保存的文件名為“yuxuejiao.wav”。語音信號(hào)的屬性為“8.000KHz,8位,單聲道 7KB/秒” ,其它選項(xiàng)為默認(rèn)。2.4語音信號(hào)的分析將“yuxuejiao.wav”語音文件復(fù)制到計(jì)算機(jī)裝有Matlab軟件的磁盤中相應(yīng)Matlab目錄中的“work”文件夾中。打開Matlab軟件,在菜單欄中選擇“File”“new”將代碼復(fù)制到空白處,并保存到“work”文件夾中,和聲音文件在同一個(gè)文件夾內(nèi)。2.4.1語音信號(hào)的截取處理和頻譜分析編寫MATLAB編碼實(shí)現(xiàn)對(duì)信號(hào)的截取處理,
4、也就是對(duì)yf.wav語音的截取處理,截取范圍為17000-57000。代碼如下:s, fs=wavread(yuxuejiao.wav);s1=s(17000:57000);figure(1);subplot(211)plot(s)title(原始語音信號(hào))subplot(212)plot(s1)title(截短語音信號(hào));wavwrite(s1,fs,s1.wav);S1=fft(s1);figure(2)subplot(311);plot(s1);title(截短預(yù)處理語音信號(hào))subplot(312)plot(abs(S1)title(預(yù)處理語音信號(hào)頻譜);subplot(313);k=
5、0:39999;plot(k(1:20000)*1,abs(S1(1:20000);title(預(yù)處理語音信號(hào)單邊帶頻譜)在m文件編輯器中輸入相應(yīng)的指令將自己的語音信號(hào)導(dǎo)入Matlab工作臺(tái),點(diǎn)擊“run”或者 “F5”運(yùn)行文件。效果如圖3和圖4圖3語音信號(hào)的截取處理圖在圖3中, 其中第一個(gè)圖為原始語音信號(hào);第二個(gè)圖是截短后的信號(hào)圖。圖4頻譜分析圖其中第二個(gè)圖是信號(hào)的FFT結(jié)果,其橫坐標(biāo)的具體值是X(k)中的序號(hào)k;第三個(gè)圖是確定濾波頻率范圍的參考圖,其橫坐標(biāo)的具體值應(yīng)當(dāng)是遵循DFT定義式和頻率分辨率求得的:當(dāng)k等于0時(shí), ,從數(shù)字角頻率上看,對(duì)應(yīng)的正好是即直流的位置,也就是說,在取濾波頻段時(shí)
6、,當(dāng)將主要能量(即紅色框的部分)保留,其余頻段部分的信號(hào)濾除。相當(dāng)于是信號(hào)的實(shí)際頻譜采樣,而又是連續(xù)時(shí)間語音信號(hào)的采樣。的每兩個(gè)相鄰取值之間的頻率間隔大小對(duì)應(yīng)到語音信號(hào)的頻譜中去,其頻率間隔大小正好是稱頻率分辨率,其中, 由于在第三個(gè)圖中是遵循DFT定義式和頻率分辨率求得,表示的是濾波頻率的范圍2.5濾波器的設(shè)計(jì)2.5.1濾波器理論參數(shù)設(shè)定(1) 從圖4可知我們已經(jīng)確定了通帶截止頻率和阻帶截止頻率和。 由于hanning窗過渡帶滿足 (2)給定所要求的頻率響應(yīng)函數(shù) (3)求單位采樣響應(yīng)(4)計(jì)算線性相位延遲,hanning窗=(5)濾波器的單位采樣響應(yīng):2.5.2濾波器的MATLAB仿真 在M
7、文件中繼續(xù)編寫代碼,把計(jì)算出來的參數(shù)帶入代碼中。代碼如下:s2=awgn(s1,30); %完成加噪!30dbwavwrite(s2,fs,s2.wav);figure(3);subplot(211);plot(s2);title(加噪后語音信號(hào));subplot(212);S2=fft(s2);plot(abs(S2);title(加噪后信號(hào)頻譜);figure(4)subplot(211);plot(s1);title(語音信號(hào));subplot(212);plot(s2);title(加噪后語音信號(hào));%濾波器完成相關(guān)參數(shù)配置wp=900*2*pi/8000;wst=1000*2*pi/
8、8000;wc=(wp+wst)/2;N=ceil(3.1*2*pi/(wst-wp);r=(N-1)/2;hn1=fir1(N-1,wc/pi,low,hanning(N);%s3=conv(s2,hn1);wavwrite(s3,fs,s3.wav);S3=fft(s3);figure(5)freqz(hn1);title(濾波器幅頻特性與相頻特性)figure(6)subplot(111)stem(hn1);title(濾波器單位采樣響應(yīng));figure(7)subplot(211)plot(s3)title(濾波器處理之后信號(hào)圖)subplot(212);plot(abs(S3); t
9、itle(濾波器處理之后頻譜);figure(8)subplot(211)plot(s2);title(加噪后語音信號(hào));subplot(212);plot(s3);title(濾波器處理之后信號(hào)圖);s4=conv(s1,hn1);p1=sum(s1.2);p2=sum(s2.2)-sum(s1.2);SNR1=10*log10(p1/p2);p3=sum(s4.2)/8000;p4=sum(s3.2)/8000-sum(s4.2)/8000;SNR2=10*log10(p3/p4);2.6 噪聲疊加圖5 語音信號(hào)與加噪聲后語音信號(hào)對(duì)比圖五為語音信號(hào)與加噪聲后語音信號(hào)對(duì)。比計(jì)算機(jī)隨即產(chǎn)生的噪
10、聲指令為:awgn() 所加的噪聲為30 。圖6 加噪后語音信號(hào)與加噪后信號(hào)頻譜圖圖7濾波器頻幅特性與相頻特性圖8濾波器單位采樣響應(yīng) 圖7是濾波器頻幅特性與相頻特性,圖8是hanning窗函數(shù)單位沖激響應(yīng)系數(shù)。圖9 濾波器處理后的信號(hào)及頻譜圖圖10 加噪后的語音信號(hào)及處理后的信號(hào)圖設(shè)計(jì)的濾波器是用單位采樣響應(yīng)h(n)表示的,可以利用帶噪聲語音與h(n)做時(shí)域卷積,即: 。在Matlab中,卷積運(yùn)算可以用函數(shù)“conv( )”實(shí)現(xiàn)。濾波前(含噪聲)的信號(hào)和濾波后信號(hào)的信噪比, 利用sum求信號(hào)的功率。在m文件中繼續(xù)編寫信噪比代碼段:s4=conv(s1,hn1);p1=sum(s1.2);p2=
11、sum(s2.2)-sum(s1.2);SNR1=10*log10(p1/p2);p3=sum(s4.2)/8000;p4=sum(s3.2)/8000-sum(s4.2)/8000;SNR2=10*log10(p3/p4);三、設(shè)計(jì)結(jié)果與分析濾波就是將外界干擾的不需要的頻率濾除掉,一面影響要測試數(shù)據(jù)的結(jié)果。采樣間隔也會(huì)對(duì)實(shí)驗(yàn)產(chǎn)生影響,取的過大會(huì)導(dǎo)致頻率混疊。3.1比較濾波前(含噪聲信號(hào)的文件)和濾波后的語音信號(hào)效果。圖10 利用hamming窗FIR低通濾波器濾波效果圖圖10為濾波前后的效果圖,第一個(gè)圖是加噪后的圖,圖中波形的毛刺部分(即干擾噪聲)被濾除,輸入的帶噪聲信號(hào)經(jīng)過濾波器濾波后,即
12、第二個(gè)圖中所示波形。圖中有明顯的濾波效果,濾波器不僅把30的噪聲濾掉了,同時(shí)也也能濾掉了原信號(hào)的聲音,因此在實(shí)際的實(shí)驗(yàn)會(huì)和理論有誤差,許多時(shí)候得不到完美的效果。3.2濾波前(含噪聲)的信號(hào)和濾波后信號(hào)的信噪比圖11相關(guān)參數(shù)的截圖,其中snr為信噪比圖11是利用“信噪比代碼段”做出來的效果圖,圖中名稱snr1為最濾波前的信噪比,snr2為濾波后的信噪比,snr2大于snr1,最終達(dá)到了效果。四、設(shè)計(jì)體會(huì)在本次課程設(shè)計(jì)中,讓我更加了解了數(shù)字信號(hào)處理在實(shí)際上的應(yīng)用,課程設(shè)計(jì)不僅要求對(duì)濾波器理論的研究,更重要的是實(shí)際設(shè)計(jì)中遇到的問題。因?yàn)橛辛诉@次課程設(shè)計(jì),我不得不對(duì)其設(shè)計(jì)原理進(jìn)行更深一層次的理解,對(duì)書
13、中原來學(xué)到的理論,僅知道了其表面,而不知其原因。在設(shè)計(jì)中也使我對(duì)一些概念有了更深刻的認(rèn)識(shí)。除此之外,對(duì)程序的編譯不是一蹴而就的,而是經(jīng)過多次的編譯與調(diào)試。以前用MATLAB就是簡單的輸入,可以說都不是自己的勞動(dòng)成果,但這次不一樣,課程設(shè)計(jì)沒有別人給你編好,而是自己寫,出錯(cuò)率就大大提高了,但這是過程,學(xué)習(xí)就是在過程中進(jìn)行的,經(jīng)過自己幾天的腦力勞動(dòng),再加上同學(xué)們的幫助,不僅對(duì)讀程序有了很大提高,而且自己的編譯水平也上了一個(gè)新臺(tái)階,更加熟系了MATLAB的應(yīng)用,也對(duì)其中的函數(shù)有了大概的了解,對(duì)其中一些函數(shù)也相當(dāng)熟練。完成整個(gè)設(shè)計(jì)過程后,學(xué)到的東西已經(jīng)不僅僅上面的那些東西,還有同學(xué)們互相幫助,共同探討和設(shè)計(jì)過程中的每一個(gè)細(xì)節(jié),也許每一個(gè)細(xì)節(jié)的錯(cuò)誤就可能導(dǎo)致結(jié)局的失敗,所以我認(rèn)為這次收獲最大的莫過于靜心,學(xué)習(xí)不能急,一定要冷靜,心無旁騖,不放過任何一個(gè)細(xì)節(jié),就能帶
溫馨提示
- 1. 本站所有資源如無特殊說明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請(qǐng)下載最新的WinRAR軟件解壓。
- 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請(qǐng)聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
- 3. 本站RAR壓縮包中若帶圖紙,網(wǎng)頁內(nèi)容里面會(huì)有圖紙預(yù)覽,若沒有圖紙預(yù)覽就沒有圖紙。
- 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
- 5. 人人文庫網(wǎng)僅提供信息存儲(chǔ)空間,僅對(duì)用戶上傳內(nèi)容的表現(xiàn)方式做保護(hù)處理,對(duì)用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對(duì)任何下載內(nèi)容負(fù)責(zé)。
- 6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容,請(qǐng)與我們聯(lián)系,我們立即糾正。
- 7. 本站不保證下載資源的準(zhǔn)確性、安全性和完整性, 同時(shí)也不承擔(dān)用戶因使用這些下載資源對(duì)自己和他人造成任何形式的傷害或損失。
最新文檔
- 微商線下引流活動(dòng)方案
- 開展聯(lián)賽活動(dòng)方案
- 徐州紅跑車會(huì)員活動(dòng)方案
- 開學(xué)禮品活動(dòng)方案
- 開展社工日宣傳活動(dòng)方案
- 開展養(yǎng)殖活動(dòng)方案
- 形體瘦身活動(dòng)方案
- 心理購買活動(dòng)方案
- 徒步龍?zhí)痘顒?dòng)方案
- 影樓元旦活動(dòng)方案
- 消防系統(tǒng)維護(hù)保養(yǎng)方案
- 骨科護(hù)理實(shí)習(xí)生小講課
- 四川省南充市2023-2024學(xué)年七年級(jí)下學(xué)期期末考試道德與法治試卷(含答案)
- 2025至2030中國汽車散熱器行業(yè)市場發(fā)展分析及商業(yè)模式與投融資發(fā)展報(bào)告
- 統(tǒng)編版語文二下園地三+單元復(fù)習(xí)課 課件
- 2025年輕人情緒消費(fèi)趨勢報(bào)告-抖音商城xsocialbeta-202506
- 培訓(xùn)中心項(xiàng)目管理制度
- 承包企業(yè)食堂管理制度
- 智能合約的自適應(yīng)優(yōu)化與動(dòng)態(tài)執(zhí)行研究-洞察闡釋
- 高中教科研課題:《新課程背景下高中語文情境教學(xué)改革研究》課題工作匯報(bào)
- 2025年寧波市奉化區(qū)紅果文體產(chǎn)業(yè)運(yùn)營管理有限公司招聘筆試參考題庫含答案解析
評(píng)論
0/150
提交評(píng)論