您好,登錄后才能下訂單哦!
今天就跟大家聊聊有關python二分法怎樣查找函數底值,可能很多人都不太了解,為了讓大家更加了解,小編給大家總結了以下內容,希望大家根據這篇文章可以有所收獲。
假設連續函數f(x)在區間(a,b)上有一個底值m,且在該底值下的函數輸出值為M,即f(m)=M,利用二分法查找該底值:(s為足夠小的數)
令t=(a+b)/2,若|f(t)-M|<=s,則m=t,若|f(t)-M|>s,如果(f(t)-M)和(f(a)-M)同號,a=t,反之b=t,繼續二分法t=(a+b)/2...直到|f(t)-M|<=s,則m=t。
例如:一項一年期投資,每個季度初投入10000元,期滿時收入44163.225,求內部收益率(已設定為0.04)。
收益函數為:
從(0,1)區間中查找該內部收益率,用二分法令t=(0+1)/2=0.5,比較發現|f(t)-44163.225|>10**(-10),并且f(t)-44163.225和f(1)-44163.225同號,則繼續令t=(0+0.5)/2=0.25,,,直到查找出內部收益率為0.03999999445689362,約等于0.04
def f(i): y = 10000 * (1 + i) ** 4 + 10000 * (1 + i) ** 3 + 10000*(1+i)**2+10000*(1+i)**1 return y def division(a,b,M,n): #a,b為自選值范圍,(f(a)-M)*(F(b)-M)<或=0 #M為f(x)輸出值,n為輸出值精度要求小數位數 while True: t = (a + b) / 2 if abs(f(t) - M) <= 10 ** (-n): m = t break if (f(t) - M) / abs(f(t) - M) == (f(a) - M) / abs(f(a) - M): a = t else: b = t return m y=division(0,1,44163.225,10) print(y)
看完上述內容,你們對python二分法怎樣查找函數底值有進一步的了解嗎?如果還想了解更多知識或者相關內容,請關注億速云行業資訊頻道,感謝大家的支持。
免責聲明:本站發布的內容(圖片、視頻和文字)以原創、轉載和分享為主,文章觀點不代表本網站立場,如果涉及侵權請聯系站長郵箱:is@yisu.com進行舉報,并提供相關證據,一經查實,將立刻刪除涉嫌侵權內容。