时域卷积定理及频域卷积定理-程序员宅基地

技术标签: 专栏03-图神经网络  频域卷积定理  时域卷积定理  

前言

在这里插入图片描述

之前的文章介绍了卷积积分的计算,作为图神经网络的前置知识点,理解卷积定理有助于理解谱图卷积

卷积定理是傅立叶变换满足的一个重要性质。维基百科中卷积定理的定义如下:

卷积定理指出,两个函数(或信号)的卷积的傅里叶变换是它们傅里叶变换的逐点乘积。

更一般地说,一个域(例如,时域)中的卷积等于另一个域(例如,频域)中的逐点乘法。具体分为时域卷积定理和频域卷积定理,时域卷积定理即时域内的卷积对应频域内的乘积;频域卷积定理即频域内的卷积对应时域内的乘积,两者具有对偶关系。

连续变量的函数

卷积定理

设两个函数 g ( t ) g(t) g(t) h ( t ) h(t) h(t) 的傅里叶变换为 G G G H H H

G ( f ) ≜ F { g } ( f ) = ∫ − ∞ ∞ g ( t ) e − i 2 π f t   d t , f ∈ R H ( f ) ≜ F { h } ( f ) = ∫ − ∞ ∞ h ( t ) e − i 2 π f t   d t , f ∈ R \begin{aligned} G(f) &\triangleq \mathcal{F}\{g\}(f) = \int_{-\infty}^{\infty}g(t) e^{-i 2 \pi f t} \, dt, \quad f \in \mathbb{R}\\ H(f) &\triangleq \mathcal{F}\{h\}(f) = \int_{-\infty}^{\infty}h(t) e^{-i 2 \pi f t} \, dt, \quad f \in \mathbb{R} \end{aligned} G(f)H(f)F{ g}(f)=g(t)ei2πftdt,fRF{ h}(f)=h(t)ei2πftdt,fR其中, ≜ \triangleq 表示定义为, F \mathcal{F} F是傅里叶变换算子。通常, F { g } ( f ) \mathcal{F}\{g\}(f) F{ g}(f) 也可以表示为 F { g ( t ) } \mathcal{F}\{g(t)\} F{ g(t)} F [ g ( t ) ] \mathcal{F}[g(t)] F[g(t)] g g g h h h 的卷积定义为:
r ( t ) = { g ∗ h } ( t ) ≜ ∫ − ∞ ∞ g ( τ ) h ( t − τ )   d τ = ∫ − ∞ ∞ g ( t − τ ) h ( τ )   d τ . r(t) = \{g*h\}(t) \triangleq \int_{-\infty}^{\infty} g(\tau) h(t-\tau)\, d\tau = \int_{-\infty}^{\infty} g(t-\tau) h(\tau)\, d\tau. r(t)={ gh}(t)g(τ)h(tτ)dτ=g(tτ)h(τ)dτ.这里, ∗ * 表示卷积。两个函数其他形式的卷积有 g ( t ) ∗ h ( t ) g(t)*h(t) g(t)h(t),博文卷积积分的计算有助于理解卷积。

卷积定理表述为:
R ( f ) ≜ F { r } ( f ) = G ( f ) H ( f ) . f ∈ R ( 1 ) R(f) \triangleq \mathcal{F}\{r\}(f) = G(f) H(f). \quad f \in \mathbb{R}\quad\quad(1) R(f)F{ r}(f)=G(f)H(f).fR(1)其中 R ( f ) R(f) R(f)也可以写成如下形式: R ( f ) ≜ F [ r ( t ) ] ≜ F [ h ( t ) ∗ g ( t ) ] R(f)\triangleq \mathcal{F}[r(t)] \triangleq \mathcal{F}[h(t)*g(t)] R(f)F[r(t)]F[h(t)g(t)]

证明

