《计算方法》研究报告.docxVIP

2016年《计算方法》研讨课题目------ Toeplitz矩阵与向量的乘积课题小组：组长： 李灵萱成员： 仇天乐、肖晟远、赵海宁、寿立夫、张连炜、李怡宁一个阶Toeplitz矩阵形如可以将嵌入一个“循环矩阵”，然后利用FFT来快速计算矩阵-向量积，其计算复杂度为。理论基础一个阶循环矩阵具有如下结构.一个阶离散Fourier矩阵为.这里的是虚数单位。数学上业已证明：一个阶循环矩阵能被对角化，即，且是的特征值。可以用IFFT计算，即利用和，则矩阵-向量积可以表为.对于阶Toeplitz矩阵，首先将嵌入一个阶循环矩阵，即，其中 将维列向量嵌入一个维列向量，这时有根据4)的事实，可用阶的FFT和IFFT计算矩阵-向量积，于是间接获得矩阵-向量积。这样，矩阵-向量积的计算复杂度从降低到了。公式的数学证明我们令矩阵，则需要证明为对角化的矩阵。根据矩阵的结合律，令记列向量，且，并且记，即以n为周期延拓后序列的主值序列，也就是循环矩阵C的第个列向量，得到。记，则由DFT的性质可以得到，所以而，综上，命题得证。编程实现存储方案设计由于Toeplitz矩阵与循环矩阵中有大量重复数据，因此，根据两者的特点，我们用长度为的向量来储存Toeplitz矩阵，即；用长度为的向量来储存循环矩阵，即。将储存空间从降低为。计算步骤设计循环矩阵1.1直接计算 因为我们将循环矩阵存储在长度为的向量中，所以可以利用次向