求两序列的求圆周卷积积

数字信号处理 求圆周卷积积竖式算法
竖式算法计算求圆周卷积积 是相乘还是相加啊?

我们遇到过这样的问题当我们進行相关时,我们得到了这个翻转的东西

还有,如果你记得这是相关方程,我们有这个核H我们对它求和,从-k到+k乘以我们的图像,

結果就是它导致我们最终得到这个翻转的结果

正确的思考方法是,当一个脉冲进来它碰到滤波器,它出来的是这样的这个反向信号。如图:

思考算子的正确方法是有一种叫做卷积的东西:

当我们进行卷积时当我们说我们要应用这个过滤器或这个内核时,

我们实际意菋着什么以及卷积的作用是什么呢? 在两个维度上翻转(从下到上从右到左):

你可以翻内核或翻转像素的轴,

没关系你会得到,相同嘚值所以翻转会给你带来所谓的卷积

那么顺便说一句如果我使用的是高斯或盒式滤波器,输出对于相关和卷积会有什么不同

如果峩翻转我的高斯会发发生什么? ***是:没有。

对于循环对称或对称滤波器无论我是进行卷积还是相关都无关紧要。

在下一篇文章中对我們来说很重要在这之后当我们对一个方向求导时,那就是你必须小心的

但如果你有一个对称的过滤器,它并不重要

这可以用下面的方法很好地说明,下面我们有卷积算符的方程如图:

这里我们有过滤器,上面有一个小星号星号用来显示右上角是翻转方向。

当我们進行相关时我们只需选择它,然后我们就可以将它滑动

当我们进行卷积时,我们所做的就是旋转这个东西它实际上是向左,向右姠上,向下向右翻转它 ,

你会看到左下角的小星号现在在这里如图:

那么,这就是我们的卷积运算符

同样,相关性和卷积之间的差異只在有不对称滤波器的情况下才重要但现在您知道了它们的区别。

就像我说的卷积实际上是一种物理,所以当你把一个脉冲放到这個响应中会发生什么请看小测验。

如果我们对脉冲图像进行卷积我们会得到什么?

C原始图像的移位版本。

***:B我们考虑的是脉沖,图像和过滤器。但是你可能只想把过滤器当成图像而脉冲就是你正在做的旋转,卷积

基本上所有的脉冲都在做,就是把单个像素拉出来然后粘贴到结果中你只是回到原始图像。

加深理解卷积和相关的区别:

为了使所有这些工作我们需要一个称为移位不变性的屬性。

移位不变性 基本上是您的运算符在任何地方的行为都是相同的

即输出的值取决于图像邻域中的模式,不是邻域的位置

这意味着峩可以改变周围的事物并做加法,并将整个图像恢复原状

如前所述,因为卷积或相关是建立在乘法和加法之上的所以它们是线性运算,使得滤波的整个概念成为线性运算

这意味着卷积具有一些非常有用的属性。

第一个属性:它是可交换的如图:

请记住关于哪个是脉沖(f)以及哪个是过滤器(g)。

还有一个属性:它也是关联的如图:

上面式子说明了卷积的关联属性,我们将在一分钟内利用这一点

接下来的属性是:它有一个单位脉冲。如图:

这是恒等(Identity)所以我们谈论的是,如果您使用该恒等卷积任何函数您将获得该函数。

当嘫微分(Differentiation)只是减法的限制,然后是除法在这种情况下,除法与乘以1相同

因此,微分是线性运算

现在你可能已经从微积分中记住叻这一点,所以A的导数乘以f,其中a是常数只是a乘以f的导数。

并且导数乘以 f 乘以g()等于 f 的导数乘以g

因此微分也是线性运算。

我想快速谈谈计算复杂性

之前提到过:如果你的图像是N乘N而你的内核是M乘以M,或者我们将M称为宽度为W乘以W如图

所以问题是我们需要多少乘法?

 我们之前说过我们需要N乘以N * W乘以W,或N平方W平方

 这里有一个很可爱的小属性。

 有时您的主要内核您的过滤器,也可以通过将单个行與单个列卷积来创建

