使用牛顿法求平方根的示意图如下:
求解某个数n的平方根可以看作是求方程的解。而方程即是上面图片中的二次曲线(蓝色),只要求出该二次曲线与x轴的交点x即求出了n的平方根。而牛顿法的思想就是通过迭代的方法不断逼近x。
首先取任意初值x0,则可以得到在二次曲线上对应于x0点的值y0和斜率,根据点(x0, y0)和斜率k可以求出过(x0, y0)的切线(红色)。可以看出相对于x0,切线与x轴的交点是x更好的一个近似。通过不断的重复该过程可以得到x值一个很好的近似,即求出了n的平方根。
很容易求出过点(x0, y0)斜率为k的直线方程为:,那该直线与x轴的交点x1为:,带入实际的值为:。现在我们就求出了从xn到xn+1的递推公式,只要设置好结束条件,即可求出x的近似值,也就是n平方根的近似值!
下面是根据牛顿法就n的平方根的C++程序:
#include <iostream>
#include <algorithm>
using namespace std;
double mySqrt(double n)
{
double x = 1.0;
while(abs(x * x - n) > 1e-9)
x = (x + n / x) / 2;
return x;
}
int main()
{
cout << mySqrt(2) << endl;
return 0;
}