根据傅里叶变换及卷积积分有:
F [ g ( t ) ∗ h ( t ) ] = ∫ − ∞ + ∞ g ( t ) ∗ h ( t ) e − i ω t d t = ∫ − ∞ + ∞ ∫ − ∞ + ∞ g ( τ ) ∗ h ( t − τ ) d τ e − i ω t d t ( 卷积展开 ) = ∫ − ∞ + ∞ ∫ − ∞ + ∞ h ( t − τ ) e − i ω t d t   g ( τ ) d τ ( 1 ) \begin{aligned} \mathcal{F}[g(t)*h(t)]&=\int_{-\infin}^{+\infin} g(t)*h(t) e^{-i\omega t}dt \\ &=\int_{-\infin}^{+\infin}\int_{-\infin}^{+\infin} g(\tau)*h(t-\tau) d\tau e^{-i\omega t}dt \qquad (卷积展开)\\ &=\int_{-\infin}^{+\infin}\int_{-\infin}^{+\infin} h(t-\tau)e^{-i\omega t}dt\, g(\tau) d\tau \qquad (1)\\\\ \end{aligned} F[g(t)h(t)]=+g(t)h(t)etdt=++g(τ)h(tτ)dτetdt(卷积展开)=++h(tτ)etdtg(τ)dτ(1) x = t − τ x=t-\tau x=tτ,则 t = x + τ t=x+\tau t=x+τ d t = d x dt=dx dt=dx,所以:
( 1 ) = ∫ − ∞ + ∞ ∫ − ∞ + ∞ h ( x ) e − i ω ( x + τ ) d x   g ( τ ) d τ = ∫ − ∞ + ∞ ∫ − ∞ + ∞ h ( x ) e − i ω x e − i ω τ d x   g ( τ ) d τ = ∫ − ∞ + ∞ h ( x ) e − i ω x ∫ − ∞ + ∞ e − i ω τ d x   g ( τ ) d τ ( x 是在 τ 的积分中是常量 , 可以移出来 ) = g ( ω ) ⋅ h ( ω ) \begin{aligned} (1)&=\int_{-\infin}^{+\infin}\int_{-\infin}^{+\infin} h(x)e^{-i\omega (x+\tau)}dx \, g(\tau) d\tau \\ &=\int_{-\infin}^{+\infin}\int_{-\infin}^{+\infin} h(x) e^{-i\omega x} e^{-i\omega \tau} dx \, g(\tau) d\tau \\ &=\int_{-\infin}^{+\infin}h(x) e^{-i\omega x} \int_{-\infin}^{+\infin} e^{-i\omega \tau} dx\, g(\tau) d\tau \qquad (x是在\tau的积分中是常量,可以移出来)\\ &=g(\omega) \cdot h(\omega) \end{aligned} (1)=++h(x)e(x+τ)dxg(τ)dτ=++h(x)exeτdxg(τ)dτ=+h(x)ex+eτdxg(τ)dτ(x是在τ的积分中是常量,可以移出来)=g(ω)h(ω)其中 ω = 2 π f \omega=2\pi f ω=2πf
证毕。 ■ \blacksquare

小结

上述为时域卷积定理的公式及证明,频域卷积定义与之类似。所以连续变量函数的卷积定理为:

时域卷积定理为: F [ g ( t ) ∗ h ( t ) ] = G ( ω ) ⋅ H ( ω ) \mathcal{F}[g(t)*h(t)]=G(\omega) \cdot H(\omega) F[g(t)h(t)]=G(ω)H(ω)

频域卷积定理为: F [ g ( t ) ⋅ h ( t ) ] = 1 2 π   G ( ω ) ∗ H ( ω ) \mathcal{F}[g(t)\cdot h(t)]={\frac1 {2\pi}} \, G(\omega) * H(\omega) F[g(t)h(t)]=2π1G(ω)H(ω)

离散变量序列

卷积定理

对于离散变量序列,有着与连续变量函数类似的定理,这里 F \mathcal{F} F 表示离散傅里叶变换算子,设两个序列 g [ n ] g[n] g[n] h [ n ] h[n] h[n] 的傅里叶变换为 G G G H H H
G ( f ) ≜ F { g } ( f ) = ∑ n = − ∞ ∞ g [ n ] ⋅ e − i 2 π f n    , f ∈ R H ( f ) ≜ F { h } ( f ) = ∑ n = − ∞ ∞ h [ n ] ⋅ e − i 2 π f n    . f ∈ R \begin{aligned} G(f) &\triangleq \mathcal{F}\{g\}(f) = \sum_{n=-\infty}^{\infty} g[n]\cdot e^{-i 2\pi f n}\;, \quad f \in \mathbb{R}\\ H(f) &\triangleq \mathcal{F}\{h\}(f) = \sum_{n=-\infty}^{\infty} h[n]\cdot e^{-i 2\pi f n}\;. \quad f \in \mathbb{R} \end{aligned} G(f)H(f)F{ g}(f)=n=g[n]ei2πfn,fRF{ h}(f)=n=h[n]ei2πfn.fR

