基本介紹
簡介,基本內涵,定義,性質,常見函式的卷積積分,卷積定理,群上卷積,套用,地震中的套用,
簡介
褶積（又名卷積）和反褶積（又名去卷積）是一種積分變換的數學方法，在許多方面得到了廣泛套用。用褶積解決試井解釋中的問題，早就取得了很好成果；而反褶積，直到最近，Schroeter、Hollaender和Gringarten等人解決了其計算方法上的穩定性問題，使反褶積方法很快引起了試井界的廣泛注意。有專家認為，反褶積的套用是試井解釋方法發展史上的又一次重大飛躍。他們預言，隨著測試新工具和新技術的增加和套用，以及與其它專業研究成果的更緊密結合，試井在油氣藏描述中的作用和重要性必將不斷增大。
基本內涵
簡單定義：設:f(x),g(x)是R1上的兩個可積函式，作積分：
可以證明，關於幾乎所有的實數x，上述積分是存在的。這樣，隨著x的不同取值，這個積分就定義了一個新函式h(x)，稱為函式f與g的卷積，記為h(x)=(f*g)(x)。
容易驗證，(f * g)(x) = (g * f)(x)，並且(f * g)(x)仍為可積函式。這就是說，把卷積代替乘法，L1（R1）空間是一個代數，甚至是巴拿赫代數。
卷積與傅立葉變換有著密切的關係。利用一點性質，即兩函式的傅立葉變換的乘積等於它們卷積後的傅立葉變換，能使傅立葉分析中許多問題的處理得到簡化。
由卷積得到的函式f*g一般要比f和g都光滑。特別當g為具有緊緻集的光滑函式，f為局部可積時，它們的卷積f * g也是光滑函式。利用這一性質，對於任意的可積函式f，都可以簡單地構造出一列逼近於f的光滑函式列fs，這種方法稱為函式的光滑化或正則化。
定義
卷積是兩個變數在某範圍內相乘後求和的結果。如果卷積的變數是序列x(n)和h(n)，則卷積的結果
，
序列的卷積和
序列的卷積和其中星號*表示卷積。當時序n=0時，序列h(-i)是h(i)的時序i取反的結果；時序取反使得h(i)以縱軸為中心翻轉180度，所以這種相乘後求和的計算法稱為卷積和，簡稱卷積。另外，n是使h(-i)位移的量，不同的n對應不同的卷積結果。
如果卷積的變數是函式x(t)和h(t)，則卷積的計算變為
，
函式的卷積積分
函式的卷積積分其中p是積分變數，積分也是求和，t是使函式h(-p)位移的量，星號*表示卷積。
性質
各種卷積運算元都滿足下列性質：
令x(t)*h(t)=y(t)，a為任意常數或復常數，則卷積有如下性質：
性質名稱 | 函式的卷積積分 | 序列的卷積和 |
交換律 | x(t)*h(t) = h(t)*x(t) | x[k]*h[k] = h[k]*x[k] |
分配律 | x(t)*[g(t)+h(t)] = x(t)*g(t)+x(t)*h(t) | x[k]*(g[k]+h[k]) = x[k]*g[k]+x[k]*h[k] |
結合律 | [x(t)*g(t)]*h(t) = x(t)*[g(t)*h(t)] | {x[k]*g[k]}*h[k] = x[k]*{g[k]*h[k]} |
數乘結合律 | a[x(t)*h(t)]=[ax(t)]*h(t)=x(t)*ah(t) |  |
平移特性 |  | x[k-n]*h[k-l] = y[k-(n+l)] |
微分特性 (差分特性) | y'(t)=x'(t)*h(t) = x(t)*h'(t) |   |
積分特性 (求和特性) |  |  |
等效特性 |  |  |
其中，
，
，
性質名稱 | 函式的卷積積分 | 序列的卷積和 |
 |  | |
延時特性 |  |  |
微分特性 (差分特性) |  |   |
積分特性 (求和特性) |  |  |
常見函式的卷積積分
在信號分析中，以下卷積積分的結果會比較常用：
 |  |  |
 | |  |
 | |  |
