Yvon Shong

人生应当如盛夏🔆

  • 主页
  • 技术
  • 随笔
所有文章 关于我 友链

Yvon Shong

人生应当如盛夏🔆

  • 主页
  • 技术
  • 随笔

深入理解卷积

2016-04-13

卷积相关知识的梳理

内积

(a1,b1,c1)·(a2,b2,c2)=a1*a2+b1*b2+c1*c2

相关

有一个和卷积很容易混淆的概念,叫做相关

首先我们来比较一下:

自相关函数 R(x)=∫f(t)f(t-x)dt 移-积-和

卷积 g(x)=∫f(x)f(t-x)dt 褶-移-积-和

相关和卷积只差一个反转的过程

相关描述的是:把一个信号 平移一段后,与原函数的相似程度

具体怎么个【褶-移-积-和】后面再讲

相关的物理意义

相关最早是用来概率论中描述随机变量之间关系的概念,如相关系数。

实际上信号一般是一个随机过程,为了实现信号的检测、识别与提取,经常要了解两个信号的相似性,或一个信号经过一段延迟后自身的相似性。 但相关系数有缺陷。

因为分子是两个信号的内积,如sinx和cosx,从波形上看只是相位不同,而相关系数为零(因为正弦和余弦正交),因此引进相关函数,将原来两函数直接内积改为一个函数和另一个函数的延迟作内积。

其实二者从概念没有任何联系。由于相关函数第一个函数和第二个函数的延迟作内积(内积就是乘积累加),所以相关函数不满足交换,而卷积可以。延迟不同结果不同,所以相关系数是个数,而相关函数是函数,是延迟的函数。

卷积

在泛函分析中,卷积是通过两个函数f和g生成第三个函数的一种数学算子,表征函数f与 经过翻转和平移的g 的重叠部分的累积。 如果将参加卷积的一个函数看作区间的指示函数,卷积还可以被看作是“滑动加权平均”的推广。

一维信号

g(x)=∫f(x)f(t-x)dt

1
2
3
4
褶:f(x) -> f(-x)
移:f(-x) -> f(t-x)
积:f(x)f(t-x)
和:∫f(x)f(t-x)dt

反转平移之后重叠部分做内积(乘积累加)

某种意义上,对任意一个t,h(t)对它做了一个相同模式的贡献/变化

一些便于理解的例子

辐射:设某核电站事故中,某工作人员每天到抢险现场工作T分钟,接受一定剂量的辐射,辐射会自然地衰减,如此工作N天,如何计算总的辐射量?

信号:一个物体时刻会受到一个周期性的固定模式的的能量的冲击。并且,物体在受到某个时刻的能量后,会按照一个关于时间的函数进行变化。 如何计算 某个时刻 之前所有冲击信号的遗留和现在受到的冲击信号的能量之和?

复利:

本金 第一年 第二年 第三年 第四年 第五年
+100 100*(1.05)^1 100*(1.05)^2 100*(1.05)^3 100*(1.05)^4 100*(1.05)^5
+200 200*(1.05)^1 200*(1.05)^2 200*(1.05)^3 200*(1.05)^4
+300 300*(1.05)^1 300*(1.05)^2 300*(1.05)^3
+300 300*(1.05)^1 300*(1.05)^2
+200 200*(1.05)^1
+100

多项式乘法:

多项式a(n1)乘多项式b(n2)

对应的系数矩阵(向量)为a[最高次的系数,...,常数项系数] 和 b[最高次的系数,...,常数项系数]

多项式相乘结果的系数为c[最高次的系数,...,常数项系数]

1
c(i)=\sum_{k=0}^i a(k)b(i-k),i=0~(n1+n2)

二维图像

对应一维信号处理过程,把时间t的定义域空间换成图像的平面矩阵形状空间

二维离散卷积公式

二维连续卷积公式

(f*g)(t1,t2)=∫∫f(x1,x2)g(t1-x1,t2-x2)dx1 dx2

一些例子的理解

将图像按照模板进行处理

对于图像而言,离散卷积的计算过程是模板翻转,然后在原图像上滑动模板,把对应位置上的元素相乘后加起来,得到最终的结果。如果不考虑翻转,这个滑动-相乘-叠加的过程就是相关操作。

