实 验 报 告
实验题目: 0.618法的MATLAB实现
学生姓名:
学 号:
实验时间:
实验名称: 0.618法求解单峰函数极小点
实验目的及要求:
1. 了解并熟悉0.618法的方法原理, 以及它的MATLAB实现.
2. 运用0.618法解单峰函数的极小点.
实验内容:
0.618法方法原理:
定理: 设是区间上的单峰函数, , 且. 如果, 则对每一个, 有; 如果,
则对每一个, 有.
根据上述定理, 只需选择两个试探点, 就可将包含极小点的区间缩短. 事实上, 必有
如果, 则;
如果, 则.
0.618 法的基本思想是, 根据上述定理, 通过取试探点使包含极小点的区间(不确定区间)不断缩短, 当区间长度小到一定程度时, 区间上各点的函数值均接近极小值, 因此任意一点都可作为极小点的近似.
0.618 法计算试探点的公式:
0.618法的算法步骤:
①置初始区间及精度要求, 计算试探点和, 计算函数值和. 计算公式是
令.
②若, 则停止计算. 否则, 当时, 转步骤③; 当时, 转步骤④.
③置, , ,, 计算函数值, 转步骤⑤.
④置, , ,, 计算函数值, 转步骤⑤.
⑤置, 返回步骤②.
实验流程图及其MATLAB实现:
1. 流程图:
代码及数值算例:
程序源代码:
function [x,k]=GSe(f,a,b,delta)
% 0.618法求解单峰函数极小点
f=inline(f);
N=10000;
for k=1:N
m=a+0.382*(b-a);
n=a+0.618*(b-a);
if f(m)>f(n)
a=m;
m=n;
else b=n;
n=m;
end
if abs(b-a)
x=0.5*(b+a);break;
end
end
数值算例:
,
初始区间, 精度.
(i) 键入命令:
>> [x,k]=GSe('2*x^2-x-1',-1,1,0.16)
(ii) 运行结果:
x =
0.2229
k =
6
总结:
0.618法(黄金分割法)适用于单峰函数, 故应先确定目标函数的单峰区间, 方可进行迭代计算. 但单峰区间不是很明显就能确定, 故可用进退法寻找并确定单峰区间.
开始
置初始区间及精度要求
计算试探点和
计算函数值和
Y
N
Y
N
置, ,
计算函数值
置, ,
计算函数值
置