基于Mathematica卷积计算.docVIP

基于Mathematica卷积计算 　　摘 要:介绍了用Mathematica计算连续信号卷积和离散信号卷积和的方法。所提出的方法可得到连续信号卷积的精确解析解,而非用Matlab得到的近似解。对于离散信号,无论是有限长序列还是无限长序列的卷积和均可求解。与传统计算方法相比,该方法避开了繁琐的做图或积分求和过程,使计算更简便、快速。关键词:Mathematica; Matlab; 延续信号卷积; 离散信号卷积和 　　中图分类号:TN911??6-34文献标识码:A 　　文章编号:1004-373X(2010)19-0081-02 　　 　　Mathematica-based Computation of Convolution 　　LI Chun-ran 　　(Department of Physics, Bohai University, Jinzhou 121000, China) 　　Abstract: The method to compute the convolution of continuous signals and the convolution sum of discrete signals by ??using?? Mathematica is introduced. The precise analytical solutions rather than approximate solutions gained by Matlab are obtained by the method mentioned above. As for the discrete signals, the convolution sum of both finite and infinite sequences can be solved. In comparison with the traditional computation, this method avoids the complicated procedures such as drawing or integral summations, and makes the computation more convenient and swift.Keywords: Mathematica; Matlab; continuous signal convolution; convolution sum of discrete signal 　　 　　0 引 言 　　卷积在分析系统的零状态响应过程中有着重要作用,是线性时不变系统时域分析的基本数学工具,广泛应用于通信系统、信号处理、电子工程和微分方程等领域[1]。许多文献讨论了卷积的计算方法,文献[2]对连续信号的卷积计算方法做了总结;文献[3]讨论了离散信号的卷积和的4种求解方法,包括解析法、列表法、利用卷积和性质求解及应用单位样值信号求解。一般参考书[4-7]都给出利用定义式、图解法和卷积性质这些通用方法进行求解连续信号的卷积或离散信号的卷积和。其中,利用定义式及卷积性质的求解,需要根据信号形式做繁琐的积分或求和;利用图解法进行的反转、平移、相乘、求和四个步骤,绘图过程麻烦。总之这些传统的手工计算方法显得十分复杂而令人气馁。 　　随着计算机技术的发展,越来越多的计算问题交由计算机处理。Matlab和Mathematica作为优秀的科学计算软件,在工程计算、信号处理与通讯、图像处理等领域均得到广泛的应用。文献[8-9]分别讨论了利用Matlab求解连续和离散系统的卷积与卷积和。?┪南?[8]用插值近似解决问题,得到了很好的结果,但Matlab只得到了近似的数值解,没有得到精确解。?┪南?[9]的方法对有限长序列很有效,对无限长序列无能为力。Mathematica具有强大的符号及数值运算能力及方便实用的绘图功能,应用Mathematica有望得到精确的解析解。遗憾的是Mathematica并没有现成计算卷积的命令,考虑到用Mathematica来计算积分、求和、画图都非常容易,可以利用这些特点,做好计算卷积的程序包或函数,使计算过程大大简化,并得到精确的数值解或解析解。 　　1 连续时间信号的卷积积分 　　信号??f(t)与g(t)的卷积积分y(t)定义为[4]:?? 　　y(t)=f(t)*g(t)=∫∞??-∞??f(τ)g(t-τ)dτ(1) 　　式中:??f(t)和g(t)为连续信号;τ??为虚设积分变量,积分的结果为一新的时间信号。对于任意两个连续信号的卷积,只要知道其函数形式,即可以通过计算这两个函数乘积的积分来计算卷积。在Matlab中可以通过对两连续函数采样,然后调用库函