畢業(yè)論文--基于dsp的語音信號的去噪方法_第1頁
已閱讀1頁,還剩25頁未讀 繼續(xù)免費(fèi)閱讀

下載本文檔

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

文檔簡介

1、<p>  畢 業(yè) 設(shè) 計(論 文)</p><p>  設(shè)計(論文)題目: 基于DSP的語音信號的去噪方法</p><p>  2013年 5月5日 </p><p>  基于DSP的語音信號的去噪方法</p><p>  摘要:語音是語言信息的載體,是人類交流信息最自然、最有效、最方便的手段。在許多情況下,人們所關(guān)心的

2、語音信號不可避免地被其他信號所污染,影響了人們的聽辨。</p><p>  小波變換在當(dāng)今信號處理領(lǐng)域中是一種十分活躍的理論,小波閾值消減法是對小波分解系數(shù)進(jìn)行閾值化的降噪技術(shù)。本文簡要介紹主要針對基于小波變換的閾值去噪方法進(jìn)行了研究,對小波閾值消噪從兩個方面進(jìn)行改進(jìn)。一方面,通過對傳統(tǒng)小波閾值消噪算法中的軟閾值函數(shù)和硬閾值函數(shù)以及L.Breiman提出的非負(fù)死區(qū)閾值函數(shù)進(jìn)行分析,提出了兩種改進(jìn)的閾值函數(shù)。另一方

3、面,受到譜減法思想的啟發(fā),提出了新的閾值估計思想。經(jīng)實驗證明,本文閾值函數(shù)和閾值估計的消噪結(jié)果保持了信號的邊緣特性,并且提高了語音信噪比。</p><p>  關(guān)鍵詞:小波消噪 閾值估計 閾值函數(shù) 高斯白噪聲</p><p>  Abstract:Speech signal is the carder of the speech information,and also is the mo

4、st nature,convenient and shortcut way of intercommunication.However in most cases,speech signal is disturbed and polluted inevitably so that we can’t catch the explicit meaning of the speech signal.</p><p> 

5、 Wavelet transform theory is one of the modem analysis methods.Wavelet threshold de-noising is a de-noising technique based on the threshold of the wavelet coefficients.In this paper,several approaches of speech enhancem

6、ent are introduced briefly.Wavelet threshold de-noising is the main part in this paper.Two aspects improvement is presented in this paper ,In one aspect,two new threshold function is presented based on the analysis of th

7、e traditional soft—threshold,hard—threshold functions and the</p><p>  Key words:wavelet de-noising threshold estimation threshold functions white Gaussi

8、an noise</p><p><b>  目 錄</b></p><p>  摘要 基于DSP的語音信號的去噪方法1</p><p>  目錄 ...............................................................................................

9、.....................................2</p><p><b>  第一章 緒論3</b></p><p>  1.1課題的研究背景及意義4</p><p>  1.2語音去噪的基礎(chǔ)知識4</p><p>  1.2.1語音信號的特點5</p><p&g

10、t;  1.2.2噪聲特性5</p><p>  1.3 TI DSP 介紹6</p><p>  1.3.1 C54x芯片6</p><p>  1.4. CCS開發(fā)環(huán)境7</p><p>  1.4.1 CCS集成開發(fā)環(huán)境8</p><p>  1.5 本文重點總結(jié)9</p><p

11、>  第二章 小波在語音消噪中的應(yīng)用9</p><p>  2.1經(jīng)典的小波算法10</p><p>  2.1.1 模極大值去噪算法10</p><p>  2.1.2相關(guān)性去噪算法10</p><p>  2.1. 3小波閾值去噪算法10</p><p>  2.2一種新的小波閾值去噪方法11&

12、lt;/p><p>  2.3 本章總結(jié)13</p><p>  第三章 實驗設(shè)計原理及過程14</p><p>  3.1 小波去噪原理分析14</p><p>  3.1.1. 小波去噪原理14</p><p>  3.1.2 小波去噪步驟15</p><p>  3.2.閾值的選

13、取與量化15</p><p>  3.2.1. 軟閾值和硬閾值15</p><p>  3.2.2. 閾值的幾種形式16</p><p>  3.2.3. 閥值的選取16</p><p>  3.3實驗步驟17</p><p>  3.3.1小波去噪驗證仿真17</p><p>  

14、3.3.2. 小波去噪與FFT去噪效果對比18</p><p>  3.4本章總結(jié)19</p><p>  第四章 實驗仿真及圖20</p><p>  4.1實驗結(jié)論20</p><p><b>  結(jié)論22</b></p><p><b>  致謝23</b>

15、</p><p>  參 考 文 獻(xiàn)24</p><p><b>  附錄25</b></p><p><b>  第一章 緒論</b></p><p>  語音信號處理是信息高速公路、多媒體技術(shù)、辦公自動化、現(xiàn)代通信及智能系統(tǒng)等新興領(lǐng)域應(yīng)用的核心技術(shù)之一。在高度發(fā)達(dá)的信息社會,用數(shù)字化的方法

16、進(jìn)行語音的傳送、存儲、識別、合成、增強(qiáng)等是整個數(shù)字化通信網(wǎng)中最重要、最基本的組成部分之一。同時,語言不僅是人類相互間進(jìn)行溝通最自然和最方便的形式,也是人與機(jī)器之間進(jìn)行通信的重要工具,它是一種理想的人機(jī)通信方式,因而可為計算機(jī)、自動化系統(tǒng)等建立良好的人機(jī)交互環(huán)境,進(jìn)一步推動計算機(jī)和其他智能機(jī)器的應(yīng)用,提高社會的信息化和自動化程度。語音處理技術(shù),其中最重要的包括語音編碼、語音合成、語音識別、說話人識別及語音增強(qiáng),它的應(yīng)用極其廣泛,包括工業(yè)、

17、軍事、交通、醫(yī)學(xué)、民用等各個領(lǐng)域。目前,語音處理技術(shù)處于蓬勃發(fā)展的時期,己有大量產(chǎn)品投放市場,并且不斷有新產(chǎn)品被開發(fā)研制,具有及其廣闊的市場需要和應(yīng)用前景。</p><p>  1.1課題的研究背景及意義</p><p>  語音是語言信息的載體,是人類交流信息最自然、最有效、最方便的手段。在許多情況下,人們所關(guān)心的語音信號不可避免地被其他信號所污染,即語音總會受到外界環(huán)境噪聲的干擾,這些

