小波分析理論及其應(yīng)用(共8頁(yè))

1、精選優(yōu)質(zhì)文檔-傾情為你奉上上海大學(xué)20102011學(xué)年冬季學(xué)期研究生課程課程名稱： 信息采集與處理技術(shù) 課程編號(hào)： 論文題目: 小波分析理論及其應(yīng)用 研究生姓名: 劉金鼎 學(xué) 號(hào): 論文評(píng)語(yǔ):成 績(jī): 任課教師: 昝鵬 評(píng)閱日期: 專心-專注-專業(yè)小波分析理論及其應(yīng)用劉金鼎（上海大學(xué) 機(jī)電工程與自動(dòng)化學(xué)院，上海 ）摘要：小波分析的理論與方法是從Fourier分析的思想方法演變而來(lái)的。就象Fourier分析分為積分Fourier變換和Fourier級(jí)數(shù)一樣，小波分析也分為(積分)小波變換和小波級(jí)數(shù)兩部分，(積分)小波變換的主體是連續(xù)小波變換，多尺度小波變換和s進(jìn)小波變換；而小波級(jí)數(shù)的主體部分是關(guān)

2、于小波框架的理論。小波分析理論深刻，應(yīng)用廣泛，并且仍在迅速發(fā)展之中。本文作者作為初學(xué)者，單單就（積分）小波變換這一理論中比較基本和初步的東西所作的一點(diǎn)歸納和整理，介紹了小波變換的定義及特點(diǎn)，以及多分辨率分析的問題，最后以一些圖像去噪應(yīng)用來(lái)形象說(shuō)明小波分析的作用。關(guān)鍵詞：傅里葉分析；小波分析；多分辨率PXI BusLIU Jin-ding(School of Mechatronics Engineering & Automation, Shanghai University, Shanghai , China)Abstract: The theory and methods of wav

3、elet analysis comes from Fourier analysis .Just as Fourier analysis is divided into Fourier transform and Fourier series, wavelet analysis is divided into the wavelet transform and wavelet series. The main body of the wavelet transform is the continuous wavelet transform, multi-scale wavelet transfo

4、rm and s-dyadic wavelet transform, while the main part of the wavelet series is wavelet frame. Wavelet analysis is a kind of profound theory, which is used widely and develops rapidly. The author of the paper is a beginner of wavelet theory; he just summarized and organized some fundamental theory o

5、f wavelet analysis. The paper introduced the definition and characteristics of wavelet analysis, and then talked about the theory of multi- resolution ratio. In the end, a few of image denoising abstract applications were used to explain the function of wavelet analysis vividly.Key words: Fourier an

6、alysis; wavelet analysis; multi- resolution ratio1 引言1.1 問題的提出Fourier變換只能告訴我們信號(hào)尺度的范圍，而無(wú)法給出信號(hào)的結(jié)構(gòu)以及它蘊(yùn)含的大小不同尺度的串級(jí)過(guò)程，即Fourier變換在時(shí)空域中沒有任何分辨率。此外，傅立葉分析無(wú)法解決信號(hào)奇異性的位置。20世紀(jì)80年代初由法國(guó)油氣工程師Morlet提出的小波分析1（wavelet Analysis，又稱子波分析)能成功地解決這些問題。因此小波分析是Fourier分析發(fā)展史上的一個(gè)里程碑。小波分析一面世，立刻成為國(guó)際研究熱點(diǎn)。目前小波分析在信號(hào)處理、圖像壓縮、語(yǔ)音編碼、模式識(shí)別、地震勘

7、探、大氣科學(xué)以及許多非線性科學(xué)領(lǐng)域內(nèi)取得了大量的研究成果。小波分析之所以廣泛得到應(yīng)用在于:它具有時(shí)域和頻域同時(shí)具有良好的局部性質(zhì)；能將信號(hào)(時(shí)間序列)分解成交織在一起的多尺度成分，從而能夠不斷地聚集到所研究對(duì)象的任意微小細(xì)節(jié)；同時(shí)具有數(shù)學(xué)上嚴(yán)格意義的突變點(diǎn)診斷能力。1.2 小波分析的形成及發(fā)展小波分析是一調(diào)和分析方法2,3，是Fourier分析發(fā)展史上的一個(gè)里程碑式的進(jìn)展，被人們譽(yù)為數(shù)學(xué)“顯微鏡”。小波分析理論及其方法的形成和應(yīng)用在科學(xué)技術(shù)界引起一場(chǎng)軒然大波并成蔓延之勢(shì)。小波理論形成經(jīng)歷了三個(gè)階段2:(1)Fourier變換(FT)階段：在信號(hào)分析中，我們對(duì)信號(hào)的基本刻化，往往采取時(shí)域和頻域兩

