jjzjj

线性卷积、循环卷积、周期卷积的定义、计算方法及三者之间的关系

持续进步中 2023-07-10 原文

文章目录


前言

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


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

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

有关线性卷积、循环卷积、周期卷积的定义、计算方法及三者之间的关系的更多相关文章

  1. ruby - 如何使用 Nokogiri 的 xpath 和 at_xpath 方法 - 2

    我正在学习如何使用Nokogiri,根据这段代码我遇到了一些问题:require'rubygems'require'mechanize'post_agent=WWW::Mechanize.newpost_page=post_agent.get('http://www.vbulletin.org/forum/showthread.php?t=230708')puts"\nabsolutepathwithtbodygivesnil"putspost_page.parser.xpath('/html/body/div/div/div/div/div/table/tbody/tr/td/div

  2. ruby - 如何从 ruby​​ 中的字符串运行任意对象方法? - 2

    总的来说,我对ruby​​还比较陌生,我正在为我正在创建的对象编写一些rspec测试用例。许多测试用例都非常基础,我只是想确保正确填充和返回值。我想知道是否有办法使用循环结构来执行此操作。不必为我要测试的每个方法都设置一个assertEquals。例如:describeitem,"TestingtheItem"doit"willhaveanullvaluetostart"doitem=Item.new#HereIcoulddotheitem.name.shouldbe_nil#thenIcoulddoitem.category.shouldbe_nilendend但我想要一些方法来使用

  3. ruby - 为什么我可以在 Ruby 中使用 Object#send 访问私有(private)/ protected 方法? - 2

    类classAprivatedeffooputs:fooendpublicdefbarputs:barendprivatedefzimputs:zimendprotecteddefdibputs:dibendendA的实例a=A.new测试a.foorescueputs:faila.barrescueputs:faila.zimrescueputs:faila.dibrescueputs:faila.gazrescueputs:fail测试输出failbarfailfailfail.发送测试[:foo,:bar,:zim,:dib,:gaz].each{|m|a.send(m)resc

  4. ruby - Facter::Util::Uptime:Module 的未定义方法 get_uptime (NoMethodError) - 2

    我正在尝试设置一个puppet节点,但ruby​​gems似乎不正常。如果我通过它自己的二进制文件(/usr/lib/ruby/gems/1.8/gems/facter-1.5.8/bin/facter)在cli上运行facter,它工作正常,但如果我通过由ruby​​gems(/usr/bin/facter)安装的二进制文件,它抛出:/usr/lib/ruby/1.8/facter/uptime.rb:11:undefinedmethod`get_uptime'forFacter::Util::Uptime:Module(NoMethodError)from/usr/lib/ruby

  5. ruby - 树顶语法无限循环 - 2

    我脑子里浮现出一些关于一种新编程语言的想法,所以我想我会尝试实现它。一位friend建议我尝试使用Treetop(Rubygem)来创建一个解析器。Treetop的文档很少,我以前从未做过这种事情。我的解析器表现得好像有一个无限循环,但没有堆栈跟踪;事实证明很难追踪到。有人可以指出入门级解析/AST指南的方向吗?我真的需要一些列出规则、常见用法等的东西来使用像Treetop这样的工具。我的语法分析器在GitHub上,以防有人希望帮助我改进它。class{initialize=lambda(name){receiver.name=name}greet=lambda{IO.puts("He

  6. ruby-on-rails - 在 Ruby 中循环遍历多个数组 - 2

    我有多个ActiveRecord子类Item的实例数组,我需要根据最早的事件循环打印。在这种情况下,我需要打印付款和维护日期,如下所示:ItemAmaintenancerequiredin5daysItemBpaymentrequiredin6daysItemApaymentrequiredin7daysItemBmaintenancerequiredin8days我目前有两个查询,用于查找maintenance和payment项目(非排他性查询),并输出如下内容:paymentrequiredin...maintenancerequiredin...有什么方法可以改善上述(丑陋的)代

  7. Ruby 方法() 方法 - 2

    我想了解Ruby方法methods()是如何工作的。我尝试使用“ruby方法”在Google上搜索,但这不是我需要的。我也看过ruby​​-doc.org,但我没有找到这种方法。你能详细解释一下它是如何工作的或者给我一个链接吗?更新我用methods()方法做了实验,得到了这样的结果:'labrat'代码classFirstdeffirst_instance_mymethodenddefself.first_class_mymethodendendclassSecond使用类#returnsavailablemethodslistforclassandancestorsputsSeco

  8. ruby-on-rails - Rails 3.2.1 中 ActionMailer 中的未定义方法 'default_content_type=' - 2

    我在我的项目中添加了一个系统来重置用户密码并通过电子邮件将密码发送给他,以防他忘记密码。昨天它运行良好(当我实现它时)。当我今天尝试启动服务器时,出现以下错误。=>BootingWEBrick=>Rails3.2.1applicationstartingindevelopmentonhttp://0.0.0.0:3000=>Callwith-dtodetach=>Ctrl-CtoshutdownserverExiting/Users/vinayshenoy/.rvm/gems/ruby-1.9.3-p0/gems/actionmailer-3.2.1/lib/action_mailer

  9. ruby-on-rails - Rails 应用程序之间的通信 - 2

    我构建了两个需要相互通信和发送文件的Rails应用程序。例如,一个Rails应用程序会发送请求以查看其他应用程序数据库中的表。然后另一个应用程序将呈现该表的json并将其发回。我还希望一个应用程序将存储在其公共(public)目录中的文本文件发送到另一个应用程序的公共(public)目录。我从来没有做过这样的事情,所以我什至不知道从哪里开始。任何帮助,将不胜感激。谢谢! 最佳答案 无论Rails是什么,几乎所有Web应用程序都有您的要求,大多数现代Web应用程序都需要相互通信。但是有一个小小的理解需要你坚持下去,网站不应直接访问彼此

  10. ruby - Highline 询问方法不会使用同一行 - 2

    设置:狂欢ruby1.9.2高线(1.6.13)描述:我已经相当习惯在其他一些项目中使用highline,但已经有几个月没有使用它了。现在,在Ruby1.9.2上全新安装时,它似乎不允许在同一行回答提示。所以以前我会看到类似的东西:require"highline/import"ask"Whatisyourfavoritecolor?"并得到:Whatisyourfavoritecolor?|现在我看到类似的东西:Whatisyourfavoritecolor?|竖线(|)符号是我的终端光标。知道为什么会发生这种变化吗? 最佳答案

随机推荐