18、噪聲包括從周圍環(huán)境、傳輸媒質(zhì)中引入的噪聲、電器設(shè)備的噪聲以及其他說話人干擾等等。環(huán)境噪聲會影響語音質(zhì)量,嚴(yán)重的情況下語音將完全淹沒到噪聲中,無法分辨。這些背景噪聲對人們得到語言信號中夾帶的有效信息產(chǎn)生了很大的干擾,會引起有用語音信息的相對失效。計算機(jī)無法準(zhǔn)確地分離出噪聲與有用信息,因此要讓它準(zhǔn)確無誤地執(zhí)行所要求執(zhí)行的操作,就必須對語音進(jìn)行降噪處理。研究如何進(jìn)行語音去噪的科學(xué)意義很大,它能幫助人們解決很多問題,同時它具有廣闊的應(yīng)用前景。較

19、好的降噪處理結(jié)果,有利于語音信號的編碼,從而也會減少數(shù)據(jù)的存儲量,利于在網(wǎng)絡(luò)中的傳送。語音增強(qiáng)技術(shù)就是從噪聲背景中提取、增強(qiáng)有用的信號,抑制、降低噪聲干擾的技術(shù)。比如,語音識別系統(tǒng)在實驗室環(huán)境中可取得相當(dāng)好的效果,但在噪聲環(huán)境中,尤其是在強(qiáng)噪聲環(huán)境中使用時,系統(tǒng)的識別率將受到嚴(yán)重影響。由于干擾通常都是隨機(jī)的,因而從帶噪語音中提取相對完全純凈的語音幾乎不可能。。語音降噪處理的意義非常重大,它可</p><p>  

20、1.2語音去噪的基礎(chǔ)知識</p><p>  語音去噪不但與信號處理技術(shù)有關(guān),而且涉及到人的聽覺感知和語音學(xué)。另外,噪聲來源不同,特性也不相同,因而很難找出一種適用于各種噪聲環(huán)境的去噪算法。對不同噪聲必須使用不同的語音去噪方法。語音去噪算法的基礎(chǔ)是對語音和噪聲特性的了解和分析。</p><p>  1.2.1語音信號的特點</p><p>  語音是時變的、非平穩(wěn)、

21、非遍歷的隨機(jī)過程。語音可以分為周期性的濁音和非周期的清音。濁音和清音經(jīng)常在一個音節(jié)中同時出現(xiàn)。濁音部分和音質(zhì)關(guān)系密切,在時域上呈現(xiàn)出明顯的周期性,在頻域上有共振峰結(jié)構(gòu),而且能量大部分集中在較低頻段內(nèi),是語音中大幅度高能量的部分,清音則沒有明顯的時域和頻域特征,類似于白噪聲,能量較小,在強(qiáng)噪聲中容易被掩蓋,但信噪比較高時能提供較多的信息。根據(jù)中心極限定理,傅里葉展開系數(shù)被認(rèn)為是獨(dú)立的高斯隨機(jī)變量,均值為零但方差是時變的。語音的短時譜幅度的

22、統(tǒng)計特性是時變的,只有當(dāng)分析幀長趨于無窮大時,才能近似認(rèn)為其具有高斯分布。這種高斯模型應(yīng)用于有限幀長時只是一種近似的描述,在寬帶噪聲污染的帶噪語音的語音去噪中,這種假設(shè)可作為分析的前提。</p><p><b>  1.2.2噪聲特性</b></p><p>  噪聲的來源很多,因此噪聲的特性也是多樣的。根據(jù)噪聲對語音頻譜的干擾方式,噪聲主要分為加性噪聲和乘性噪聲。當(dāng)

23、噪聲對語音的干擾表現(xiàn)為兩者信號在時域進(jìn)行相加時,該噪聲被稱為加性噪聲。顯然噪聲和語音在頻域也為相加關(guān)系。在實際環(huán)境中,背景噪聲可以看作加性噪聲,加性噪聲是對背景噪聲一種比較貼切的表述。乘性噪聲是指噪聲和語音在頻譜上是相乘的關(guān)系,在時域則是卷積的關(guān)系,因此也稱為卷積噪聲。</p><p>  我們關(guān)注的噪聲通常分為沖激噪聲、周期噪聲、寬帶噪聲、語音干擾噪聲等。非加性噪聲主要是傳送網(wǎng)絡(luò)的電路噪聲等。</p>

24、;<p>  (1) 周期噪聲:往往來源于電動機(jī)、風(fēng)扇等周期運(yùn)轉(zhuǎn)的機(jī)械,50或60Hz交流電源聲也是周期噪聲。它們的頻譜圖特點為有許多離散的窄譜峰,通常可以通過功率譜發(fā)現(xiàn),用濾波或變換技術(shù)予以濾除。</p><p>  (2) 沖激噪聲:通常是放電時會產(chǎn)生。它的時域波形是類似于沖激函數(shù)的窄脈沖。沖激噪聲可以在時域內(nèi)進(jìn)行消除。影響通常有兩種方法:對帶噪語音信號的幅度求均值,將該均值作為判斷閾值,凡

25、是超過該閾值的均判為沖激噪聲,把它進(jìn)行衰減或者完全置零:當(dāng)沖激脈沖不太密集時,也可以通過某些點內(nèi)插的方法避開或者平滑掉沖激點,從而能在重建語音信號中去掉沖激噪聲。</p><p>  (3) 寬帶噪聲:呼吸噪聲、風(fēng)、隨機(jī)噪聲源產(chǎn)生的噪聲、量化噪聲等均可視為寬帶噪聲,應(yīng)用中常近似為高斯噪聲或白噪聲。量化噪聲通常作為白噪聲來處理,也可以當(dāng)作寬帶噪聲處理。其顯著特點是噪聲頻譜遍布于語音信號域頻譜中,導(dǎo)致消除噪聲較為困

26、難。一般需要采用非線性處理方法。</p><p>  (4) 語音干擾:可能是受到其他語音源的干擾造成的。干擾語音信號和待傳語音信號同時在一個信道中傳輸所造成的干擾稱為語音干擾??紤]到一般情況下兩種語音的基音不同,也不成整數(shù)倍,利用它們的基音差別可以用梳狀濾波器提取基音和各次諧波,再恢復(fù)出有用語音信號。</p><p>  (5) 傳輸噪聲:這是傳輸系統(tǒng)的電路噪聲。與背景噪聲不同的是,

27、語音和噪聲需要通過卷積變換,轉(zhuǎn)化為加性噪聲來處理。</p><p>  1.3 TI DSP 介紹</p><p>  TI公司DSP種類多,品種齊全,適合各種需要。目前,使用較為廣泛的主要有三個系列:C2000,C5000和C6000。其他系列產(chǎn)品應(yīng)用較少或已經(jīng)淘汰。每個系列又有多種DSP可供選擇。同一系列的DSP具有相同的內(nèi)核、相同或兼容的匯編指令集;它們之間的差別是具有不同大小的片

