快速傅立葉變換的意義及應用

快速傅立葉變換的意義及應用

 1、爲什麼要進行傅里葉變換,其物理意義是什麼?

 

傅立葉變換是數字信號處理領域一種很重要的算法。要知道傅立葉變換算法的意義,首先要了解傅立葉原理的意義。傅立葉原理表明:任何連續測量的時序或信號,都可以表示爲不同頻率的正弦波信號的無限疊加。而根據該原理創立的傅立葉變換算法利用直接測量到的原始信號,以累加方式來計算該信號中不同正弦波信號的頻率、振幅和相位。

 

和傅立葉變換算法對應的是反傅立葉變換算法。該反變換從本質上說也是一種累加處理,這樣就可以將單獨改變的正弦波信號轉換成一個信號。

 

因此,可以說,傅立葉變換將原來難以處理的時域信號轉換成了易於分析的頻域信號(信號的頻譜),可以利用一些工具對這些頻域信號進行處理、加工。最後還可以利用傅立葉反變換將這些頻域信號轉換成時域信號。

 

從現代數學的眼光來看,傅里葉變換是一種特殊的積分變換。它能將滿足一定條件的某個函數表示成正弦基函數的線性組合或者積分。在不同的研究領域,傅里葉變換具有多種不同的變體形式,如連續傅里葉變換和離散傅里葉變換。

 

在數學領域,儘管最初傅立葉分析是作爲熱過程的解析分析的工具,但是其思想方法仍然具有典型的還原論和分析主義的特徵。"任意"的函數通過一定的分解,都能夠表示爲正弦函數的線性組合的形式,而正弦函數在物理上是被充分研究而相對簡單的函數類:1. 傅立葉變換是線性算子,若賦予適當的範數,它還是酉算子;2. 傅立葉變換的逆變換容易求出,而且形式與正變換非常類似;3. 正弦基函數是微分運算的本徵函數,從而使得線性微分方程的求解可以轉化爲常係數的代數方程的求解.在線性時不變雜的卷積運算爲簡單的乘積運算,從而提供了計算卷積的一種簡單手段;5.  離散形式的傅立葉的物理系統內,頻率是個不變的性質,從而系統對於複雜激勵的響應可以通過組合其對不同頻率正弦信號的響應來獲取;4. 著名的卷積定理指出:傅立葉變換可以化復變換可以利用數字計算機快速的算出(其算法稱爲快速傅立葉變換算法(FFT))。

 

正是由於上述的良好性質,傅里葉變換在物理學、數論、組合數學、信號處理、概率、統計、密碼學、聲學、光學等領域都有着廣泛的應用。

 

2、圖像傅立葉變換的物理意義

 

圖像的頻率是表徵圖像中灰度變化劇烈程度的指標,是灰度在平面空間上的梯度。如:大面積的沙漠在圖像中是一片灰度變化緩慢的區域,對應的頻率值很低;而對於地表屬性變換劇烈的邊緣區域在圖像中是一片灰度變化劇烈的區域,對應的頻率值較高。傅立葉變換在實際中有非常明顯的物理意義,設f是一個能量有限的模擬信號,則其傅立葉變換就表示f的譜。從純粹的數學意義上看,傅立葉變換是將一個函數轉換爲一系列周期函數來處理的。從物理效果看,傅立葉變換是將圖像從空間域轉換到頻率域,其逆變換是將圖像從頻率域轉換到空間域。換句話說,傅立葉變換的物理意義是將圖像的灰度分佈函數變換爲圖像的頻率分佈函數,傅立葉逆變換是將圖像的頻率分佈函數變換爲灰度分佈函數

 

