python二分法查找函数底值
更新时间：2021年12月30日 16:32:30 作者：一大盘洋芋
大家好，本篇文章主要讲的是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二分法查找函数底值的文章就介绍到这了,更多相关python函数底值内容请搜索脚本之家以前的文章或继续浏览下面的相关文章希望大家以后多多支持脚本之家！
相关文章
python3.4用循环往mysql5.7中写数据并输出的实现方法
下面小编就为大家带来一篇python3.4用循环往mysql5.7中写数据并输出的实现方法。小编觉得挺不错的，现在就分享给大家，也给大家做个参考。一起跟随小编过来看看吧2017-06-06
今天小编就为大家分享一篇python 实现查询Neo4j多节点的多层关系，具有很好的参考价值，希望对大家有所帮助。一起跟随小编过来看看吧2019-12-12
jupyter notebook是一个以网页形式来使用的python编辑器,很多小伙伴在第一次安装它的时候选择的都是默认安装，那么jupyter默认工作目录如何更改，本文就来介绍一下2023-08-08
这篇文章主要介绍了Django 对象关系映射(ORM)源码详解,文中通过示例代码介绍的非常详细，对大家的学习或者工作具有一定的参考学习价值,需要的朋友可以参考下2019-08-08
在日常使用jupyter-notebook时,可能会碰到需要切换不同虚拟环境的场景,下面这篇文章主要给大家介绍了关于jupyter notebook更换环境的方法,需要的朋友可以参考下2023-05-05
Google Bard 是一种开源数据可视化和探索工具,可为 开发人员 提供支持,本文主要为大家介绍了Python调用Google Bard的方法,需要的可以参考下2023-08-08
这篇文章主要介绍了Python迭代器模块itertools使用原理解析,文中通过示例代码介绍的非常详细，对大家的学习或者工作具有一定的参考学习价值,需要的朋友可以参考下2019-12-12
这篇文章主要介绍了Python实现人脸识别，首选抓取多张图片，从中获取特征数据集和平均特征值然后写入 csv 文件 - 计算特征数据集的欧式距离作对比，下面一起来看具体得实现过程吧2022-01-01
这篇文章主要介绍了python实现文件按组复制到不同目录的例子,需要的朋友可以参考下2014-06-06
这篇文章主要为大家详细介绍了python使用pygame实现笑脸乒乓球弹珠球游戏，文中示例代码介绍的非常详细，具有一定的参考价值，感兴趣的小伙伴们可以参考一下2019-11-11
最新评论