卷积,toeplitz矩阵与傅里叶变换FT

news/2025/2/22 23:07:07/文章来源:https://www.cnblogs.com/eien/p/18719137

相关与卷积

一维相关和卷积运算是两个向量共同作用得到一个新的向量。
相关的分量形式:

\[\begin{aligned} 映射(f: \mathbb{Z} \to \mathbb{R}),对于(w = (w(-a), \ldots, w(0), \ldots, w(a))), \\ 一维相关: g(x) = \sum_{s = -a}^{a} w(s)f(x + s), x \in \mathbb{Z}. \end{aligned} \]

卷积的分量形式:

\[\begin{aligned} 映射(f: \mathbb{Z} \to \mathbb{R}),对于(w = (w(-a), \ldots, w(0), \ldots, w(a))), \\ 一维卷积(f \star w : \mathbb{Z} \to \mathbb{R}): g(x) := (f \star w)(x) = \sum_{s = -a}^{a} w(- s)f(x + s), x \in \mathbb{Z}. \end{aligned} \]

卷积是一种线性变换。
对于长度为M的卷积核\(w = h[0]...h[M - 1]\),卷积\(g(x) := (f \star w)(x)\) 等价于对向量\(x\)左乘了一个矩阵H。其中H是一个循环矩阵(也就是toeplitz矩阵).

\[H = \begin{bmatrix} h[0] & 0 & \cdots & 0 \\ h[1] & h[0] & \cdots & 0 \\ \vdots & \vdots & \ddots & \vdots \\ h[M - 1] & h[M - 2] & \cdots & h[0] \\ 0 & h[M - 1] & \cdots & h[1] \\ \vdots & \vdots & \ddots & \vdots \\ 0 & 0 & \cdots & h[M - 1] \end{bmatrix} \]


二维卷积示例:
这里卷积核为

\[\begin{bmatrix} 0 & 2 & 0\\ 1 & 2 & 1\\ 2 & 0 & 2\\ \end{bmatrix} \]

image


二维卷积和相关和一维情况类似,都是左乘一个卷积核形成的一个循环矩阵。一般在机器学习里面卷积和相关相同,因为一般用到的卷积核是对称的。

傅里叶变换(FT)

相关概念:数学分析中的傅里叶展开/数值分析里的最佳平方逼近。
即选择一簇线性无关的函数基(三角函数是正交的,自然线性无关)逼近要拟合的函数。
image

卷积,toeplitz矩阵与FT的关系

通过上述分析可以知道,卷积对应的线性变换的矩阵表示就是toeplitz矩阵,对一个向量做卷积就是左乘一个toeplitz矩阵。同时根据傅里叶变换中的性质,有在时间域的卷积操作等价于是频率域中的乘积。所以左乘一个toeplitz矩阵可以通过傅里叶变换来进行计算。
也就是$$g=f\star h=fft^{-1}(fft(f)fft(h))=Ah$$

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.hqwc.cn/news/884923.html

如若内容造成侵权/违法违规/事实不符,请联系编程知识网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

记录一次修复 JetBrains Rider 控制台输出乱码

在使用 JetBrains Rider 调试程序时,控制台输出日志出现了乱码。歪打正着结果困扰许久的问题得到了解决,于是记录下了这个小短文。 具体的修复建议如下:将终端编码设置为 GB2312 具体操作:设置->编辑器->常规->控制台->设置为 GB2312 ->保存之前在网上找各…

联合体union