28、內(nèi)存儲器、不同的片內(nèi)外設(shè)和外部接口等,工作電壓和速度也有所區(qū)別。以上3大系列DSP實現(xiàn)功能的側(cè)重點不同,也就是說應(yīng)用領(lǐng)域有所不同:</p><p>  (1)C2000系列是16位定點DSP。它是一個控制器系列,主要應(yīng)用于工業(yè)控制領(lǐng)域,它除了具有一個DSP內(nèi)核外,還有大量的片內(nèi)外設(shè)資源,如A/D、定時器、各種同步和異步串口、看門狗、CAN總線接口等;加上其價格低廉,速度更高,可靠性更強(qiáng),可以取代傳統(tǒng)單片機(jī)。內(nèi)部具

29、有Flash,方便固化程序,而其他系列DSP都沒有內(nèi)部Flash。</p><p>  (2)C5000系列DSP是16位定點低功耗DSP,性價比極高,主要應(yīng)用于無線通信系統(tǒng)及手持式通訊產(chǎn)品,如手機(jī),PDA和GPS等。C5000又分為C54X和C55X兩個系列。相對C2000系列來說,其內(nèi)部存儲更大,運(yùn)行速度更快,更適合執(zhí)行較為復(fù)雜的數(shù)字信號處理任務(wù),但控制功能相對較弱。一般由核心電壓和I/O電壓兩種電壓供電。核

30、心電壓較低,所以功耗很低,且體積很小,方便集成。</p><p>  (3)C6000系列是32位的DSP系列。在TI的所有系列DSP中運(yùn)行速度最快。其中C62XX是定點DSP,而C64XX和C67XX是浮點DSP,它主要應(yīng)用于需要大量快速運(yùn)算的場合,如數(shù)字視頻處理、無線基站等。由于速度很高,所以功耗也很大。</p><p>  1.3.1 C54x芯片</p><p&

31、gt;  在本設(shè)計中使用的DSP是TMS320VC5410。它屬于TI C5000系列中的C54X系列,正如前面所說,同一系列的DSP具有相同的內(nèi)核、相同或兼容的匯編指令集,差別僅在于內(nèi)存儲器的大小,片內(nèi)外設(shè)等等,所以就首先介紹一下C54系列DSP普遍具有的特點和性能。</p><p>  54X具有改進(jìn)的哈佛結(jié)構(gòu),使其處理能力達(dá)到最大。分開的程序空間和地址空間提供了高度的并行性,可以同時訪問程序指令和數(shù)據(jù),例如

32、三次讀操作和一次寫操作可以在一個周期內(nèi)完成。帶并行存儲的指令和具有特殊應(yīng)用的指令充分利用了這種結(jié)構(gòu)。這種并行性支持一套強(qiáng)大的算術(shù)運(yùn)算、邏輯運(yùn)算和位操作運(yùn)算,所以使得這些運(yùn)算可以在單個機(jī)器周期內(nèi)完成。而且54X的運(yùn)行機(jī)制還支持中斷處理、重復(fù)操作和函數(shù)調(diào)用等等。</p><p>  對于任意通用可編程芯片來說,一般都具備以下幾個部分:</p><p>  (1)內(nèi)部存儲器:主要用來存儲程序、執(zhí)

33、行程序、存儲數(shù)據(jù)等;</p><p>  (2)中央處理單元(CPU):用來實現(xiàn)各種運(yùn)算功能;</p><p>  (3)片內(nèi)外設(shè):用來實現(xiàn)一些特定功能,如時鐘發(fā)生器、硬件定時器等等;</p><p>  (4)外部總線接口:用來和其他芯片接口,協(xié)同工作;</p><p>  (5)通信接口:用來從外圍芯片獲得數(shù)據(jù)或者將處理完的數(shù)據(jù)傳輸出去;

34、</p><p>  (6)內(nèi)部總線:用來連接芯片中不同的單元。也就是說,以上各個部分的通信是通過內(nèi)部總縣來完成的。</p><p>  1.4. CCS開發(fā)環(huán)境</p><p>  本節(jié)將介紹CCS Code Composer Studio 的基本開發(fā)環(huán)境、軟件開發(fā)過程、CCS組件。CCS 提供了配置、建立、調(diào)試、跟蹤和分析程序的工具,它便于實時、嵌入式信號處理程

35、序的編制和測試,能夠加速開發(fā)進(jìn)程,提高工作效率。</p><p>  CCS 全稱是 Code Composer Studio 它提供了基本的代碼生成工具,具有一定的調(diào)試、分析能力,在CCS下的程序開發(fā)過程如圖2.4示。</p><p>  圖2.4程序開發(fā)過程</p><p><b>  CCS 包括:</b></p><

36、p>  1、CCS代碼生成工具</p><p>  2、CCS集成開發(fā)環(huán)境IDE</p><p>  3、DSP/BIOS插件程序和API RTDX插件,主機(jī)接口和API等。</p><p>  下面將著重介紹前兩項。</p><p>  1.4.1 CCS集成開發(fā)環(huán)境 </p><p>  調(diào)試DSP目標(biāo)程序,

37、它主要由幾個主要的窗口組成:工程組顯示窗口、程序內(nèi)容顯示窗口、編輯信息提示窗口和主要工具欄。另外,在編輯過程中 還可以顯示諸如存儲器觀察窗口、變量監(jiān)視框、圖形顯示框等調(diào)試界面,他們?yōu)槌绦蚓帉懻{(diào)試提供多種手段為軟件開發(fā)提供了極大的方便。圖 2.5 就是基本編輯界面。</p><p>  圖2.5 CCS基本編輯界面</p><p>  在利用CCS編程的過程中,不可避免的要遇到如何將程序變量

38、分配到內(nèi)存中去的問題。因為,即使 DSP 的存儲空間比較大,速度也十分快,但是如果內(nèi)存空間分配不當(dāng)?shù)脑挘€是會出現(xiàn)空間不夠行速度下降,甚至程序跑飛的情況,這樣會時程序調(diào)試起來十分麻煩。所以,一定要分配好內(nèi)存空間CCS提供了兩種分配空間的方法:利用cmd文件或rcp文件分配空間,其中cmd 文件是純文本格式的描述性的空間分配方式, 它的優(yōu)點是程序員對空間的可控制性較高,可以將不同的塊分配到指定的地址,并規(guī)定長度。但它對于初學(xué)者來說,要求對

