简介
基本内涵
可以证明,关于几乎所有的实数
x
,上述积分是存在的。这样,随着
x
的不同取值,这个积分就定义了一个新函数
,称为函数
与
的卷积,记为
。
由卷积得到的函数
一般要比
和
都光滑。特别当
为具有
紧致
集的
光滑函数
,
为局部可积时,它们的卷积
也是光滑函数。利用这一性质,对于任意的
可积函数
,都可以简单地构造出一列逼近于
的
光滑函数
列
,这种方法称为函数的光滑化或
正则化
。
定义
性质
性质名称
|
函数的卷积积分
|
序列的卷积和
|
|
|
|
延时特性
|
|
|
微分特性
(差分特性)
|
|
|
积分特性
(求和特性)
|
|
|
常见函数的卷积积分
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
卷积定理
这一定理对
拉普拉斯变换
、
双边拉普拉斯变换
、
Z变换
、
Mellin变换
和
Hartley变换
(参见
Mellin inversion theorem
)等各种傅里叶变换的变体同样成立。在
调和分析
中还可以推广到在局部
紧致
的
阿贝尔群
上定义的傅里叶变换。
利用卷积定理可以简化卷积的运算量。对于长度为
n
的序列,按照卷积的定义进行计算,需要做(2
n
- 1)组对位
乘法
,其
计算复杂度
为;而利用
傅里叶变换
将序列变换到频域上后,只需要一组对位乘法,利用傅里叶变换的
快速算法
之后,总的计算复杂度为。这一结果可以在快速
乘法
计算中得到应用。
群上卷积
对于这些群上定义的卷积同样可以给出诸如卷积定理等性质,但是这需要对这些群的
表示理论
以及调和分析的
Peter-Weyl定理
。
应用
统计学
中,加权的滑动平均是一种卷积。
概率论
中,两个统计独立变量X与Y的和的
概率密度函数
是X与Y的概率密度函数的卷积。
光学
中,反射光可以用光源与一个反映各种
反射效应
的函数的卷积表示。
电子工程
与信号处理中,任一个
线性系统
的输出都可以通过将输入信号与系统函数(系统的
冲激响应
)做卷积获得。
物理学
中,任何一个线性系统(符合
叠加原理
)都存在卷积。
卷积应用
(1张)
其中 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自选
首先,在提到卷积之前,必须提到卷积出现的背景。卷积是在
信号与线性系统
的基础上或背景中出现的,脱离这个背景单独谈卷积是没有任何意义的,除了那个所谓卷反
公式
上的数学意义和积分(或求和,离散情况下)。
信号与线性系统,讨论的就是信号经过一个线性系统以后发生的变化(就是输入 输出 和所经过的所谓系统,这三者之间的数学关系)。所谓线性系统的含义,就是,这个所谓的系统,带来的输出信号与输入信号的数学关系式之间是线性的运算关系。
因此,实际上,都是要根据我们需要待处理的信号形式,来设计所谓的系统传递函数,那么这个系统的传递函数和输入信号,在数学上的形式就是所谓的卷积关系。
卷积关系最重要的一种情况,就是在信号与线性系统或
数字信号处理
中的卷积定理。利用该定理,可以将
时间域
或空间域中的卷积运算等价为
频率域
的相乘运算,从而利用FFT等快速算法,实现有效的计算,节省运算代价。
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];