关于开方的迭代算法

去年某司有一道面试题是不用 开根号计算开方，当时第一反应就是迭代逼近，但是对方还要求能不能把速度加快，先稍作总结。

迭代逼近意思就是 我先随便取一个数，如果平方大于给定 N，那么就取一个小一点的；如果小了，那就取个大一点，最后达到精度要求即可。

实现代码：

class solution{
public： float my_sqrt(float N){
       float x=num/2;
float x2=num,x1=0;//要是x取大了就去小的x1，反正取x2

float delta=0.0001;//精度

while(!(fabs(x-x1)<delta || fabs(x-x2)<delta))
{

if (x*x>num)
{
x2=x;
x=(x+x1)/2;
}
else(x*x<num)
{
x1=x;
x=(x+x2)/2;
}

}
return x;
}

快速算法就是牛顿迭代法：

f(x)=pow(x,2)-N;找到f(x)=0的x即可。函数梯度为2x，所以y(k+1)-y(k)=g(x(k))(x(k+1)-x(k)),其中g(x(k)) 为在x(k)处的梯度。

推导得到x(k+1)=(x(k)+N/x(k))/2;

class solution{
public： float my_sqrt(float N){
        float x=num/2，x1;
float delta=0.0001;//精度

while(!(fabs(x-x1)<delta))
{


x1=x;
x=(x+N/x)/2;


}
return x;
}

因为x的初值选择会对迭代次数有影响，所以还有一些关于初值选择的研究。