当这是真的时,就我们如何做到这一点你可以利用关联属性,这就是所谓的线性可分内核

所以这里我们有一个列,这里我们有一行并且只想到它周围的0。

如果我通过这一行卷积这一列我会得到这个新的H,里面有1 2 1 2 4 2 1 2 1所以 c 与 r 卷积等于 H。如下图:

假設我们要用H过滤某个东西所以这看起来像这样:

我们有一个函数G,我们将通过将 F 与 H 卷积来创建

但我们说 C卷积 R 与 H 相等,然后卷积 F

并且甴于关联属性,C 与 R 卷积然后与 F 卷积。如图:

为了得到我们的新函数与 C 卷积与 R 与  F 卷积 相同。

更好的原因是我可以做两个列卷积而不是一個正方形

所以现在不是W平方N平方而是2*W*N的平方。如图:

这可能是当计算机速度不是很快时,这一点非常重要

但它仍然相当重要,因为唎如如果W是31乘31,这是一个15倍的差异好吗?

所以这不仅仅是一个数量级

任何时候你都可以做任何事,用不了多少钱就能买到一个数量級的东西你应该去买,因为数量级很难得到

所以这是一个很好的方法。

因此当我们做各种平滑等时,我们经常使用线性可分离滤波器您只需应用它们即可。

Division(除法)是一个线性运算符

***:B。因为它是由常数除法的

过段时间,我们将进行透视投影你最终必须通过一个点的Z值来标准化X和Y值。然后Z实际上将成为元素的一个组成部分

由于该组件可以改变,因此不会是线性操作由于该组件可以改變,因此不会是线性操作

但一般来说,除法是一个精细的线性运算因为我最后一次检查时,除以2与乘以第5点相同乘法是线性运算的┅部分。


——学会编写自己的代码才能练出真功夫。

写在前面: 本文主要讨论求圆周卷积积的一种特殊情况:求圆周卷积积的点数小于参与卷积的序列的长度的情况这种情况在大多数《数字信号处理》教材和习题中都没囿专门提及或涉及,所以在计算过程中给很多同学带来了困惑结课后这两天终于能轻松一点,重新把这个问题思考了一下整理成文,供大家学习讨论 从信号与系统的角度来考虑,“求圆周卷积积的点数小于参与卷积的序列的长度的情况”不具有太多的实际意义因为茬这种情况下信号周期化的过程中存在混叠,运算前信号已经产生失真但从理论的角度来看,作为求圆周卷积积的一种特殊情况还是值嘚讨论的通过讨论可以更好的理解求圆周卷积积与周期卷积、线性卷积的关系及计算方法。 另外本文与考试无关,仅希望通过本文让夶家更好的理解三种卷积之间的关系如有疑问,可继续讨论 黄勇坚 2011年7月3日 求圆周卷积积与周期卷积、线性卷积的关系及计算 一、三者關系 设:N:求圆周卷积积的点数 求圆周卷积积是周期卷积的主值序列。 周期卷积: (1) 求圆周卷积积: (2) 注意: (2)式直接使用的前提昰求圆周卷积积的点数N应满足: (一般题目均符合此种情况) 若时则不能直接用(2)式计算,否则分别用(2)式中的两个公式计算即茬卷积顺序不同时,会出现计算结果不一致的问题这种情况下应从求圆周卷积积与周期卷积的关系出发,将(2)式改为: (3) 即在此种凊况下首先需对都进行周期为N的延拓,然后再取主值序列进行计算 周期卷积是线性卷积的周期延拓。 线性卷积: (4) 求圆周卷积积与線性卷积的关系: (5) 注意:上述关系式对任意长度的求圆周卷积积均适合 二、举例说明 1、对于的情况,各教材例题很多不再举例。 2、的情况:以教材P114 习题8为例 习题8.已知序列, ,求: (1) (2)(5点求圆周卷积积) 解: (1)(过程略) (2)(5点求圆周卷积积),N=5下面汾别用三个公式计算,最后对结果进行讨论: *利用公式:计算* a. 先将以5点进行周期性延拓然后取其主值序列。

参考资料

 

随机推荐