傅立葉分析和小波分析之間的關(guān)系之通俗終極版

1、傅立葉分析和小波分析之間的關(guān)系之通俗終極版轉(zhuǎn)載請(qǐng)注明出處和作者知乎作者:咚懂咚懂咚從傅里葉變換到小波變換，并不是一個(gè)完全抽象的東西，完全可以講得很形象。小波變換有著明確的物理意義，如果我們從它的提出時(shí)所面對(duì)的問(wèn)題看起，可以整理出非常清晰的思路。下面我就按照傅里葉-短時(shí)傅里葉變換-小波變換的順序，講一下為什么會(huì)出現(xiàn)小波這個(gè)東西、小波究竟是怎樣的思路。（反正題主要求的是通俗形象，沒(méi)說(shuō)簡(jiǎn)短，希望不會(huì)太長(zhǎng)不看。）一、傅里葉變換關(guān)于傅里葉變換的基本概念在此我就不再贅述了，默認(rèn)大家現(xiàn)在正處在理解了傅里葉但還沒(méi)理解小波的道路上。（在第三節(jié)小波變換的地方我會(huì)再形象地講一下傅里葉變換）下面我們主要將傅里葉變換的

2、不足。即我們知道傅里葉變化可以分析信號(hào)的頻譜，那么為什么還要提出小波變換？答案就是方沁園所說(shuō)的，“對(duì)非平穩(wěn)過(guò)程，傅里葉變換有局限性”。看如下一個(gè)簡(jiǎn)單的信號(hào)：做完FFT（快速傅里葉變換）后，可以在頻譜上看到清晰的四條線，信號(hào)包含四個(gè)頻率成分。一切沒(méi)有問(wèn)題。但是，如果是非平穩(wěn)信號(hào)呢？如上圖，最上邊的是頻率始終不變的平穩(wěn)信號(hào)。而下邊兩個(gè)則是頻率隨著時(shí)間改變的非平穩(wěn)信號(hào)，它們同樣包含和最上信號(hào)相同頻率的四個(gè)成分。做FFT后，我們發(fā)現(xiàn)這三個(gè)時(shí)域上有巨大差異的信號(hào)，頻譜卻非常一致。尤其是下邊兩個(gè)非平穩(wěn)信號(hào)，我們從頻域上無(wú)法區(qū)分它們，因?yàn)樗鼈儼乃膫€(gè)頻率的信號(hào)的成分確實(shí)是一樣的，只是出現(xiàn)的先后順序不同?？?/p>

3、見(jiàn)，傅里葉變換處理非平穩(wěn)信號(hào)有天生缺陷。它只能獲取一段信號(hào)總體上包含哪些頻率的成分，但是對(duì)各成分出現(xiàn)的時(shí)刻并無(wú)所知。因此時(shí)域相差很大的兩個(gè)信號(hào)，可能頻譜圖一樣。然而平穩(wěn)信號(hào)大多是人為制造出來(lái)的，自然界的大量信號(hào)幾乎都是非平穩(wěn)的，所以在比如生物醫(yī)學(xué)信號(hào)分析等領(lǐng)域的papers中，基本看不到單純傅里葉變換這樣naive的方法。上圖所示的是一個(gè)正常人的事件相關(guān)電位。對(duì)于這樣的非平穩(wěn)信號(hào)，只知道包含哪些頻率成分是不夠的，我們還想知道各個(gè)成分出現(xiàn)的時(shí)間。知道信號(hào)頻率隨時(shí)間變化的關(guān)系，各個(gè)時(shí)刻的瞬時(shí)頻率及其幅值這也就是時(shí)頻分析。二、短時(shí)傅里葉變換（Short-time Fourier Transform,

