牛顿迭代法(最优化理论与方法-牛顿迭代法)

目录:
简介牛顿迭代法
牛顿迭代法公式
牛顿迭代法的收敛性
牛顿法的改进
牛顿法和梯度下降法的区别
一、简介牛顿法
迭代法也称为辗转法,是一种不断用变量的旧值递推新值的过程,跟迭代法相对应的是直接法,即一次性解决问题。但多数方程不存在求根公式,因此求解根非常困难,甚至不可能,从而寻找方程的近似根就显得特别重要。
迭代算法是用计算机解决问题的一种基本方法。它利用计算机运算速度快、适合做重复性操作的特点,让计算机对一组指令或一定步骤进行重复执行,在每次执行这组指令或这些步骤时,都从变量的原值推出它的一个新值。
利用迭代算法解决问题,需要做好以下3方面的工作:
确定迭代变量。在可以用迭代算法解决的问题中,至少存在一个直接或间接地不断由旧值递推出新值的变量,这个变量就是迭代变量。
建立迭代关系式。所谓迭代关系式,是指如何从变量的前一个值推出下一个值的公式或关系。迭代关系式的建立是解决迭代问题的关键,通常可以使用递推或倒推的方法来完成。
对迭代过程进行控制。在什么时候结束迭代过程?这是编写迭代程序必须考虑的问题。不能让迭代过程无休止地重复下去。迭代过程的控制通常可以分为两种情况:(一)所需的迭代次数是个确定的值,可以计算出来。(二)所需的迭代次数无法确定,需要程序员进步一分析出用来结束迭代过程结束的条件。?
牛顿迭代法(Newton’smethod)又称为牛顿-拉夫逊(拉夫森)方法。它是牛顿在17世纪提出的一种在实数域和复数域上近似求解方程的方法。多数方程不存在求根公式,因此求精确根非常困难,甚至不可能,从而寻找方程的近似根就显得特别重要。该方法使用函数f(x)的泰勒级数的前面几项来寻找方程f(x)=0的根。牛顿迭代法是求方程根的重要方法之一,其最大优点是在方程f(x)=0的单根附近具有平方收敛,而且该法还可以用来求方程的重根、复根。此时一定线性收敛,但是可通过一些方法变成超线性收敛。另外该方法广泛用于计算机编程中。
二、牛顿迭代法公式
用迭代法解非线性方程时,如何构造迭代函数是非常重要的,那么怎样构造的迭代函数才能保证迭代法收敛呢?牛顿迭代法就是常用的方法之一,其迭代格式的来源或者理解方式大概有两种方式:
1.设,对在点作泰勒展开:
略去二次项,得到的线性近似式:。由此得到方程0的近似根(假定0),,即可构造出迭代格式(假定0):。这就是牛顿迭代公式,若得到的序列{}收敛于,则就是非线性方程的根。
2.牛顿迭代法也称为牛顿切线法,设ξ是=0的根,选取作为ξ的初始近似值,过点做曲线的切线L,L的方程为,求出L与x轴交点的横坐标,称x1为ξ的一次近似值。过点做曲线的切线,并求该切线与x轴交点的横坐标,称为ξ的二次近似值。重复以上过程,得ξ的近似值序列,其中称为ξ的次近似值,上式称为牛顿迭代公式,如图1左图所示,这就是牛顿切线法的几何解释。实际上,牛顿迭代法也可以从几何意义上推出。利用牛顿迭代公式,由得到,从图1右图上看,就是过点作函数的切线,切线与轴的交点就是,所以有,整理后也能得出牛顿迭代公式:。
三、牛顿迭代法的收敛性
这里我就不在证明了,直接给出结论:
牛顿迭代法具有较高的收敛速度,它的收敛阶数=2;而牛顿迭代法的局部收敛性较强,只有初值充分地接近ξ,才能确保迭代序列的收敛性。如果是连续的,并且待求的零点是孤立的,那么在零点周围存在一个区域,只要初始值位于这个邻近区域内,那么牛顿法必定收敛。并且,如果0,?那么牛顿法将具有平方收敛的性能.?粗略的说,这意味着每迭代一次,牛顿法结果的有效数字将增加一倍。
四、牛顿法的改进
用牛顿法解方程,虽然在单根附近具有较快的收敛速度,但它有个明显的缺点,就是每次都要计算导数,当比较复杂时,计算可能很困难。
从计算导数困难方面,改进的算法有:简化牛顿法。
为了避免频繁地计算导数值,可将它取为固定值,比如在牛顿迭代公式中用代,即在迭代过程中始终保持分母不变,则有简化牛顿迭代公式(或固定斜率切线法):
从初值选择方面,改进的算法有:牛顿下山法。
由牛顿法的收敛性定理知,牛顿法对初始值的选取要求是很高的。一般地说,牛顿法只有局部收敛性。当初始值取得离根太远时,迭代将不收敛,而一旦初始值进入收敛域内,牛顿法就有平方收敛的速度,为了扬长避短,扩大初始值选取的范围。这就是牛顿下山法的思想。
这里,我就不在展开来讲解了。需要了解详细过程的同学,请自行学习。
五、牛顿法与梯度下降法的区别
1.?牛顿法优点:最优化问题中,牛顿法比梯度下降法求解需要的迭代次数更少。
原因:牛顿法是二阶收敛,梯度下降是一阶收敛,所以牛顿法就更快。如果更通俗地说的话,比如你想找一条最短的路径走到一个盆地的最底部,梯度下降法每次只从你当前所处位置选一个坡度最大的方向走一步,牛顿法在选择方向时,不仅会考虑坡度是否够大,还会考虑你走了一步之后,坡度是否会变得更大。所以牛顿法比梯度下降法看得更远一点,能更快地走到最底部。根据wiki上的解释如图2所示,从几何上说,牛顿法就是用一个二次曲面去拟合你当前所处位置的局部曲面,而梯度下降法是用一个平面去拟合当前的局部曲面,通常情况下,二次曲面的拟合会比平面更好,所以牛顿法选择的下降路径会更符合真实的最优下降路径。红色的牛顿法的迭代路径,绿色的是梯度下降法的迭代路径。
图2:牛顿法与梯度下降法的区别
2.牛顿法缺点:
(1)?对目标函数有严格的要求,必须有连续的一、二阶偏导数,海森矩阵必须是正定的。
(2)?计算量大,除计算梯度外,还需要计算二阶偏导矩阵及其逆矩阵。
Reference:
最优化问题中,牛顿法为什么比梯度下降法求解需要的迭代次数更少? - 大饼土博的回答 - 知乎 https://www.zhihu.com/question/19723347/answer/14636244