二分法开方sqrt函数编写

约定不等于承诺〃 2023-01-06 11:36 141阅读 0赞

#include "stdio.h"
    
    /************************************************************************* * 函数名称：abs_float * 功能说明：求绝对值 * 参数说明：dataIn 输入数据 * 函数返回：绝对值 * 修改时间：2021年01月08日 * 备 注： *************************************************************************/
    float abs_float(float dataIn)
    { 
        return (dataIn>=0)?dataIn:(-dataIn);
    }
    
    /************************************************************************* * 函数名称：my_sqrt * 功能说明：开方 * 参数说明：dataIn 输入数据 * 参数说明：acc 精度 * 函数返回：绝对值 * 修改时间：2021年01月14日 * 备 注： *************************************************************************/
    float my_sqrt(float DataIn)
    { 
        #define X2_DataIn(x) (x*x-DataIn)
        #define my_sqrt_acc 1e-6
    
        float start;
        float end;
        float mid;
    
        if(DataIn<0) return -1;
        if(DataIn==0) return 0;
        if(DataIn==1) return 1;
        if(DataIn<1)
        { 
            start = DataIn;
            end = 1;
        }
        else if(DataIn>1)
        { 
            start = 1;
            end = DataIn;
        }
    
        mid = start + (end - start) / 2.0;
        while(abs_float(start - mid) > my_sqrt_acc)
        { 
            if(X2_DataIn(mid)==0)
                return mid;
            else if(X2_DataIn(mid)*X2_DataIn(start)<0)
                end = mid;
            else
                start = mid;
            mid = start + (end - start) / 2.0;
        }
        return mid;
    }
    
    /************************************************************************* * 函数名称：main * 功能说明：主函数 * 参数说明：无 * 函数返回：0 * 修改时间：2021年01月14日 * 备 注： *************************************************************************/
    int main(void)
    { 
    	float a[10] = { 0, 1, 0.25, 2, 3, 5, 7, 0.18, 0.04, 25};
    	int i;
    	for(i=0; i<10; i++)
    	{ 
    		printf("%f\n", my_sqrt(a[i]));
    	}
    	return 0;
    }

![在这里插入图片描述][watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3dlaXhpbl80Mzk2NDk5Mw_size_16_color_FFFFFF_t_70]

--------------------

**相关文章：**

*  [sqrt函数实现（多种算法比较）][sqrt]
 *  [一个Sqrt函数引发的血案][Sqrt]
 *  [sqrt函数实现（神奇的算法）][sqrt 1]
 *  [sqrt函数实现（神奇的算法）][sqrt 2]
 *  [Best Square Root Method - Algorithm - Function (Precision VS Speed)][Best Square Root Method - Algorithm - Function _Precision VS Speed]
 *  [金庸老先生笔下的数学谜题][Link 1]
 *  [面试题：实现一个开方函数][Link 2]
 *  [【算法 | Python】二分法求函数在区间\[a,b\]中的根][_ Python_a_b]

[watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3dlaXhpbl80Mzk2NDk5Mw_size_16_color_FFFFFF_t_70]: /images/20221119/78df85feb74e48048de8f6bc1babc7ff.png
[sqrt]: https://blog.csdn.net/weixin_43964993/article/details/112618194
[Sqrt]: https://diducoder.com/sotry-about-sqrt.html
[sqrt 1]: https://blog.csdn.net/qq_26499321/article/details/73724763
[sqrt 2]: https://blog.csdn.net/weixin_43964993/article/details/112604518
[Best Square Root Method - Algorithm - Function _Precision VS Speed]: https://www.codeproject.com/Articles/69941/Best-Square-Root-Method-Algorithm-Function-Precisi
[Link 1]: https://k.sina.cn/article_2897328623_acb1b9ef01900cfkm.html
[Link 2]: https://blog.csdn.net/mynameisbill2/article/details/100185542
[_ Python_a_b]: https://blog.csdn.net/weixin_43964993/article/details/108351554