24小时热门版块排行榜    

CyRhmU.jpeg
查看: 2651  |  回复: 9
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

s20090281

银虫 (小有名气)

[求助] 遗传算法的一段程序请教解释

各位高手我在做遗传算法的东西,但是自己刚开始研究,高手给我解决一下下面的程序是什么意思吧:
while generation<=maxgen
   ind=randperm(N-2)+2; % 随机配对交叉
   A=fatherrand(:,ind(1N-2)/2));
   B=fatherrand(:,ind((N-2)/2+1:end));
%     多点交叉
   rnd=rand(num,(N-2)/2);
   ind=rnd   tmp=A(ind);
   A(ind)=B(ind);
   B(ind)=tmp;
% % 两点交叉
%     for kk=1N-2)/2
%         rndtmp=randint(1,1,num)+1;
%         tmp=A(1:rndtmp,kk);
%         A(1:rndtmp,kk)=B(1:rndtmp,kk);
%         B(1:rndtmp,kk)=tmp;
%     end
   fatherrand=[fatherrand(:,1:2),A,B];
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangyuewu08

木虫 (著名写手)

【答案】应助回帖

★ ★
xiegangmai(金币+2): 鼓励讨论交流! 2011-05-15 22:10:42
s20090281(金币+2): 2011-05-16 08:15:05
楼主,你给的程序里面好像有很多的自定义函数,比如fatherrand(  ),A(  ),B(  ),randint(  ),randtemp(  ),遗传算法的流程并不困难,把这些自定义的函数的功能理解清楚就应该没问题了。
我也是刚开始学遗传算法。
往事并不如烟
2楼2011-05-14 08:59:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangyuewu08

木虫 (著名写手)

交叉的意思就是在两个染色体相同的位置上进行互换,这个在高中的生物中讲的很多。

楼主,请教一下细节:
1.图标“”那里有没有运算符?
2.在你的遗传算法中,个体是怎样表示的,染色体如何编码,如何存储,感觉这段代码中的A和B就是染色体的矩阵,可否详细告诉它们的维数,行、列的意义?
3.N代表什么?

有了这些,对理解会有帮助
往事并不如烟
4楼2011-05-16 09:06:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangyuewu08

木虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★
xiegangmai(金币+4): 谢谢应助! 2011-05-18 22:56:13
s20090281(金币+6): 很感谢你,以后还要多多请教呢 2011-05-19 08:28:59
因为楼主和管理员都给了金币,今天又花时间把这代码看了几遍,组要是揣测那些没有交代清楚的条件(在其它子函数中定义的变量)。
我推测fatherrand是num行N列的数组,其中每一列元素代表种群中的一个个体,fatherrand是一个种群,num应该是染色体代码的长度。
先说说这一段:
ind=randperm(N-2)+2; % 随机配对交叉
   A=fatherrand(:,ind(1N-2)/2));
   B=fatherrand(:,ind((N-2)/2+1:end));

ind=randperm(N-2)+2;      作用是将( N - 2 )进行全排,排列之后加上2是为了能够选取最后两列的个体,但是这样的话fatherrand的前2列始终不参加交叉,这个从代码的最后一句 fatherrand=[fatherrand(:,1:2),A,B]可以看出,是将前2列与交叉之后的A和B重新组成新的种群。
A=fatherrand(:,ind(1N-2)/2));       作用是将( N - 2 )个个体的一半放在数组A中。
B=fatherrand(:,ind((N-2)/2+1:end));    将另一半放在B中。

多点交叉我还没有看懂,在这里只说说我对两点交叉的理解:
% % 两点交叉
%     for kk=1: ( N-2)/2
%         rndtmp=randint(1,1,num)+1;
%         tmp=A(1:rndtmp,kk);
%         A(1:rndtmp,kk)=B(1:rndtmp,kk);
%         B(1:rndtmp,kk)=tmp;
%     end

for kk=1: ( N-2)/2;   作用是能够遍历A和B中的( N - 2 ) / 2个个体(从第一列到最后一列)。
rndtmp=randint(1,1,num)+1;      作用是产生一个[ 1, num ]的随机数,注意是闭区间,这个随机数确定了交叉点在染色体中的位置。由于只确定了一个交叉点,所以我觉得这应该叫单点交叉(我看到的一个GA算法的例子中是这样定义单点交叉的)。

%         tmp=A(1:rndtmp,kk);
%         A(1:rndtmp,kk)=B(1:rndtmp,kk);
%         B(1:rndtmp,kk)=tmp;
这3行的作用是将A中某一列的( 1 : rndtmp )个元素与B中相应列的相应元素进行交换,也就是交叉。循环( N - 2 ) / 2次,将A和B中的相应个体均进行交叉。

最后一句的作用是将交叉后的个体与没有参与交叉的前2个个体进行组合,生成新的种群。

有空再帮楼主看看多点交叉的代码,写了这么大一坨字,希望楼主别嫌我太啰嗦
往事并不如烟
7楼2011-05-18 20:13:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 s20090281 的主题更新
信息提示
请填处理意见