求解 x 2 = C x^2 = C x2=C等价于求解 f ( x ) = x 2 − C f(x)=x^2 - C f(x)=x2−C的零点。注意牛顿迭代法不同于牛顿法,牛顿法进行一阶展开而牛顿迭代法进行二阶展开。
- 牛顿法,根据当前点
x
n
x_n
xn求解下一个点
x
n
+
1
x_{n+1}
xn+1:
对 f ( x ) f(x) f(x)进行二阶展开 f ( x ) = f ( x n ) + f ′ ( x n ) ( x − x n ) + f ′ ′ ( x n ) 2 ( x − x n ) 2 f(x) = f(x_n) + f'(x_n)(x-x_n) +\frac{f''(x_n)}{2}(x-x_n)^2 f(x)=f(xn)+f′(xn)(x−xn)+2f′′(xn)(x−xn)2,进行求导的到 f ′ ( x ) = f ′ ( x n ) + f ′ ′ ( x n ) ( x − x n ) = 0 f'(x) = f'(x_n)+f''(x_n)(x-x_n)=0 f′(x)=f′(xn)+f′′(xn)(x−xn)=0,进一步得到 x = x n − f ′ ( x n ) f ′ ′ ( x n ) x=x_n - \frac{f'(x_n)}{f''(x_n)} x=xn−f′′(xn)f′(xn)。这里其实求的是极值点。 - 牛顿迭代法,进行一节展开,求近似零点。
对 f ( x ) f(x) f(x)进行一阶展开 f ( x ) = f ( x n ) + f ′ ( x n ) ( x − x n ) = 0 f(x) = f(x_n) + f'(x_n)(x-x_n) = 0 f(x)=f(xn)+f′(xn)(x−xn)=0,进一步得到 x = x n − f ( x n ) f ′ ( x n ) x=x_n - \frac{f(x_n)}{f'(x_n)} x=xn−f′(xn)f(xn)。这里其实求的是极值点。带入 f ( x ) = x 2 − C f(x)=x^2 - C f(x)=x2−C,得到 x = 1 2 ( x n + c x n ) x=\frac{1}{2}(x_n + \frac{c}{x_n}) x=21(xn+xnc)