python实现黄金分割法的示例代码
作者：jtwty 发布时间：2022-08-27 19:41:17
标签：python,黄金分割法
一.问题
使用黄金分割法来计算
二.代码
#黄金分割法python求解PPT上第一个例题
#因为函数要求解最大值而这个方法一般求解最小值所以把函数取负
import numpy as np
import matplotlib.pyplot as plt
rate = 0.618034
def f(x):
#求解体积函数公式，乘1.0将结果变为浮点数
return -1.0*x*(350-2*x)*(260-2*x)
def tarceback(f,a0,b0,accuracy):
a = a0
b = b0
x2 = a+rate*(b-a)
x1 = b-rate*(b-a)
f1 = f(x1)
f2 = f(x2)
print(x1,x2)
arr = search(f,a,b,x1,x2,f1,f2,accuracy)
printFunc(f,a,b,arr[0],arr[1])
def search(f,a,b,x1,x2,f1,f2,accuracy):
if f1<=f2:
if x2-a<accuracy:
print(x1,f1)
return (x1,f1)
else:
b = x2
x2 = x1
f2 = f1
x1 = a+b-x2
f1 = f(x1)
print(x1,x2)
return search(f,a,b,x1,x2,f1,f2,accuracy)
else:
if b-x1<accuracy:
print(x2,f2)
return (x2,f2)
else:
a = x1
x1 = x2
f1 = f2
x2 = a+b-x1
f2 = f(x2)
print(x1,x2)
return search(f,a,b,x1,x2,f1,f2,accuracy)
def printFunc(f,a,b,x,y):
t = np.arange(a,b,0.01)
s = f(t)
plt.plot(t,s)
plt.plot([x],[y],'ro')
plt.plot([x,x],[y,0],'k--')
plt.plot([0,x],[y,y],'k--')
# plt.annotate(r'$(x,y)$',xy=(x,y))
plt.show()
tarceback(f,0,130,0.05)
三.结果
来源：https://blog.csdn.net/weixin_46308081/article/details/116164984
0
投稿猜你喜欢
- 一、前言mysql5.7版本开始支持JSON类型字段，本文详细介绍json_extract函数如何获取mysql中的JSON类型数据json
- <?php$dsn = "mysql:host=localhost;dbname=group";$db = new
这是一个特别简单的用原生js实现的一个小鸟游戏，比较简单，适合新手练习。html结构<div id="game"&- 前言在awesomes上寻找移动端框架的时候意外发现了vux的页面切换效果，后面由于其他考虑没有选用vuex但是这个切换效果确实感觉很有新意
- 如下所示：安装kafka支持库pip install kafka-pythonfrom kafka import KafkaProducer
- 本文实例为大家分享了Python模拟用户登录验证的具体代码，供大家参考，具体内容如下1.功能简介此程序模拟用户登录验证的过程，实现用户名输入
- Ajax 是异步的JavaScript和XML的简称，是一种更新页面某部分的机制。它赋予了你从服务器获取数据后，更新页面某部分的权力，从而避
- Shapely是一个Python库，用于操作和分析笛卡尔坐标系中的几何对象。引入包from shapely.geometry import
- 这篇文章是入门级别的应用Python + Selenium进行自动化测试，包括环境搭建及简单的实例。基本思想是用Firefox Seleni
- 例如下面这段代码 { var temp = "12"; } alert(temp); //输出 12 如果按照通常的编程
- 本文实例讲述了python单元测试unittest用法。分享给大家供大家参考。具体分析如下：单元测试作为任何语言的开发者都应该是必要的，因为
- 很多人会把Primary Key和聚集索引搞混起来，或者认为这是同一个东西。这个概念是非常错误的。 主键是一个约束（constraint），
- Python 代码库之Tuple如何append元素tuple不像array给我们提供了append函数，我们可以通过下面的方式添加t=[1
- 举例吧，这样更容易看的懂。例1：from django.core.exceptions import ObjectDoesNotExistt
- 1.sys模块sys模块的常见函数列表：sys.argv: 实现从程序外部向程序传递参数。sys.exit([arg]): 程序中间的退出，
- 在使用ros的时候经常会用到rosbag来录制或者回放算法，是个非常有用的工具。rosbag 命令列表命令作用record录制一个包，并且指
- 引子之前clubot使用的pyxmpp2的默认mainloop也就是一个poll的主循环,但是clubot上线后资源占用非常厉害,使用str
- 1. Python模块和包：一切从基础开始Python模块是一个Python文件，包含一些相关的函数、类或变量的定义，可以通过 i
- 在前面的文章中，我分别介绍了基于列表元素柱状图和基于表格元素柱状图的实现方法，虽然方法比较简单，但是它却包含了基本的实现原理。在了解了前面两
- 今天我们会讲到一个[装饰器]注记：链接“装饰器”指Python3教程中的装饰器教程。可以在这里快速了解什么是装饰器。@functools.l