8、種基本形式。時(shí)域分析無(wú)法得到關(guān)于信號(hào)變化的更多信息(如采樣、周期等)。1822年Fourier提出的頻域分析法Fourier變換（），能揭示信號(hào)f(t)的能量在各個(gè)頻率成分中的分布情況。設(shè)信號(hào)為，其Fourier變換為：許多時(shí)域上看不清的問題，通過(guò)就顯得清晰了。Fourier變換將信號(hào)的時(shí)域特征和頻率特征聯(lián)系起來(lái)，能分別從時(shí)域和頻域上觀察信號(hào)，但不能把二者有機(jī)結(jié)合起來(lái)。另外，F(xiàn)ourier變換是整個(gè)時(shí)間域內(nèi)的積分，識(shí)別出的頻率在什么時(shí)候產(chǎn)生并不知道，因此不能反映某一局部時(shí)間內(nèi)信號(hào)的頻譜特性，即在時(shí)間域上沒有任何分辨率。這樣在信號(hào)分析中就面臨一對(duì)矛盾：時(shí)域和頻域的局部化矛盾。Fourier變換對(duì)

9、具有突變的信號(hào)，如地震波、暴雨、洪水等的分析帶來(lái)諸多不便和困難。這就促使尋求一種信號(hào)時(shí)頻局部分析新方法。 (2)短時(shí)Fourier變換(SFT)階段1946年Gabor提出SFT。短時(shí)Fourier變換又稱加窗Fourier變換，由Gabor1946年提出。其基本思想是:把信號(hào)劃分成許多小的時(shí)間間隔，用Fourier變換分析每一個(gè)時(shí)間間隔，以確定該間隔存在的頻率，以達(dá)到時(shí)頻局部化之目的。短時(shí)Fourier變換的表達(dá)式為：SFT能實(shí)現(xiàn)信號(hào)時(shí)頻局部化分析，但窗函數(shù)一選定，其窗口的大小和形狀固定不變，其分辨率是有限的。由于頻率與周期成反比，反映信號(hào)高頻成分需要較高的時(shí)間分辨率(窄的時(shí)間窗)，反映低頻

10、成分需要較低的時(shí)間分辨率(寬的時(shí)間窗)。因此，加窗Fourier變換對(duì)研究高頻率信號(hào)和低頻率信號(hào)都不是有效的。(3)小波分析階段小波分析是一種窗口的大小固定、形狀可變的時(shí)頻局部化信號(hào)分析方法，即在低頻部分具有較高的頻率分辨率和較低的時(shí)間分辨率，在高頻部分具有較高的時(shí)間分辨率和較低頻率分辨率。小波在繼承SFT的基礎(chǔ)上，Morlet提出了小波變換法(WT)。WT可研究信號(hào)在各個(gè)時(shí)刻或各空間位置在不同尺度上的演變情況，實(shí)現(xiàn)了時(shí)頻局部化分析。小波理論的思想源于信號(hào)分析的伸縮與平移。1980年由Morlet首創(chuàng)。1984年他與Grossman共同提出連續(xù)小波變換的幾何體系，成為小波分析發(fā)展的里程碑。19

11、85年，法國(guó)數(shù)學(xué)家Meyer創(chuàng)造性構(gòu)造了規(guī)范正交基，提出了多分辨率概念和框架理論。小波熱由此興起。1986年Battle和Lemarie記又分別獨(dú)立地給出了具有指數(shù)衰減的小波函數(shù)；同年，Mallat創(chuàng)造性地發(fā)展了多分辨分析概念和理論并提出子決速小波變換算法Mallat算法。Daubechies(1988)構(gòu)造了具有有限緊支集的正交小波基，Chui和王建忠(1990)構(gòu)造了基于樣條函數(shù)的正交小波。至此，小波分析的系統(tǒng)理論得以建立。最近有人又提出了小波包理論，它是小波理論的進(jìn)一步發(fā)展。 2 小波變換的基本理論小波即小區(qū)域的波，是一種特殊的長(zhǎng)度有限、平均值為零的波形。它有兩個(gè)特點(diǎn)：一是“小”，即在時(shí)

12、域具有緊支集或近似緊支集；二是正負(fù)交替的“波動(dòng)性”，也即支流分量為零。2.1連續(xù)小波變換4,52.1.1 連續(xù)小波基函數(shù)所謂小波(Wavelet)，即存在于一個(gè)較小區(qū)域的波。小波函數(shù)的數(shù)學(xué)定義是：設(shè)為一平方可積函數(shù)，即，若其傅立葉變換滿足： 時(shí)，則稱為一個(gè)基本小波或小波母函數(shù)，并稱上式是小波函數(shù)的可容許條件。根據(jù)小波函數(shù)的定義，小波函數(shù)一般在時(shí)域具有緊支集或近似緊支集，即函數(shù)的非零值定義域具有有限的范圍，這即所謂“小”的特點(diǎn);另一方面，根據(jù)可容許性條件可知，即直流分量為零，因此小波又具有正負(fù)交替的波動(dòng)性。 將小波母函數(shù)進(jìn)行伸縮和平移，設(shè)其伸縮因子(亦稱尺度因子)為，平移因子為，并記平移伸縮后的

