论文部分内容阅读
求解非线性方程的根在数值分析中很重要,它广泛应用于工程和其它应用领域的科学计算中。迭代算法是求方程根的众多方法中应用得最为广泛的一类方法,它从某个初始点出发,由迭代格式生成一组收敛于方程根的序列。众所周知,利用目标函数f(x)的高阶导数信息可以更加方便地构造出高阶收敛的迭代格式。然而,在很多情况下计算高阶导数是很费计算时间的,甚至是不可能的。因此,这类方法在实际应用中受到很多限制。对于求解方程的单根,目前已经有了相当多的高阶迭代方法,相应的构造技术手段也比较丰富和成熟。然而,这些技术手段一旦用来构造求解重根的迭代算法时,就显得是相当的复杂甚至无效。另一方面,大量的理论和数值实验都表明,求解单根的高阶迭代算法如果不加修正地直接用来求解方程的重根,其收敛阶数将大为降低(通常只有一阶)。一个简单明了的例子就是经典的牛顿迭代法。它二阶收敛于方程的单根,但只能线性收敛于方程的重根,收敛速度变慢。因此,如何构造出求解方程重根的高阶,尤其最优阶的迭代格式是一项具有挑战性的工作。到目前为止,这方面的成果还不是很丰富。目前绝大多数求解重根的最优阶迭代算法在构造过程中,都利用了方程重根的重数信息。另外,对单根的情形,有许多工作是讨论方法的收敛半径问题。然而对重根的情形,相应的工作却少之又少。这篇论文主要研究求解非线性方程重根的迭代算法。我们首先给出了一类最优的、四阶收敛的迭代格式,该类迭代格式几乎包含了所有已知的最优阶迭代算法。通过对求解单根的多步迭代算法的观察分析,我们发现,它们的共同点是第一步都借助牛顿迭代法或其它具有二阶收敛速度的迭代格式。然而,这种预迭代方法似乎在构造求解重根的迭代算法时并没有充分利用起来。因此,我们尝试给出两类校正步为修正牛顿迭代步的迭代算法。此外,我们还给出了两个不需要知道根的重数的迭代算法。而且,它们也不需要计算目标函数的导数以及“精挑细选”地选取初始值。尽管这两种方法的收敛阶并非最优,但它们的优势是不言而喻的。对迭代算法而言,收敛半径的问题是至关重要。对于求解重根的迭代算法而言,这方面的工作很少。直到最近,Ren和Argyros基于高阶差分和多重积分,给出了修正的牛顿迭代法的收敛半径估计。我们首先利用他们的方法,估算出Osad;迭代算法的收敛半径。进一步,我们提出了一种基于带有积分余项的Taylor展开式的方法,改进了上述收敛半径。这种方法相比于Ren和Argyros给出的方法简单而有效。因此,我们又重新估算了修正牛顿迭代法的收敛半径,同样得到了更好的结果。本文的最后一部分在复平面区域上,讨论了一些求解重根的迭代算法的动力学行为。我们绘制了大量的吸引域分形图并加以分析。结合已有的结果,我们可以发现,本文给出的一些迭代方法在方程根的附近具有较好的动力学行为,很少出现混沌现象。尽管Halley迭代算法表现最好,但它在迭代过程中需要计算函数的二阶导数,这在实际应用中很不现实。需要说明的是,为了验证我们的理论分析结果,在每一节最后,我们都将给出几个数值试验来增强说服力。