滤波

比如最简单的高斯模板,就是把模板内像素乘以不同的权值然后加起来作为模板的中心像素值,如果模板取值全为1,就是滑动平均;如果模板取值为高斯,就是加权滑动平均,权重是中间高,四周低,在频率上理解就是低通滤波器。

投影

因为当前模板内部图像和模板的相乘累加操作就是图像局部片段和模板的内积操作,这个过程可以看作图像的片段的空间 向 模板空间上的投影,一幅图像和一个模板卷积,得到的结果就是图像各个片段在这个空间上的投影累加和。如果这样的模板有一组,我们可以把这一组看成一组基,常见的如用一组Garbor滤波器提取图像的特征。

总的来说

在泛函分析中,卷积是通过两个函数f 和g 生成第三个函数的一种数学算子,表征函数f 与经过翻转和平移的g的重叠部分的累积。如果将参加卷积的一个函数看作区间的指示函数,卷积还可以被看作是“滑动加权平均”的推广。

(赋予不同的权值的想法)(加权累加)

二者在时间上不一定同步。卷积关心的正是二者在时间不同步的情况下(各个时间)激励和响应作用的结果。并且,这种响应是非因果的,即是不仅要考虑某一时刻输入信号的响应,还要全部包含这之前所有输入的响应。

离散卷积:

图像本身是离散的 所以权值也是离散的,模板是格子般的

忽略远距离的像素点的低权值影响,所以模板是有限的格子

权值分布的按照不同的函数分布,便是不同的卷积核,只不过受像素格的影响,是取离散的值。

整体就是个:褶 -> 移 -> 积 -> 和 的过程

卷积的类型

离散卷积和连续函数的卷积

线性卷积和循环卷积

通常说卷积时指的是线性卷积,在信号处理领域是跟线性时不变系统有关的一种运算

循环卷积我理解是使用DFT(FFT)计算线性卷积时的衍生品。

一些问题:

为什么要反转?

卷积的其中一方参与者是冲击响应,它所描述的的曲线方向与时间流逝一致。而卷积的输出等于以前的信号效果累加,这个累加必然从当前时间点逆时间流逝方向进行的。很显然,离当前时间越近,那个输入残留在系统中的回响就越大。

也可以从之前的多项式乘法的解释上理解。

卷积核可不可以是改变的?

某种意义上,对任意一个t,h(τ)对它做了一个相同模式的贡献/变化

h(t)能否是变化的?

g(x)=∫f(x)(t-x)dt -----> g(x)=∫f(x)f(x,t-x)dt

比如考虑边缘情况,边缘时模板就不相同了。

模板的大小是有限的

二维卷积时 什么情况才能忽略远距离的像素值呢?万一远距离的像素的权值更高呢?

(收敛的想法 h(t)得是收敛的)

线性时不变系统

连续时间系统与离散时间系统

线性系统有两种定义:
(1)根据系统的输入和输出关系是否具有线性来定义满足叠加原理的系统具有线性特性。即若对两个激励x1(n)和x2(n),有T[ax1(n)+bx2(n)]=aT[x1(n)]+bT[x2(n)],式中a、b为任意常数。不满足上述关系的为非线性系统。
(2)根据组成系统的元件特性来定义由线性元件和独立电源组成的系统。

时不变系统

系统的参数不随时间而变化,即不管输入信号作用的时间先后,输出信号响应的形状均相同,仅是从出现的时间不同。用数学表示为T[x(n)]=y[n]则 T[x(n-n0)]=y[n-n0],这说明序列x(n)先移位后进行变换与它先进行变换后再移位是等效的。 

线性时不变系统

既满足叠加原理又具有时不变特性,它可以用单位脉冲响应来表示。单位脉冲响应是输入端为单位脉冲序列时的系统输出,一般表示为h(n),即h(n)=T[δ(n)]。 
任一输入序列x(n)的相应y(n)=T[x(n)]=T[ δ(n-k)]; 
由于系统是线性的,所以上式可以写成y(n)=T[δ(n-k)];
又由于系统是时不变的,即有T[δ(n-k)]=h(n-k);
从而得y(n)=h(n-k)=x(n)*h(n);   
这个公式称为离散卷积,用“*”表示。