39、C5410的內(nèi)存空間的分配有較清楚的認(rèn)識,否則容易將數(shù)據(jù)分配到不該分配的地方,引起程序運(yùn)行的沖突,甚至程序跑飛。所以,一般在編程過程中,使用 rcp 文件對程序進(jìn)行內(nèi)存分rcp文件全稱是 recipe 文件,這是一種圖形化界面的內(nèi)存分配文件。對于一個新生成的工程組rcp文件的生成方法是:選擇CCS菜單上的Tools項,選擇linker configuration項,將分配方式改為use the visual linker之后直接編譯,編

40、譯信息提示窗口中會顯示出錯信息,提示找不到rcp文件,雙擊提示C</p><p>  當(dāng)程序中新增加了變量后,在rcp文件中會出現(xiàn)Not yet placed項,只要將其下的文件夾,根據(jù)類型拖動到run_view下的相應(yīng)的數(shù)據(jù)存儲器或程序存儲器即可。還可以選擇用何種類型的存儲器空間裝載,通過觀察存儲空間以使用的狀況自行分配空間,所以十分靈活,并且不會出現(xiàn)空間重疊的現(xiàn)象,避免了程序跑飛。</p>&l

41、t;p>  1.5 本文重點總結(jié)</p><p>  人們研究了許多種語音降噪算法。目前有代表性的語音去噪算法主要有自相關(guān)相減法、自適應(yīng)噪聲濾波法、諧波增強(qiáng)法、基于語音生成模型的降噪算法、短時譜幅度估計算法、譜相減算法、小波變換法、人耳掩蔽效應(yīng)法等。本文所研究的重點在于通過小波分析方法對含高斯白噪聲的語音信號的去噪。論文中的實驗中采用的語音信號有兩種,一種是在安靜的環(huán)境下,用普通計算機(jī)錄制的語音信號。用Ma

42、tlab軟件實現(xiàn)仿真,用函數(shù)生成噪聲信號并對錄制的語音信號進(jìn)行軟件加噪,然后用提出的改進(jìn)思想進(jìn)行降噪處理。另一種就是錄制的帶有機(jī)箱噪聲的語音信號。</p><p>  本文研究的主要目的就是適應(yīng)目前信息社會對語音信號處理技術(shù)的要求,結(jié)合DSP理論在小波信號去噪方面的應(yīng)用,對小波語音去噪的處理過程提出改進(jìn)。主要在以下幾個方面進(jìn)行研究:</p><p>  (1)分析已有的語音降噪方法,包括短

43、時譜估計算法、譜減法等。</p><p>  (2)對小波去噪方法進(jìn)行研究,介紹小波去噪基本知識及原理,重點研究和討論小波閾值函數(shù)和閾值兩個方面。</p><p>  第二章 小波在語音消噪中的應(yīng)用</p><p>  小波分析理論是一種新興的信號處理理論,它在時間上和頻率上都有很好的局部性,這使得小波分析非常適合于時-頻分析,借助時- 頻局部分析特性,小波分析理

44、論已經(jīng)成為信號去噪中的一種重要的工具。利用小波方法去噪,是小波分析應(yīng)用于實際的重要方面。小波去噪的關(guān)鍵是如何選擇閾值和如何利用閾值來處理小波系數(shù),通過對小波閾值化去噪的原理介紹,運(yùn)用MATLAB 中的小波工具箱,對一個含噪信號進(jìn)行閾值去噪,實例驗證理論的實際效果,證實了理論的可靠性。本文簡述了幾種小波去噪方法,其中的閾值去噪的方法是一種實現(xiàn)簡單、效果較好的小波去噪方法。</p><p>  2.1經(jīng)典的小波算法&

45、lt;/p><p>  2.1.1 模極大值去噪算法</p><p>  模極大值去噪算法是根據(jù)信號和噪聲在多尺度空間上小波變換系數(shù)的模極值傳播規(guī)律的不同而發(fā)展起來的一種去噪算法。理論上只要信號與噪聲的奇異性有差異,就能產(chǎn)生很好的去噪效果。一般信號小波系數(shù)的模極大值將隨著小波分解層數(shù)的增大而增大;而對于白噪聲信號,其模值隨著分解層數(shù)的增大而減小。因此,觀察不同尺度間小波變換模極大值變化的規(guī)律,

46、去除幅度隨尺度的增加而減小的點(對應(yīng)噪聲的極值點),保留幅度隨尺度增加而增大的點(對應(yīng)于有用信號的極值點)。然后再由保留的模極大值點用交替投影法進(jìn)行重建,即可以達(dá)到去噪的目的。但是,交替投影法算法復(fù)雜,容易造成投影信號的偏差,難以在實際應(yīng)用中對信號進(jìn)行實時處理。</p><p>  2.1.2相關(guān)性去噪算法</p><p>  相關(guān)性去噪算法是根據(jù)信號經(jīng)小波變換后,其小波系數(shù)在各尺度上有較

47、強(qiáng)的相關(guān)性,尤其是在信號的邊緣附近,其相關(guān)性更加明顯,而噪聲對應(yīng)的小波系數(shù)在各尺度間卻沒有這種明顯的相關(guān)性來去噪的。在尺度空間上的相關(guān)運(yùn)算能使噪聲的幅值大為減小,從而抑制了噪聲和小的邊緣,增強(qiáng)了信號的主要邊緣,更好地刻畫了原始信號。并且在小尺度上,這種作用明顯大于在大尺度上的作用。由于噪聲能量主要是分布在小尺度上,因而這種隨尺度增大而作用強(qiáng)度遞減的性質(zhì),恰好濾除了噪聲,很好的保留了有用信號。</p><p>  

48、2.1. 3小波閾值去噪算法</p><p>  小波變換具有一種“集中”的能力,能將信號的能量集中到少數(shù)的小波系數(shù)上,而白噪聲在任何正交基上的變換仍然是白噪聲,其分量分布在大多數(shù)展開系數(shù)上。相對來說,有用信號所對應(yīng)的小波系數(shù)幅值較大,但數(shù)目較少,而噪聲對應(yīng)的小波系數(shù)是一致分布的,個數(shù)較多,但幅值小。基于這一思想,Donoho等人提出硬閾值和軟閾值去噪方法,即在眾多小波系數(shù)中,把絕對值較小的系數(shù)置為零,而讓絕對值

49、較大的系數(shù)保留或收縮,得到估計小波系數(shù),然后利用估計小波系數(shù)直接進(jìn)行信號重構(gòu),即可達(dá)到去噪的目的。</p><p>  小波閾值去噪算法主要分為4個步驟:</p><p>  第1步:選擇一個小波基函數(shù),確定小波分解層數(shù)并對信號進(jìn)行小波分解。常用于去噪的小波函數(shù)有dbN小波、symN小波和coifN小波,層數(shù)一般為3-5層。</p><p><b>  第

