遗传算法工具箱的编码问题(请教)

skyflowing · 发表于 2006-6-20 21:40

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要登录才可以下载或查看，没有账号？我要加入

x

目标函数代码如下: function [sol,eval]=f553(sol,options) m(1)=sol(1); m(2)=sol(2); m(3)=sol(3);
q=[0.01 0.05 0.10 0.18;0.08 0.02 0.15 0.12;0.04 0.05 0.20 0.10]; g1=51-(m(1)+3).^2+m(2).^2+m(3).^2; g2=20*sum(m+exp(-m))-120; g3=20*sum(m.*exp(-m/4))-65;
if((g1>=0)&(g2>=0)&(g3>=0)) multi=1; for i=1:3 summ=0; for j=2:4 summ=summ+q(i,j).^(m(i)+1); end multi=multi*(1-(1-(1-q(i,1)).^(m(i)+1))-summ); end eval=multi; else eval=-500; end 求解部分代码如下: bounds=[1 4;1 7;1 7]; options=[1e-6 1]; initPop=initializega(80,bounds,'f553',[],options); [p,endPop,bestSols,trace]=ga(bounds,'f553',[],initPop,[1e-6 1 1],'maxGenTerm',100,'normGeomSelect',[0.08],['arithXover'],[2],'nonUnifMutation',[2 5 3]) %[p,endPop,bestSols,trace]=ga(bounds,'f553') 为什么options里面的1设为2 运行会出错而并不会得到二进制的初始种群呢.请指教~,非常感谢!

skyflowing · 发表于 2006-6-20 21:44

还有一个疑问如果我要得到一个最优向量,向量中每一位的值为0或者1,应该在程序中怎么表示? 我是把initializega子程序中的options(2)=1部分改为 xZomeLength = numVars+1; %Length of string is numVar + fit pop = zeros(num,xZomeLength); %Allocate the new population pop(:,1:numVars)=(ones(num,1)*rng).*(rand(num,numVars))+... (ones(num,1)*bounds(:,1)'); for i=1:1:numVars for j=1:1:num if pop(j,i)<=0.5 pop(j,i)=0; else pop(j,i)=1; end end end 但是不知道在ga程序中其他部分如交叉,变异等部分是不是不能保证每一位的值为0或者1. 应该怎么改啊[em06][em06]

账号		自动登录	找回密码
密码			我要加入

[综合讨论] 遗传算法工具箱的编码问题(请教)

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。