【学习笔记】牛顿迭代法求立方根

float CubeRoot(float num)
{
    float x = num;
    float error = 1e-5;
    
    while (fabs(num - (x * x * x)) >= error)
    {
        x = (2 * x + num / (x * x)) / 3.0;
    }
    
    return x;
}

代码很简单，就是使用牛顿迭代法计算，然后判断是否达到想要的精度，达到精度后退出。

这里精度设置是0.00001，可以根据自己实际情况调整。

公式推导

这里说明代码x = (2 * x + num / (x * x)) / 3.0;是如何得到的。

牛顿迭代法公式如下。

$x_{n+1}=x_n?\frac{f(x_n)}{f^{\prime }(x_n)}$

我们计算立方根的公式：

$y=\sqrt[3]{x}=x^{\frac{1}{3}}$

所以

$y^3=x$

构建以y为自变量的函数方程为

$f(y)=y^3?x=0$

$f^{\prime }(y)=3y^2=0$

将$f(y)$和$f^{\prime }(y)$带入

$y?\frac{f(y)}{f^{\prime }(y)}=y?\frac{y^3?x}{3y^2}$

? $=\frac{1}{3}(2y+\frac{x}{y^2})$