Python math庫 ln(x)運算的實現(xiàn)及原理
這個是很有用的一個運算,除了本身可以求自然對數(shù),還是求指數(shù)函數(shù)需要用到的基礎(chǔ)函數(shù)。
實現(xiàn)原理就是泰勒展開,最簡單是在x=1處進(jìn)行泰勒展開:
但該函數(shù)離1越遠(yuǎn)越難收斂,同時大于2時無法收斂,所以需要進(jìn)行換元,然后重新展開:
但是該換元在接近0時或者接近無窮大時收斂困難,處在1到10范圍內(nèi)收斂快且精度高,所以對大于10或小于1的值進(jìn)行分解如下:
ln(55000)=ln(5.5)+4ln10
ln(0.0015)=ln(1.5)-4ln10
ln10為算好的值,可直接由ln_h1(10)得到
Epsilon 為精度控制
輸出的i可以檢測收斂次數(shù)。
Epsilon = 10e-16 ln10 = 2.30258509299404568401 def ln_h(x): ''' ln函數(shù)泰勒換元展開 :param x: 0<x :return:ln(x) ''' def ln_h1(x): s2 = 0.0 delta = x = (x - 1.0) / (x + 1.0) i = 0 while fab_h(delta * 2) / (i * 2 + 1) > Epsilon: s2 += delta / (i * 2 + 1) delta *= x * x i += 1 print(i) return 2 * s2 coef = 0 if x > 10: while x / 10 > 1: coef += 1 x /= 10 return ln_h1(x) + coef*ln10 elif x < 1: while x * 10 < 10: coef += 1 x *= 10 return ln_h1(x) - coef*ln10 else: return ln_h1(x)
以上就是本文的全部內(nèi)容,希望對大家的學(xué)習(xí)有所幫助,也希望大家多多支持腳本之家。
相關(guān)文章
Python selenium實現(xiàn)斷言3種方法解析
這篇文章主要介紹了Python selenium實現(xiàn)斷言3種方法解析,文中通過示例代碼介紹的非常詳細(xì),對大家的學(xué)習(xí)或者工作具有一定的參考學(xué)習(xí)價值,需要的朋友可以參考下2020-09-09python dataframe astype 字段類型轉(zhuǎn)換方法
下面小編就為大家分享一篇python dataframe astype 字段類型轉(zhuǎn)換方法,具有很好的參考價值,希望對大家有所幫助。一起跟隨小編過來看看吧2018-04-04使用Protocol Buffers的C語言拓展提速Python程序的示例
這篇文章主要介紹了使用Protocol Buffers的C語言拓展提速Python程序的示例,使用C拓展Python是Python編程進(jìn)階中的重要技巧,需要的朋友可以參考下2015-04-04python+appium+yaml移動端自動化測試框架實現(xiàn)詳解
這篇文章主要介紹了python+appium+yaml移動端自動化測試框架實現(xiàn)詳解,文中通過示例代碼介紹的非常詳細(xì),對大家的學(xué)習(xí)或者工作具有一定的參考學(xué)習(xí)價值,需要的朋友們下面隨著小編來一起學(xué)習(xí)學(xué)習(xí)吧2020-11-11