4、 STFT）一個(gè)簡(jiǎn)單可行的方法就是加窗。我又要套用方沁園同學(xué)的描述了，“把整個(gè)時(shí)域過(guò)程分解成無(wú)數(shù)個(gè)等長(zhǎng)的小過(guò)程，每個(gè)小過(guò)程近似平穩(wěn)，再傅里葉變換，就知道在哪個(gè)時(shí)間點(diǎn)上出現(xiàn)了什么頻率了。”這就是短時(shí)傅里葉變換。看圖：時(shí)域上分成一段一段做FFT，不就知道頻率成分隨著時(shí)間的變化情況了嗎！用這樣的方法，可以得到一個(gè)信號(hào)的時(shí)頻圖了：圖上既能看到300Hz, 200 Hz, 100 Hz, 50 Hz四個(gè)頻域成分，還能看到出現(xiàn)的時(shí)間。兩排峰是對(duì)稱(chēng)的，所以大家只用看一排就行了。是不是棒棒噠？但是STFT依然有缺陷。使用STFT存在一個(gè)問(wèn)題，我們應(yīng)該用多寬的窗函數(shù)？窗太寬太窄都有問(wèn)題：窗太窄，窗內(nèi)的信號(hào)太短，

5、會(huì)導(dǎo)致頻率分析不夠精準(zhǔn)，頻率分辨率差。窗太寬，時(shí)域上又不夠精細(xì)，時(shí)間分辨率低。（這里插一句，這個(gè)道理可以用海森堡不確定性原理來(lái)解釋。類(lèi)似于我們不能同時(shí)獲取一個(gè)粒子的動(dòng)量和位置，我們也不能同時(shí)獲取信號(hào)絕對(duì)精準(zhǔn)的時(shí)刻和頻率。這也是一對(duì)不可兼得的矛盾體。我們不知道在某個(gè)瞬間哪個(gè)頻率分量存在，我們知道的只能是在一個(gè)時(shí)間段內(nèi)某個(gè)頻帶的分量存在。 所以絕對(duì)意義的瞬時(shí)頻率是不存在的。）看看實(shí)例效果吧：上圖對(duì)同一個(gè)信號(hào)（4個(gè)頻率成分）采用不同寬度的窗做STFT，結(jié)果如右圖。用窄窗，時(shí)頻圖在時(shí)間軸上分辨率很高，幾個(gè)峰基本成矩形，而用寬窗則變成了綿延的矮山。但是頻率軸上，窄窗明顯不如下邊兩個(gè)寬窗精確。所以窄窗口時(shí)

6、間分辨率高、頻率分辨率低，寬窗口時(shí)間分辨率低、頻率分辨率高。對(duì)于時(shí)變的非穩(wěn)態(tài)信號(hào)，高頻適合小窗口，低頻適合大窗口。然而STFT的窗口是固定的，在一次STFT中寬度不會(huì)變化，所以STFT還是無(wú)法滿足非穩(wěn)態(tài)信號(hào)變化的頻率的需求。三、小波變換那么你可能會(huì)想到，讓窗口大小變起來(lái)，多做幾次STFT不就可以了嗎？！沒(méi)錯(cuò)，小波變換就有著這樣的思路。但事實(shí)上小波并不是這么做的（關(guān)于這一點(diǎn)，方沁園同學(xué)的表述“小波變換就是根據(jù)算法，加不等長(zhǎng)的窗，對(duì)每一小部分進(jìn)行傅里葉變換”就不準(zhǔn)確了。小波變換并沒(méi)有采用窗的思想，更沒(méi)有做傅里葉變換。）至于為什么不采用可變窗的STFT呢，我認(rèn)為是因?yàn)檫@樣做冗余會(huì)太嚴(yán)重，STFT做不

7、到正交化，這也是它的一大缺陷。于是小波變換的出發(fā)點(diǎn)和STFT還是不同的。STFT是給信號(hào)加窗，分段做FFT；而小波直接把傅里葉變換的基給換了將無(wú)限長(zhǎng)的三角函數(shù)基換成了有限長(zhǎng)的會(huì)衰減的小波基。這樣不僅能夠獲取頻率，還可以定位到時(shí)間了【解釋】來(lái)我們?cè)倩仡櫼幌赂道锶~變換吧，沒(méi)弄清傅里葉變換為什么能得到信號(hào)各個(gè)頻率成分的同學(xué)也可以再借我的圖理解一下。傅里葉變換把無(wú)限長(zhǎng)的三角函數(shù)作為基函數(shù)：這個(gè)基函數(shù)會(huì)伸縮、會(huì)平移?？s得窄，對(duì)應(yīng)高頻；伸得寬，對(duì)應(yīng)低頻。然后這個(gè)基函數(shù)不斷和信號(hào)做相乘。某一個(gè)尺度（寬窄）下乘出來(lái)的結(jié)果，就可以理解成信號(hào)所包含的當(dāng)前尺度對(duì)應(yīng)頻率成分有多少。于是，基函數(shù)會(huì)在某些尺度下，與信號(hào)相

