數學之美:牛頓-拉夫遜迭代法原理及其實現

[導讀] 前面剛轉了一篇文章提到了牛頓-拉夫遜（拉弗森）（Newton-Raphson method）方法，感覺這個數學方法很有必要相對深入寫一篇文章來總結分享印證一下自己的理解。這是寫本文的由來，如果發現文章中有錯誤之處，請留言交流討論。

什麼是牛頓-拉夫遜方法？ 牛頓其人：Isaac Newton（1642年12月25日– 1727年3月20日）是一位英國數學家，物理學家，天文學家，神學家和作家，被公認為有史以來最有影響力的科學家之一，並且是科學革命的關鍵人物。他的書《自然哲學的數學原理》於1687年首次出版，奠定了古典力學的基礎。牛頓還為光學做出了開創性的貢獻，並與戈特弗裡德·威廉·萊布尼茲（Gottfried Wilhelm Leibniz）發展了無窮微積分的學科。
牛頓
拉弗森Joseph Raphson 生卒不詳，其最著名的著作是1690年出版的《通用分析方程》。它包含一種方法，現在稱其為牛頓-拉夫森方法，用於近似方程式的求根。艾薩克·牛頓（Isaac Newton）在1671年寫的《通量法》中開發了一個非常相似的公式，但是這項工作要到1736年才出版，這是拉夫森分析之後近50年。但是，該方法的Raphson版本比Newton方法更簡單，因此通常被認為是更好的方法。所以，牛頓迭代法（簡寫）就是一種近似求解實數域與複數域求解方程的數學方法。那麼這個方法是具體是什麼原理呢？
牛頓迭代如何迭代？ 直接看數學公式描述如何迭代不直觀，先來看動圖就很容易理解牛頓迭代法為什麼叫迭代法以及怎樣迭代的：牛頓迭代法是原理是根據一個初始點由前面描述知道，牛頓迭代法是用來近似求解方程的，這裡有兩個點需要說明：
迭代法非常適合計算機編程實現，實際上計算機編程對於牛頓迭代法廣為應用知道一個直線的一個坐標點

#include <stdio.h>
#include <math.h>
#include <stdlib.h>

/*假定待求根函數如下*/
#define    F(x)    (2*(x)*(x)-10*cos(x)+(x)-80)

/*其一階導數為*/
#define   DF(x)    (4*(x)+10*sin(x)+1)

float newton_rooting(float x0,float precision,float min_deltax,int max_iterations)
{
     float xn,xn1,fn,fn1,dfn;
     float deltax;
     int step = 0;
     xn  = x0;
     xn1 = x0;
     do{
       xn  = xn1;
       fn  = F(xn);
       dfn = DF(xn);
       /*判0*/
       if( fabs(dfn) <1e-6 )
       {
            if( fabs(fn)>precision )
                return NAN;
            else
                return fn;
       }

       xn1 = xn - fn/dfn;
       fn1 = F(xn1);
       deltax = fabs(xn1-xn);
       
       step++;
       if( step>max_iterations )
       {
           if( fabs(fn1)<precision )
               return xn;
           else
               return NAN;
       }
     }
     while( fabs(fn1)>precision || deltax>min_deltax );

     return xn1;
}

void main()
{
     float root_guess = 23.0f;
     float precision = 0.00001f;
     float min_deltax = 0.001f;
     float root;
     int step = 7;

     root = newton_rooting( root_guess,precision,min_deltax,step );
     printf("根為: %f,函數值為：%f\n", root,F(root));

     root_guess = -23;
     root = newton_rooting( root_guess,precision,min_deltax,step );
     printf("根為: %f,函數值為：%f\n", root,F(root));
}

根為: 6.457232, 函數值為：0.000004
根為: -6.894969,函數值為：-0.000008

函數值已經很接近於0了，如果還需要更為精確的值，則可以選擇在此基礎上進一步求解，比如利用二分法逼近。需要注意些啥？ 求斜率可能為0，如為0時，則可能找到了函數的極值，比如：如果選擇的初始猜測根的接近方程f(x)=0中函數f(x)的拐點 ，Newton-Raphson方法可能開始偏離根。然後，它可能會又收斂回到根。例如

來源：嵌入式客棧

作者：逸珺

聲明：本文由作者原創，文章內容系作者個人觀點，電子發燒友網轉載僅作為傳達一種不同的觀點，不代表電子發燒友網對該觀點的贊同或支持，如有異議，歡迎聯繫電子發燒友網。