13、函數(shù)為，則：并稱 為參數(shù) 和 小波基函數(shù)。由于 和 均取連續(xù)變換的值，因此又稱為連續(xù)小波基函數(shù)，它們是由同一母函數(shù) 經(jīng)伸縮和平移后得到的一組函數(shù)系列。定義小波母函數(shù)的窗口寬度為，窗口中心為，則可以求得連續(xù)小波基函數(shù)的窗口中心及窗口寬度分別為：設(shè)是的傅立葉變換，頻域窗口中心為，窗口寬度為，的傅立葉變換為，則有:所以此時(shí)頻域窗口中心及窗口寬度分別為： 由此可見，連續(xù)小波的時(shí)、頻窗口中心和寬度均是尺度因子的函數(shù)，均隨著的變化而伸縮，并且還有即連續(xù)小波基函數(shù)的窗口面積是不變的，這正是Heisenberg測(cè)不準(zhǔn)原理。將不同a、b值下的時(shí)頻窗口繪在同一個(gè)圖上，就得到小波基函數(shù)的相平面（如圖1所示）。圖1小

14、波基函數(shù)的相平面對(duì)不同的頻率成分，在時(shí)域上的取樣步長(zhǎng)是可調(diào)的，高頻者(對(duì)應(yīng)小的m值)采樣步長(zhǎng)小，低頻者(對(duì)應(yīng)大的m值)采樣步長(zhǎng)大。也就是說(shuō)，小波變換能實(shí)現(xiàn)了窗口的大小固定，形狀可變的時(shí)頻局部化，見圖1。正是這個(gè)意義上小波變換被譽(yù)為數(shù)學(xué)“顯微鏡”。2.1.2 連續(xù)小波變換將空間的任意函數(shù)在小波基下進(jìn)行展開，稱其為函數(shù)的連續(xù)小波變換CWT，變換式為： 當(dāng)小波的容許性條件成立時(shí)，其逆變換為：其中為的容許性條件我們可以這樣理解，傅立葉分析是將信號(hào)分解成一系列不同頻率的正弦波的疊加，同樣小波分析是將信號(hào)分解為一系列小波函數(shù)的疊加，而這些小波函數(shù)都是由一個(gè)母小波函數(shù)經(jīng)過(guò)平移和尺度伸縮得來(lái)的。小波分析優(yōu)于傅

15、立葉分析的地方是，它在時(shí)域和頻域同時(shí)具有良好的局部化性質(zhì)。而且由于對(duì)高頻成分采用逐漸精細(xì)的時(shí)域或頻域取樣步長(zhǎng)，從而可以聚焦到對(duì)象的任何細(xì)節(jié)，所以被稱為“數(shù)學(xué)顯微鏡”。小波分析廣泛應(yīng)用與信號(hào)處理、圖像處理、語(yǔ)音識(shí)別等領(lǐng)域?？梢赃@樣理解小波變換的含義：打個(gè)比喻，我們用鏡頭觀察目標(biāo)信號(hào)f (t)， (t)代表鏡頭所起的所用。b 相當(dāng)于使鏡頭相對(duì)于目標(biāo)平行移動(dòng)，a的所用相當(dāng)于鏡頭向目標(biāo)推進(jìn)或遠(yuǎn)離。由此可見，小波變換有以下特點(diǎn)：Ø 多尺度/多分辨的特點(diǎn)，可以由粗及細(xì)地處理信號(hào)；Ø 可以看成用基本頻率特性為()的帶通濾波器在不同尺度a下對(duì)信號(hào)做濾波。Ø 適當(dāng)?shù)剡x擇小波，使(t

16、)在時(shí)域上為有限支撐,Y()在頻域上也比較集中，就可以使WT在時(shí)、頻域都具有表征信號(hào)局部特征的能力。2.2離散小波變換6計(jì)算機(jī)中的圖像信息是以離散信號(hào)形式存放的，所以需要將連續(xù)小波變換離散化。而最基本的離散化方法就是二進(jìn)制離散，一般將這種經(jīng)過(guò)離散化的小波及其變換叫做二進(jìn)小波和二進(jìn)變換。需要注意的是這里的離散化都是針對(duì)連續(xù)的尺度因a和連續(xù)平移因子b的，而不是針對(duì)時(shí)間t的。這兒限制尺度因子a總是正數(shù)。（1）尺度與位移的離散化對(duì)連續(xù)小波基函數(shù)尺度因子和平移因子進(jìn)行離散化可以得到離散小波變換，從而減少小波變換系數(shù)的冗余度。在離散化時(shí)通常對(duì)尺度因子和平移因子按冪級(jí)數(shù)進(jìn)行離散化，即?。檎麛?shù)，但一般都假定

