4
$\开始组$

我经常在材料计算中使用Newton-Raphson方法,在那里我必须解一组小的非线性方程(大小=1.5)。在大多数情况下,它是有效的。然而,在极端情况下,如振荡和发散,常常会遇到收敛失败。问题是,给定一组方程和一个初始起点,有没有办法分析解的存在性?

$\端组$
  • $\开始组$ 我相信你想问的是,是否有可能提前知道该方法是否收敛,对吗?或者你真的在问如何确定解决方案何时存在? $\端组$
    – 参宿七
    评论 5月22日12:12
  • $\开始组$ @Rigel好问题。我会说后者,如果解决方案存在的话。 $\端组$
    – 卡什坦
    评论 5月22日12:22
  • $\开始组$ 关于牛顿方法的一个很好的参考是P.Deufhard的《非线性问题的牛顿方法》一书。 $\端组$ 评论 5月28日15:35

3个答案

重置为默认值
7
$\开始组$

牛顿方法的收敛性意味着根的存在,但函数和初始猜测必须满足一些条件才能保证牛顿方法的局部收敛。康托洛维奇定理保证了牛顿方法的局部收敛性(这意味着根的存在)。这些假设有点技术性,但在牛顿方法的分析中是标准的。假设的形式和定理的精确表述因来源而异,这里我陈述的内容与1995年凯利发现的内容类似。

假设1

存在常量$\测试版$,美元\eta$,$\bar{r}$、和$\伽马$具有$\beta\eta\gamma\leq\frac{1}{2}$$x_0\in\mathbb{R}^n$这样的话

  1. $F(美元)$在处可微分$x_0美元$、和$$\|F'(x_0)^{-1}\|\leq\beta,\\text{和}\|F'(x_0)^{-1}华氏度(x_0)\|\leq\eta。$$
  2. $F'(美元)$Lipschitz连续且Lipschit常数为$\伽马$在一个半径为$\bar{r}\geqr_-$关于$x_0美元$哪里$$r_-=\frac{1-\sqrt{1-2\β\β\γ}}{β\γ}。$$

在英语中,这些假设大致如下:

  1. 线性近似的根$F(美元)$蚂蚁$x_0美元$是唯一的,并且线性近似的根位于$x_0美元$距离不太远$x_0美元$

  2. 高阶导数$F(美元)$在大约$x_0美元$所以当我们在线性近似的根上做同样的程序时,第一点仍然为真。

事实上,值得注意的是,牛顿方法和各种修正的根和精确误差估计的存在可以通过简单的估计来保证$F(美元)$及其最初的几个衍生产品。

定理的表述是

定理(坎托罗维奇)

让假设1成立并让$\beta\eta\gamma\leq\frac{1}{2}$并定义$\mathcal美元{B} _0(0)=\{x\in\mathbb{R}^n:\|x-x_0\|\leq R_-\}$.然后有一个唯一的根x美元^*$属于$F(美元)$在里面$\mathcal美元{B} _0(0)$,牛顿迭代$x_0美元$当初始迭代收敛到$x美元^*$、和迭代$x_k\in\mathcal美元{B} _0(0)$为所有人千美元$此外,x美元^*$是的唯一根$F(美元)$在半径范围内$$r=\min\left\{\bar{r},\frac{1+\sqrt{1-2\beta\eta\gamma}}{beta\gamma{right\}$$关于$x_0美元$和错误$e_k(_k)$满足估计$$\|e_k\|leq\frac{(2\beta\eta\gamma)^{2^k}}{2^k\beta\gamma}。$$也就是说,牛顿的方法是二次收敛的。

实际上,你很少能保证先验的这些假设都得到了满足,但它们可以定量地衡量“局部”对于牛顿方法的局部收敛意味着什么,以及这种局部性是如何受到诸如$F'(美元)$或大型二阶导数。最大化收敛机会的最实际方法是通过缩放或其他方式减少问题的条件反射,或者简单地进行更好的初始猜测。显然,“猜测更接近正确答案”是一个循环建议,但如果方程组对应的问题有一些现有的领域知识,那么可以利用这一点进行更好的初始化。例如,根据文献中的参数值或简化模型或缺少物理特性的模型预测的值进行初始化可能会非常有效。有一些领域不赞成这种做法(CFD浮现在脑海中),因为它会使结果产生偏差,但如果你在一开始就遇到了问题,那么它是有效的。

凯利,C.T。,线性和非线性方程的迭代方法,应用数学前沿。16.宾夕法尼亚州费城:SIAM,工业和应用数学学会。xiii,165页(1995年)。ZBL0832.65046号

$\端组$
$\开始组$

除了阻尼牛顿法或信赖域法之外,对于小问题,你可能会发现斯梅尔阿尔法理论或值得研究的数字认证的一般主题。Smale所做的是表明,在某些情况下,可以根据函数及其高阶导数的大小,判断牛顿方法在特定点是否处于二次收敛域。这给出了一个充分但不是必要的条件,所以它可以工作,但不是完全等价的。

$\端组$
2
$\开始组$

你可能想试试阻尼牛顿法:而不是采取完整的步骤($x{k+1}=xk+d_k$,其中$d_k(千美元)$是牛顿步),找到$\alpha\in(0,1])$这样的话$x_k+\αd_k$充分降低了方程的范数。请参见幻灯片11:https://www.stat.cmu.edu/~ryantibs/convexopt-S15/讲座/14-newton.pdf

$\端组$

你的答案

单击“发布您的答案”,表示您同意我们的服务条款并确认您已阅读我们的隐私政策

不是你想要的答案吗?浏览标记的其他问题问你自己的问题