python二分法查找函數(shù)底值

更新時(shí)間：2021年12月30日 16:32:30 作者：一大盤洋芋

大家好，本篇文章主要講的是python二分法查找函數(shù)底值，感興趣的同學(xué)趕快來看一看吧，對你有用的話記得收藏一下，方便下次瀏覽

假設(shè)連續(xù)函數(shù)f(x)在區(qū)間(a,b)上有一個(gè)底值m，且在該底值下的函數(shù)輸出值為M，即f(m)=M，利用二分法查找該底值：(s為足夠小的數(shù)）

令t=(a+b)/2，若|f(t)-M|<=s，則m=t，若|f(t)-M|>s，如果(f(t)-M)和(f(a)-M)同號，a=t，反之b=t，繼續(xù)二分法t=(a+b)/2...直到|f(t)-M|<=s，則m=t。

例如：一項(xiàng)一年期投資，每個(gè)季度初投入10000元，期滿時(shí)收入44163.225，求內(nèi)部收益率(已設(shè)定為0.04)。

收益函數(shù)為：

$f(x)=10000\sum_{k=1}^{4}(1+i)^{k}$

從(0,1)區(qū)間中查找該內(nèi)部收益率，用二分法令t=(0+1)/2=0.5，比較發(fā)現(xiàn)|f(t)-44163.225|>10**(-10)，并且f(t)-44163.225和f(1)-44163.225同號，則繼續(xù)令t=(0+0.5)/2=0.25，，，直到查找出內(nèi)部收益率為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為輸出值精度要求小數(shù)位數(shù)
    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)

到此這篇關(guān)于python二分法查找函數(shù)底值的文章就介紹到這了,更多相關(guān)python函數(shù)底值內(nèi)容請搜索腳本之家以前的文章或繼續(xù)瀏覽下面的相關(guān)文章希望大家以后多多支持腳本之家！

您可能感興趣的文章:

相關(guān)文章

python3.4用循環(huán)往mysql5.7中寫數(shù)據(jù)并輸出的實(shí)現(xiàn)方法
下面小編就為大家?guī)硪黄猵ython3.4用循環(huán)往mysql5.7中寫數(shù)據(jù)并輸出的實(shí)現(xiàn)方法。小編覺得挺不錯的，現(xiàn)在就分享給大家，也給大家做個(gè)參考。一起跟隨小編過來看看吧
2017-06-06
python 實(shí)現(xiàn)查詢Neo4j多節(jié)點(diǎn)的多層關(guān)系
今天小編就為大家分享一篇python 實(shí)現(xiàn)查詢Neo4j多節(jié)點(diǎn)的多層關(guān)系，具有很好的參考價(jià)值，希望對大家有所幫助。一起跟隨小編過來看看吧
2019-12-12
jupyter默認(rèn)工作目錄的更改方法
jupyter notebook是一個(gè)以網(wǎng)頁形式來使用的python編輯器,很多小伙伴在第一次安裝它的時(shí)候選擇的都是默認(rèn)安裝，那么jupyter默認(rèn)工作目錄如何更改，本文就來介紹一下
2023-08-08
Django 對象關(guān)系映射(ORM)源碼詳解
這篇文章主要介紹了Django 對象關(guān)系映射(ORM)源碼詳解,文中通過示例代碼介紹的非常詳細(xì)，對大家的學(xué)習(xí)或者工作具有一定的參考學(xué)習(xí)價(jià)值,需要的朋友可以參考下
2019-08-08
手把手教你jupyter?notebook更換環(huán)境的方法
在日常使用jupyter-notebook時(shí),可能會碰到需要切換不同虛擬環(huán)境的場景,下面這篇文章主要給大家介紹了關(guān)于jupyter?notebook更換環(huán)境的方法,需要的朋友可以參考下
2023-05-05
Python調(diào)用Google?Bard的圖文詳解
Google?Bard?是一種開源數(shù)據(jù)可視化和探索工具,可為?開發(fā)人員?提供支持,本文主要為大家介紹了Python調(diào)用Google?Bard的方法,需要的可以參考下
2023-08-08
Python迭代器模塊itertools使用原理解析
這篇文章主要介紹了Python迭代器模塊itertools使用原理解析,文中通過示例代碼介紹的非常詳細(xì)，對大家的學(xué)習(xí)或者工作具有一定的參考學(xué)習(xí)價(jià)值,需要的朋友可以參考下
2019-12-12
Python實(shí)現(xiàn)人臉識別
這篇文章主要介紹了Python實(shí)現(xiàn)人臉識別，首選抓取多張圖片，從中獲取特征數(shù)據(jù)集和平均特征值然后寫入?csv?文件?-?計(jì)算特征數(shù)據(jù)集的歐式距離作對比，下面一起來看具體得實(shí)現(xiàn)過程吧
2022-01-01
python實(shí)現(xiàn)文件分組復(fù)制到不同目錄的例子
這篇文章主要介紹了python實(shí)現(xiàn)文件按組復(fù)制到不同目錄的例子,需要的朋友可以參考下
2014-06-06
python使用pygame實(shí)現(xiàn)笑臉乒乓球彈珠球游戲
這篇文章主要為大家詳細(xì)介紹了python使用pygame實(shí)現(xiàn)笑臉乒乓球彈珠球游戲，文中示例代碼介紹的非常詳細(xì)，具有一定的參考價(jià)值，感興趣的小伙伴們可以參考一下
2019-11-11