这篇文章主要讲解了“怎么改进Python二分法和牛顿迭代法求算术平方根”,文中的讲解内容简单清晰,易于学习与理解,下面请大家跟着小编的思路慢慢深入,一起来研究和学习“怎么改进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二分法和牛顿迭代法求算术平方根这一问题有了更深刻的体会,具体使用情况还需要大家实践验证。这里是亿速云,小编将为大家推送更多相关知识点的文章,欢迎关注!
亿速云「云服务器」,即开即用、新一代英特尔至强铂金CPU、三副本存储NVMe SSD云盘,价格低至29元/月。点击查看>>
免责声明:本站发布的内容(图片、视频和文字)以原创、转载和分享为主,文章观点不代表本网站立场,如果涉及侵权请联系站长邮箱:is@yisu.com进行举报,并提供相关证据,一经查实,将立刻删除涉嫌侵权内容。