function r=ff(t,h)%适应度函数
r=t-h;
%%%%%%%%%%%%%%%%%%%%%function k=g(t,h)%约束函数
f1=-5t/12+h-35/12;
f2=-5t/2-h+35/2;
f3=t+5;
f4=5-h;
k=max(0,f1)+max(0,f2)+abs(f3)+abs(f4);
%%%%%%%%%%%%%%%%%%%%%
function xnext=Prey(tt,hh)%觅食行为
Trynumber=5;%可输入（改变）
visual=1;%可输入（改变）
step=0.01;
p=0;
for i=1:Trynumber
xt=[tt hh]+rand*step*visual;
if ff(xt(1),xt(2))>ff(tt,hh)
xnext=[tt hh]+rand*step*(xt-[tt hh]);
p=1;
break;
end
end
if p==0
xnext=[tt hh]+rand*step*visual;
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function q=ASF()
N=20;
visual=1;
step=0.01;
maxnum=50;
x=10*rand(N,2);
%计算当前适应度
m=zeros(1,N);
u=zeros(1,N);
for i=1:Nm(i)=ff(x(i,1),x(i,2));
n(i)=g(x(i,1),x(i,2));
end
%创建公告板 g,root
[ymax,index]=max(m);
xmax=x(index,:);
%开始迭代
num=0;
e=0.01;%可行域
