卷積
維基百科,自由的百科全書
圖示兩個方形脈衝波的卷積。其中函數 "g" 首先對 反射,接着平移 "t" ,成爲 。那麼重疊部份的面積就相當於 "t" 處的卷積,其中橫座標代表待積變量 以及新函數 的自變量 "t" 。
圖示方形脈衝波和指數衰退的脈衝波的卷積(後者可能出現於 RC電路中),同樣地重疊部份面積就相當於 "t" 處的卷積。注意到因爲 "g" 是對稱的,所以在這兩張圖中,反射並不會改變它的形狀。
在泛函分析中,卷積(捲積)、旋積或摺積,是通過兩個函數f 和g 生成第三個函數的一種數學算子,表徵函數f與經過翻轉和平移的g 的重疊部分的累積。如果將參加卷積的一個函數看作區間的指示函數,卷積還可以被看作是“滑動平均”的推廣。
目錄
簡單介紹
卷積是分析數學中一種重要的運算。設: ,是上的兩個可積函數,作積分:
可以證明,關於幾乎所有的 ,上述積分是存在的。這樣,隨着 的不同取值,這個積分就定義了一個新函數,稱爲函數 與 的卷積,記爲。我們可以輕易驗證:,並且 仍爲可積函數。這就是說,把卷積代替乘法, 空間是一個代數,甚至是巴拿赫代數。
卷積與傅里葉變換有着密切的關係。例如兩函數的傅里葉變換的乘積等於它們卷積後的傅里葉變換,利用此一性質,能簡化傅里葉分析中的許多問題。
由卷積得到的函數 一般要比 和 都光滑。特別當 爲具有緊支集的光滑函數, 爲局部可積時,它們的卷積 也是光滑函數。利用這一性質,對於任意的可積函數 ,都可以簡單地構造出一列逼近於 的光滑函數列 ,這種方法稱爲函數的光滑化或正則化。
卷積的概念還可以推廣到數列、測度以及廣義函數上去。
定義
函數f 與g 的卷積記作,它是其中一個函數翻轉並平移後與另一個函數的乘積的積分,是一個對平移量的函數。
積分區間取決於f 與g 的定義域。
對於定義在離散域的函數,卷積定義爲
圖解卷積
- 首先將兩個函數都用來表示。
- 對其中一個函數做水平翻轉: →
- 加上一個時間偏移量,讓 能沿着 軸滑動。
- 讓t從-∞滑動到+∞。兩函數交會時,計算交會範圍中兩函數乘積的積分值。 換句話說,我們是在計算一個滑動的的加權平均值。也就是使用當做加權函數,來對取加權平均值。
- 最後得到的波形(未包含在此圖中)就是f和g的卷積。
如果f(t)是一個單位脈衝,我們得到的乘積就是g(t)本身,稱爲衝激響應 。
計算卷積的方法
當 爲有限長度 , 爲有限長度 的信號,計算卷積 有三種主要的方法,分別爲 1.直接計算(Direct Method) 2.快速傅里葉轉換(FFT) 和 3.分段卷積 (sectioned convolution)。方法1是直接利用定義來計算卷積,而方法2和3都是用到了FFT來快速計算卷積。也有不需要用到FFT的作法,如使用數論轉換。
方法1 直接計算
- 作法: 利用卷積的定義
- 若 和 皆爲實數信號,則需要 個乘法。
- 若 和 皆爲更一般性的複數信號,不使用複數乘法的快速算法,會需要 個乘法;但若使用複數乘法的快速算法,則可簡化至 個乘法。
- 因此,使用定義直接計算卷積的複雜度爲 。
方法2 快速傅里葉轉換(FFT)
- 概念:由於兩個離散信號在時域(time domain)做卷積相當於這兩個信號的離散傅里葉轉換在頻域(frequency domain)做相乘:
- ,可以看出在頻域的計算較簡單。
- 作法: 因此這個方法即是先將信號從時域轉成頻域:
- ,於是
- ,最後再將頻域信號轉回時域,就完成了卷積的計算:
- 總共做了 2 次 DFT 和 1 次 IDFT。
- 特別注意 DFT 和 IDFT 的點數 要滿足 。
- 由於 DFT 有快速算法 FFT,所以運算量爲
- 假設 點 DFT 的乘法量爲 , 和 爲一般性的複數信號,並使用複數乘法的快速算法,則共需要 個乘法。
方法3 分段卷積(sectioned convolution)
- 概念: 將 切成好幾段,每一段分別和 做卷積後,再將結果相加。
- 作法: 先將 切成每段長度爲 的區段 (),假設共切成S段:
-
- Section 1:
-
- Section 2:
-
- Section r:
-
- Section S:
- ,爲各個section的和
。
- 因此,
,
- 每一小段作卷積則是採用方法2,先將時域信號轉到頻域相乘,再轉回時域:
。
- 總共只需要做 點 FFT 次,因爲 只需要做一次FFT。
- 假設 點 DFT 的乘法量爲 , 和 爲一般性的複數信號,並使用複數乘法的快速算法,則共需要 個乘法。
- 運算量:
- 運算複雜度: ,和 呈線性,較方法2小。
應用時機
以上三種方法皆可用來計算卷積,其差別在於所需總體乘法量不同。基於運算量以及效率的考量,在計算卷積時,通常會選擇所需總體乘法量較少的方法。
以下根據 和 的長度( )分成5類,並列出適合使用的方法:
- 爲一非常小的整數 - 直接計算
- - 快速傅里葉轉換
- - 分段卷積
- - 快速傅里葉轉換
- 爲一非常小的整數 - 直接計算
基本上,以上只是粗略的分類。在實際應用時,最好還是算出三種方法所需的總乘法量,再選擇其中最有效率的方法來計算卷積。
例子
Q1: 當 ,適合用哪種方法計算卷積?
Ans:
- 方法1: 所需乘法量爲
- 方法2: ,而2016點的DFT最少乘法數 ,所以總乘法量爲
- 方法3:
- 若切成 8 塊(),則。選,則總乘法量爲 ,比方法1和2少了很多。
- 但是若要找到最少的乘法量,必須依照以下步驟
- (1)先找出 : 解 :
- (2)由 算出點數在 附近的DFT所需最少的乘法量,選擇DFT的點數
- (3)最後由 算出
- 因此,
- (1)由運算量對 的偏微分爲0而求出
- (2),所以選擇101點 DFT 附近點數乘法量最少的點數 或 。
- (3-1)當 ,總乘法量爲 。
- (3-2)當 ,總乘法量爲 。
- 由此可知,切成 20 塊會有較好的效率,而所需總乘法量爲 21480。
- 因此,當,所需總乘法量: 分段卷積 < 快速傅里葉轉換 < 直接計算。故,此時選擇使用分段卷積來計算卷積最適合。
Q2: 當 ,適合用哪種方法計算卷積?
Ans:
- 方法1: 所需乘法量爲
- 方法2: ,選擇1026點
DFT 附近點數乘法量最少的點數,。
-
- 因此,所需乘法量爲
-
- 方法3:
-
- (1)由運算量對 的偏微分爲0而求出
- (2),所以選擇7點 DFT 附近點數乘法量最少的點數 或 或 。
- (3-1)當 ,總乘法量爲 。
- (3-2)當 ,總乘法量爲 。
- (3-3)當 ,總乘法量爲 。
- 由此可知,切成 171 塊會有較好的效率,而所需總乘法量爲 5476。
-
- 因此,當,所需總乘法量: 分段卷積 < 直接計算 < 快速傅里葉轉換。故,此時選擇使用分段卷積來計算卷積最適合。
- 雖然當 是個很小的正整數時,大致上適合使用直接計算。但實際上還是將3個方法所需的乘法量都算出來,才能知道用哪種方法可以達到最高的效率。
Q3: 當 ,適合用哪種方法計算卷積?
Ans:
- 方法1: 所需乘法量爲
- 方法2: ,選擇1026點
DFT 附近點數乘法量最少的點數,。
-
- 因此,所需乘法量爲
-
- 方法3:
-
- (1)由運算量對 的偏微分爲0而求出
- (2),所以選擇1623點 DFT 附近點數乘法量最少的點數 。
- (3)當 ,總乘法量爲 。
- 由此可知,此時切成一段,就跟方法2一樣,所需總乘法量爲 44232。
-
- 因此,當,所需總乘法量: 快速傅里葉轉換 = 分段卷積 < 直接計算。故,此時選擇使用分段卷積來計算卷積最適合。
多元函數卷積
按照翻轉、平移、積分的定義,還可以類似的定義多元函數上的積分:
性質
各種卷積算子都滿足下列性質:
其中Df 表示f的微分,如果在離散域中則是指差分算子,包括前向差分與後向差分兩種:
- 前向差分:
- 後向差分:
卷積定理
卷積定理指出,函數卷積的傅里葉變換是函數傅里葉變換的乘積。即,一個域中的卷積相當於另一個域中的乘積,例如時域中的卷積就對應於頻域中的乘積。
其中表示f 的傅里葉變換。
這一定理對拉普拉斯變換、雙邊拉普拉斯變換、Z變換、Mellin變換和Hartley變換(參見Mellin inversion theorem)等各種傅里葉變換的變體同樣成立。在調和分析中還可以推廣到在局部緊緻的阿貝爾羣上定義的傅里葉變換。
利用卷積定理可以簡化卷積的運算量。對於長度爲的序列,按照卷積的定義進行計算,需要做組對位乘法,其計算複雜度爲;而利用傅里葉變換將序列變換到頻域上後,只需要一組對位乘法,利用傅里葉變換的快速算法之後,總的計算複雜度爲。這一結果可以在快速乘法計算中得到應用。
在羣上的卷積
若 G 是有某 m 測度的羣(例如豪斯多夫空間上哈爾測度下局部緊緻的拓撲羣),對於G 上 m-勒貝格可積的實數或複數函數f 和g,可定義它們的卷積:
對於這些羣上定義的卷積同樣可以給出諸如卷積定理等性質,但是這需要對這些羣的表示理論以及調和分析的彼得-外爾定理。
應用
卷積在工程和數學上都有很多應用:
- 統計學中,加權的滑動平均是一種卷積。
- 概率論中,兩個統計獨立變量X與Y的和的概率密度函數是X與Y的概率密度函數的卷積。
- 聲學中,回聲可以用源聲與一個反映各種反射效應的函數的卷積表示。
- 電子工程與信號處理中,任一個線性系統的輸出都可以通過將輸入信號與系統函數(系統的衝激響應)做卷積獲得。
- 物理學中,任何一個線性系統(符合疊加原理)都存在卷積。
參見
外部鏈接
- PlanetMath上Convolution的資料。
- Visual convolution Java Applet
- Lecture notes, Jian-Jiun Ding (2013), Advanced Digital Signal Processing