几个结构矩阵乘积的Strassen算法.pdfVIP

第22卷 第 4期 五邑大学学报 (自然科学版 ) Vb1．22 No．4 2008年 11月 JOURNAL OF WUYIUNIVERSITY (Natural NOV． 2008 文章编号：1006—7302 (2008)04·0028-04 几个结构矩阵乘积的Strassen算法 曹寒冬，曹文胜 (五 邑大学 数理 系，广东 江门 529020) 摘要：运用广义中心对称矩阵和广义中心Hermitian矩阵的约化性质得到 了计算此类矩阵乘积 的Strassen算法．此算法和传统算法相比，大约是传统算法计算量的一半． 关键词 ：广义中心对称矩阵；广义中心Hermitian矩阵；矩阵乘积；Strassen算法 中图分类号：O241．6 文献标识码：A SeveralStrassenAlgorithmsofStructuredM atrixMultiplication CA0 Han—dong， CAO Wen—sheng (DepartmentofMathematicsandPhysics，WuyiUniversity，Jiangmen529020，China) Abstract：Inthispaper，wetakeadvantageofthereducibilityofthegeneralizedcentr0symmetric matrixandthegeneralizedcentrohermitianmatrixtoproposetwoStrassenalgorithmsforcomputing theproductsofsuchmatrices．It’Sshownthattheamountofcomputationofouralgorithm ishalfof thatoftheconventionalalgorithm． Keywords：generalizedcentr0symmetricmatrix；generalizedcentrohermitian matrix；matrix-matrix multiplicati0n；Strassenalgorithm 1 引言及准备知识 在计算信号复原问题时，离散方程所表现的大多是一个大规模的线性方程组．对于许多已在信 号复原问题 中出现的结构矩阵，如应用小波技术处理信号复原问题时，就会出现诸如 (广义 )中心 对称矩阵、(广义 )中心 Hermitian矩阵等结构矩阵，至今还没有找到有效的快速算法．用传统的算 法计算矩阵乘积w-AP，其 中A、P都是nXn矩阵，需要，z次乘法运算 以及 ，l一Ft次加法运算， 即大约需要2n浮点运算．因此，当讨论的结构矩阵为非奇异矩阵时，为了减少计算量和节省内存， 设计快速的算法显得尤其重要．本文运用广义中心对称矩阵和广义中心 Hermitian矩阵的约化性质 得到计算此类矩阵乘积的Strassen算法，有效地减少 了计算量． 文中，C 表示所有n×m复矩阵集，R 表示所有 tlxm实矩阵集，，表示反对角线上元素为 l其余元素为 0的置换矩阵，C／表示不相交对换的固定积，特别地 表示对应于 的n阶置换矩阵． 定义 屯” 设矩阵A∈C ，如果 I，A，=A成立 ，则称 A为中心对称矩阵；如果 ．，A．，=一A 成立，则称A为斜中心对称矩阵；如果．，A，‘=A成立，则称A为中心 Hermitian矩阵，其中A表 示矩阵A的共轭． 收稿 日期 ：2008—07—16 作者简介 ：曹寒 冬 (1980一 )，男 ，湖 南郴 州人 ，助教 ，硕士 ，研 究方 向 ：数值代数 ．E-mail ahan love1314@126．com． — 第22卷 第 4期 曹寒冬等：儿个结构矩阵乘积的Strasse