序列 g g g h h h 的离散卷积的定义为:

r [ n ] ≜ ( g ∗ h ) [ n ] = ∑ m = − ∞ ∞ g [ m ] ⋅ h [ n − m ] = ∑ m = − ∞ ∞ g [ n − m ] ⋅ h [ m ] . r[n] \triangleq (g * h)[n] = \sum_{m=-\infty}^\infty g[m]\cdot h[n - m] = \sum_{m=-\infty}^\infty g[n-m]\cdot h[m]. r[n](gh)[n]=m=g[m]h[nm]=m=g[nm]h[m]. ( g ∗ h ) [ n ] (g * h)[n] (gh)[n]也可以表示为 g [ n ] ∗ h [ n ] g[n]*h[n] g[n]h[n]

离散序列的卷积定理为:

R ( f ) = F { g ∗ h } ( f ) =   G ( f ) H ( f ) R(f) = \mathcal{F}\{g * h\}(f) =\ G(f) H(f) R(f)=F{ gh}(f)= G(f)H(f) \quad
R ( ω ) = F { g ∗ h } ( ω ) =   G ( ω ) H ( ω ) R(\omega) = \mathcal{F}\{g * h\}(\omega) =\ G(\omega) H(\omega) R(ω)=F{ gh}(ω)= G(ω)H(ω)

证明

本证明来自劉奕文教授的Lecture《The Discrete-Time Fourier Transform and
Convolution Theorems: A Brief Tutorial》(pdf)

Y ( ω ) = ∑ n = − ∞ ∞ y [ n ] e − j ω n = ∑ n = − ∞ ∞ ( ∑ m = − ∞ ∞ g [ m ] ⋅ h [ n − m ] ) e − j ω n = ∑ m = − ∞ ∞ g [ m ] ( ∑ m = − ∞ ∞ h [ n − m ] e − j ω n ) = ∑ m = − ∞ ∞ g [ m ] ( ∑ m = − ∞ ∞ h [ n − m ] e − j ω ( n − m )   e − j ω ( m ) ) = ∑ m = − ∞ ∞ g [ m ] e − j ω ( m ) ( ∑ m = − ∞ ∞ h [ n − m ] e − j ω ( n − m ) ) = ∑ m = − ∞ ∞ g [ m ] e − j ω ( m ) H ( ω ) = G ( ω ) H ( ω ) \begin{aligned} Y(\omega) &= {\sum_{n=-\infty}^\infty }y[n]e^{-j\omega n} \\ &={\sum_{n=-\infty}^\infty } \Bigg( {\sum_{m=-\infty}^\infty } g[m] \cdot h[n-m] \Bigg)e^{-j\omega n} \\ &={\sum_{m=-\infty}^\infty } g[m] \Bigg( {\sum_{m=-\infty}^\infty } h[n-m] e^{-j\omega n} \Bigg) \\ &={\sum_{m=-\infty}^\infty } g[m] \Bigg( {\sum_{m=-\infty}^\infty } h[n-m] e^{-j\omega (n-m)}\, e^{-j\omega (m)}\Bigg) \\ &={\sum_{m=-\infty}^\infty } g[m] e^{-j\omega (m)} \Bigg( {\sum_{m=-\infty}^\infty } h[n-m] e^{-j\omega (n-m)} \Bigg) \\ &={\sum_{m=-\infty}^\infty } g[m] e^{-j\omega (m)} H(\omega) \\ &=G(\omega) H(\omega) \\ \end{aligned} Y(ω)=n=y[n]ejωn=n=(m=g[m]h[nm])ejωn=m=g[m](m=h[nm]ejωn)=m=g[m](m=h[nm]e(nm)e(m))=m=g[m]e(m)(m=h[nm]e(nm))=m=g[m]e(m)H(ω)=G(ω)H(ω) ■ \blacksquare