|  |  |
| |  |
| |  |
其中
, 
在卷積運算中，如果能運用以上結論，將大大簡化積分運算，節省解題時間。
卷積定理
perfect spaces卷積混響F[g(x)*f(x)] = F[g(x)]F[f(x)]
其中F表示的是傅立葉變換。
這一定理對拉普拉斯變換、雙邊拉普拉斯變換、Z變換、Mellin變換和Hartley變換（參見Mellin inversion theorem）等各種傅立葉變換的變體同樣成立。在調和分析中還可以推廣到在局部緊緻的阿貝爾群上定義的傅立葉變換。
利用卷積定理可以簡化卷積的運算量。對於長度為n的序列，按照卷積的定義進行計算，需要做2n- 1組對位乘法，其計算複雜度為；而利用傅立葉變換將序列變換到頻域上後，只需要一組對位乘法，利用傅立葉變換的快速算法之後，總的計算複雜度為。這一結果可以在快速乘法計算中得到套用。
群上卷積
卷積與相關分析對於這些群上定義的卷積同樣可以給出諸如卷積定理等性質，但是這需要對這些群的表示理論以及調和分析的Peter-Weyl定理。
套用
卷積在工程和數學上都有很多套用：
統計學中，加權的滑動平均是一種卷積。機率論中，兩個統計獨立變數X與Y的和的機率密度函式是X與Y的機率密度函式的卷積。聲學中，回聲可以用源聲與一個反映各種反射效應的函式的卷積表示。電子工程與信號處理中，任一個線性系統的輸出都可以通過將輸入信號與系統函式（系統的衝激回響）做卷積獲得。物理學中，任何一個線性系統（符合疊加原理）都存在卷積。
介紹一個實際的機率學套用例子。假設需求到位時間的到達率為poisson(λ)分布，需求的大小的分布函式為D（.)，則單位時間的需求量的分布函式為 F(x)：
卷積套用
卷積套用其中 D(k)(x)為k階卷積。
卷積是一種線性運算,圖像處理中常見的mask運算都是卷積，廣泛套用於圖像濾波。castlman的書對卷積講得很詳細。
高斯變換就是用高斯函式對圖像進行卷積。高斯運算元可以直接從離散高斯函式得到：
for(i=0; i<N; i++)
{
for(j=0; j<N; j++)
{
g[i*N+j]=exp(-((i-(N-1)/2)^2+(j-(N-1)/2)^2))/(2*delta^2));
sum += g[i*N+j];
}
}
再除以 sum 得到歸一化運算元
N是濾波器的大小，delta自選
信號與線性系統，討論的就是信號經過一個線性系統以後發生的變化（就是輸入 輸出 和所經過的所謂系統，這三者之間的數學關係）。所謂線性系統的含義，就是，這個所謂的系統，帶來的輸出信號與輸入信號的數學關係式之間是線性的運算關係。
因此，實際上，都是要根據我們需要待處理的信號形式，來設計所謂的系統傳遞函式，那么這個系統的傳遞函式和輸入信號，在數學上的形式就是所謂的卷積關係。
C++語言代碼：
void convolution(double *input1, double *input2, double *output, int mm, int nn)
{
double *xx = new double[mm + nn - 1];
// do convolution
for (int i = 0; i < mm + nn - 1; i++)
{
xx[i] = 0.0;
for (int j = 0; j < mm; j++)
{
if (i - j >= 0 && i - j < nn)
xx[i] += input1[j] * input2[i - j];
}
}
// set value to the output array
for (int i = 0; i < mm + nn - 1; i++)
output[i] = xx[i];
delete[] xx;
}
地震中的套用
地震勘探中,在地表激發點激發的地震子波(seismic wavelet)向地下傳播,當遇到地下波阻抗界面時,一部分能量就會作為反射地震波向上反射回地表,被地面的感測器接收,隨著地震波不斷向下傳播、反射、接收,就會記錄一系列時間延遲的地震波(大地濾波後的地震子波),稱為地震記錄.
這一過程或地震記錄可以用數學模型描述.如果假設地下介質為古皮奧(Goupilaud)的水平層狀介質模型,子波為雷克(Ricker)子波,地震記錄可以看作是由震源子波與地下反射率函式、多次反射、儀器等諸多因素的相褶
積的過程,令x(t),w(t)和n(t)分別表示地震記錄,地震子波及噪聲,褶積過程數學模型描述為：
公式
公式長期以來,褶積模型廣泛用於描述地震信號.顧名思義,反褶積就是褶積的逆過程,從地震記錄x(t)中恢復出反射率函式r(t).