傅立葉變換以前,圖像(未壓縮的位圖)是由對在連續空間(現實空間)上的採樣得到一系列點的集合,我們習慣用一個二維矩陣表示空間上各點,則圖像可由z=f(x,y)來表示。由於空間是三維的,圖像是二維的,因此空間中物體在另一個維度上的關係就由梯度來表示,這樣我們可以通過觀察圖像得知物體在三維空間中的對應關係。爲什麼要提梯度?因爲實際上對圖像進行二維傅立葉變換得到頻譜圖,就是圖像梯度的分佈圖,當然頻譜圖上的各點與圖像上各點並不存在一一對應的關係,即使在不移頻的情況下也是沒有。傅立葉頻譜圖上我們看到的明暗不一的亮點,實際上圖像上某一點與鄰域點差異的強弱,即梯度的大小,也即該點的頻率的大小(可以這麼理解,圖像中的低頻部分指低梯度的點,高頻部分相反)。一般來講,梯度大則該點的亮度強,否則該點亮度弱。這樣通過觀察傅立葉變換後的頻譜圖,也叫功率圖,我們首先就可以看出,圖像的能量分佈,如果頻譜圖中暗的點數更多,那麼實際圖像是比較柔和的(因爲各點與鄰域差異都不大,梯度相對較小),反之,如果頻譜圖中亮的點數多,那麼實際圖像一定是尖銳的,邊界分明且邊界兩邊像素差異較大的。對頻譜移頻到原點以後,可以看出圖像的頻率分佈是以原點爲圓心,對稱分佈的。將頻譜移頻到圓心除了可以清晰地看出圖像頻率分佈以外,還有一個好處,它可以分離出有週期性規律的干擾信號,比如正弦干擾,一副帶有正弦干擾,移頻到原點的頻譜圖上可以看出除了中心以外還存在以某一點爲中心,對稱分佈的亮點集合,這個集合就是干擾噪音產生的,這時可以很直觀的通過在該位置放置帶阻濾波器消除干擾

 

 

快速傅氏變換 英文名是fast Fourier transform

快速傅氏變換(FFT)是離散傅氏變換(DFT)的快速算法,它是根據離散傅氏變換的奇、偶、虛、實等特性,對離散傅立葉變換的算法進行改進獲得的。它對傅氏變換的理論並沒有新的發現,但是對於在計算機系統或者說數字系統中應用離散傅立葉變換,可以說是進了一大步。

設x(n)爲N項的複數序列,由DFT變換,任一X(m)的計算都需要N次複數乘法和N-1次複數加法,而一次複數乘法等於四次實數乘法和兩次實數加法,一次複數加法等於兩次實數加法,即使把一次複數乘法和一次複數加法定義成一次“運算”(四次實數乘法和四次實數加法),那麼求出N項複數序列的X(m),即N點DFT變換大約就需要N2次運算。當N=1024點甚至更多的時候,需要N2=1048576次運算,在FFT中,利用WN的週期性和對稱性,把一個N項序列(設N=2k,k爲正整數),分爲兩個N/2項的子序列,每個N/2點DFT變換需要(N/2)2次運算,再用N次運算把兩個N/2點的DFT變換組合成一個N點的DFT變換。這樣變換以後,總的運算次數就變成N+2(N/2)2=N+N2/2。繼續上面的例子,N=1024時,總的運算次數就變成了525312次,節省了大約50%的運算量。而如果我們將這種“一分爲二”的思想不斷進行下去,直到分成兩兩一組的DFT運算單元,那麼N點的DFT變換就只需要Nlog2N次的運算,N在1024點時,運算量僅有10240次,是先前的直接算法的1%,點數越多,運算量的節約就越大,這就是FFT的優越性。

 

小波分析 (Wavelet)

小波分析是當前數學中一個迅速發展的新領域,它同時具有理論深刻和應用十分廣泛的雙重意義。

小波變換的概念是由法國從事石油信號處理的工程師J.Morlet在1974年首先提出的,通過物理的直觀和信號處理的實際需要經驗的建立了反演公式,當時未能得到數學家的認可。正如1807年法國的熱學工程師J.B.J.Fourier提出任一函數都能展開成三角函數的無窮級數的創新概念未能得到著名數學家J.L.Lagrange,P.S.Laplace以及A.M.Legendre的認可一樣。幸運的是,早在七十年代,A.Calderon表示定理的發現、Hardy空間的原子分解和無條件基的深入研究爲小波變換的誕生做了理論上的準備,而且J.O.Stromberg還構造了歷史上非常類似於現在的小波基;1986年著名數學家Y.Meyer偶然構造出一個真正的小波基,並與S.Mallat合作建立了構造小波基的同意方法棗多尺度分析之後,小波分析纔開始蓬勃發展起來,其中比利時女數學家I.Daubechies撰寫的《小波十講(Ten Lectures on Wavelets)》對小波的普及起了重要的推動作用。它與Fourier變換、窗口Fourier變換(Gabor變換)相比,這是一個時間和頻率的局域變換,因而能有效的從信號中提取信息,通過伸縮和平移等運算功能對函數或信號進行多尺度細化分析(Multiscale Analysis),解決了Fourier變換不能解決的許多困難問題,從而小波變化被譽爲“數學顯微鏡”,它是調和分析發展史上里程碑式的進展。