50、2步:閾值的確定</b></p><p>  小波閾值在去噪過程中起到?jīng)Q定性的作用。如果閾值太小,那么閾值處理后的小波系數(shù)中包含了過多的噪聲分量;如果閾值太大,那么將會丟失信號的一部分有用信息, 從而造成小波系數(shù)重構(gòu)后的信號失真。常用的閾值選擇方案有四種:</p><p> ?。ㄒ唬┗跓o偏似然估計的軟閾值估計(rigsure)</p><p> ?。ǘ?/p>

51、)長度對數(shù)閾值(sqtwolog)計算公式為:</p><p>  (三)啟發(fā)式SURE閾值(heursure)是前兩種閾值的綜合。</p><p> ?。ㄋ模┳钚O大方差閾值(minimaxi)產(chǎn)生一個最小均方誤差的極值,在給定的函數(shù)集中實現(xiàn)最大均方誤差最小化。</p><p>  在實際應(yīng)用中,斯坦無偏似然估計和最大最小值法比較保守,而其他兩種方法產(chǎn)生的閾值則過

52、大。</p><p>  第3步:選擇合適的閾值函數(shù)對小波系數(shù)進(jìn)行閾值處理。</p><p>  常用的閾值函數(shù)有以下兩種:</p><p>  硬閾值法: </p><p><b> ?。?)</b></p><p>  軟閾值法: (

53、2)</p><p>  第4步:小波重構(gòu)。根據(jù)閾值化處理后的高頻小波系數(shù)以及未處理的低頻小波系數(shù)進(jìn)行離散小波反變換重構(gòu)信號。</p><p>  2.2一種新的小波閾值去噪方法</p><p>  小波閾值算法比較簡單,運(yùn)算量小,因此在信號去噪方面得到了廣泛的應(yīng)用。但是這種方法還存在以下兩種不足:①白噪聲信號的模極大值隨著分解層數(shù)的增大而減小,對不同的分解層數(shù)采用

54、相同的閾值進(jìn)行處理,會在低頻系數(shù)中濾除過多有用信息而在高頻系數(shù)保留了一些過多的噪聲;②在進(jìn)行閾值處理時,硬閾值處理能更多的保留真實信號的尖峰等特征,但由于其本身的不連續(xù)性,去噪的的結(jié)果會出現(xiàn)震蕩,軟閾值是一種更平滑的方式,在去噪后能產(chǎn)生更光滑的結(jié)果,但估計的小波系數(shù)與原小波系數(shù)之間存在恒定的偏差。本文采用的改進(jìn)方法如下:</p><p>  一、各層采用不同的閾值</p><p>  由文

55、獻(xiàn)[4]可知: (3)</p><p>  其中K為一個常數(shù),j為分解層數(shù),為Lip指數(shù),為第j層的小波系數(shù)。Lip指數(shù)與信號的奇異性有關(guān),越大,信號越平滑。對于一般信號,即有用信號對應(yīng)小波系數(shù)隨分解尺度j的增大而變大。而白噪聲的Lip指數(shù)則為負(fù)值, ,即噪聲對應(yīng)的小波系數(shù)隨分解尺度j的增大而減小。</p><p>  有文獻(xiàn)[8]

56、可知,白噪聲的Lip指數(shù)滿足(4)式:</p><p><b> ?。?)</b></p><p>  由(3)式和(4)式可知</p><p><b> ?。?)</b></p><p>  其中為噪聲對應(yīng)的第j層小波系數(shù),由式 可知噪聲對應(yīng)的第j+1層小波系數(shù)的最大值小于第j層小波系數(shù)的最大值

57、,因此,本論文在閾值處理時每層系數(shù)采用不相同的閾值,用前面所述的四種閾值計算方法確定第一層閾值,以后各層閾值為前一層閾值的倍,即。</p><p>  二、采用新的閾值函數(shù)</p><p>  鑒于軟閾值和硬閾值的缺點,本論文采用一種新的閾值函數(shù)。該閾值函數(shù)既要保證其連續(xù)性,又要盡可能的消除軟閾值函數(shù)中的恒定偏差。閾值函數(shù)如下:

58、 (6)</p><p><b>  其中,</b></p><p>  式中N為一個正常數(shù),為處理前的小波系數(shù),為處理后的小波系數(shù)。該閾值函數(shù)和軟閾值函數(shù)一樣具有連續(xù)性,而且當(dāng)時, 函數(shù)是高階可導(dǎo)的,并且隨著系數(shù)的增大,t的值逐漸減小,使得處理前后的小波系數(shù)不變。當(dāng)N取值很大時, 新閾值函數(shù)相似于軟閾值函數(shù); 當(dāng)N 趨近于0時, 新閾值函數(shù)相似于硬閾值

59、函數(shù)。通過調(diào)節(jié)N的大小可以改變新閾值函數(shù)的類型,與經(jīng)典的閾值函數(shù)相比更具有靈活性。</p><p>  圖1硬閾值函數(shù)、軟閾值函數(shù)及改進(jìn)法制函數(shù)的曲線</p><p>  圖1為硬閾值函數(shù)、軟閾值函數(shù)及改進(jìn)法制函數(shù)的曲線圖,其中閾值。從圖中可以看出,硬閾值函數(shù)在閾值點處不連續(xù),軟閾值函數(shù)存在恒定的偏差。改進(jìn)閾值函數(shù)在閾值點處連續(xù),并且以硬閾值曲線為漸近線。</p><p

60、>  本章首先介紹小波去噪的定義及其特點,然后了解小波變換進(jìn)行降噪的原理,對閾值去噪方法做重點研究。下一章將用仿真實驗來說明消噪效果。</p><p>  2.3 本章總結(jié)</p><p>  本章首先介紹了小波閾值去噪的基本原理以及~般方法,并且消噪時兩個重要的方面進(jìn)行重點研究,提出了一些改進(jìn)。</p><p>  第一方面是對閾值函數(shù)進(jìn)行改進(jìn)。由于傳統(tǒng)

61、的軟閾值函數(shù)和硬閾值函數(shù)各有各的優(yōu)缺點。因此針對它們的缺點,提出了兩種改進(jìn)的閾值函數(shù),一種是介于軟閾值和硬閾值函數(shù)。另一個是基于非負(fù)死區(qū)閾值函數(shù)改進(jìn)的函數(shù)。本章就其單調(diào)性、連續(xù)性方面給出的證明。</p><p>  第二個改進(jìn)的方面是閾值的估計。受到了譜減法思想的啟發(fā),對噪聲的估計也是基于無音段的噪聲分析。由于信號的小波系數(shù)與噪聲的小波系數(shù)的不同表現(xiàn),所以把噪聲經(jīng)過小波分析后的模最大值作為閾值。這樣就可以較好地恢

