主題
Search

拉蓋爾方法


一種求演算法,它從任何起始位置收斂到一個複數。 為了解釋這個公式,考慮一個n階多項式及其導數,

P_n(x)=(x-x_1)(x-x_2)...(x-x_n)
(1)
P_n^'(x)=(x-x_2)...(x-x_n)+(x-x_1)...(x-x_n)+...
(2)
(3)
=P_n(x)(1/(x-x_1)+...+1/(x-x_n)).
(4)

現在考慮 P_n(x) 的對數和對數導數

ln|P_n(x)|=ln|x-x_1|+ln|x-x_2|+...+ln|x-x_n|
(5)
(dln|P_n(x)|)/(dx)=1/(x-x_1)+1/(x-x_2)+...+1/(x-x_n)
(6)
=(P_n^'(x))/(P_n(x))
(7)
=G(x)
(8)
-(d^2ln|P_n(x)|)/(dx^2)=1/((x-x_1)^2)+1/((x-x_2)^2)+...+1/((x-x_n)^2)
(9)
=[(P_n^'(x))/(P_n(x))]^2-(P_n^('')(x))/(P_n(x))=H(x).
(10)

現在做出“一組相當極端的假設”,即正在尋找的根 x_1 與當前最佳猜測相距 a,因此

 a=x-x_1,
(11)

而所有其他根都位於相同的距離 b,因此

 b=x-x_i
(12)

對於 i=2, 3, ..., n (Acton 1990; Press et al. 1992, p. 365)。 這使得 GH 可以用 ab 表示為

G=1/a+(n-1)/b
(13)
H=1/(a^2)+(n-1)/(b^2),
(14)

同時求解這些方程得到 a

 a=n/(max[G+/-sqrt((n-1)(nH-G^2))]),
(15)

其中符號的選擇是為了使分母的幅度最大。

要應用此方法,請為試探值 a 計算 x,然後使用 x-a 作為下一個試探值,並迭代直到 a 變得足夠小。 例如,對於多項式 4x^3+3x^2+2x+1,起始點為 x_0=-1.0,演算法非常快速地收斂到實根,如(-1.0-0.58113883008419-0.60582958618827)。

設定 n=2 得到 哈雷的無理公式


另請參閱

哈雷的無理公式, 哈雷方法, 牛頓法,

使用 探索

參考文獻

Acton, F. S. Numerical Methods That Work, 2nd printing. Washington, DC: Math. Assoc. Amer., 1990.Adams, D. A. "A Stopping Criterion for Polynomial Root Finding." Comm. ACM 10, 655-658, 1967.Press, W. H.; Flannery, B. P.; Teukolsky, S. A.; and Vetterling, W. T. Numerical Recipes in FORTRAN: The Art of Scientific Computing, 2nd ed. Cambridge, England: Cambridge University Press, pp. 365-366, 1992.Ralston, A. and Rabinowitz, P. §8.9-8.13 in A First Course in Numerical Analysis, 2nd ed. New York: McGraw-Hill, 1978.

在 中被引用

拉蓋爾方法

請引用本文為

Weisstein, Eric W. “拉蓋爾方法。” 來自 --一個 資源。 https://mathworld.tw/LaguerresMethod.html

主題分類