小波變換通俗解釋

2021-01-14 數學與人工智慧
從傅立葉變換到小波變換,並不是一個完全抽象的東西,可以講得很形象。小波變換有著明確的物理意義,如果我們從它的提出時所面對的問題看起,可以整理出非常清晰的思路。

下面就按照傅立葉-->短時傅立葉變換-->小波變換的順序,講一下為什麼會出現小波這個東西、小波究竟是怎樣的思路。


關於傅立葉變換的基本概念在此我就不再贅述了,默認大家現在正處在理解了傅立葉但還沒理解小波的道路上。

下面我們主要將傅立葉變換的不足。即我們知道傅立葉變化可以分析信號的頻譜,那麼為什麼還要提出小波變換?答案「對非平穩過程,傅立葉變換有局限性」。看如下一個簡單的信號:

做完FFT(快速傅立葉變換)後,可以在頻譜上看到清晰的四條線,信號包含四個頻率成分。

一切沒有問題。但是,如果是頻率隨著時間變化的非平穩信號呢?

如上圖,最上邊的是頻率始終不變的平穩信號。而下邊兩個則是頻率隨著時間改變的非平穩信號,它們同樣包含和最上信號相同頻率的四個成分。做FFT後,我們發現這三個時域上有巨大差異的信號,頻譜(幅值譜)卻非常一致。尤其是下邊兩個非平穩信號,我們從頻譜上無法區分它們,因為它們包含的四個頻率的信號的成分確實是一樣的,只是出現的先後順序不同。

可見,傅立葉變換處理非平穩信號有天生缺陷。它只能獲取一段信號總體上包含哪些頻率的成分,但是對各成分出現的時刻並無所知。因此時域相差很大的兩個信號,可能頻譜圖一樣。

然而平穩信號大多是人為製造出來的,自然界的大量信號幾乎都是非平穩的,所以在比如生物醫學信號分析等領域的論文中,基本看不到單純傅立葉變換這樣naive的方法。

上圖所示的是一個正常人的事件相關電位。對於這樣的非平穩信號,只知道包含哪些頻率成分是不夠的,我們還想知道各個成分出現的時間。知道信號頻率隨時間變化的情況,各個時刻的瞬時頻率及其幅值——這也就是時頻分析。

短時傅立葉變換

(Short-time Fourier Transform,STFT)

一個簡單可行的方法就是——加窗。「把整個時域過程分解成無數個等長的小過程,每個小過程近似平穩,再傅立葉變換,就知道在哪個時間點上出現了什麼頻率了。」這就是短時傅立葉變換。

看圖:

時域上分成一段一段做FFT,不就知道頻率成分隨著時間的變化情況了嗎!

用這樣的方法,可以得到一個信號的時頻圖了:

圖上既能看到10Hz, 25 Hz, 50 Hz, 100 Hz四個頻域成分,還能看到出現的時間。兩排峰是對稱的,所以大家只用看一排就行了。

是不是棒棒的?時頻分析結果到手。但是STFT依然有缺陷。
使用STFT存在一個問題,我們應該用多寬的窗函數?
窗太寬太窄都有問題:

窗太窄,窗內的信號太短,會導致頻率分析不夠精準,頻率解析度差。窗太寬,時域上又不夠精細,時間解析度低。

(這裡插一句,這個道理可以用海森堡不確定性原理來解釋。類似於我們不能同時獲取一個粒子的動量和位置,我們也不能同時獲取信號絕對精準的時刻和頻率。這也是一對不可兼得的矛盾體。我們不知道在某個瞬間哪個頻率分量存在,我們知道的只能是在一個時間段內某個頻帶的分量存在。所以絕對意義的瞬時頻率是不存在的。)

看看實例效果吧:

上圖對同一個信號(4個頻率成分)採用不同寬度的窗做STFT,結果如右圖。用窄窗,時頻圖在時間軸上解析度很高,幾個峰基本成矩形,而用寬窗則變成了綿延的矮山。但是頻率軸上,窄窗明顯不如下邊兩個寬窗精確。

所以窄窗口時間解析度高、頻率解析度低,寬窗口時間解析度低、頻率解析度高。對於時變的非穩態信號,高頻適合小窗口,低頻適合大窗口。然而STFT的窗口是固定的,在一次STFT中寬度不會變化,所以STFT還是無法滿足非穩態信號變化的頻率的需求。



