一类非线性矩阵方程的迭代算法.pdfVIP

第27卷第2期 贵州大学学报(自然科学版) VoL 27N0．2 2010年 4月 JournalofGuizhouUniversity(NaturalSciences) Apr．2010 文章编号 1000—5269(2010)02一OOO4—03 一 类非线性矩阵方程的迭代算法 龙建辉h，何佑梅 ，胡锡炎2 (1．福建工程学院数理系，福建 福州350108；2．湖南大学数学与计量经济学院，湖南 长沙410082) 摘 要：研究了一类在电路设计，信号处理中有广泛应用的非线性矩阵方程。给出了求该方程解 的迭代算法，证明了其收敛性。数值例子说明所给算法是有效的。 关键词：矩阵方程 ；迭代算法；收敛性 中图分类号：0241．5，0241．7 文献标识码 ：A 考虑矩阵方程 Newton法解方程组(3)的算法如下， R=( + ) ( — ) (1) 算法 (1)(Newton法) 其中M是已知的凡×n阶M一矩阵， 是未知 步 1：给定充分好的初始向量 。‘= ( ’， (0) (0) (0) (0) 的凡×n阶且对角元为0的矩阵，是未知的n×n … ， n ， r12 ， … ， r】R ， … ， rnl ， … ， 阶正对角矩阵。在实际应用中 往往是需要的。 r 。‘))； ． 一 l 方程[1]在电路设计，信号处理等领域有广泛的应 步2：解关于 的线性方程组 用，详见文献 [1]．令 m ，分别表示矩阵 ，R的 G( “’)^ =c(x“)； (3) (iJ)元素。 是矩阵 的第 i个对角元，则矩阵方 步3： “¨= “+h． 程(1)可化为如下的非线性代数方程组： 其中G( “)=(g)ER ， r —I IT~ilrli+ —mii 0 (2) ( )，_l’…n2． 【∑ r+m“+f)r一m ：0 Newton法具有局部二次收敛性，它在理论上 是吸引人的．但是，在它的每轮迭代中都需要求 方程组(2)中有n个方程和n个未知数 ( 。阶的线性方程组(3)的解，这可能是困难的工 ≠ )， ．对于方程组 (2)的求解方法一般是用多 作．此外每轮迭代不仅要G( “’)的儿个分量，也 变量的Newton法，割线法等 ]。然而Newton法是 需要计算 G( “)．初始向量若选取不好，则 局部二次收敛，它对初始值的要求较高。所以在应 Newton法的二次收敛得不到保证，甚至不收敛。 用中常用其它的办法求得较好的初始值，然后用 因此Newton在实际应用中有较大的困难． Newton法作进一步的校正。文[3]提供了一个较