占用大小按字节数最大的变量来,下面的联合体的大小为4字节typedef union MyUnion {char c; //1byteint i; //4byte } _MyUnion; int main() {_MyUnion v;printf("%x, %x\n", v.c, v.i); //未初始化时的值printf("%x, %x, %x, %x\n", &v, &v.c, &a…

使用docker搭建php开发环境

本文年代久远,可能已经不适用。仅供参考! 新的php镜像参考 多容器 apache+php+mariadb+redis主要配置 构建php镜像 基于官方php镜像进行定制,增加所需的php扩展,修改php配置,以及创建虚拟主机。 Dockerfile 文件php/Dockerfile FROM php:apache# php configure COPY php.i…

基于SVPWM矢量控制的无速度传感器电机控制系统simulink建模与仿真

1.课题概述基于SVPWM矢量控制的无速度传感器电机控制系统simulink建模与仿真,包括电机,SVPWM模块,矢量控制器模块等。2.系统仿真结果3.核心程序与模型 版本:MATLAB2022a 4.系统原理简介基于空间电压矢量控制(Space Vector PWM,SVPWM)的无速度传感器电机控制系统,特别是…

如何为 .NET 在线客服系统的 Open Api 开放接口实现 QPS 限流功能

在 .NET 中实现 QPS 限流有很多方案,对于小流量的应用,内存限流和基于中间件的限流方式可能已经足够,而对于高并发的分布式系统,可能需要 Redis 或者更复杂的算法,如令牌桶或滑动窗口。我在业余时间开发了一款自己的独立产品:升讯威在线客服与营销系统。陆陆续续开发了几…

C# 部署 tensorRT resnet分类模型 deploy TensorRT model with C#

基于生态的优势,一般都是用pytorch写模型;但是工业现场的部署,目前比较方便的还是onnx;当下为了追求效率,使用了TensorRT C#: .net8 GPU: A2000 / 4070Ti OS: windows Cuda: 11.8 预处理:Emgu.CVinstall v11.8 Cuda 默认安装,会生成以下路径 C:\Program Files\NVIDIA GP…

自然语言处理入门【第6章】:自然语言处理的应用

搜广推业务 信息检索 搜索、广告和推荐是很多互联网大厂赖以生存的命脉。比如,百度的搜索引擎、京东淘宝的推荐系统、腾讯的广告业务等。这一节我们就看看搜广推业务当中的自然语言处理方法。 1. 什么是信息检索 信息检索(Information Retrieval,IR)是指从大量的非结构化或…

提取东京食尸鬼主题曲

钢琴的声音实在是太好听了,这演奏者的指法真的逆天了! 利用IE缓存来提取视频文件是屡试不爽啊,只要视频不是经过加密的,均能提取成功。打开IE浏览器,关闭所有标签,仅留下一个空白标签,并打开Internet选项删除浏览历史记录点设置,查看文件,打开缓存文件夹,删除文件夹内…

树梅派XBMC解决方案比较

本文是Raspberry Pi XBMC Solutions Compared: Raspbmc vs OpenELEC vs XBian的译文,可能已经过时。XMBC是我们最喜欢的媒体中心解决方案之一。Raspberry Pi使得打造一台豪华的XMBC机器仅需35$。有三个为Pi定制的版本,但是,你该选择哪一个呢?这些版本略有差异,各有长短。本…

deepseek破限的一种方法

起因是https://www.south-plus.net/read.php?tid=2455666这篇帖子的作者在B22F提到他的Deepseek的破甲关键词和使用流程,咱直接复制粘贴直接被deepseek拒绝了呜呜呜 失败案例: 我大致提取了一段不会被办的话,试了一下可以和其他的混合使用破甲效果极佳(但是荤段子得靠其他…

实体机双系统

本文年代久远,可能已经不再适用。仅供参考首先说一下,双系统并不会影响两个系统各自的速度。两个系统都是独立的,分别位于硬盘的不同位置,独立运行。也就是说,当一个系统在运行时,另一个系统只是储存在硬盘上的文件,这自然不会对系统的速度有影响。 另外,安装双系统是有…

2024.1.17web HTTP协议基础

2.1 HTTP的基本概念与交互模型 HTTP是Hyper Text Transfer Protocol(超文本传输协议)的缩写。它是一个应用层协议,由请求和响应构成,是一个标准的客户端服务器模型。HTTP是一个无状态的协议。 HTTP 协议的特点 1.无连接—限制每次连接只处理一个请求 2.无状态—协议对于事务…