那麼你可能會想到,讓窗口大小變起來,多做幾次STFT不就可以了嗎?!沒錯,小波變換就有著這樣的思路。

但事實上小波並不是這麼做的(有人認為「小波變換就是根據算法,加不等長的窗,對每一小部分進行傅立葉變換」,這是不準確的。小波變換並沒有採用窗的思想,更沒有做傅立葉變換。)

至於為什麼不採用可變窗的STFT呢,我認為是因為這樣做冗餘會太嚴重,STFT做不到正交化,這也是它的一大缺陷。

於是小波變換的出發點和STFT還是不同的。STFT是給信號加窗,分段做FFT;而小波直接把傅立葉變換的基給換了——將無限長的三角函數基換成了有限長的會衰減的小波基。這樣不僅能夠獲取頻率,還可以定位到時間了~

【解釋】
來我們再回顧一下傅立葉變換吧,沒弄清傅立葉變換為什麼能得到信號各個頻率成分的同學也可以再借我的圖理解一下。
傅立葉變換把無限長的三角函數作為基函數:

這個基函數會伸縮、會平移(其實是兩個正交基的分解)。縮得窄,對應高頻;伸得寬,對應低頻。然後這個基函數不斷和信號做相乘。某一個尺度(寬窄)下乘出來的結果,就可以理解成信號所包含的當前尺度對應頻率成分有多少。於是,基函數會在某些尺度下,與信號相乘得到一個很大的值,因為此時二者有一種重合關係。那麼我們就知道信號包含該頻率的成分的多少。

仔細體會可以發現,這一步其實是在計算信號和三角函數的相關性。

看,這兩種尺度能乘出一個大的值(相關度高),所以信號包含較多的這兩個頻率成分,在頻譜上這兩個頻率會出現兩個峰。

以上,就是粗淺意義上傅立葉變換的原理。
如前邊所說,小波做的改變就在於,將無限長的三角函數基換成了有限長的會衰減的小波基。

這就是為什麼它叫「小波」,因為是很小的一個波嘛~

從公式可以看出,不同於傅立葉變換,變量只有頻率ω,小波變換有兩個變量:尺度a(scale)和平移量 τ(translation)。尺度a控制小波函數的伸縮,平移量 τ控制小波函數的平移。尺度就對應於頻率(反比),平移量 τ就對應於時間。

當伸縮、平移到這麼一種重合情況時,也會相乘得到一個大的值。這時候和傅立葉變換不同的是,這不僅可以知道信號有這樣頻率的成分,而且知道它在時域上存在的具體位置。

而當我們在每個尺度下都平移著和信號乘過一遍後,我們就知道信號在每個位置都包含哪些頻率成分。

看到了嗎?有了小波,我們從此再也不害怕非穩定信號啦!從此可以做時頻分析啦!
做傅立葉變換隻能得到一個頻譜,做小波變換卻可以得到一個時頻譜!

↑:時域信號

↑:傅立葉變換結果

↑:小波變換結果
小波還有一些好處:
1. 我們知道對於突變信號,傅立葉變換存在吉布斯效應,我們用無限長的三角函數怎麼也擬合不好突變信號:

然而衰減的小波就不一樣了:

2. 小波可以實現正交化,短時傅立葉變換不能。
以上,就是小波的意義。
--
以上只是用形象地給大家展示了一下小波的思想,希望能對大家的入門帶來一些幫助。畢竟如果對小波一無所知,直接去看那些堆砌公式、照搬論文語言的教材,一定會痛苦不堪。

在這裡推薦幾篇入門讀物,都是以感性介紹為主,易懂但並不深入,對大家初步理解小波會很有幫助。文中有的思路和圖也選自於其中:

