matlab黄金分割法求解.docx

文档介绍：------------- 精选文档 -----------------
4 次作业 用 matlab 求解 f(t)=(t^2-1)^2+(t-1)^2+3 [-10,10]最小
植（黄金分割法）。
x 为最优的 t 值，minf 为最优的函数值， f 为函数，t 为函数 f 的变量， a,b 为下单峰区间 [a,b],eps 为精确度，
函数文件：
function [x,minf] = minHJ(f,a,b,eps)
format long;
if nargin == 3
eps = -6;
End
l = a + *(b-a);
u = a + *(b-a);a=-10
k=1;
tol = b-a;
while tol>eps && k<100000
fl= subs(f ,findsym(f),l);%subs 是置换函数
if fl > fu
a = l;
l = u;
u = a + *(b - a);
else
b = u;
可编辑
------------- 精选文档 -----------------
u = l;
l = a + *(b-a);
end
k = k+1;
tol = abs(b - a);
end
if k == 100000
disp(' 找不到最小值！ ');
x= NaN;
minf = NaN;
return;
end
x= (a+b)/2;
minf = subs(f, findsym(f),x);
format short;
窗口 M 文件调用：
syms t