工程数学上机实验（一）

一维寻优法（0.618 法）程序设计

1、基本要求

通过一维寻优黄金分割法的程序设计，培养学生计算机的应用能力，并为今后无约束优化方法

的学习和编程，奠定基础；掌握缩小寻优区间的黄金分割法。

2、主要内容

（1）请用 0.618 法求解优化问题：

        在区间[0,1]上的极小点和极小值；

（2）根据 0.618 法算法步骤编写 Matlab 的程序实现 0.618 搜索法；

（3）要求输出内容包括：极小点、极小值、每次迭代的 a、b、al、ak 的值；

（4）按照模板撰写实验报告，要求规范整洁。

3、操作要点

（1）学习 MATLAB 基础编程知识；

（2）认真学习 0.618 法算法步骤，编写 MATLAB 程序 golds 函数；

（3）调用 golds 函数，输入题目参数；

（4）撰写实验报告。

4、主要仪器设备

微机及 Matlab 软件。

代码:
golds.m：

function [xmin, fmin, iter] = golds(fun, a, b, tol)
 
rho = (sqrt(5) - 1) / 2;
 
x1 = a + (1 - rho) * (b - a);
x2 = a + rho * (b - a);
 
f1 = fun(x1);
f2 = fun(x2);
 
iter = [a, b, x1, x2; zeros(floor(log2((b - a) / tol)), 4)];
 
for k = 1:size(iter, 1) - 1
    if f1 < f2
        b = x2;
        x2 = x1;
        x1 = a + (1 - rho) * (b - a);
        f2 = f1;
        f1 = fun(x1);
    else
        a = x1;
        x1 = x2;
        x2 = a + rho * (b - a);
        f1 = f2;
        f2 = fun(x2);
    end
    iter(k + 1, :) = [a, b, x1, x2];
    
    if abs(b - a) < tol
        break;
    end
end
 
if f1 < f2
    xmin = x1;
    fmin = f1;
else
    xmin = x2;
    fmin = f2;
end
 
iter = iter(1:k + 1, :);

main.m

[umin, fmin, iter] = golds(@(x) x^2 - sin(x), 0, pi, 1e-6);
disp(['极小点 = ' num2str(umin)]);
disp(['极小值 = ' num2str(fmin)]);
disp(['每次迭代的 a、b、al、ak 的值:'])
disp(iter);

运行结果: