在
泛函分析中,卷积、旋积或褶积(英语:Convolution)是通过两个
函数f和g生成第三个函数的一种
数学运算,其本质是一种特殊的积分变换,表征函数f与g经过翻转和平移的重叠部分函数值乘积对重叠长度的积分。
简介
卷积(又名褶积)和反卷积(又名反褶积)是一种积分变换的数学方法,在许多方面得到了广泛应用。用卷积解决试井解释中的问题,早就取得了很好成果;而反卷积,直到最近,Schroeter、Hollaender和Gringarten等人解决了其计算方法上的稳定性问题,使反卷积方法很快引起了数学界的广泛注意。有专家认为,反卷积的应用是试井解释方法发展史上的又一次重大飞跃。他们预言,随着测试新工具和新技术的增加和应用,以及与其它专业研究成果的更紧密结合,试井在油气藏描述中的作用和重要性必将不断增大。
基本内涵
简单定义:设:f(x),g(x)是R1上的两个可积函数,作
积分:
可以证明,关于几乎所有的实数x,上述积分是存在的。这样,随着x的不同取值,这个积分就定义了一个新函数,称为函数与的卷积,记为。
容易验证,,并且仍为
可积函数。这就是说,把卷积代替乘法,L1(R1)空间是一个代数,甚至是巴拿赫代数。
卷积与傅里叶变换有着密切的关系。利用一点性质,即两函数的傅里叶变换的乘积等于它们卷积后的傅里叶变换,能使
傅里叶分析中许多问题的处理得到简化。
由卷积得到的函数一般要比和都光滑。特别当为具有
紧致集的
光滑函数,为局部可积时,它们的卷积也是光滑函数。利用这一性质,对于任意的
可积函数,都可以简单地构造出一列逼近于的
光滑函数列,这种方法称为函数的光滑化或
正则化。
卷积的概念还可以推广到
数列、
测度以及
广义函数上去。
定义
卷积是两个变量在某范围内相乘后求和的结果。如果卷积的变量是序列和,则卷积的结果
其中星号*表示卷积。当时序时,序列是的时序i取反的结果;时序取反使得以纵轴为中心翻转180度,所以这种相乘后求和的计算法称为卷积和,简称卷积。另外,是使位移的量,不同的对应不同的卷积结果。
如果卷积的变量是函数和,则卷积的计算变为
其中p是积分变量,积分也是求和,t是使函数h(-p)位移的量,星号*表示卷积。
性质
各种卷积算子都满足下列性质:
令,为任意常数或复常数,则卷积有如下性质:
其中,,,
单位阶跃序列与单位冲激序列而言,卷积还具有下列性质:
常见函数的卷积积分
在信号分析中,以下卷积积分的结果会比较常用:
其中,
在卷积运算中,如果能运用以上结论,将大大简化积分运算,节省解题时间。
卷积定理
卷积定理指出,函数卷积的
傅里叶变换是函数傅里叶变换的乘积。即,一个域中的卷积相当于另一个域中的乘积,例如
时域中的卷积就对应于
频域中的乘积。
这一定理对
拉普拉斯变换、
双边拉普拉斯变换、
Z变换、
Mellin变换和Hartley变换(参见Mellin inversion theorem)等各种傅里叶变换的变体同样成立。在
调和分析中还可以推广到在局部
紧致的
阿贝尔群上定义的傅里叶变换。
利用卷积定理可以简化卷积的运算量。对于长度为n的序列,按照卷积的定义进行计算,需要做(2n- 1)组对位
乘法,其计算复杂度为;而利用
傅里叶变换将序列变换到频域上后,只需要一组对位乘法,利用傅里叶变换的
快速算法之后,总的计算复杂度为。这一结果可以在快速
乘法计算中得到应用。
群上卷积
若G是有某m
测度的
群(例如
豪斯多夫空间上Harr测度下局部紧致的
拓扑群),对于G上m-勒贝格可积的
实数或
复数函数f和g,可定义它们的卷积:
对于这些群上定义的卷积同样可以给出诸如卷积定理等性质,但是这需要对这些群的表示理论以及调和分析的Peter-Weyl定理。
应用
卷积在工程和数学上都有很多应用:
统计学中,加权的滑动平均是一种卷积。
概率论中,两个统计独立变量X与Y的和的
概率密度函数是X与Y的概率密度函数的卷积。
光学中,反射光可以用光源与一个反映各种
反射效应的函数的卷积表示。
电子工程与信号处理中,任一个
线性系统的输出都可以通过将输入信号与系统函数(系统的
冲激响应)做卷积获得。
物理学中,任何一个线性系统(符合
叠加原理)都存在卷积。
介绍一个实际的概率学应用例子。假设需求到位时间的到达率为poisson(λ)分布,需求的大小的分布函数为D(.),则单位时间的需求量的分布函数为 F(x):
其中 D(k)(x)为k阶卷积。
卷积是一种线性运算,图像处理中常见的mask运算都是卷积,广泛应用于图像滤波。castlman的书对卷积讲得很详细。
高斯变换就是用高斯函数对图像进行卷积。高斯算子可以直接从离散
高斯函数得到:
再除以 sum 得到归一化算子
N是滤波器的大小,delta自选
首先,在提到卷积之前,必须提到卷积出现的背景。卷积是在
信号与线性系统的基础上或背景中出现的,脱离这个背景单独谈卷积是没有任何意义的,除了那个所谓卷反
公式上的数学意义和积分(或求和,离散情况下)。
信号与线性系统,讨论的就是信号经过一个线性系统以后发生的变化(就是输入 输出 和所经过的所谓系统,这三者之间的数学关系)。所谓线性系统的含义,就是,这个所谓的系统,带来的输出信号与输入信号的数学关系式之间是线性的运算关系。
因此,实际上,都是要根据我们需要待处理的信号形式,来设计所谓的系统传递函数,那么这个系统的传递函数和输入信号,在数学上的形式就是所谓的卷积关系。
卷积关系最重要的一种情况,就是在信号与线性系统或
数字信号处理中的卷积定理。利用该定理,可以将
时间域或空间域中的卷积运算等价为
频率域的相乘运算,从而利用FFT等快速算法,实现有效的计算,节省运算代价。
C++语言代码:
地震中的应用
地震勘探中,在地表激发点激发的地震子波(seismic wavelet)向地下传播,当遇到地下波阻抗界面时,一部分能量就会作为反射地震波向上反射回地表,被地面的传感器接收,随着地震波不断向下传播、反射、接收,就会记录一系列时间延迟的地震波(大地滤波后的地震子波),称为地震记录.
这一过程或地震记录可以用数学模型描述.如果假设地下介质为古皮奥(Goupilaud)的水平层状介质模型,子波为雷克(Ricker)子波,地震记录可以看作是由震源子波与地下反射率函数、多次反射、仪器等诸多因素的相卷积的过程,令x(t),w(t)和n(t)分别表示地震记录,地震子波及杂波,卷积过程数学模型描述为:
长期以来,卷积模型广泛用于描述地震信号.顾名思义,反卷积就是卷积的逆过程,从地震记录x(t)中恢复出反射率函数r(t).