62、復(fù)出純凈語音,當(dāng)然這個過程中語音信號也會受到一定程度的損害。</p><p><b>  實驗設(shè)計原理及過程</b></p><p>  3.1 小波去噪原理分析</p><p>  3.1.1. 小波去噪原理</p><p>  疊加性高斯白噪聲是最常見的噪聲模型,受到疊加性高斯白噪聲“污染”的觀測信號可以表示為:&l

63、t;/p><p>  (7) </p><p>  其中yi為含噪信號,為“純凈”采樣信號,zi為獨(dú)立同分布的高斯白噪聲,為噪聲水平,信號長度為n. 為了從含噪信號yi中還原出真實信號,可以利用信號和噪聲在小波變換下的不同的特性,通過對小波分解系數(shù)進(jìn)行處理來達(dá)到信號和噪聲分離的目的。在實際工程應(yīng)用中,有用信號通常表現(xiàn)為低頻信號或是一些比較平穩(wěn)的信號,而噪聲信號則通常表現(xiàn)為高頻信號,

64、所以我們可以先對含噪信號進(jìn)行小波分解(如進(jìn)行三層分解):</p><p><b>  (8)</b></p><p>  圖1 三層小波分解示意圖</p><p>  其中為分解的近似部分, 為分解的細(xì)節(jié)部分,,則噪聲部分通常包含在,,中,用門限閾值對小波系數(shù)進(jìn)行處理,重構(gòu)信號即可達(dá)到去噪的目的。</p><p>  3

65、.1.2 小波去噪步驟</p><p>  總結(jié)去噪過程,可以分成以下三個步驟:</p><p>  1) 對觀測數(shù)據(jù)作小波分解變化:</p><p><b>  (9)</b></p><p>  其中y表示觀測數(shù)據(jù)向量y1,y2,…y,f是真實信號向量f1,f2,…fn,z是高斯隨 機(jī)向量z1,z2,…zn ,其中

66、用到了小波分解變換是線性變換的性質(zhì)。</p><p>  2)對小波系數(shù)W0作門限閾值處理(根據(jù)具體情況可以使用軟閾值處理或硬閾值 處理,而且可以選擇不同的閾值形式,這將在后面作詳細(xì)討論),比如選取最著 名的閾值形式:</p><p><b>  (10)</b></p><p>  門限閾值處理可以表示為,可以證明當(dāng)n 趨于無窮大時使

67、用閾值公式(4)對 小波系數(shù)作軟閾值處理可以幾乎完全去除觀測數(shù)據(jù)中的噪聲。</p><p>  3) 對處理過的小波系數(shù)作逆變換重構(gòu)信號:</p><p><b>  (11)</b></p><p>  即可得到受污染采樣信號去噪后的信號。</p><p>  3.2.閾值的選取與量化</p>&

68、lt;p>  Donoho-Johnstone小波收縮去噪方法的關(guān)鍵步驟是如何選擇閾值和如何進(jìn)行門限閾值處理,在這將作較為詳細(xì)的討論。</p><p>  3.2.1. 軟閾值和硬閾值</p><p>  在對小波系數(shù)作門限閾值處理操作時,可以使用軟閾值處理方法或硬閾值處理方法,硬閾值處理只保留較大的小波系數(shù)并將較小的小波系數(shù)置零:</p><p><b

69、>  (12)</b></p><p>  軟閾值處理將較小的小波系數(shù)置零但對較大的小波系數(shù)向零作了收縮:</p><p><b>  (13)</b></p><p>  直觀形式見圖2(圖中取t=1)從圖上我們可以看出軟閾值處理是一種更為平滑的形式,在去噪后能產(chǎn)生更為光滑的結(jié)果,而硬閾值處理能夠更多的保留真實信號中的尖峰等

70、特征軟閾值處理實質(zhì)上是對小波分解系數(shù)作了收縮,從而Donoho-Johnstone將這種去噪技術(shù)稱之為小波收縮。</p><p>  圖2 硬閥值和軟閥值</p><p>  3.2.2. 閾值的幾種形式</p><p>  閾值的選取有多種形式,選取規(guī)則都是基于含噪信號模型式(1.1)中信號水平為1 的情況,對于噪聲水平未知或非白噪聲的情況可以在去噪時重新調(diào)整得

71、到的閾值。</p><p>  在MATLAB中有4種閾值函數(shù)形式可以選用:</p><p>  (1) sqtwolog:采用固定的閾值形式,如式(10),因為這種閾值形式在軟門限閾值處理中能夠得到直觀意義上很好的去噪效果。</p><p>  (2) minimaxi采用極大極小原理選擇的閾值,和sqtwolog一樣也是一種固定的閾值,它產(chǎn)生一個最小均方誤差的極

72、值,計算公式為:</p><p><b>  (14)</b></p><p>  (3) rigrsure:采用史坦的無偏似然估計原理進(jìn)行閾值選擇,首先得到一個給定閾值的風(fēng)險估計,選擇風(fēng)險最小的閾值 作為最終選擇。</p><p>  (4) heursure:選擇啟發(fā)式閾值它是sqtwolog和rigrsure 的綜合,當(dāng)信噪比很小時,估計

73、有很大的噪聲,這時heursure, 采用固定閾值sqtwolog。</p><p>  3.2.3. 閥值的選取</p><p>  閾值化處理的關(guān)鍵問題是選擇合適的閾值如果閾值(門限) 太小,去噪后的信號仍然有噪聲存在;相反,如果太大,重要信號特征將被濾掉,引起偏差。從直觀上,對于給定小波系數(shù),噪聲越大,閾值就越大。大多數(shù)閾值選擇過程是針對一組小波系數(shù),即根據(jù)本組小波系數(shù)的統(tǒng)計特性,計

74、算出一個閾值。</p><p>  Donoho 等提出了一種典型閾值選取方法,從理論上給出并證明閾值與噪聲的方差成正比,其大小為:</p><p><b>  3.3實驗步驟</b></p><p>  3.3.1小波去噪驗證仿真</p><p>  實驗信號是由wnoise()函數(shù)產(chǎn)生的含標(biāo)準(zhǔn)的高斯白噪聲信噪比為3

75、的heavy sine信號,用wden()函數(shù)進(jìn)行去噪處理。</p><p>  1) 首先產(chǎn)生一個長度為210點,包含高斯噪聲的heavy sine信號及heavy sine含噪信號, 其噪聲標(biāo)準(zhǔn)差為3 , 如圖3 a及b所示。</p><p>  2) 利用‘sym8’小波對信號分解,在分解的第5層上,利用軟閾值法去噪,結(jié)果如圖3 c所示</p><p>  3