8、乘得到一個(gè)很大的值，因?yàn)榇藭r(shí)二者有一種重合關(guān)系。那么我們就知道信號(hào)包含多少該頻率的成分。（看，這兩種尺度能乘出一個(gè)大的值，所以信號(hào)包含較多的這兩個(gè)頻率成分，在頻譜上這兩個(gè)頻率會(huì)出現(xiàn)兩個(gè)峰）以上，就是傅里葉變換的原理。如前邊所說(shuō)，小波做的改變就在于，將無(wú)限長(zhǎng)的三角函數(shù)基換成了有限長(zhǎng)的會(huì)衰減的小波基。這就是為什么它叫“小波”，因?yàn)槭呛苄〉囊粋€(gè)波嘛從公式可以看出，不同于傅里葉變換，變量只有頻率，小波變換有兩個(gè)變量：尺度a（scale）和平移量 （translation）。尺度a控制小波函數(shù)的伸縮，平移量控制小波函數(shù)的平移。尺度就對(duì)應(yīng)于頻率（反比），平移量就對(duì)應(yīng)于時(shí)間。當(dāng)伸縮、平移到這么一種重合情況時(shí)

9、，也會(huì)相乘得到一個(gè)大的值。這時(shí)候和傅里葉變換不同的是，這不僅可以知道信號(hào)有這樣頻率的成分，而且知道它在時(shí)域上存在的具體位置。而當(dāng)我們?cè)诿總€(gè)尺度下都平移著和信號(hào)乘過(guò)一遍后，我們就知道信號(hào)在每個(gè)位置都包含哪些頻率成分?？吹搅藛幔坑辛诵〔?，我們從此再也不害怕非穩(wěn)定信號(hào)啦！從此可以做時(shí)頻分析啦！做傅里葉變換只能得到一個(gè)頻譜，做小波變換卻可以得到一個(gè)時(shí)頻譜！：時(shí)域信號(hào)：傅里葉變換結(jié)果：小波變換結(jié)果小波還有一些好處：1. 我們知道對(duì)于突變信號(hào)，傅里葉變換存在吉布斯效應(yīng)，我們用無(wú)限長(zhǎng)的三角函數(shù)怎么也擬合不好突變信號(hào)：然而衰減的小波就不一樣了：2. 小波可以實(shí)現(xiàn)正交化，短時(shí)傅里葉變換不能。以上，就是小波的意義

10、。-以上只是用形象地給大家展示了一下小波的思想，希望能對(duì)大家的入門(mén)帶來(lái)一些幫助。畢竟如果對(duì)小波一無(wú)所知，直接去看那些堆砌公式、照搬論文語(yǔ)言的教材，一定會(huì)痛苦不堪。在這里推薦幾篇入門(mén)讀物，都是以感性介紹為主，易懂但并不深入，對(duì)大家初步理解小波會(huì)很有幫助。文中有的思路和圖也選自于其中：1. THE WAVELET TUTORIAL （強(qiáng)烈推薦，點(diǎn)擊鏈接：INDEX TO SERIES OF TUTORIALS TO WAVELET TRANSFORM BY ROBI POLIKAR）2. WAVELETS：SEEING THE FOREST AND THE TREES3. A Really Friendly Guide to Wavelets4. Conceptual wavelets但是真正理解透小波變換，這些還差得很遠(yuǎn)。比如你至少還要知道有一個(gè)“尺度函數(shù)”的存在，它是構(gòu)造“小波函數(shù)”的