非线性方程组一个迭代解法.pdfVIP

维普资讯 第33卷第3期 大 连 理 工 大 学 学 报 Vo1．33．NO．3 1993年 5月 JournalofDallanUniversityofTechnology M By 1993 l )非线性方程组的一个迭代解法 罗远诠 刘 悦 (应用数学系 ) (电磁工程系 ) 摘 要 给 出了一个解 阶非线性方程组 的具有三阶收敛速度的迭代法 ，它可 看成解单个非线性方程的抛物线迭代法的推广．其一次迭代所需工作量是 牛顿迭代法的 1~2／n倍．当一阶导数阵奇异时计算也可进行． √ 蓑分类号：：方02程41移‘性；兰 三竖竺兰度一 才1；我 掂丝 ． 7 l＼ 、 ’ … 0 前 言 考虑非线性方程组 ，(z)= 0 (1) 其中；向量 z=( ， 。，…， )，向量函数 ，(z)一(fl )，^ )，…， ))．至今求解 (1)的各种数值解法中，牛顿迭代法及在它基础上发展起来 的各种迭代法仍是最常用的方法 ； 这些方法实际上都用到且只用到 ，(z)的一阶导数．对于 =1即单个方程的情形 ，使用 ，(z) 的二阶导数的一些迭代法的收敛性及收敛速度都比牛顿法好．很 自然不少作者都想把一维 的 这些方法推广到 ”维的情形 “．但是对 维情形 ， )的二阶导数的计算工作量太大，因此 目 前这些方法还不实用．本文提出一个方法 ；它只用到牛顿方 向上的二阶导数 ，因而计算量增 加不多．但是收敛速度可达到三阶，而且能适应 ， )存在奇异点的情形．数值例子中，收敛 域也比牛顿法有明显的扩大．当 =1时，本文方法退化为Muller法。，因此本方法可看成是 一 维情形的抛物线法 (Muller法)的一种推广． 1 一个迭代算法 设z是式 (1)的一个解， 是第 次迭代值 ，(m，毗，…， )一一，() f(x)／ fl (x) f(x) ；又设 工是使 I乱I一 I4I的整数 ·作变换 收稿 日期：1992—012O；修订 日期 ；1993—03一O1 维普资讯 大 连 理 工 大 学 学 报 33卷 (2) 其中矩阵 dL … I~L-- 乱· aL+ dL十 一 aL 式 (2)的逆变换 一 一 “ (3) 其中：口=(，)．当 jA睁工时，—maj／(aL~a；)1／乱；当 ≠ Vi=L时 口= 岛／ (aL= )．令 F()一f(x+vt)，则式 (1)变为 Fn)一 0 (4) 显然式 (4)有解 t= ’一 )．将式 (4)按泰勒公式展开 ，考虑到 F(。] ，( )， (O)= ， ()口，得到 F0)一 ，( )+ ， ( )vt上 1／Z ·a2F(O)／at[·畦+ R0) (5) 其中：tL是t的第 ￡个分量 ；R()是余项 ，它包括除畦外的其余二次项及更高阶的