基于MATLAB完成卷积运算,基于Matlab的离散卷积

基于Matlab 的离散卷积

刘国良

(洛阳理工学院　河南洛阳　471000)

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

关键词:Matlab ;数字信号处理;离散卷积运算;快速傅里叶运算

中图分类号:TP311 文献标识码:B 文章编号:10042373X (2009)052125202

Discreted Convolution B ased on Matlab

L IU Guoliang

(L uoyang Institute of Science and Technology ,L uoyang ,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 convolu 2tion emulation operation can be adopted in Matlab.Linear convolution is the foundation of project application ,but circular con 2volution and FF T 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 ,FF T is a kind of most suitable method ,it is more rapidly in opera 2tional 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.

K eywords :Matlab ;DSP ;discreted convolution operation ;FF T

收稿日期:2008207202

0　引　言

在数字信号处理(DSP )系统中,无论在时域或频域都离不开卷积运算和快速傅里叶(FF T )运算。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 )的卷积。就是说频谱函数的乘积相当于时

间函数的卷积,反之亦然。即:

y (t )=x (t )3h (t )=

∫∞

-∞

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

(2

)

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

由于DSP 主要依靠计算机完成,而计算机无论在

时域或频域只能处理有限长的离散信号,因此上述关系

《现代电子技术》2009年第5期总第292期 测试?试量?自动化