;
因果序列卷积的ListConvolve算法程序设计
[摘要] 采用调用Mathematica的内部函数ListConvolve进行因果序列卷积算法的程序设计，并编写因果序列卷积的ListConvolve算法和其它两种算法的程序，以及与其它两种算法的比较。并从实例展示其程序用法和高效性。
[关键字] Mathematica程序;卷积;ListConvolve算法;因果序列
Causal sequence convolution algorithm design ListConvolve
Abstract Invoke the Mathematica ListConvolve internal function is adopted to improve the program L causal sequence convolution algorithm design, and write a causal sequence of convolution ListConvolve algorithm program and other two kinds of algorithm, and compared with the other two algorithms. From examples show its program usage and efficiency.
Keywords Mathematica program;Convolution;ListConvolve algorithm;Causal sequence；
目录
引言 1
1因果序列卷积的意义 2
2常见有限长因果序列卷积的计算方法 2
2.1ListConvolve算法 2
2.2矩阵乘法算法 2
2.3矢量乘法算法 2
3各种卷积计算方法程序实现的编程思路 2
3.1各种算法编程前的准备 2
3.1.1构造滤波器 2
3.1.2计算滤波器时域系统函数的采样数据 3
3.1.3构造信号 4
3.2因果序列卷积的各种算法编程实现 6
3.2.1 ListConvolve算法程序设计 6
3.2.2矩阵乘法算法程序设计 7
3.2.3矢量乘法算法程序设计 8
4各种卷积算法的比较 9
4.1 ListConvolve算法本身的比较 9
4.2卷积的三种算法比较 9
5应用实例 9
6结语 11
致谢 11
参考文献 12
引言
Mathematica是由英国科学家斯蒂芬·沃尔夫勒姆领导的著名公司沃尔夫勒姆研究公司开发的一款广泛使用的科学计算软件。它拥有强大的数值计算和符号运算能力，是使用最广泛的数学软件之一。Mathematica的发布标志着现代科技计算的开始[1]。当然Mathematica是世界上通用计算系统中最强大的系统。自从1988发布以来，它已经对如何在科技和其它领域运用计算机产生了深刻的影响。就像人们常说，Mathematica的发布标志着现代科技计算的开始。自从上世纪六十年代以来，在数值、代数、图形、以及其它方面一直有个别的软件包存在。但是，Mathematica的基本概念是用一个连贯的和统一的方法创造一个能适用于科技计算各个方面的软件系统。至于如何实现这一点的关键之处是发明了一种新的计算机符号语言。这种语言可以仅仅用很少量的基本元素制造出广泛的物体，而满足科技计算的广泛性。这在人类历史上还是第一次出现。而且长期以来，Mathematica核心设计的普遍性使得它涉及的领域不断增长，扩大。从刚开始的是一个主要用于数学和科技计算的系统，到而今发展成为许多计算领域的主要力量，Mathematica已经成为世界上最强大的通用计算系统之一。Mathematica 应用广泛，Mathematica的用户群中最主要的是科技工作者和其它专业研究人士。不仅如此，Mathematica还被广泛地用于教学中。数学中的许多计算是非常繁琐的，特别是函数的作图费时又费力，而且所画的图形很不规范，所以现在流行用Mathematica符号计算系统进行学习，从高中到研究生院的数以百计的课程都是在使用它，很多问题便可以迎刃而解。此外，随着学生版的出现，Mathematica已经在全世界的学生中流行起来，成为了一个著名的学习，研究工具。在运算领域里，Mathematica作为当前国际控制界最流行的面向工程与科学计算的高级语言，近年来得到了业界的一致认可[2-7]。而且在控制各种系统的分析、设计算法方面得到了非常广泛的应用，其自身也因此得到了迅速的发展，功能也在不断扩充。便于对通信系统的物理概念和运行过程的直观理解，所以近年来在通信工程专业中得到了广大师生的特别重视和广泛的应用，在理论教学、课程实践环节，以及理论和技术前沿的研究中发挥了非常重要作用。而传统的因果序列卷积运算方法主要可以归纳为图解法, 解析法, 位移法和乘式法，矩阵乘法和矢量乘法等等。至于图解法直观明了但烦琐,不便于编程计算;而解析法适合于短序列编程, 当序列较长时, 表达式冗长, 不便于编程计算; 乘式法与解析法有相同特点, 不适合较长的两个序列的卷积和运算; 移位法有效的解决了较长序列的卷积和运算,是一种较好的计算方法, 特别对于长序列更显示出其计算能力, 但是需要考虑数据存储结构问题,而且也不便于理解。而本次毕业设计中主要提出了另一种计算方法, 即 ListConvolve算法，易于编程，算法高速[8]。
ListConvolve算法适用于因果序列卷积的运算,而且算法高速，占据内存小，因此是一种有效的计算方法. 在工程实际中具有较大的应用价值。
在此调用Mathematica软件中的软件ListConvolve算法进行程序设计，使其能够实现因果序
列卷积的ListConvolve算法的运算，使卷积运算更加简便。
许多文献讨论了卷积的计算方法以及对各种算法的总结，离散信号的卷积的多种求解方法，包括解析法、列表法、利用卷积和性质求解，应用单位样值信号求解，ListConvolve算法，矩阵算法及矢量空间算法。一般利用定义式、图解法和卷积性质这些通用方法进行求解连续信号的卷积或离散信号的卷积和。其中，利用定义式及卷积性质的求解，需要根据信号形式做繁琐的积分或求和；利用图解法进行的反转、平移、相乘、求和四个步骤，绘图过程麻烦[9]。而利用ListConvolve算法则简单高速，但算法不够透明，而矩阵算法则属于算法过程透明，但却占内存过大，矢量空间算法则运算时间过长。总之各种算法都有利有弊。本文将利用矩阵乘法和矢量乘法与ListConvolve算法相比较来突显出ListConvolve算法的优越性。
对于本文来说，主要研究了ListConvolve算法的如何调用，矩阵乘法算法和矢量乘法算法的计算方法是如何实现的。以及如何各种算法的程序编写[10]。
因为卷积在信号处理过程中是很常见的，所以研究如何实现卷积滤波是非常只有用的。相信卷积滤波会得到广泛应用。
1因果序列卷积的意义
任意两个序列的卷积并无意义的，只有一个代表离散信号的序列f(k)与描述离散LTI系统的序列g(k)（g(k)是离散LTI系统在单位脉冲序列k在它激励下的零状态响应。）卷积y(k)=f(k)*g(k) 才有真实的物理意义，它就是离散LTI系统在输入离散信号f(k)激励下的零状态响应，即系统的输出，能够代表真是的离散信号的序列，只能是因果序列，所以，本文只讨论一个描述输入信号的有限长序列与一个描述系统滤波特性的有限长序列之间的卷积，其结果是输入信号经过滤波器滤波后的结果。
因为在信号处理过程中会常遇到卷积计算问题，所以对于研究因果序列卷积的算法是非常有用的。
2常见有限长因果序列卷积的计算方法
2.1ListConvolve算法
ListConvolve算法，对于学者而言，只知道经过编程实践发现，它的算法很优越，省时省空间，目前来说它的时间复杂度和空间复杂度都是最低的。它是由世界著名软件公司（Mathematica软件）的精英们设计的，所以也只有他们清楚其内部算法。这种算法的透明度不高，至于为什么还要对其进行研究，是因为要学习如何熟练的使用它，如何熟练地调用其内部函数，以充分的利用其高速性与占据内存小的特点。
2.2矩阵乘法算法
矩阵乘法是一种高效的算法可以把一些一维递推优化到log（ n )，还可以求路径方案等，所以更是一种应用性极强的算法。矩阵，是线性代数中的基本概念之一。一个m×n的矩阵就是m×n个数排成m行n列的一个数阵。由于它把许多数据紧凑的集中到了一起，所以有时候可以简便地表示一些复杂的模型。例如：