卷积定理

两个函数卷积后的傅里叶变换=两个函数的傅里叶变换后的乘积

卷积定理指出,函数卷积的傅里叶变换是函数傅里叶变换的乘积。即一个域中的卷积对应于另一个域中的乘积,例如时域中的卷积对应于频域中的乘积。 这一定理对拉普拉斯变换、双边拉普拉斯变换、Z变换、Mellin变换和Hartley变换等各种傅里叶变换的变体同样成立。在调和分析中还可以推广到在局部紧致的阿贝尔群上定义的傅里叶变换。 利用卷积定理可以简化卷积的运算量。对于长度为n的序列,按照卷积的定义进行计算,需要做2n-1组对位乘法,其计算复杂度为O(n^2);而利用傅里叶变换将序列变换到频域上后,只需要一组对位乘法,利用傅里叶变换的快速算法之后,总的计算复杂度为O(n*log n)。这一结果可以在快速乘法计算中得到应用。

赏

人民币版的赞

微信
  • matlab
  • tech

扫一扫,分享到微信

微信分享二维码
MatlabLens
深入理解快速傅里叶变换
  1. 1. 内积
  2. 2. 相关
    1. 2.1. 相关的物理意义
  3. 3. 卷积
    1. 3.1. 一维信号
    2. 3.2. 一些便于理解的例子
    3. 3.3. 二维图像
      1. 3.3.1. 二维离散卷积公式
      2. 3.3.2. 二维连续卷积公式
      3. 3.3.3. 一些例子的理解
        1. 3.3.3.1. 将图像按照模板进行处理
        2. 3.3.3.2. 滤波
        3. 3.3.3.3. 投影
    4. 3.4. 总的来说
      1. 3.4.1. 离散卷积:
      2. 3.4.2. 卷积的类型
        1. 3.4.2.1. 离散卷积和连续函数的卷积
        2. 3.4.2.2. 线性卷积和循环卷积
    5. 3.5. 一些问题:
    6. 3.6. 为什么要反转?
    7. 3.7. 卷积核可不可以是改变的?
    8. 3.8. 模板的大小是有限的
      1. 3.8.1. 线性时不变系统
  4. 4. 卷积定理
© 2025 Yvon Shong
Hexo Theme Yilia by Litten
  • 所有文章
  • 关于我
  • 友链

