你知道怎么改进Python 二分法和牛顿迭代法求算术平方根吗

更新时间：2021年08月31日 18:10:52 作者：Hann Yang

这篇文章主要介绍了Python编程实现二分法和牛顿迭代法求平方根代码的改进，具有一定参考价值，需要的朋友可以了解下，希望能够给你带来帮助

二分法

def sqrtb(n):
    if n<0: raise ValueError('n>=0')
    left,right,x=0,n,n/2
    while not -1e-15<x*x-n<1e-15:
        if x*x>n:
            right,x = x,left+(x-left)/2
        else:
            left,x = x,right-(right-x)/2
    return x

求最接近算术平方根的整数

def sqrtB(x):
    if x==0: return 0
    #y,x=x,round(x)
    left,right,ret = 1,x,0
    while left<=right:
        mid = left + (right-left)//2
        if mid<x/mid:
            left = mid+1
            ret = mid
        elif mid==x/mid:
            ret = mid
            break
        else:
            right = mid-1
    return ret

>>> sqrtB(9)
3
>>> sqrtB(8)
2
>>> sqrtB(9.2)
3.0
>>> sqrtB(7.8)
2.0
>>> sqrtB(4)
2
>>>

二分法原理

牛顿迭代法

def sqrtn(n):
    if n<0: raise ValueError('n>=0')
    x = n/2
    while not -1e-15<x*x-n<1e-15:
        x = (x+n/x)/2
    return x

一点小改进：不用1e-15来比较

def sqrt2(n):
    x = n
    while x*x>n:
        x = (x+n/x)/2
    return x

缺点：碰到n=7,13,...等，会进入死循环

增加判断跳出循环：

def sqrt(n):
    x = n
    while x*x>n:
        y,x = x,(x+n/x)/2
        if y==x: break
    return x

# sqrt(n) n=1~25的精度测试：

0.0
-2.220446049250313e-16
0.0
0.0
0.0
0.0
0.0
-4.440892098500626e-16
0.0
-4.440892098500626e-16
0.0
0.0
4.440892098500626e-16
0.0
0.0
0.0
0.0
8.881784197001252e-16
-8.881784197001252e-16
0.0
0.0
0.0
0.0
0.0
0.0
>>>

牛顿迭代法原理

从函数意义上理解：要求函数f(x)=x²，使f(x)=num的近似解，即x²-num=0的近似解。

从几何意义上理解：要求抛物线g(x)=x²-num与x轴交点（g(x)=0）最接近的点。

假设g(x0)=0，即x0是正解，让近似解x不断逼近x0，x0 ~ x - f(x)/f'(x)

def cubeN(n):
    x,y = n/3,0
    while not -1e-15<x-y<1e-15:
        y,x = x,(2/3)*x+n/(3*x*x)
    return x
'''
>>> cubeN(27)
3.0
>>> cubeN(9)
2.080083823051904
>>>
'''

总结

本篇文章就到这里了，希望能够给你带来帮助，也希望您能够多多关注脚本之家的更多内容！

您可能感兴趣的文章:

自己使用总结Python程序代码片段
这篇文章主要介绍了自己使用总结Python程序代码片段,本文收集了如反向读取文件、往文件中所有添加指定的前缀、匿名函数作为返回值、将二进制数转为10进制数等实用代码片段,需要的朋友可以参考下
2015-06-06
Windows下Pycharm远程连接虚拟机中Centos下的Python环境(图文教程详解)
由于最近学习tensorflow的需要，tensorflow是在Linux环境下，使用的是Python。为了方便程序的调试，尝试在Windows下的Pycharm远程连接到虚拟机中Centos下的Python环境，感兴趣的朋友跟随小编看看吧
2020-03-03
python 基于空间相似度的K-means轨迹聚类的实现
这篇文章主要介绍了python 基于空间相似度的K-means轨迹聚类的实现，文中通过示例代码介绍的非常详细，对大家的学习或者工作具有一定的参考学习价值，需要的朋友们下面随着小编来一起学习学习吧
2021-03-03
Python中的通函数numpy.ufunc详解
这篇文章主要介绍了什么是通函数numpy.ufunc,简单说就是numpy的函数，因为numpy针对的是数组张量，因此，几乎每一个函数都是ufunc。本文针对ufunc的属性进行研究，需要的朋友可以参考下
2023-04-04
5个Python中实现文字转语音模块的使用讲解
这篇文章主要为大家详细介绍了5个Python中实现文字转语音模块的使用,文中的示例代码讲解详细,对我们深入掌握Python有一定的帮助,需要的可以参考下
2023-11-11
python字符串替换re.sub()方法解析
这篇文章主要介绍了python字符串替换re.sub()方法解析,文中通过示例代码介绍的非常详细，对大家的学习或者工作具有一定的参考学习价值,需要的朋友可以参考下
2019-09-09
Python supervisor强大的进程管理工具的使用
这篇文章主要介绍了Python supervisor强大的进程管理工具的使用，本文主要跟大家分享在类unix操作系统下supervisor的使用以及一些关于进程的知识，具有一定的参考价值，感兴趣的小伙伴们可以参考一下
2019-04-04
python引用DLL文件的方法
这篇文章主要介绍了python引用DLL文件的方法,涉及Python调用dll文件的相关技巧,需要的朋友可以参考下
2015-05-05
python中break、continue 、exit() 、pass终止循环的区别详解
这篇文章主要介绍了python中break、continue 、exit() 、pass终止循环的区别详解，文中通过示例代码介绍的非常详细，对大家的学习或者工作具有一定的参考学习价值，需要的朋友们下面随着小编来一起学习学习吧
2019-07-07
pytorch中的scatter_add_函数的使用解读
这篇文章主要介绍了pytorch中的scatter_add_函数的使用解读，具有很好的参考价值，希望对大家有所帮助。如有错误或未考虑完全的地方，望不吝赐教
2023-06-06