线性卷积、循环卷积、周期卷积的定义、计算方法及三者之间的关系-程序员宅基地

技术标签: matlab  线性卷积  算法  周期卷积  循环卷积  c语言  


前言

本文首先介绍了卷积的物理意义及卷积的性质(交换律、结合律、分配律),并简单阐述了线性时不变系统;之后分别就定义、适用范围及计算方法对线性卷积、循环卷积、周期卷积进行了介绍;最后总结了线性卷积、循环卷积、周期卷积之间的关系


一、卷积的物理意义及性质

1. 物理意义

卷积物理意义:一个函数(如:单位响应)在另一个函数(如:输入信号)上的加权叠加

对于线性时不变系统,如果知道该系统的单位响应,那么将单位响应和输入信号求卷积,就相当于把输入信号的各个时间点的单位响应加权叠加,就直接得到了输出信号。

线性时不变系统:既满足叠加原理又具有时不变特性
1)叠加定理 f ( a + b ) = f ( a ) + f ( b ) f(a+b) = f(a) + f(b) f(a+b)=f(a)+f(b)
2)时不变特性:特定输入的输出结果不会随时间变化。

2. 卷积性质

卷积积分满足交换律、结合律以及分配律:
1)交换律 f ( t ) ∗ g ( t ) = g ( t ) ∗ f ( t ) f(t)∗g(t)=g(t)∗f(t) f(t)g(t)=g(t)f(t);
2)结合律 f ( t ) ∗ g ( t ) ∗ h ( t ) = f ( t ) ∗ ( g ( t ) ∗ h ( t ) ) f(t)∗g(t)∗h(t)=f(t)∗(g(t)∗h(t)) f(t)g(t)h(t)=f(t)(g(t)h(t));
3)分配律 f ( t ) ∗ ( g ( t ) + h ( t ) ) = f ( t ) ∗ g ( t ) + f ( t ) ∗ h ( t ) f(t)∗(g(t)+h(t))=f(t)∗g(t)+f(t)∗h(t) f(t)(g(t)+h(t))=f(t)g(t)+f(t)h(t);


二、线性卷积定义及计算方法

1. 定义公式

y ( n ) = x ( n ) ∗ h ( n ) = ∑ x ( k ) h ( n − k ) y(n)=x(n)∗h(n)= ∑x(k)h(n−k) y(n)=x(n)h(n)=x(k)h(nk)

线性时不变离散系统中,若序列 x ( n ) x(n) x(n)是系统的输入信号 h ( n ) h(n) h(n)是系统在单位脉冲作用下的单位脉冲响应,由于输入离散时间序列 x ( n ) x(n) x(n)可表示为一系列脉冲的线性组合,根据线性系统的齐次性可加性 x ( n ) x(n) x(n)作用于系统所引起的零状态响应 y ( n ) y(n) y(n)就是序列 x ( n ) x(n) x(n) h ( n ) h(n) h(n)的卷积和。

2. 适用范围

线性卷积表示一个信号通过一个系统的输出,这个信号可以是无限长的,也可以是有限长的,可以的离散的也可以是连续的。

3. 计算方法

被卷积序列 x [ n ] x[n] x[n] h [ n ] h[n] h[n]的序列长度分别为 M M M N N N,则卷积得到的 y [ n ] y[n] y[n]序列长为 L = M + N − 1 L=M+N-1 L=M+N1,计算线性卷积的简单方法为进位保留法

在这里插入图片描述
C语言实现代码:

/*****************************************************************
*	线性卷积函数
*****************************************************************/

void convolution(double input1[], double input2[], double output[], int n, int m)
{
    
    int k = 0;
	int i = 0;
    int j = 0;

	for (k = 0; k < m + n - 1; k++) 
	{
    
		output[k] = 0;
	}

	//开始卷积
	//利用时延效果,记录所有乘积后,时间位置一样的相加
	for (i = 0; i < m; i++) 
	{
    
		for (j = 0; j < n; j++) 
		{
    
			output[i + j] += input1[i] * input2[j];
		}
	}
}

三、循环卷积定义及计算方法

1. 定义公式

设序列 h ( n ) h(n) h(n) x ( n ) x(n) x(n)的长度分别为 N N N M M M h ( n ) h(n) h(n) x ( n ) x(n) x(n) L L L点循环卷积定义为:
在这里插入图片描述
式中, L L L称为循环卷积区间长度, L ≥ m a x [ N , M ] L≥max[N,M] Lmax[N,M] x ( ( n − m ) ) L x\big((n-m)\big)_L x((nm))L是以 L L L为周期的周期信号, n n n m m m的变化区间均是 [ 0 , L − 1 ] [0,L-1] [0L1]

