- 积分
- 0
- 注册时间
- 2009-1-5
- 仿真币
-
- 最后登录
- 1970-1-1
|
楼主 |
发表于 2009-5-16 16:06:10
|
显示全部楼层
来自 河北石家庄
是自带的ga!!
function[sol,eval]=fit(sol,options)
x1=sol(1);
x2=sol(2);
r1=0.1;
r2=0.8;
g1=x1-2*x2+1;
g2=x1.^2/4-x2.^2+1;
if(g1>=0)&(g2>=0)
eval=(x1-2).^2+(x2-1).^2;
else
eval=(x1-2).^2+(x2-1).^2+r1*g1+r2*g2;
end
eval=-eval;
bounds=ones(2,1)*[-1,1];
initPop=initializega(100,bounds,'fit');
[p endpop bestSols trace]=ga(bounds,'Min');
p
plot(trace(:,1),trace(:,3),'r-');
hold on
plot(trace(:,1),trace(:,2),'b*');
xlable('Generation');ylable('Fitness');
Legend('解的变化','种群平均值的变化');
运行出来
??? Input argument "sol" is undefined. |
|