基于Matlab的离散卷积

龙源期刊网 http://www.qikan.com.cn

基于Ｍａｔｌａｂ的离散卷积

作者：刘国良

来源：《现代电子技术》2009年第05期

摘 要：卷积运算广泛用于通讯、电子、自动化等领域的线性系统的仿真、分析及数字信号处理等方面。在Matlab中可以使用线性卷积、圆周卷积和快速傅里叶运算实现离散卷积。线性卷积是工程应用的基础，但圆周卷积和快速傅里叶运算实现线性离散卷积具有速度快等优势，圆周卷积采用循环移位，在Matlab中没有专用函数，需要根据圆周卷积的运算过程编制程序代码；快速傅里叶运算(FFT)是DSP的核心算法，在序列比较长时FFT是一种最合适的方法，运算速度快、程序简单，序列越长其优势越明显。以同一个例子介绍了进行离散卷积仿真运算的两种方法与特点。

关键词：Matlab；数字信号处理；离散卷积运算；快速傅里叶运算 中图分类号：TP311文献标识码：B 文章编号：1004-373X(2009)05-125-02 Discreted Convolution Based on Matlab LIU Guoliang

(Luoyang Institute of Science and Technology,Luoyang,471000,China)

Abstract：The convolution operation is widely applied in the field of linear system

simulations,analysis and DSP design of communication,electronics,automatic system and so on,linear convolution,circular convolution and FFT for discreted convolution emulation operation can be

adopted in Matlab.Linear convolution is the foundation of project application,but circular convolution and FFT in linear convolution have advantages of fast speed.Circular convolution adopts circulating shift,there is not privately function in Matlab,it needs to draw up program code according to the operational course of circular convolution.FFT is the major operation of DSP,when sequence grows fairly,FFT is a kind of most suitable method,it is more rapidly in operational speed and program simple,the longer of the sequence the more of its advantage.With a same example,to carry out characteristic and two kinds of method of emulation operation for discreted convolution are introduced.

Keywords：Matlab;DSP;discreted convolution operation;FFT 0 引 言

龙源期刊网 http://www.qikan.com.cn

在数字信号处理(DSP)系统中，无论在时域或频域都离不开卷积运算和快速傅里叶(FFT)运算。Matlab具有强大的矩阵运算能力，方便实用的绘图功能和语言的高度集成性，在DSP开发中，使用Matlab可以快速对系统进行仿真运算。

设输入信号为x(t)，其频谱函数为X(jΩ)，该信号通过滤波器h(t)后，其输出信号y(t)的频谱函数Y(jΩ)是频谱函数X(jΩ)与滤波器的频谱函数H(jΩ)的乘积，即: Y(jΩ)=X(jΩ)H(jΩ)(1)

该信号通过低通滤波器的示意图如图1所示。

而在时域，输出信号y(t)实际是输入信号x(t)与滤波器h(t)的卷积。就是说频谱函数的乘积相当于时间函数的卷积，反之亦然。即：

-∞x(τ)h(t-τ)dτ(2)

图1 信号通过低通滤波器的示意图

由于DSP主要依靠计算机完成，而计算机无论在时域或频域只能处理有限长的离散信号，因此上述关系可表示为： Y(ω)=X(ω)H(ω)(3)

y(n)=∑∞m=-∞x(m)h(n-m)=x(n)*h(n)(4) 式中：n表示离散时间；数字频率频率和周期。

1 实现离散卷积的方法

可以使用线性卷积、圆周卷积实现离散卷积，也可以根据时域信号函数与频谱函数的关系，使用快速傅里叶变换(FFT，IFFT)实现离散卷积。 1.1 圆周卷积

线性卷积是求离散系统响应的主要方法之一，许多重要应用都建立在这一理论基础上,如卷积滤波等，专用函数conv(x,h)可完成线性卷积过程。

圆周卷积的计算速度要远远快于线性卷积，如果选择圆周卷积的长度以用圆周卷积取代线性卷积。方法如下: (1) 定义圆周卷积的长度：选择

-1。

-1,则可

；f为输入信号频率；

为抽样

本文来源：https://www.bwwdw.com/article/41o.html