牛頓-拉夫森算法

牛頓-拉夫森(Newton-Raphson)算法是一種非線性方程數值求根的疊代算法。設非線性方程為f(x) =0，設x₀為閾值，由泰勒公式近似地有: f(x)=f(x₀) +f'(x₀) (x-x₀)， 由此得到求根的一般疊代公式x_(k+1)=x_k-[f(x_k)/f'(x_k)]其中f'(x_k)是函式f(x)在x_k處的導數。在對時間序列MA模型的參數進行矩估計時可以採用Newton-Raphson算法進行疊代。

基本介紹

• 中文名：牛頓-拉夫森算法
• 外文名：Newton-Raphson Algorithm
• 所屬學科：數學
• 簡介：非線性方程數值求根的疊代算法

基本介紹

牛頓-拉夫森算法是對非線性方程作泰勒展開並取一次近似的結果。

在討論非線性最小二乘時，目標函式J(β)對β是非線性的， 如果是第i個疊代值，考慮在附近將J(β)展開成二次函式。令

其中是由J的二階導數組成的矩陣，它的第k行第l列的元素是(是的第k個分量)：

是在附近對J的二階逼近。

現在來找出的穩定點，令

如果是滿秩的，則有解

用(1)作為第i次疊代所定義的算法稱為牛頓-拉夫森算法，它相當於在一般的疊代格式中取。

如果J(β)是參數的二次函式(二次型)，即J(β)與一致，則是J的穩定點，而且當是正定時它是一個最小值點。這時算法是可接受的，並且一次疊代就收斂。由於牛頓-拉夫森算法具有這種性質，所以稱它為二階收斂的，如果是負定的或不定的，則算法是不可接受的。

當J不是二次型，一般與J的穩定點不一致，因此也就不可能一次疊代就收斂，但只要是正定的， 則算法是可接受的。

在實際套用時， 也可以在(1)中加上一個標量因子以加快其收斂速度， 即(1)修改為

適當選擇ρ＞0即可取得較好的效果。

牛頓-拉夫森方法的優缺點

牛頓-拉夫森方法有收斂快的優點，但是它也存在著缺點， 主要是H(β)的正定性不一定能得到保證，同時求H時需要計算二階導數這是很複雜的。為了克服這些困難，提出了各種改進方法。為了克服H的不定型問題提出了麥夸特(Marquardt)方法。為了克服求二階導數的困難提出高斯(Gauss)方法和變尺度方法。