主題
Search

牛頓法


牛頓法,也稱為牛頓-拉夫森方法,是一種求根演算法,它使用函式 f(x) 在疑似附近的泰勒級數的前幾項。牛頓法有時也稱為 牛頓迭代法,儘管在本文中,後一個術語保留用於應用牛頓法計算平方根

對於f(x),一個多項式,牛頓法本質上與霍納法相同。

泰勒級數 f(x) 關於點 x=x_0+epsilon 由下式給出

 f(x_0+epsilon)=f(x_0)+f^'(x_0)epsilon+1/2f^('')(x_0)epsilon^2+....
(1)

僅保留到一階項,

 f(x_0+epsilon) approx f(x_0)+f^'(x_0)epsilon.
(2)

方程 (2) 是曲線在 (x_0,f(x_0)) 處的切線方程,所以 (x_1,0) 是該切線x 軸的交點。因此,圖表可以很好地直觀地說明為什麼牛頓法在良好選擇的起始點處有效,以及為什麼它可能在選擇不當的起始點處發散。

上面的表示式可以用來估計從初始猜測 x_0 開始,為了更接近根所需的偏移量 epsilon。令 f(x_0+epsilon)=0 並求解 (2) 中的 epsilon=epsilon_0 得到

 epsilon_0=-(f(x_0))/(f^'(x_0)),
(3)

這是一階調整到的位置。透過令 x_1=x_0+epsilon_0,計算新的 epsilon_1,等等,可以重複該過程,直到它使用下式收斂到不動點(這恰好是一個根)

 epsilon_n=-(f(x_n))/(f^'(x_n)).
(4)

不幸的是,此過程在水平漸近線區域性極值附近可能不穩定。但是,如果的位置的初始選擇良好,則可以迭代應用該演算法以獲得

 x_(n+1)=x_n-(f(x_n))/(f^'(x_n))
(5)

對於 n=1, 2, 3, .... 提供牛頓法安全收斂的初始點 x_0 稱為近似零點

牛頓法可以在 Wolfram 語言中實現為

  NewtonsMethodList[f_, {x_, x0_}, n_] :=
    NestList[# - Function[x, f][#]/
      Derivative[1][Function[x, f]][#]& , x0, n]

假設牛頓迭代 x_(k+1) 收斂於 x^*f^'(x^*)!=0,並定義第 k 步後的誤差為

 x_k=x^*+epsilon_k.
(6)

f(x_k)x^* 附近展開得到

f(x_k)=f(x^*)+f^'(x^*)epsilon_k+1/2f^('')(x^*)epsilon_k^2+...
(7)
=f^'(x^*)epsilon_k+1/2f^('')(x^*)epsilon_k^2+...
(8)
f^'(x_k)=f^'(x^*)+f^('')(x^*)epsilon_k+....
(9)

但是

epsilon_(k+1)=epsilon_k+(x_(k+1)-x_k)
(10)
=epsilon_k-(f(x_k))/(f^'(x_k))
(11)
 approx epsilon_k-(f^'(x^*)epsilon_k+1/2f^('')(x^*)epsilon_k^2)/(f^'(x^*)+f^('')(x^*)epsilon_k).
(12)

取二階展開式

 (aepsilon+1/2bepsilon^2+cepsilon^3)/(a+bepsilon+depsilon^2) approx epsilon-b/(2a)epsilon^2
(13)

得到

 epsilon_(k+1) approx (f^('')(x^*))/(2f^'(x^*))epsilon_k^2.
(14)

因此,當該方法收斂時,它是二次收斂的。

NewtonsMethodBasins

將牛頓法應用於任何二次或更高次多項式的根會產生一個 C 的有理對映,並且當存在三個或更多不同的根時,此對映的 Julia 集是一個分形。迭代方法求解 z^n-1=0 的根,起始點為 z_0 得到

 z_(i+1)=z_i-(z_i^n-1)/(nz_i^(n-1))
(15)

(Mandelbrot 1983, Gleick 1988, Peitgen and Saupe 1988, Press et al. 1992, Dickau 1997)。由於這是一個 n多項式,因此存在 n 個演算法可以收斂到的。為每個吸引盆(收斂到相同的初始點 z_0 的集合)著色不同的顏色,然後給出上面的圖。

NewtonsMethodCross

分形通常也從非多項式對映中產生。上面的圖顯示了函式 z^2-2^z (D. Cross, 私人通訊, 1月 10, 2005) 和 z^3-3^z 的牛頓法收斂所需的迭代次數。


另請參閱

Alpha-Test, 近似零點, 不動點, Halley 的無理公式, Halley 法, 霍納法, Householder 法, Laguerre 法, 牛頓迭代法, 牛頓向量場, 求根演算法 在 課堂中探索此主題

使用 探索

參考文獻

Abramowitz, M. and Stegun, I. A. (Eds.). Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables, 9th printing. New York: Dover, p. 18, 1972.Acton, F. S. Ch. 2 in Numerical Methods That Work. Washington, DC: Math. Assoc. Amer., 1990.Arfken, G. Mathematical Methods for Physicists, 3rd ed. Orlando, FL: Academic Press, pp. 963-964, 1985.Boyer, C. B. and Merzbacher, U. C. A History of Mathematics, 2nd ed. New York: Wiley, 1991.Dickau, R. M. "Basins of Attraction for z^5=1 Using Newton's Method in the Complex Plane." http://mathforum.org/advanced/robertd/newtons.html.Dickau, R. M. "Variations on Newton's Method." http://mathforum.org/advanced/robertd/newnewton.html.Dickau, R. M. "Compilation of Iterative and List Operations." Mathematica J. 7, 14-15, 1997.Gleick, J. Chaos: Making a New Science. New York: Penguin Books, plate 6 (following p. 114) and p. 220, 1988.Gourdon, X. and Sebah, P. "Newton's Iteration." http://numbers.computation.free.fr/Constants/Algorithms/newton.html.Householder, A. S. Principles of Numerical Analysis. New York: McGraw-Hill, pp. 135-138, 1953.Mandelbrot, B. B. The Fractal Geometry of Nature. San Francisco, CA: W. H. Freeman, 1983.Newton, I. Methodus fluxionum et serierum infinitarum. 1664-1671.Ortega, J. M. and Rheinboldt, W. C. Iterative Solution of Nonlinear Equations in Several Variables. Philadelphia, PA: SIAM, 2000.Peitgen, H.-O. and Saupe, D. The Science of Fractal Images. New York: Springer-Verlag, 1988.Press, W. H.; Flannery, B. P.; Teukolsky, S. A.; and Vetterling, W. T. "Newton-Raphson Method Using Derivatives" and "Newton-Raphson Methods for Nonlinear Systems of Equations." §9.4 and 9.6 in Numerical Recipes in FORTRAN: The Art of Scientific Computing, 2nd ed. Cambridge, England: Cambridge University Press, pp. 355-362 and 372-375, 1992.Ralston, A. and Rabinowitz, P. §8.4 in A First Course in Numerical Analysis, 2nd ed. New York: McGraw-Hill, 1978.Raphson, J. Analysis aequationum universalis. London, 1690.Smale, S. "On the Efficiency of Algorithms of Analysis." Bull. Amer. Math. Soc. 13, 87-121, 1985.Varona, J. L. "Graphic and Numerical Comparison Between Iterative Methods." Math. Intell. 24, 37-46, 2002.Whittaker, E. T. and Robinson, G. "The Newton-Raphson Method." §44 in The Calculus of Observations: A Treatise on Numerical Mathematics, 4th ed. New York: Dover, pp. 84-87, 1967.

在 上被引用

牛頓法

請引用本文為

Weisstein, Eric W. "牛頓法。" 來自 --一個 資源。 https://mathworld.tw/NewtonsMethod.html

學科分類