1. THE WAVELET TUTORIAL (強烈推薦,搜索:INDEX TO SERIES OFTUTORIALS TO WAVELET TRANSFORM BY ROBI POLIKAR)(http://users.rowan.edu/~polikar/WAVELETS/WTtutorial.html)
2. WAVELETS:SEEING THE FOREST AND THE TREES
3. A Really Friendly Guide to Wavelets
4. Conceptual wavelets

但是真正理解透小波變換,這些還差得很遠。比如你至少還要知道有一個「尺度函數」的存在,它是構造「小波函數」的關鍵,並且是它和小波函數一起才構成了小波多解析度分析,理解了它才有可能利用小波做一些數位訊號處理;你還要理解離散小波變換、正交小波變換、二維小波變換、小波包……這些內容國內教材上講得也很糟糕,大家就一點一點啃吧~

相關焦點

  • 小波變換通俗解釋版
    小波變換有著明確的物理意義,如果我們從它的提出時所面對的問題看起,可以整理出非常清晰的思路。下面就按照傅立葉-->短時傅立葉變換-->小波變換的順序,講一下為什麼會出現小波這個東西、小波究竟是怎樣的思路。
  • 小波變換 完美通俗解讀
    我並不claim我會把整個小波變換講清楚,這是不可能的事,我只能盡力去圍繞要點展開,比如小波變換相對傅立葉變換的好處,這些好處的原因是什麼,小波變換的幾個根本性質是什麼,背後的推導是什麼。我希望達到的目的就是一個小波變換的初學者在看完這個系列之後,就能用matlab或者別的工具對信號做小波變換的基本分析並且知道這個分析大概是怎麼回事。
  • 完美通俗解讀小波變換,終於搞懂小波是什麼東東了
    恩,以上就是通常情況下你能在國內網站上搜到的小波變換文章告訴你的。但為什麼呢?這是我希望在這個系列文章中講清楚的。不過在這篇文章裡,我先點到為止,把小波變換的重要特性以及優點cover了,在下一篇文章中再具體推導這些特性。小波變換的本質和傅立葉變換類似,也是用精心挑選的basis來表示信號方程。
  • 小波變換原理與應用_小波變換的基本原理_小波變換的應用
    打開APP 小波變換原理與應用_小波變換的基本原理_小波變換的應用 發表於 2018-01-13 09:42:37      小波變換的基本原理   傳統的信號理論,是建立在Fourier分析基礎上的,而Fourier變換作為一種全局性的變化,其有一定的局限性。在實際應用中人們開始對Fourier變換進行各種改進,小波分析由此產生了。
  • 小波變換和motion信號處理(二)
    這是《小波變換和motion信號處理》系列的第二篇,深入小波。第一篇我進行了基礎知識的鋪墊,第三篇主要講解應用。是父小波。需要提醒一點的是,這個正交純粹是為了小波分析的方便而引入的特性,並不是說小波變換的基就一定必須是正交的。但大部分小波變換的基確實是正交的,所以本文就直接默認正交為小波變換的主要性質之一了。
  • 小波變換教程(一):為什麼需要小波變換
    本文講解ROBI POLIKAR編寫的小波變換教程《THE WAVELET TUTORIAL》的第一部分:為什麼需要小波變換。
  • Matlab傅立葉變換、餘弦變換和小波變換
    圖像小波變換的 Matlab 實現函數3.1 一維小波變換的 Matlab 實現(1) dwt 函數 Matlab功能:一維離散小波變換格式:[cA,cD]=dwt(X,'wname')[cA,cD]=dwt(X,Lo_D,Hi_D)別可以實現一維、二維和 N 維 DFT說明:[cA,cD]=dwt(X,'wname') 使用指定的小波基函數 'wname
  • 小波的秘密1-小波變換概況與綜述
    變換就是把一個信號分解成一系列的小波(短時間),也許就會問,小波變換的小波是什麼啊,定義中就是告訴我們小波,因為這個小波實在是太多,一個是種類多,還有就是同一種小波還可以尺度變換,但是小波在整個時間範圍的幅度平均值是0,具有有限的持續時間和突變的頻率和振幅,可以是不規則,也可以是不對稱,很明顯正弦波就不是小波,什麼的是呢,看下面幾個圖就是:
  • 小波變換教程(1):基本原理
    小波變換是一個相對較新的概念(其出現大約是在20世紀80年代),但是有關於它的文章和書籍卻不少。這其中大部分都是由數學專業人士寫給其他同行看的,不過,仍然有大量數學專家不知道其他同行們討論的是什麼(我的一個數學教授就承認過)。換言之,大多數介紹小波變換的文獻對那些小波新手們來說用處不大(此為個人觀點)。
  • 小波變換和motion信號處理:第二篇
    這是《小波變換和motion信號處理》系列的第二篇,深入小波。第一篇我進行了基礎知識的鋪墊,第三篇主要講解應用。需要提醒一點的是,這個正交純粹是為了小波分析的方便而引入的特性,並不是說小波變換的基就一定必須是正交的。但大部分小波變換的基確實是正交的,所以本文就直接默認正交為小波變換的主要性質之一了。引入這個父小波呢,主要是為了方便做多解析度分析(multiresolution analysis, MRA)。說到這裡,你的問題可能會井噴了:好好的為什麼出來一個父小波呢?
  • 基於FPGA實現多種小波變換
    引 言 基於提升框架的小波變換方法,利用FPGA 可編程特性可實現多種小波變換。提升框架(LS :Lifting Scheme) 是由Sweldens 等人在近幾年提出的一種小波變換方法,用它的框架結構能有效地計算DWT。對於較長的濾波器,LS 的操作次數比濾波器組的操作方式減少將近一半,更適合硬體實現。作者根據提升小波變換的框架式結構,利用FPGA 可完全重構的特點構造不同的小波變換核,以滿足不同應用場合的要求。
  • 圖像的二維提升小波變換的FPGA實現
    小波分析理論以其良好的時頻區域性和多解析度分析能力,開闢了圖像處理的嶄新領域。小波變換是一種很好的圖像分解方法,非常適合於分析突變信號而用於靜止圖像邊緣的提取和壓縮。
  • 基於SI濾波器的一種小波變換的實現
    摘要:文中在應用對數域電路的基礎上,提出了一種新型的連續小波變換方法,它通過對母小波的一種數值逼近得到小波函數的有理公式,並以Marr小波為例來模擬這個逼近過程,並用Matlab對逼近過程進行仿真。
  • 基於提升小波變換的模糊圖像融合算法研究
    圖像的小波分解是一種像素級圖像融合普遍採用的方法,利用小波非冗餘性使圖像經小波分解後數據量不會增大;同時小波分解具有方向性的特性,可針對人眼對不同方向的高頻分量具有不同解析度這一視覺特性,在圖像融合時獲得效果更佳的融合圖像。基於提升方法的小波變換又稱第2代小波變換,它既保持傳統的小波時頻局部化等特性,又克服其局限性。
  • 基於FPGA 的多用途提升小波變換核
    引 言本文引用地址:http://www.eepw.com.cn/article/189802.htm基於提升框架的小波變換方法,利用FPGA 可編程特性可實現多種小波變換。
  • 高斯類小波變換的開關電流頻域法實現
    本文提出了一個以頻域中的高斯函數單元為核心的共享結構系統實現3種高斯類小波變換。在此頻域共享結構實現方案中,復用頻域高斯函數單元採用開關電流電路實現頻域高斯類小波變換系統,不同尺度上的高斯類小波變換可通過調節開關電流電路的時鐘頻率獲得,所提出的頻域共享結構高斯類小波變換系統適合於製成通用型小波變換晶片。
  • 小波變換和motion信號處理(一)
    我們實驗室主要是搞圖像的,實力在全國也是很強的,進去後和師兄師姐聊,大家都在搞什麼小波變換,H264之類的。當時的我心思都不在這方面,盡搞什麼作業系統移植,ARM+FPGA這些東西了。對小波變換的認識也就停留在神秘的「圖像視頻壓縮算法之王」上面。後來我才發現,在別的很廣泛的領域中,小波也逐漸開始流行。比如話說很早以前,我們接觸的信號頻域處理基本都是傅立葉和拉普拉斯的天下。
  • 完全搞懂傅立葉變換和小波(1)——總綱
    本文引用地址:http://www.eepw.com.cn/article/201702/344594.htm  完全搞懂傅立葉變換和小波,你至少需要知道哪些預備知識?主頁君從今天開始就將通過一些列文章告訴你他們之間的來龍去脈!本節是全部系列文章的第一節——總綱,日後我們也將按照這個思路一點一點講述所有的知識。
  • 二維9/7小波變換VLSI設計
    1 引言 JPEG2000標準與JEPG標準最顯著的不同點是選用以小波變換為主的多解析度編碼方法。JPEG2000標準給出兩種雙正交小波濾波器,即有損壓縮和無損壓縮,前者採用CDF9/7小波,後者採用5/3小波。CDF9/7小波是圖像壓縮的首選濾波器,自然圖像壓縮性能好於5/3小波。 目前圖像二維離散小波變換的硬體實現方法有兩種:一種是採用分離的一維小波變換,另一種是設計非離散的二維小波變換結構。
  • 基於小波變換的圖像壓縮算法改進研究
    1 小波變換1.1 小波變換的產生及原理  儘管傅立葉(Fourier)變換可以確切地告訴人們某個信號是否包含特定的頻率分量,但它無法說明該頻率分量發生在哪個時間段。因此,它僅適用於處理平穩信號,而不適用於處理非平穩信號。