小波(Wavelet)這一術語,顧名思義,“小波”就是小的波形。所謂“小”是指它具有衰減性;而稱之爲“波”則是指它的波動性,其振幅正負相間的震盪形式。與Fourier變換相比,小波變換是時間(空間)頻率的局部化分析,它通過伸縮平移運算對信號(函數)逐步進行多尺度細化,最終達到高頻處時間細分,低頻處頻率細分,能自動適應時頻信號分析的要求,從而可聚焦到信號的任意細節,解決了Fourier變換的困難問題,成爲繼Fourier變換以來在科學方法上的重大突破。有人把小波變換稱爲“數學顯微鏡”。

小波分析的應用是與小波分析的理論研究緊密地結合在一起地。現在,它已經在科技信息產業領域取得了令人矚目的成就。 電子信息技術是六大高新技術中重要的一個領域,它的重要方面是圖像和信號處理。現今,信號處理已經成爲當代科學技術工作的重要部分,信號處理的目的就是:準確的分析、診斷、編碼壓縮和量化、快速傳遞或存儲、精確地重構(或恢復)。從數學地角度來看,信號與圖像處理可以統一看作是信號處理(圖像可以看作是二維信號),在小波分析地許多分析的許多應用中,都可以歸結爲信號處理問題。現在,對於其性質隨實踐是穩定不變的信號,處理的理想工具仍然是傅立葉分析。但是在實際應用中的絕大多數信號是非穩定的,而特別適用於非穩定信號的工具就是小波分析。

小波分析是當前應用數學和工程學科中一個迅速發展的新領域,經過近10年的探索研究,重要的數學形式化體系已經建立,理論基礎更加紮實。與Fourier變換相比,小波變換是空間(時間)和頻率的局部變換,因而能有效地從信號中提取信息。通過伸縮和平移等運算功能可對函數或信號進行多尺度的細化分析,解決了Fourier變換不能解決的許多困難問題。小波變換聯繫了應用數學、物理學、計算機科學、信號與信息處理、圖像處理、地震勘探等多個學科。數學家認爲,小波分析是一個新的數學分支,它是泛函分析、Fourier分析、樣調分析、數值分析的完美結晶;信號和信息處理專家認爲,小波分析是時間—尺度分析和多分辨分析的一種新技術,它在信號分析、語音合成、圖像識別、計算機視覺、數據壓縮、地震勘探、大氣與海洋波分析等方面的研究都取得了有科學意義和應用價值的成果。

事實上小波分析的應用領域十分廣泛,它包括:數學領域的許多學科;信號分析、圖像處理;量子力學、理論物理;軍事電子對抗與武器的智能化;計算機分類與識別;音樂與語言的人工合成;醫學成像與診斷;地震勘探數據處理;大型機械的故障診斷等方面;例如,在數學方面,它已用於數值分析、構造快速數值方法、曲線曲面構造、微分方程求解、控制論等。在信號分析方面的濾波、去噪聲、壓縮、傳遞等。在圖像處理方面的圖像壓縮、分類、識別與診斷,去污等。在醫學成像方面的減少B超、CT、核磁共振成像的時間,提高分辨率等。

(1)小波分析用於信號與圖像壓縮是小波分析應用的一個重要方面。它的特點是壓縮比高,壓縮速度快,壓縮後能保持信號與圖像的特徵不變,且在傳遞中可以抗干擾。基於小波分析的壓縮方法很多,比較成功的有小波包最好基方法,小波域紋理模型方法,小波變換零樹壓縮,小波變換向量壓縮等。

(2)小波在信號分析中的應用也十分廣泛。它可以用於邊界的處理與濾波、時頻分析、信噪分離與提取弱信號、求分形指數、信號的識別與診斷以及多尺度邊緣檢測等。

(3)在工程技術等方面的應用。包括計算機視覺、計算機圖形學、曲線設計、湍流、遠程宇宙的研究與生物醫學方面。

 

 

 

 

傅立葉變換的基函數是sin cos函數,也就是說它只是用這2個去逼近原信號,對於原信號不一定都與正弦和餘弦函數相似,因此小波基函數的多樣性能使逼近的權重因子更加準確.正是因爲FFT的頻域局部性不強,而且時頻分開,所以纔有了以後的方法――當然包括了小波!

發表評論
所有評論
還沒有人評論,想成為第一個評論的人麼? 請在上方評論欄輸入並且點擊發布.
相關文章