怎么改進(jìn)Python二分法和牛頓迭代法求算術(shù)平方根

這篇文章主要講解了“怎么改進(jìn)Python二分法和牛頓迭代法求算術(shù)平方根”，文中的講解內(nèi)容簡(jiǎn)單清晰，易于學(xué)習(xí)與理解，下面請(qǐng)大家跟著小編的思路慢慢深入，一起來(lái)研究和學(xué)習(xí)“怎么改進(jìn)Python二分法和牛頓迭代法求算術(shù)平方根”吧！

目錄

• 二分法

• 二分法原理

• 牛頓迭代法

• 牛頓迭代法原理

  
  

二分法

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

求最接近算術(shù)平方根的整數(shù)

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

一點(diǎn)小改進(jìn)：不用1e-15來(lái)比較

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

缺點(diǎn)：碰到n=7,13,...等，會(huì)進(jìn)入死循環(huán)

增加判斷跳出循環(huán)：

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的精度測(cè)試：

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
>>>

牛頓迭代法原理

從函數(shù)意義上理解：要求函數(shù)f(x)=x²，使f(x)=num的近似解，即x²-num=0的近似解。

從幾何意義上理解：要求拋物線g(x)=x²-num與x軸交點(diǎn)（g(x)=0）最接近的點(diǎn)。

假設(shè)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
>>>
'''

感謝各位的閱讀，以上就是“怎么改進(jìn)Python二分法和牛頓迭代法求算術(shù)平方根”的內(nèi)容了，經(jīng)過(guò)本文的學(xué)習(xí)后，相信大家對(duì)怎么改進(jìn)Python二分法和牛頓迭代法求算術(shù)平方根這一問(wèn)題有了更深刻的體會(huì)，具體使用情況還需要大家實(shí)踐驗(yàn)證。這里是億速云，小編將為大家推送更多相關(guān)知識(shí)點(diǎn)的文章，歡迎關(guān)注！