小结

上述为时域卷积定理的公式及证明,频域卷积定义与之类似。所以离散变量序列的卷积定理为:

时域卷积定理为: F [ g ∗ h ] = G ( ω ) ⋅ H ( ω ) \mathcal{F}[g*h]=G(\omega) \cdot H(\omega) F[gh]=G(ω)H(ω)

频域卷积定理为: F [ g ⋅ h ] = 1 2 π   G ( ω ) ∗ H ( ω ) \mathcal{F}[g\cdot h]={\frac1 {2\pi}} \, G(\omega) * H(\omega) F[gh]=2π1G(ω)H(ω)

版权声明:本文为博主原创文章,遵循 CC 4.0 BY-SA 版权协议,转载请附上原文出处链接和本声明。
本文链接:https://blog.csdn.net/u012762410/article/details/127104625

智能推荐

hive使用适用场景_大数据入门:Hive应用场景-程序员宅基地

文章浏览阅读5.8k次。在大数据的发展当中,大数据技术生态的组件,也在不断地拓展开来,而其中的Hive组件,作为Hadoop的数据仓库工具,可以实现对Hadoop集群当中的大规模数据进行相应的数据处理。今天我们的大数据入门分享,就主要来讲讲,Hive应用场景。关于Hive,首先需要明确的一点就是,Hive并非数据库,Hive所提供的数据存储、查询和分析功能,本质上来说,并非传统数据库所提供的存储、查询、分析功能。Hive..._hive应用场景

zblog采集-织梦全自动采集插件-织梦免费采集插件_zblog 网页采集插件-程序员宅基地

文章浏览阅读496次。Zblog是由Zblog开发团队开发的一款小巧而强大的基于Asp和PHP平台的开源程序,但是插件市场上的Zblog采集插件,没有一款能打的,要么就是没有SEO文章内容处理,要么就是功能单一。很少有适合SEO站长的Zblog采集。人们都知道Zblog采集接口都是对Zblog采集不熟悉的人做的,很多人采取模拟登陆的方法进行发布文章,也有很多人直接操作数据库发布文章,然而这些都或多或少的产生各种问题,发布速度慢、文章内容未经严格过滤,导致安全性问题、不能发Tag、不能自动创建分类等。但是使用Zblog采._zblog 网页采集插件

Flink学习四:提交Flink运行job_flink定时运行job-程序员宅基地

文章浏览阅读2.4k次,点赞2次,收藏2次。restUI页面提交1.1 添加上传jar包1.2 提交任务job1.3 查看提交的任务2. 命令行提交./flink-1.9.3/bin/flink run -c com.qu.wc.StreamWordCount -p 2 FlinkTutorial-1.0-SNAPSHOT.jar3. 命令行查看正在运行的job./flink-1.9.3/bin/flink list4. 命令行查看所有job./flink-1.9.3/bin/flink list --all._flink定时运行job

STM32-LED闪烁项目总结_嵌入式stm32闪烁led实验总结-程序员宅基地

文章浏览阅读1k次,点赞2次,收藏6次。这个项目是基于STM32的LED闪烁项目,主要目的是让学习者熟悉STM32的基本操作和编程方法。在这个项目中,我们将使用STM32作为控制器,通过对GPIO口的控制实现LED灯的闪烁。这个STM32 LED闪烁的项目是一个非常简单的入门项目,但它可以帮助学习者熟悉STM32的编程方法和GPIO口的使用。在这个项目中,我们通过对GPIO口的控制实现了LED灯的闪烁。LED闪烁是STM32入门课程的基础操作之一,它旨在教学生如何使用STM32开发板控制LED灯的闪烁。_嵌入式stm32闪烁led实验总结

Debezium安装部署和将服务托管到systemctl-程序员宅基地

