用二
4224min ()f t t t t =--[,.]t ∈内的极小值点,要求准
1.
function [t d]=erfenfa(a,b)
k=1; %记录循环次数 while abs(a-b)>0.0005
c=(a+b)/2;
C(k)=c; %存储每次循环中点c 的值
if ff(c)<0
a=c;
end
if ff(c)==0
t1=c;
break ;
end
if ff(c)>0
b=c;
end
k=k+1;
end
t=(a+b)/2; %最终符合要求的值
d=f(t); %最优解
C
k
function y=f(t)
y=t^4-2*t^2-4*t;
function y=ff(t)
y=4*t^3-4*t-4;
运行结果
>> [t d]=erfenfa(1,1.5)
C =
Columns 1 through 9
1.2500 1.3750 1.3125 1.3438 1.3281 1.3203 1.3242 1.3262
1.3252
Column 10
1.3247
k =
11