tag:

  • event
  • matlab
  • obsidian
  • onedrive
  • zotero
  • hexo
  • web
  • openwrt
  • IME
  • photograph
  • cdqz
  • slam
  • sensor
  • gps
  • gnss
  • lidar
  • radar
  • camera
  • photography

    缺失模块。
    1、请确保node版本大于6.2
    2、在博客根目录(注意不是yilia根目录)执行以下命令:
    npm i hexo-generator-json-content --save

    3、在根目录_config.yml里添加配置:

      jsonContent:
        meta: false
        pages: false
        posts:
          title: true
          date: true
          path: true
          text: false
          raw: false
          content: false
          slug: false
          updated: false
          comments: false
          link: false
          permalink: false
          excerpt: false
          categories: false
          tags: true
    

  • 输入法杂谈

    2023-06-15

    #IME

  • 我们把往事的记忆幻化成光影纸质想使之永垂不朽

    2023-06-07

    #cdqz

  • PDF 文献库解决方案——Zotero

    2023-05-29

    #onedrive#zotero

  • 从零开始 NanoPi R2S

    2023-05-16

    #openwrt

  • 聚焦激光雷达(五)——处理器

    2023-05-11

    #sensor#lidar

  • 观星指南

    2023-05-07

    #photograph

  • 聚焦激光雷达(四)——光学系统

    2023-04-19

    #sensor#lidar

  • 聚焦激光雷达(三)——接收器

    2023-04-01

    #sensor#lidar

  • 聚焦激光雷达(二)——激光器

    2023-03-28

    #sensor#lidar

  • 欧拉角、四元数、旋转矩阵、旋转向量的相互转换(C++、python)

    2023-03-23

    #slam

  • 聚焦激光雷达(一)——扫描器

    2023-03-12

    #sensor#lidar

  • Get started with Obsidian

    2022-01-28

    #obsidian#onedrive

  • 当时的月亮

    2022-01-01

  • 键盘摄影(十)——滤镜的种类和作用

    2021-09-21

    #sensor#camera#photography

  • 走进自动驾驶传感器(三)——卫星导航系统

    2021-09-21

    #sensor#gps#gnss

  • 与世界连接

    2021-02-21

  • 走进自动驾驶传感器(二)——毫米波雷达

    2021-01-24

    #sensor#radar

  • 计算摄影——“看得见”的未来

    2020-10-19

    #camera

  • 键盘摄影(九)——光学像差

    2020-10-09

    #sensor#camera#photography

  • 键盘摄影(八)——相机结构和卡口

    2020-10-08

    #sensor#camera#photography

  • 键盘摄影(七)——深入理解图像信号处理器 ISP

    2020-06-11

    #sensor#camera#photography

  • 键盘摄影(六)——相机成像模型

    2020-05-22

    #sensor#camera#photography

  • 键盘摄影(五)——相机成像元件:CMOS/CCD

    2020-05-12

    #sensor#camera#photography

  • 键盘摄影(四)——相机成像元件:胶片与彩色暗房

    2020-05-10

    #sensor#camera#photography

  • 走进自动驾驶传感器(一)——激光雷达

    2020-05-09

    #sensor#lidar

  • 键盘摄影(三)——等效焦距和等效光圈

    2020-05-09

    #sensor#camera#photography

  • 键盘摄影(二)——曝光三要素

    2020-05-07

    #sensor#camera#photography

  • 键盘摄影(一)——相机成像基本元件

    2020-05-06

    #sensor#camera#photography

  • 于是我便战胜了这寂寞的命运

    2020-01-26

  • 见过花火

    2019-02-08

  • 充实人生的常态是成长的阵痛

    2018-06-17

  • 当我们约会时我们会玩什么

    2017-10-22

    #event

  • 富士山下/爱情转移/春夏秋冬

    2017-03-12

  • Hexo 进阶与优化

    2017-02-20

    #hexo

  • Hexo 利用 AppVeyor 持续集成

    2017-02-20

    #hexo

  • Hello world on Hexo, my new personal blog.

    2017-02-18

    #hexo

  • 2016 年 MSP 的创新之旅

    2017-01-24

  • 在世界中心呼唤爱的野兽

    2017-01-02

  • 那么多的人,你要去哪里

    2016-09-19

  • 窗外日光弹指过·席间花影坐前移

    2016-08-25

  • Hospital Website based on VUE.js+PHP+MySQL

    2016-08-01

    #web

  • MatlabLens

    2016-06-11

    #matlab

  • 深入理解卷积

    2016-04-13

    #matlab

  • 深入理解快速傅里叶变换

    2016-04-09

    #matlab

  • 兰因絮果

    2016-02-25

  • 随轻轨呼啸而过的夏天 -- 记重庆社会实践之消失的棒棒

    2015-08-31

  • 一个人活

    2014-08-30

  • 我所领略的禅与佛光因缘

    2014-07-30

  • 写在前路未知时·爱在心口难开

    2014-05-24

  • 后高考时代

    2013-10-13

  • 青春荒唐我不负你

    2013-06-11

  • 我不怕千万人阻挡

    2012-08-24

  • 西西弗续

    2012-08-17

  • 玫瑰色

    2012-04-30

  • 17

    2012-04-11

  • 光年外自我拉扯

    2012-03-25

  • 像不变星空陪着我

    2011-12-31

  • 转头一切成苍茫

    2011-10-29

  • 却是因为爱太轻

    2011-08-29

  • 一样的夏了夏天

    2011-06-12

  • 知乎专栏 | 自律走行
  • 禾赛科技 | Hesai Tech
  • 图森未来 | TuSimple
  • 早稻田大学 | Waseda University
  • 东南大学 | Southeast University
  • 成都七中 | Chengdu No.7 Middle High School
  • 泡泡机器人 | PaoPaoRobot
  • 微软校园精英 | Microsoft Student Partner
自动驾驶从业者🤖 | Hesai Tech
早稻田修士🌸 | Waseda University
vlogger🎥 | B站id yvonshong
少年热忱🙋🏻‍♂️