使用math中的sqrt函數(shù)
創(chuàng)新互聯(lián)專業(yè)為企業(yè)提供故城網(wǎng)站建設(shè)、故城做網(wǎng)站、故城網(wǎng)站設(shè)計(jì)、故城網(wǎng)站制作等企業(yè)網(wǎng)站建設(shè)、網(wǎng)頁設(shè)計(jì)與制作、故城企業(yè)網(wǎng)站模板建站服務(wù),十多年故城做網(wǎng)站經(jīng)驗(yàn),不只是建網(wǎng)站,更提供有價(jià)值的思路和整體網(wǎng)絡(luò)服務(wù)。
1、示例代碼
import math
amk = math.sqrt(100)
print(amk)
2、示例結(jié)果
10.0
pow(2,0.5)是2開根號
20/pow(2,0.5)是20除以(2開根號)
m-=20/pow(2,0.5) 等效于 m = m - 20/pow(2,0.5)
不得不說這是一個(gè)形勢所趨,
現(xiàn)在高校的很多老師,及學(xué)生做科學(xué)計(jì)算相關(guān)的項(xiàng)目都是python。原因大約有以下幾點(diǎn):
1.
python的語法簡單,這對很少接觸編程的搞學(xué)術(shù)老師的福音。
2.
python相較于其他語言有更豐富的模塊,比如科學(xué)計(jì)算的numpy。
3.
python越來越流行。
1:二分法
求根號5
a:折半:?????? 5/2=2.5
b:平方校驗(yàn):? 2.5*2.5=6.255,并且得到當(dāng)前上限2.5
c:再次向下折半:2.5/2=1.25
d:平方校驗(yàn):1.25*1.25=1.56255,得到當(dāng)前下限1.25
e:再次折半:2.5-(2.5-1.25)/2=1.875
f:平方校驗(yàn):1.875*1.875=3.5156255,得到當(dāng)前下限1.875
每次得到當(dāng)前值和5進(jìn)行比較,并且記下下下限和上限,依次迭代,逐漸逼近平方根:
代碼如下:
import math
from math import sqrt
def sqrt_binary(num):
x=sqrt(num)
y=num/2.0
low=0.0
up=num*1.0
count=1
while abs(y-x)0.00000001:
print count,y
count+=1
if (y*ynum):
up=y
y=low+(y-low)/2
else:
low=y
y=up-(up-y)/2
return y
print(sqrt_binary(5))
print(sqrt(5))
2:牛頓迭代
仔細(xì)思考一下就能發(fā)現(xiàn),我們需要解決的問題可以簡單化理解。
從函數(shù)意義上理解:我們是要求函數(shù)f(x) = x2,使f(x) = num的近似解,即x2 - num = 0的近似解。
從幾何意義上理解:我們是要求拋物線g(x) = x2 - num與x軸交點(diǎn)(g(x) = 0)最接近的點(diǎn)。
我們假設(shè)g(x0)=0,即x0是正解,那么我們要做的就是讓近似解x不斷逼近x0,這是函數(shù)導(dǎo)數(shù)的定義:
從幾何圖形上看,因?yàn)閷?dǎo)數(shù)是切線,通過不斷迭代,導(dǎo)數(shù)與x軸的交點(diǎn)會不斷逼近x0。
可以使用math庫
import matha = 4print math.sqrt(4) # 2
也可以直接利用python的**運(yùn)算符
a = 8a**(1/3) # 開3次方相當(dāng)于1/3次乘方 結(jié)果是2 math中其他常用的數(shù)學(xué)函數(shù):ceil(x) 取頂floor(x) 取底fabs(x) 取絕對值factorial (x) 階乘hypot(x,y) sqrt(x*x+y*y)pow(x,y) x的y次方sqrt(x) 開平方log(x)log10(x)trunc(x) 截?cái)嗳≌麛?shù)部分isnan (x) 判斷是否NaN(not a number)degree (x) 弧度轉(zhuǎn)角度radians(x) 角度轉(zhuǎn)弧度