g-r循环矩阵求逆的Euclid算法.pdfVIP

第 13卷 第 5期 宁波教育学院学报 Vo1．13 No．5 2011年 10月 JOURNALOFNINGBO INSTITUTEOFEDUCATION 0ct．2011 g-r循环矩阵求逆的Euclid算法 郑多闻，翁晓春 ，张荣娥 (宁波大学理学院，浙江 宁波 315211) 摘 要 ：利用交换环 的同构理论 ，结合多项式最大公因式的EuC11d算法 ，给出了求 g—r循环矩阵逆矩阵的一种 新算法，并结合数值例子给出了该算法的应用 ． 关键词 ：r--循环矩阵；g—r循环矩阵；逆矩阵；算法 中图分类号 ：0151．2 文献标识码 ：A 文章编号 ：1009—2560(2011)05一o091—04 在矩阵理论 中，循环矩阵是一类在数学和物理上有着广泛应用的特殊矩阵，在编码理论、计算机 时序分析、信号和数字图像处理、分子振动和动力系统中，许多实际问题都与周期性有关 ，导致出现了 许多特殊的循环矩阵，诸如r一循环矩阵、g一循环矩阵、对角因子循环矩阵、分块循环矩阵等一些特殊循 环矩阵，且已有许多研究成果 l-【6]．文3[]提出了一类新的循环矩阵称为 g—r循环矩 阵，研究了它的若干特 征性质 ，给出了g—r循环矩阵求逆的快速富里叶算法 ，本文利用多项式的Euclid算法 ，得到了g—r循环 矩阵求逆的另一种快速算法． 1概念与引理 以下设 C…为复数域 C上所有 n阶方阵组成的集合． 定义 1m：设矩阵A=(aij)∈C… ，若满足(1)当j≥i时，aij=_-i；(2)当ji时，aij=ra+i_i，其中i,j=0，1， … ， n一1且 r为常数 ，则称 A为 n阶r一循环矩阵，简记为A=C (a0，a -．，an一 )． 定义2囝：设矩阵A=(aj)∈C…，若满足(1)当i=0时，aoi=a~；(2)当i#0时，ra~=ai4 ，其中j--0，1，…，n-1 且 g为非负整数，j—g转化为模 n的最小非负剩余，则称A为n阶g-循环矩阵，记为A=C(ao，a --，a )． 定义 3 ：设矩阵A=(aij)∈C… ，若满足 (1)当i=0时，aoj=aj；(2)当i≠O时，aij--ai (j≥g)或 aij= rani(jg)，其 中j=0，1，…，n一1且 r为常数，g为非负整数，j—g转化为模 n的最小非负剩余，则称A为 n阶g—r循环矩阵，记为A=C：(a0，a 一，an一)． 特别 ，当g=r=1时为通常 的循环矩阵，此时 r，g可忽略．当r=1时为 g一循环矩阵，此时r可忽 略；当g=1时为r一循环矩阵，此时g可忽略． 从上述定义可知，g—r循环矩阵是r一循环矩阵和g一循环矩阵的结合，是分别对 r一循环矩阵和 g一 循环矩阵进行拓广．它是这样形成 的矩阵：首行元素为 a。，a 一，a ，从第二行起 ，每一行元素都是 由 前一行元素经过从左至右循环地移动g个元素，而且每次从右边移到左边的g个元素须乘上参数 r． 收稿 日期 ：2011-09—06 基金项 目：浙江省教育厅科研项 目，课题编号：Y200803551；宁波大学科研项 目。课题编号：XK1059。 作者简介：郑多1~(1988一)，女，浙江温州人，宁波大学理学院2008级本科生。 通讯作者 ：张荣~(1965一)，女 ，浙江奉化人，宁波大学理学院副教授，主要从事矩阵代数研究。 91 宁波教育学 院学报 2011年第5期 0 l 0 … 0 O O 1 … 0 _ ● ● ● ● ● 记 一= = C (0，1