文章浏览阅读63次。本文介绍了安装和部署Debezium的详细步骤,并演示了如何将Debezium服务托管到systemctl以进行方便的管理。本文将详细介绍如何安装和部署Debezium,并将其服务托管到systemctl。解压缩后,将得到一个名为"debezium"的目录,其中包含Debezium的二进制文件和其他必要的资源。注意替换"ExecStart"中的"/path/to/debezium"为实际的Debezium目录路径。接下来,需要下载Debezium的压缩包,并将其解压到所需的目录。

Android 控制屏幕唤醒常亮或熄灭_android实现拿起手机亮屏-程序员宅基地

文章浏览阅读4.4k次。需求:在诗词曲文项目中,诗词整篇朗读的时候,文章没有读完会因为屏幕熄灭停止朗读。要求:在文章没有朗读完毕之前屏幕常亮,读完以后屏幕常亮关闭;1.权限配置:设置电源管理的权限。

随便推点

目标检测简介-程序员宅基地

文章浏览阅读2.3k次。目标检测简介、评估标准、经典算法_目标检测

记SQL server安装后无法连接127.0.0.1解决方法_sqlserver 127 0 01 无法连接-程序员宅基地

文章浏览阅读6.3k次,点赞4次,收藏9次。实训时需要安装SQL server2008 R所以我上网上找了一个.exe 的安装包链接:https://pan.baidu.com/s/1_FkhB8XJy3Js_rFADhdtmA提取码:ztki注:解压后1.04G安装时Microsoft需下载.NET,更新安装后会自动安装如下:点击第一个傻瓜式安装,唯一注意的是在修改路径的时候如下不可修改:到安装实例的时候就可以修改啦数据..._sqlserver 127 0 01 无法连接

js 获取对象的所有key值,用来遍历_js 遍历对象的key-程序员宅基地

文章浏览阅读7.4k次。1. Object.keys(item); 获取到了key之后就可以遍历的时候直接使用这个进行遍历所有的key跟valuevar infoItem={ name:'xiaowu', age:'18',}//的出来的keys就是[name,age]var keys=Object.keys(infoItem);2. 通常用于以下实力中 <div *ngFor="let item of keys"> <div>{{item}}.._js 遍历对象的key

粒子群算法(PSO)求解路径规划_粒子群算法路径规划-程序员宅基地

文章浏览阅读2.2w次,点赞51次,收藏310次。粒子群算法求解路径规划路径规划问题描述    给定环境信息,如果该环境内有障碍物,寻求起始点到目标点的最短路径, 并且路径不能与障碍物相交,如图 1.1.1 所示。1.2 粒子群算法求解1.2.1 求解思路    粒子群优化算法(PSO),粒子群中的每一个粒子都代表一个问题的可能解, 通过粒子个体的简单行为,群体内的信息交互实现问题求解的智能性。    在路径规划中,我们将每一条路径规划为一个粒子,每个粒子群群有 n 个粒 子,即有 n 条路径,同时,每个粒子又有 m 个染色体,即中间过渡点的_粒子群算法路径规划

量化评价:稳健的业绩评价指标_rar 海龟-程序员宅基地

文章浏览阅读353次。所谓稳健的评估指标,是指在评估的过程中数据的轻微变化并不会显著的影响一个统计指标。而不稳健的评估指标则相反,在对交易系统进行回测时,参数值的轻微变化会带来不稳健指标的大幅变化。对于不稳健的评估指标,任何对数据有影响的因素都会对测试结果产生过大的影响,这很容易导致数据过拟合。_rar 海龟

IAP在ARM Cortex-M3微控制器实现原理_value line devices connectivity line devices-程序员宅基地

文章浏览阅读607次,点赞2次,收藏7次。–基于STM32F103ZET6的UART通讯实现一、什么是IAP,为什么要IAPIAP即为In Application Programming(在应用中编程),一般情况下,以STM32F10x系列芯片为主控制器的设备在出厂时就已经使用J-Link仿真器将应用代码烧录了,如果在设备使用过程中需要进行应用代码的更换、升级等操作的话,则可能需要将设备返回原厂并拆解出来再使用J-Link重新烧录代码,这就增加了很多不必要的麻烦。站在用户的角度来说,就是能让用户自己来更换设备里边的代码程序而厂家这边只需要提供给_value line devices connectivity line devices