Matlab学习手记——黄金分割法求极小值

最新推荐文章于 2023-11-14 18:51:50 发布

程序猿老王。

最新推荐文章于 2023-11-14 18:51:50 发布

阅读量1.1w

点赞数 5

分类专栏：算法系列 Matlab学习手记文章标签：黄金分割法函数极小值最优化

本文链接：https://blog.csdn.net/u012366767/article/details/81631098

版权

Matlab学习手记同时被 2 个专栏收录

80 篇文章 75 订阅

订阅专栏

算法系列

29 篇文章 9 订阅

订阅专栏

功能：黄金分割法求函数极小值点。

源码

function [xo, fo] = MinValue_Gold(func, a, b, eps, MaxIter)
%{
函数功能：黄金分割法确定最小值点，func在[a, b]上为单峰；
输入：
  func：目标函数句柄；
  a：左端点值；
  b：右端点值；
  eps：精度，默认精度1e-6；
  MaxIter：迭代次数；
输出：
  xo：最小值点；
  fo：最小值；
示例：
func = @(x) x+6*sin(4*x)+9*cos(5*x);
[x, f] = MinValue_Gold(func, 0, 9, 1e-6, 1000)
fplot(func, [0, 9]);
hold on 
plot(x, f, '*')
hold off
%}
% =======================================================================
if nargin < 3
   error('输入参数不足！');
end
if nargin == 3
   eps = 1e-6;
   MaxIter = 1000;
end
if nargin == 4
   MaxIter = 1000;
end
% 初始情况
a1 = a + 0.382*(b - a);
a2 = a + 0.618*(b - a);
f1 = func(a1);
f2 = func(a2);
ite = 0;
while abs(b - a) >= eps && ite < MaxIter
   if f1 < f2
      b = a2;
      a2 = a1;
      f2 = f1;
      a1 = a + 0.382*(b - a);
      f1 = func(a1);
   else
      a = a1;
      a1 = a2;
      f1 = f2;
      a2 = a + 0.618*(b - a);
      f2 = func(a2);
   end
   ite = ite + 1;
   if ite == MaxIter
      disp('没有找到最小值！');
      xo = NaN;
      fo = NaN;
      return;
   end
end;
xo = 0.5*(a + b);
fo = func(xo);