黄金分割算法（算法分析+python代码解释）_云模型黄金分割率法计算(1)

2401_84562529

于 2024-05-04 00:29:06 发布

阅读量449

点赞数 5

分类专栏：程序员文章标签：算法 python java

本文链接：https://blog.csdn.net/2401_84562529/article/details/138428720

版权

程序员专栏收录该内容

191 篇文章 1 订阅

订阅专栏

λk = $a_{k}$ + （1-α）（ $b_{k}$ - $a_{k}$ ）

若令α^2 = 1 - α，则有μk+1 = λk，*因此μ*k+1不必重新计算

$\Rightarrow$ α = （-1+ $\sqrt{5}$ ）/2 （α $>$ 0）

另当f（λk） $>$ f（μk）时，同理，最终可知0.618法计算试探点的公式：

λk = $a_{k}$ + 0.382（ $b_{k}$ - $a_{k}$ ）

μk = $a_{k}$ + 0.618（ $b_{k}$ - $a_{k}$ ）

如下图所示，为黄金分割法的示意图：

黄金分割法示意图

算法步骤：

步骤1：选取a $<$ b,及精度ε> 0，

步骤2：计算 $x_{1}$ :=a + 0.382(b-a), $x_{2}$ :=a + 0.618(b-a)。

步骤3：计算 $f_{1}$ ：= $f(x_{1})$ , $f_{2}$ ：= $f(x_{2})$ 。

步骤4：如果 $f_{1}$ $>$ $f_{2}$ ，则令a:= $x_{1}$ ,若|b - a| $<$ ε，则转步骤5，否则令 $f_{1}$ ：= $f_{2}$ ， $x_{1}$ ：= $x_{2}$ ， $x_{2}$ ：=a + 0.618（b-a）， $f_{2}$ ：= $f(x_{2})$ ，转步骤4。否则令b:= $x_{2}$ ,若|b - a| $<$ ε，则转步骤5，否则令 $f_{2}$ ：= $f_{1}$ ， $x_{2}$ ：= $x_{1}$ ， $x_{1}$ ：=a + 0.382（b-a）， $f_{1}$ ：= $f(x_{1})$ ，转步骤4。

步骤5：停止，输出 $\bar{x}$ = （a+b)/2 。

缺点：收敛速度慢，无理数取小数点后三位，有一定误差。优点：不要求函数可微，且每次迭代只需计算一个函数值，计算量小。

例题分析：

例题：试用0.618法求目标函数 $f(x)=x^{3}-2x+1$ 的最优解。给定初始区间[0，2],收敛精度ε=0.002.

解：第一次区间缩短计算过程：计算两点及对应函数值：

a =0，b = 2,

$x_{1}$ =a + 0.382(b-a) = 0.764, $x_{2}$ =a + 0.618(b-a) = 1.236

$f(x_{1})$ = -0.0821， $f(x_{2})$ = 0.4162

可见： $f(x_{1})$ $<$ $f(x_{2})$ ，进行置换，b = $x_{2}$ = 1.236，[a，b] = [0，1.236],|b-a| = 1.236 $>$ ε.

第二次区间缩短计算过程：计算两点及对应函数值：

a =0，b = 1.236,

$x_{1}$ =a + 0.382(b-a) = 0.472, $x_{2}$ =a + 0.618(b-a) = 0.764

$f(x_{1})$ = 0.1612， $f(x_{2})$ =-0.0821

可见： $f(x_{1})$ $>$ $f(x_{2})$ ，进行置换，a = $x_{1}$ = 0.472，[a，b] = [0.472，1.236],|b-a| = 0.764 $>$ ε.

接着进行第三次，第四次迭代，直至|b-a| $<$ ε.结束迭代。

算法代码：

'''
黄金分割算法（0.618法）
2023.10.9
'''

def fun(x):
    return x**3-2*x+1
	
