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

技术标签: 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

智能推荐

JWT(Json Web Token)实现无状态登录_无状态token登录-程序员宅基地

文章浏览阅读685次。1.1.什么是有状态?有状态服务,即服务端需要记录每次会话的客户端信息,从而识别客户端身份,根据用户身份进行请求的处理,典型的设计如tomcat中的session。例如登录:用户登录后,我们把登录者的信息保存在服务端session中,并且给用户一个cookie值,记录对应的session。然后下次请求,用户携带cookie值来,我们就能识别到对应session,从而找到用户的信息。缺点是什么?服务端保存大量数据,增加服务端压力 服务端保存用户状态,无法进行水平扩展 客户端请求依赖服务.._无状态token登录

SDUT OJ逆置正整数-程序员宅基地

文章浏览阅读293次。SDUT OnlineJudge#include<iostream>using namespace std;int main(){int a,b,c,d;cin>>a;b=a%10;c=a/10%10;d=a/100%10;int key[3];key[0]=b;key[1]=c;key[2]=d;for(int i = 0;i<3;i++){ if(key[i]!=0) { cout<<key[i.

年终奖盲区_年终奖盲区表-程序员宅基地

文章浏览阅读2.2k次。年终奖采用的平均每月的收入来评定缴税级数的,速算扣除数也按照月份计算出来,但是最终减去的也是一个月的速算扣除数。为什么这么做呢,这样的收的税更多啊,年终也是一个月的收入,凭什么减去12*速算扣除数了?这个霸道(不要脸)的说法,我们只能合理避免的这些跨级的区域了,那具体是那些区域呢?可以参考下面的表格:年终奖一列标红的一对便是盲区的上下线,发放年终奖的数额一定一定要避免这个区域,不然公司多花了钱..._年终奖盲区表

matlab 提取struct结构体中某个字段所有变量的值_matlab读取struct类型数据中的值-程序员宅基地

文章浏览阅读7.5k次,点赞5次,收藏19次。matlab结构体struct字段变量值提取_matlab读取struct类型数据中的值

Android fragment的用法_android reader fragment-程序员宅基地

文章浏览阅读4.8k次。1,什么情况下使用fragment通常用来作为一个activity的用户界面的一部分例如, 一个新闻应用可以在屏幕左侧使用一个fragment来展示一个文章的列表,然后在屏幕右侧使用另一个fragment来展示一篇文章 – 2个fragment并排显示在相同的一个activity中,并且每一个fragment拥有它自己的一套生命周期回调方法,并且处理它们自己的用户输_android reader fragment

FFT of waveIn audio signals-程序员宅基地

文章浏览阅读2.8k次。FFT of waveIn audio signalsBy Aqiruse An article on using the Fast Fourier Transform on audio signals. IntroductionThe Fast Fourier Transform (FFT) allows users to view the spectrum content of _fft of wavein audio signals

随便推点

Awesome Mac:收集的非常全面好用的Mac应用程序、软件以及工具_awesomemac-程序员宅基地

文章浏览阅读5.9k次。https://jaywcjlove.github.io/awesome-mac/ 这个仓库主要是收集非常好用的Mac应用程序、软件以及工具,主要面向开发者和设计师。有这个想法是因为我最近发了一篇较为火爆的涨粉儿微信公众号文章《工具武装的前端开发工程师》,于是建了这么一个仓库,持续更新作为补充,搜集更多好用的软件工具。请Star、Pull Request或者使劲搓它 issu_awesomemac

java前端技术---jquery基础详解_简介java中jquery技术-程序员宅基地

文章浏览阅读616次。一.jquery简介 jQuery是一个快速的,简洁的javaScript库,使用户能更方便地处理HTML documents、events、实现动画效果,并且方便地为网站提供AJAX交互 jQuery 的功能概括1、html 的元素选取2、html的元素操作3、html dom遍历和修改4、js特效和动画效果5、css操作6、html事件操作7、ajax_简介java中jquery技术

Ant Design Table换滚动条的样式_ant design ::-webkit-scrollbar-corner-程序员宅基地

文章浏览阅读1.6w次,点赞5次,收藏19次。我修改的是表格的固定列滚动而产生的滚动条引用Table的组件的css文件中加入下面的样式:.ant-table-body{ &amp;amp;::-webkit-scrollbar { height: 5px; } &amp;amp;::-webkit-scrollbar-thumb { border-radius: 5px; -webkit-box..._ant design ::-webkit-scrollbar-corner

javaWeb毕设分享 健身俱乐部会员管理系统【源码+论文】-程序员宅基地

文章浏览阅读269次。基于JSP的健身俱乐部会员管理系统项目分享:见文末!

论文开题报告怎么写?_开题报告研究难点-程序员宅基地

文章浏览阅读1.8k次,点赞2次,收藏15次。同学们,是不是又到了一年一度写开题报告的时候呀?是不是还在为不知道论文的开题报告怎么写而苦恼?Take it easy!我带着倾尽我所有开题报告写作经验总结出来的最强保姆级开题报告解说来啦,一定让你脱胎换骨,顺利拿下开题报告这个高塔,你确定还不赶快点赞收藏学起来吗?_开题报告研究难点

原生JS 与 VUE获取父级、子级、兄弟节点的方法 及一些DOM对象的获取_获取子节点的路径 vue-程序员宅基地

文章浏览阅读6k次,点赞4次,收藏17次。原生先获取对象var a = document.getElementById("dom");vue先添加ref <div class="" ref="divBox">获取对象let a = this.$refs.divBox获取父、子、兄弟节点方法var b = a.childNodes; 获取a的全部子节点 var c = a.parentNode; 获取a的父节点var d = a.nextSbiling; 获取a的下一个兄弟节点 var e = a.previ_获取子节点的路径 vue

推荐文章

热门文章

相关标签