17、），得到離散小波函數(shù)為：其對(duì)應(yīng)系數(shù)為：（2）二進(jìn)制小波變換二進(jìn)小波變換是一種特殊的離散小波變換，特別地令參數(shù)，則有。該二進(jìn)尺度分解的原理在二十世紀(jì)三十年代由 Littlewood 和 Paley 在數(shù)學(xué)上進(jìn)行了研究證明。離散小波變換為：離散二進(jìn)小波變換為：2.3 多分辨率分析7Mallat在構(gòu)造正交小波基時(shí)提出了多分辨率分析（Multi-Resolution Analysis）的概念，從空間概念上形象地說(shuō)明了小波的多分辨率特性，并將在此之前的所有正交小波基的構(gòu)造法統(tǒng)一起來(lái)，給出了正交小波的構(gòu)造方法以及正交小波的快速算法Mallat算法。Mallat算法在小波分析中的地位相當(dāng)于快速傅立葉變換在經(jīng)

18、典傅立葉分析中的地位。多分辨率分析可形象地表示為一組嵌套的多分辨率子空間（如圖2所示）。W1W2W3V3圖2嵌套的多分辨率子空間假設(shè)原信號(hào)的頻率空間為，經(jīng)第一級(jí)分解后被分解成兩個(gè)子空間：低頻的和高頻的；經(jīng)第二級(jí)分解后被分解成低頻的和高頻的。這種子空間的分解過(guò)程可以記為：其中符號(hào)表示兩個(gè)子空間的“正交和”；代表與分辨率對(duì)應(yīng)的多分辨率分析子空間；與尺度函數(shù)相對(duì)應(yīng)的小波函數(shù)的伸縮和平移構(gòu)成的矢量空間是的正交補(bǔ)空間；各是反映空間信號(hào)細(xì)節(jié)的高頻子空間，是反映空間信號(hào)概貌的低頻子空間。由離散小波框架可得到子空間的以下特性：這一結(jié)果表明：分辨率為20=1的多分辨率分析子空間可以用有限個(gè)子空間來(lái)逼近。3 小波

19、分析的應(yīng)用3.1 利用小波對(duì)信號(hào)進(jìn)行處理的一般步驟小波的應(yīng)用主要是信號(hào)的處理，其中最典型的應(yīng)用是小波圖象壓縮。另外，小波在諸如信號(hào)去噪、特征提取等多方面均有成功的應(yīng)用。下面以圖象去噪為例說(shuō)明小波應(yīng)用策略。小波的各種應(yīng)用均可分為以下三步7：1）取樣：這是一個(gè)預(yù)處理過(guò)程。取樣方法應(yīng)遵循取樣定理8。1）對(duì)原始信號(hào)作小波變換，將信號(hào)由空域變換到頻域；2）對(duì)小波系數(shù)做相應(yīng)處理；3）對(duì)處理后的小波系數(shù)做小波逆變換，重構(gòu)還原原信號(hào)。3.2小波圖像去噪因?yàn)樵肼曅盘?hào)多包含在具有較高頻率的細(xì)節(jié)中，所以小波去噪首先對(duì)圖像信號(hào)進(jìn)行小波分解，可利用門限閾值對(duì)所分解的小波系數(shù)進(jìn)行處理，然后對(duì)圖像信號(hào)進(jìn)行小波重構(gòu)，抑制圖像信號(hào)中的無(wú)用部分，恢復(fù)圖像信號(hào)中的有用部分。如圖3所示，具體步驟為9：1）圖像信號(hào)的小波分解：選擇合適的小波及恰當(dāng)?shù)姆纸鈱哟蜰，對(duì)目標(biāo)圖像進(jìn)行N層的小波分解；2）對(duì)分解后的高頻系數(shù)進(jìn)行閾值量化：對(duì)于分解的每一層，選擇恰當(dāng)?shù)拈撝?，?duì)該層高頻系數(shù)進(jìn)行閾值量化處理。利用軟閾值或硬閾值門限處理相應(yīng)的小波系數(shù), 獲得新的被壓縮的小波系數(shù)；3）重構(gòu)圖像：根據(jù)小波分解后的第N層近似的低頻系數(shù)和經(jīng)過(guò)閾值量化處理后的細(xì)節(jié)高頻系數(shù)