2. 适用范围

循环卷积(圆周卷积)针对离散信号,是一个有限长序列通过一个系统后的输出序列,在计算这个序列之前,必须先定义卷积运算的点数,不然这个运算就无法确定。

3. 计算方法

循环卷积相较于线性卷积区别在于卷积结果的序列长度与被卷积序列长度(输入信号)相同。由上面公式可知,直接计算比较复杂,通常采用线性卷积、矩阵相乘或FFT的方法计算循环卷积

1)利用线性卷积计算循环卷积
在这里插入图片描述
2)利用循环卷积矩阵计算循环卷积
在这里插入图片描述
L = m a x ( M , N ) L = max(M,N) L=max(M,N),则 x ( n ) x(n) x(n) L L L点“循环卷积矩阵”,其特点是:

① 第一行是序列 { x ( 0 ) , x ( 1 ) , x ( 2 ) , ⋯ , x ( L − 1 ) x(0),x(1),x(2),⋯,x(L−1) x(0),x(1),x(2),,x(L1)} 的循环倒相序列。
注:如果 x ( n ) x(n) x(n)的长度 M < L M<L M<L,则需要在 x ( n ) x(n) x(n)末尾补 L-M 个零后,再形成第一行的循环倒相序列;
② 第一行以后的各行均是前一行向右循环移1位形成的;
③ 矩阵的各主对角线上的序列值均相等。

3)利用FFT计算循环卷积
在这里插入图片描述
① 4点DFT

x1=[1,2,3,4];
x2=[1,2,3];

X1=fft(x1,4);
X2=fft(x2,4);
Y=X1.*X2;
y=ifft(Y);

在这里插入图片描述
② 6点DFT

x1=[1,2,3,4];
x2=[1,2,3];

X1=fft(x1,6);
X2=fft(x2,6);
Y=X1.*X2;
y=ifft(Y);

在这里插入图片描述


四、周期卷积定义及适用范围

1. 定义公式

在这里插入图片描述

2. 适用范围

周期卷积针对离散信号,是无限长周期离散信号通过一个离散系统后的输出。


总结

设卷积序列 x [ n ] x[n] x[n] h [ n ] h[n] h[n]的序列长度分别为 M M M N N N,则卷积得到的 y [ n ] y[n] y[n]序列长为 L L L

1. 线性卷积和循环卷积的关系

1)当L<N+M-1时,循环卷积是线性卷积长度为L的混叠;
2)当L=N+M-1时,循环卷积=线性卷积;
3)当L>N+M-1时,循环卷积是线性卷积末尾补L-(N+M-1)个零;

2. 三者之间的关系

1)周期卷积是线性卷积以一定序列长度(设为L)为周期的周期延拓。
2)对周期卷积取主值序列得到循环卷积。


参考

https://blog.csdn.net/abc123mma/article/details/120251384
https://blog.csdn.net/weixin_45810064/article/details/121322948?
https://blog.csdn.net/qq_39408570/article/details/110429928
https://blog.csdn.net/qq_45732223/article/details/110526379

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

智能推荐

微软发布WF教程及大量示例_微软wf-程序员宅基地

文章浏览阅读648次。继前面微软公司发布Windows Communication Foundation (WCF)和Windows CardSpace的示例程序之后,微软今天又发布了WF的教程和大量示例,对于学习WF的朋友来说,是一份很好的学习资料。下载包括教程和示例两个压缩包,在示例部分包括如下的Application:l OrderingStateMachinel Ou_微软wf

第一个错误的版本(Java算法每日一题)二分-程序员宅基地