76、) 同樣的條件下,利用固定閾值選擇算法對信號去噪,結(jié)果如圖3 d 所示。</p><p><b>  驗證仿真程序如下:</b></p><p>  x=wnoise(3,10);</p><p>  ind=linspace(0,1,2^10);</p><p>  subplot(4,1,1);</p>

77、<p><b>  plot(x);</b></p><p>  title('(a)');</p><p>  [x,noisyx]=wnoise(3,10,3,2^10);</p><p>  subplot(4,1,2);</p><p>  plot(noisyx);</p>

78、;<p>  title('(b)');</p><p>  xd=wden(x,'rigrsure','s','sln',5,'sym8');</p><p>  subplot(4,1,3);</p><p><b>  plot(xd);</b>

79、</p><p>  title('(c)')</p><p>  xd=wden(x,'sqtwolog','h','sln',5,'sym8');</p><p>  subplot(4,1,4);</p><p><b>  plot(xd);&l

80、t;/b></p><p>  title('(d)');</p><p>  3.3.2. 小波去噪與FFT去噪效果對比</p><p>  選擇MATLAB6.5 中含有噪聲的仿真信號noisbloc 作為原始信號,分別使用FFT 和小波分析方法對信號進(jìn)行去噪處理,采用的小波是sym8,分解層數(shù)為5,對比結(jié)果如圖4所示。</p>

81、<p><b>  驗證仿真程序如下:</b></p><p>  load noisbloc;</p><p>  x=noisbloc;</p><p>  subplot(2,2,1);</p><p>  plot(x);title('a')</p><p>

82、  xd=wden(x,'rigrsure','s','sln',5,'sym8');</p><p>  subplot(2,2,2);</p><p>  plot(xd);title('b')</p><p>  p1=1/length(x)*norm(x)^2; </p&g

83、t;<p>  p2=1/length(x)*norm(x-xd)^2; </p><p>  snr1=10*log(p1/p2) </p><p>  RMSE1=sqrtm(p2)</p><p>  xd=wden(x,'sqtwolog','h','sln',5,'sym8');

84、</p><p>  subplot(2,2,3);</p><p>  plot(xd);title('c')</p><p>  p1=1/length(x)*norm(x)^2; </p><p>  p2=1/length(x)*norm(x-xd)^2; </p><p>  snr2=10*

85、log(p1/p2)</p><p>  RMSE2=sqrtm(p2)</p><p><b>  wc=0.3;</b></p><p><b>  N=5;</b></p><p>  [b,a]=butter(N,wc);</p><p>  xd=filter(b,

86、a,x);</p><p>  subplot(2,2,4);plot(xd);title('d');</p><p>  p1=1/length(x)*norm(x)^2; </p><p>  p2=1/length(x)*norm(x-xd)^2; </p><p>  snr3=10*log(p1/p2)</p&

87、gt;<p>  RMSE3=sqrtm(p2)</p><p><b>  3.4本章總結(jié)</b></p><p>  基于DSP的算法 實現(xiàn)算法實現(xiàn)主要分為三個 部分 :小波分解部分、小波系數(shù)閾值 量化部分和小波重構(gòu)部分。具體程序?qū)崿F(xiàn)流程為圖1。</p><p><b>  實驗仿真及圖</b></

88、p><p><b>  4.1實驗結(jié)論</b></p><p>  圖3 小波去噪驗證仿真圖</p><p>  (a)為原始信號(b)為含噪信號(c)為軟閥值去噪信號(d)為硬閥值去噪信號</p><p>  通過對圖3的分析,可以看出對原始信號添加噪聲后得到含噪信號,利用MATLAB 中的小波工具箱對含噪信號分別進(jìn)行軟閾

89、值化和硬閾值化去噪處理,得出的去噪結(jié)果與原始信號效果非常接近, 由此可以看出利用MATLAB 中的小波變換工具箱對信號進(jìn)行去噪處理是非常理想的。</p><p>  由圖4可以看出,利用小波分析去噪的結(jié)果明顯優(yōu)于Fourier 變換,這是由于Fourier 變換只能在頻域范圍內(nèi)表述,對系數(shù)進(jìn)行處理的方法也相對單一,而利用小波對信號進(jìn)行分解后,可以采用多種計算閾值和處理閾值的方法. 對信號的噪聲成分進(jìn)行抑制,手段更

90、加靈活。為了更加精確地表示去噪結(jié)果,可以計算去噪后信號的信噪比(RSN) 和均方根誤差(RMSE)。</p><p>  圖4 小波去噪和FFT去噪效果對比圖</p><p>  a為含噪信號圖,b為軟閥值去噪信號圖,c為硬閥值去噪圖,d為FFT去噪圖</p><p>  表4 幾種方法去噪后的RSN 和RMSE</p><p>  信號

91、的信噪比越高,原始信號和去噪信號的均方根誤差越小,去噪信號就越接近原信號,去噪的效果也就越好。表4.給出了3種方法去噪后信噪比和均方根誤差的比較,可以看出,小波分析去噪結(jié)果的信噪比和均方根誤差指標(biāo)均優(yōu)于FFT。</p><p><b>  結(jié)論</b></p><p>  數(shù)模極大值之間的關(guān)系,提出對各層系數(shù)采用不同的閾值進(jìn)行處理,并針對軟閾值函數(shù)和本文利用噪聲與信號

92、在小波變換各層系數(shù)的傳輸特性的不同預(yù)計噪聲對應(yīng)的小波系硬閾值函數(shù)各自的缺陷,提出了一種新的閾值函數(shù),更增加了小波去噪的靈活性。通過最后的仿真結(jié)果表明,新的小波閾值去噪方法可大幅度提升去噪效果,具有較好的工程應(yīng)用價值。小波去噪是一個正在研究的課題,新的方法在不斷地提出。小波變換是一種信號的時頻分析方法,它具有多分辨率分析的特點,很適合探測正常信號中夾帶的瞬態(tài)反?,F(xiàn)象并展示其成分,有效區(qū)分信號中的突變部分和噪聲。因此利用小波變換能有效的對信

93、號進(jìn)行消噪的同時提取含噪信號。用傳統(tǒng)的傅立葉變換分析,顯得無能為力,因為傅立葉分析是將信號完全在頻率域中進(jìn)行分析,它不能給出信號在某個時間點的變化情況,使得信號在時間軸上的任何一個突變,都會影響信號的整個頻譜。小波變換正廣泛的應(yīng)用于各個領(lǐng)域,MATLAB給我們提供了一個很方便的工作平臺,通過MATLAB編制程序給定信號的噪聲抑制和非平穩(wěn)信號噪聲的消除。通過實例證明:基于小波變換的消噪方法是一種提取有用信號、展示噪聲和突變信號的優(yōu)越方法,