def run(a,b,ε):
	count = 0
	while True:
	    if abs(b-a) > ε:
	        x1 = a + 0.382 * (b - a)
	        x2 = a + 0.618 * (b - a)
	        count = count + 1
	        print("第{}次迭代".format(count))
	        if fun(x1) > fun(x2):
	            #print("置换x1")
	            a = x1
	            print("得到的新的区间为：[%f,%f]" %(a,b))
	            #print(abs(b-a))
	            continue
	        else:
	            #print("置换x2")
	            b = x2
	            print("得到的新的区间为：[%f,%f]" %(a,b))
	            #print(abs(b - a))
	            continue
	    elif abs(b-a) <= ε:
	        print("x的最优解是","{:.3}".format((b + a)/2))
	        print("函数在该精度上的最小值为",fun((b + a)/2))
	        break
	return a,b
run(a=0,b=2,ε=0.002)

用该代码所运算的结果为:

第1次迭代
学好 Python 不论是就业还是做副业赚钱都不错，但要学会 Python 还是要有一个学习规划。最后大家分享一份全套的 Python 学习资料，给那些想学习 Python 的小伙伴们一点帮助！



### 一、Python所有方向的学习路线



Python所有方向路线就是把Python常用的技术点做整理，形成各个领域的知识点汇总，它的用处就在于，你可以按照上面的知识点去找对应的学习资源，保证自己学得较为全面。



![](https://img-blog.csdnimg.cn/img_convert/9f49b566129f47b8a67243c1008edf79.png)



### 二、学习软件

工欲善其事必先利其器。学习Python常用的开发软件都在这里了，给大家节省了很多时间。



![](https://img-blog.csdnimg.cn/img_convert/8c4513c1a906b72cbf93031e6781512b.png)



### 三、全套PDF电子书



书籍的好处就在于权威和体系健全，刚开始学习的时候你可以只看视频或者听某个人讲课，但等你学完之后，你觉得你掌握了，这时候建议还是得去看一下书籍，看权威技术书籍也是每个程序员必经之路。

![](https://img-blog.csdnimg.cn/img_convert/46506ae54be168b93cf63939786134ca.png)



### 四、入门学习视频

我们在看视频学习的时候，不能光动眼动脑不动手，比较科学的学习方法是在理解之后运用它们，这时候练手项目就很适合了。



![](https://img-blog.csdnimg.cn/afc935d834c5452090670f48eda180e0.png?x-oss-process=image/watermark,type_d3F5LXplbmhlaQ,shadow_50,text_Q1NETiBA56iL5bqP5aqb56eD56eD,size_20,color_FFFFFF,t_70,g_se,x_16#pic_center)



### 五、实战案例



光学理论是没用的，要学会跟着一起敲，要动手实操，才能将自己的所学运用到实际当中去，这时候可以搞点实战案例来学习。



![](https://img-blog.csdnimg.cn/img_convert/252731a671c1fb70aad5355a2c5eeff0.png)



### 六、面试资料

我们学习Python必然是为了找到高薪的工作，下面这些面试题是来自阿里、腾讯、字节等一线互联网大厂最新的面试资料，并且有阿里大佬给出了权威的解答，刷完这一套面试资料相信大家都能找到满意的工作。



![](https://img-blog.csdnimg.cn/img_convert/6c361282296f86381401c05e862fe4e9.png)  

![](https://img-blog.csdnimg.cn/img_convert/d2d978bb523c810abca3abe69e09bc1a.png)




**网上学习资料一大堆，但如果学到的知识不成体系，遇到问题时只是浅尝辄止，不再深入研究，那么很难做到真正的技术提升。**

**[需要这份系统化学习资料的朋友，可以戳这里无偿获取](https://bbs.csdn.net/topics/618317507)**

**一个人可以走的很快，但一群人才能走的更远！不论你是正从事IT行业的老鸟或是对IT行业感兴趣的新人，都欢迎加入我们的的圈子（技术交流、学习资源、职场吐槽、大厂内推、面试辅导），让我们一起学习成长！**