文章浏览阅读1.1k次。问:第一个错误的版本答:/* The isBadVersion API is defined in the parent class VersionControl. boolean isBadVersion(int version); */public class Solution extends VersionControl { public int firstBadVersion(int n) { int left = 1; int righ_第一个错误的版本

opencv3.0 "ACCESS_MASK" ambiguous symbol_opencv ambiuous-程序员宅基地

文章浏览阅读1.2w次,点赞7次,收藏4次。问题如下,出现一大串“ACCESS_MASK”: 不明确的符号1>C:\Program Files (x86)\Microsoft SDKs\Windows\v7.0A\include\winspool.h(1187): error C2872: “ACCESS_MASK”: 不明确的符号1> 可能是“C:\Program Files (x86)\Microsoft SD_opencv ambiuous

bpmn流程图_【专栏丨顾问老包】泳道图OUT了?全息流程图来袭!-程序员宅基地

文章浏览阅读1.1k次,点赞2次,收藏3次。感谢包立南 老师的编写分享!微信公众号:顾问老包(ID:ourbpm)原标题:流程的语言:告别泳道图流程的表达方式就是流程图。它是一种语言,为我们描述那些业务活动是怎样发生的。流程的表达方式也经历了一个逐步演化和发展的过程,按照年代大致可以分成四种类型。业务逻辑图80年代,早期的流程图很简单,我们用来表达一个业务的基本逻辑,目的是很直观地呈现一系列动作之间的逻辑关系,我们可以叫它业务逻辑..._全息流程图

数据挖掘概念与技术pdf_数据挖掘概念与技术第三版部分习题答案-程序员宅基地

文章浏览阅读818次。反馈答案错、无答案请添加QQ:243631108关于我们大学课后答案公众号初心为:建立一个完整答案的公众号体系,致力做成全网最全的正版答案分享中心,节省大家搜索答案的宝贵时间。在未来,我们将会覆盖完大学九成的课后习题答案,还要推出一系列的复习资料,各高校的历年期末真题。这包括:哲学、经济学、法学、教育学、文学、历史学、理学、工学、农学、医学、管理学、艺术学等12个学科门类。还会包含大一基..._数据挖掘概念与技术第三版pdf

ncs v2.6.0安装教程 开发环境搭建_nrf connect sdk 2.6-程序员宅基地

文章浏览阅读720次,点赞20次,收藏27次。最简单的ncs v2.6.0开发环境搭建教程最实用的ncs v2.6.0开发环境搭建教程_nrf connect sdk 2.6

随便推点

Pytorch loss: SmoothL1Loss_-smooth 1 用pytorch实现-程序员宅基地

文章浏览阅读4.8k次。Huber loss也就是通常所说的SmoothL1 loss:SmoothL1对于异常点的敏感性不如MSE,而且,在某些情况下防止了梯度爆炸。在Pytorch中实现的SmoothL1损失是torch.nn.SmoothL1Loss, x和y可以是任何包含n个元素的Tensor,默认求均值。..._-smooth 1 用pytorch实现

ZYNQ UltraScale+ MPSoC Linux + ThreadX AMP玩法_zynqmpsoc 多核amp-程序员宅基地

文章浏览阅读2.7k次,点赞2次,收藏19次。ZYNQ UltraScale+ MPSoC Linux + ThreadX AMP玩法处理器架构ZYNQ 7000双核ARM-Cortex A9 CPUZYNQ UltraScale+ MPSoC四核ARM-Cortex A53 CPU、双核Cortex-R5 RPU、Mali-400 GPU(一个Geom_zynqmpsoc 多核amp

dumpcpp 生成COM对象源码列表_dumpcpp 导出 com组件 源文件-程序员宅基地

文章浏览阅读305次。办公软件类:版本:office_ProfessionalEnhancement2016https://download.csdn.net/download/lmhuanying1012/12132421版本:office_2010https://download.csdn.net/download/lmhuanying1012/12132514版本:WPS_2019..._dumpcpp 导出 com组件 源文件

正在写pytorch cpp前端的同行博客_py前端-程序员宅基地

文章浏览阅读175次。https://oldpan.me/别人说的深度学习圣经 https://www.jeremyjordan.me/ 找个时间一起填了_py前端

C#与Halcon联合(9)自制多功能halcon窗体_c#联合halcon做模板区域掩膜功能-程序员宅基地

文章浏览阅读5.3k次,点赞4次,收藏68次。目录1.窗体功能2.使用windows窗体控件库,制作该halcon多功能窗体3.效果展示--图像与操作说明4.如何使用自制的halcon多功能窗体5.视频展示链接与测试代码链接1.窗体功能①按比例显示图像,不填充拉伸②可鼠标拖动图像移动,滚轮缩放③可显示图像十字叉辅助显示④可以显示鼠标位于窗体位置的坐标及图像灰度值⑤可以保存当前窗体中显示的图像⑥可以对窗体进行截图⑦可以插入halcon窗体中的obj 与 msg ,并显示2.使用windows窗体控件库,制作该halcon多功能窗体①_c#联合halcon做模板区域掩膜功能

计算机二级选择题考word基础知识吗,计算机二级选择题难吗 ms office考试内容-程序员宅基地

文章浏览阅读150次。计算机二级是不少小伙伴们比较头疼的一门考试,也是一门比较重要的考试。有没有比较好的复习技巧呢?计算机二级选择题难吗?做题技巧是什么?计算机二级选择题难吗计算机二级选择题不难首先需要准备一个有计算机二级知识点的文档或者是书。文档的话一般百度一下就有,要是觉得网上的不太全面也可以通过某宝获取详细的资料。安装一个名为计算机二级宝典的APP。这个APP有点像考驾照必用的驾考宝典,关于里面二级的题很全面并且..._计算机二级选择题难吗

推荐文章

热门文章

相关标签