94、具有廣闊的實用價值。</p><p><b>  致謝</b></p><p>  短暫而充實的大學(xué)四年時光即將接近尾聲,而該畢業(yè)論文的設(shè)計是大學(xué)最后一個學(xué)期的學(xué)習(xí)目標(biāo)和任務(wù),在論文完成之際,特向我的指導(dǎo)老師xx以誠摯的謝意。</p><p>  xx老師為人誠懇,待人和藹可親,這種可貴的待人精神為我的畢業(yè)論文的設(shè)計營造了良好的氛圍。在本次設(shè)計

95、中,xx老師給予我很大的關(guān)懷和幫助,特別是在論文的選題、前期課題指導(dǎo)及各段落內(nèi)容銜接等各環(huán)節(jié)都細(xì)心的指導(dǎo)和講解,由于我在以前的課程學(xué)習(xí)中沒有扎實的數(shù)字信號處理基礎(chǔ)及Matlab程序設(shè)計經(jīng)驗,xx老師為我提供了諸多方面的資料,讓我在課題設(shè)計中更加得心應(yīng)手,最終才能圓滿的完成畢業(yè)論文的設(shè)計。</p><p>  在這次畢業(yè)論文設(shè)計中,使我深刻的體會到只有理論知識還是不夠的,要能在原有的理論基礎(chǔ)上加以動手實踐,使理論和

96、實踐很好的結(jié)合起來,這樣才能對該知識點有更加深刻的認(rèn)識,希望這種經(jīng)驗體會能夠為我以后的學(xué)習(xí)工作中不斷成長不斷進(jìn)步。</p><p>  最后祝愿我所有的老師和同學(xué)工作順利,闔家歡樂,一生平安。</p><p><b>  參 考 文 獻(xiàn)</b></p><p>  【1】胡昌華,張軍波,夏軍.基于MATLAB的系統(tǒng)分析與設(shè)計[M],

97、 西安:西安電子科技【J】大學(xué)出版社,2001,9</p><p>  【2】朱云芳,戴朝華,陳維榮.小波消噪閾值選取的一種改進(jìn)方 法[J].電測與儀表,2005,42(475),4-6 </p><p>  【3】Y Zhou,S.Lai,L.Liu,P.Lv;J.Yu.A Novel Approach To Image Denoising in Multi—wavelet

98、 Domain[C].International Conference on Communications,Circuits and Systems Proceedings,2006:463—466</p><p>  【4】L.Su,G Zhao.De-Noising of ECG Signal Using Translation-Invariant Wavelet De-Noising Method Impr

99、oved Thresholding[C].27th IEEE EMBS(Engineering in Medicine and Biology Society)Annual International Conference,2005,Shanghai,China:5946—5949</p><p>  【5】高成.Matlab小波分析及應(yīng)用.國防工業(yè)出版社.2007</p><p>  【

100、6】劉濤.使用小波分析入門.國防工業(yè)出版社.2006</p><p>  【7】 金寶龍、李輝、趙乃杰、何海峰.一種新的小波閾值去噪算法[M].彈箭與制導(dǎo)學(xué)報.2011</p><p>  【8】劉宗昂、楊莘元、王麗安.一種新的小波去噪算法[M]. 彈箭與制導(dǎo)學(xué)報.2009</p><p>  【9】吳偉、蔡培升.基于Matlab的小波去噪仿真[M]. 信息與電子工

101、程.2008</p><p>  【10】李祥兵、肖合林. 基于Matlab的小波閾值折中算法研究[M]. 電腦開發(fā)與應(yīng)用.2009</p><p>  【11】杜浩潘、叢爽. 基于Matlab小波去噪方法研究[M]. 計算機(jī)仿真.2003</p><p>  【12】張兆寧、董肖紅、潘云峰. 基于小波變換模極大值去噪方法的改進(jìn)[M]. 電力系統(tǒng)及其自動化學(xué)報.20

102、05</p><p><b>  附錄</b></p><p>  閾值去噪的驗證仿真程序:</p><p>  x=wnoise(3,10);</p><p>  ind=linspace(0,1,2^10);</p><p>  subplot(4,1,1);</p><p

103、><b>  plot(x);</b></p><p>  title('(a)');</p><p>  [x,noisyx]=wnoise(3,10,3,2^10);</p><p>  subplot(4,1,2);</p><p>  plot(noisyx);</p><

104、;p>  title('(b)');</p><p>  xd=wden(x,'rigrsure','s','sln',5,'sym8');</p><p>  subplot(4,1,3);</p><p><b>  plot(xd);</b></p

105、><p>  title('(c)')</p><p>  xd=wden(x,'sqtwolog','h','sln',5,'sym8');</p><p>  subplot(4,1,4);</p><p><b>  plot(xd);</b&g

106、t;</p><p>  title('(d)');</p><p>  小波與FFT對比驗證仿真程序:</p><p>  load noisbloc;</p><p>  x=noisbloc;</p><p>  subplot(2,2,1);</p><p>  plot

107、(x);title('a')</p><p>  xd=wden(x,'rigrsure','s','sln',5,'sym8');</p><p>  subplot(2,2,2);</p><p>  plot(xd);title('b')</p>&l

108、t;p>  p1=1/length(x)*norm(x)^2; </p><p>  p2=1/length(x)*norm(x-xd)^2; </p><p>  snr1=10*log(p1/p2) </p><p>  RMSE1=sqrtm(p2)</p><p>  xd=wden(x,'sqtwolog',&

109、#39;h','sln',5,'sym8');</p><p>  subplot(2,2,3);</p><p>  plot(xd);title('c')</p><p>  p1=1/length(x)*norm(x)^2; </p><p>  p2=1/length(x)*no

110、rm(x-xd)^2; </p><p>  snr2=10*log(p1/p2)</p><p>  RMSE2=sqrtm(p2)</p><p><b>  wc=0.3;</b></p><p><b>  N=5;</b></p><p>  [b,a]=butte

111、r(N,wc);</p><p>  xd=filter(b,a,x);</p><p>  subplot(2,2,4);plot(xd);title('d');</p><p>  p1=1/length(x)*norm(x)^2; </p><p>  p2=1/length(x)*norm(x-xd)^2; </

溫馨提示

  • 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)確性、安全性和完整性, 同時也不承擔(dān)用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。

評論

0/150

提交評論