你知道怎么改进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的三种可变参数

    这篇文章主要介绍了Python的三种可变参数,文中通过示例代码介绍的非常详细,对大家的学习或者工作具有一定的参考学习价值,需要的朋友们下面随着小编来一起学习学习吧
    2019-05-05
  • Python自动打印被调用函数变量名及对应值 

    Python自动打印被调用函数变量名及对应值 

    这篇文章主要介绍了Python自动打印被调用函数的变量名及对应的值,文章围绕主题展开详细的内容介绍,具有一定的参考价,需要的朋友可以参考一下 
    2022-07-07
  • python使用redis模块来跟redis实现交互

    python使用redis模块来跟redis实现交互

    这篇文章主要介绍了python使用redis模块来跟redis实现交互,文章围绕主题展开详细的内容介绍,具有一定的参考价值,需要的小伙伴可以参考一下
    2022-06-06
  • pandas之分组统计列联表pd.crosstab()问题

    pandas之分组统计列联表pd.crosstab()问题

    这篇文章主要介绍了pandas之分组统计列联表pd.crosstab()问题,具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑完全的地方,望不吝赐教
    2023-09-09
  • Python实现视频下载与合成的示例代码

    Python实现视频下载与合成的示例代码

    这篇文章主要为大家详细介绍了Python是如何实现视频的下载以及合成的,文中的实现步骤讲解详细,感兴趣的小伙伴快跟随小编一起学习一下吧
    2022-04-04
  • 利用python爬取斗鱼app中照片方法实例

    利用python爬取斗鱼app中照片方法实例

    最近在学习python,通过实践是学习的一个好办法,下面这篇文章就来给大家介绍了关于利用python爬取斗鱼app中照片的相关资料,文中通过示例代码介绍的非常详细,对大家的学习或者工作具有一定的参考学习价值,需要的朋友下面来一起看看吧。
    2017-12-12
  • python Matplotlib画图之调整字体大小的示例

    python Matplotlib画图之调整字体大小的示例

    本篇文章主要介绍了python Matplotlib画图之调整字体大小的示例,小编觉得挺不错的,现在分享给大家,也给大家做个参考。一起跟随小编过来看看吧
    2017-11-11
  • Python Setuptools的 setup.py实例详解

    Python Setuptools的 setup.py实例详解

    setup.py是一个 python 文件,它的存在表明您要安装的模块/包可能已经用 Setuptools 打包和分发,这是分发 Python 模块的标准。 它的目的是正确安装软件,本文给大家讲解Python Setuptools的 setup.py感兴趣的朋友跟随小编一起看看吧
    2022-12-12
  • Python实现SMTP发送邮件详细教程

    Python实现SMTP发送邮件详细教程

    这篇文章主要为大家详细介绍了Python实现SMTP发送邮件详细教程,附SMTP邮件发送代码,文中示例代码介绍的非常详细,具有一定的参考价值,感兴趣的小伙伴们可以参考一下
    2016-08-08
  • Python 抓取数据存储到Redis中的操作

    Python 抓取数据存储到Redis中的操作

    这篇文章主要介绍了Python 抓取数据存储到Redis中的操作,具有很好的参考价值,希望对大家有所帮助。一起跟随小编过来看看吧
    2020-07-07

最新评论