小波變換超清晰的理解

2021-09-25 10:02:24 字數 3083 閱讀 2803

小波變換超清晰的理解

**:從傅利葉變換到小波變換,並不是乙個完全抽象的東西,可以講得很形象。小波變換有著明確的物理意義,如果我們從它的提出時所面對的問題看起,可以整理出非常清晰的思路。

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

一、傅利葉變換

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

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

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

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

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

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

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

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

二、短時傅利葉變換(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. 我們知道對於突變訊號,傅利葉變換存在吉布斯效應,我們用無限長的三角函式怎麼也擬合不好突變訊號:

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

小波可以實現正交化,短時傅利葉變換不能。

以上,就是小波的意義。

小波變換 小波變換入門 haar小波

小波可以認為是乙個帶通濾波器,只允許頻率和小波基函式頻率相近的訊號通過。小波變換的基本思想是用一組小波函式和基函式表示乙個函式或者訊號。首先,以haar小波變換過程為例來理解小波變換。例 求只有4個畫素 9 7 3 5 的影象的哈爾小波變換係數。計算步驟如下 步驟2 求差值 differencing...

小波變換應用

改變塊變換編碼中所使用變換的種類 dft與dct 結果圖如下 dct變換結果優於dft。離散傅利葉變換dft是傅利葉變換在時域和頻域上都呈離散的形式。實際應用中通常採用快速傅利葉變換計算dft。離散余弦變換dct是與傅利葉變換相關的一種變換,dct相當於乙個長度大概是它兩倍的dft,但dct是對乙個...

小波變換與傅利葉變換

from 小波變換與傅利葉變換 分類 預設欄目 如果有人問我,如果傅利葉變換沒有學好 深入理解概念 是否能學好小波。答案是否定的。如果有人還問我,如果第一代小波變換沒學好,能否學好第二代小波變換。答案依然是否定的。但若你問我,沒學好傅利葉變換,能否操作 程式設計 小波變換,或是沒學好第一代小波,能否...