MATLAB遗传算法之 无放回式随机余数选择算子（remainder stochastic sampling with replacement）最新版
MATLAB遗传算法之 无放回式随机余数选择算子（remainder stochastic sampling with replacement）最新版
function pop = SelectionFunction(pop,fitness_value,ObjV)
[n,m]=size(pop);
num_expectation = zeros(n,1);
for i=1:1:n
num_expectation(i)=n*fitness_value(i)/sum(fitness_value);%计算群体中每个个体在下一代群体中的生存期望数目。
end
fitness_value_new = zeros(n,1);%存放新的适应度举矩阵
num_expectation_int = floor(num_expectation);
%取生存期望数目整数部分作为对应个体在下一代群体中的生存数目。sum(num_expectation_int)可以确定出下一代群体的个数
for i=1:1:n
fitness_value_new(i)=fitness_value(i) - num_expectation_int(i)*sum(num_expectation_int)/n;
%确定下一代中还未确定的的 n - sum(num_expectation_int)个个体的新适应度
end
pop_new = zeros(n,m);%用于存放新的个体。
a=zeros(1,m);%用于标识群体中一行，即一个个体。
s=0;%计数，用于标记新的群体中个数。
%以下方法为按照生存期望选择的个体。
for i=1:1:n
if num_expectation_int(i)~=0
for j=1:1:num_expectation_int(i)
if j==1
a=pop(i,:);
else
a=[a;pop(i,:)];
end
end
pop_new(s+1:s+num_expectation_int(i),:)=a;
s=s+num_expectation_int(i);
end
end
%按照生成期望选择个体结束
%以下为按照轮盘赌选择剩余的个体。
%Tips:当V为全零向量时，该算法无效，将随机选择个体;否则算法将从重要性指标不为0的个体中选择。
%index 选的m个体的位置索引
m = n - s ;%计算剩下的需要轮盘赌的个数
if m~=0
[pop_size,chromo_size] = size(pop);
fitness_table = zeros(pop_size,1);
for i=1:pop_size
if i==1
fitness_table(i) = fitness_table(i) + fitness_value_new(i);
else
fitness_table(i) = fitness_table(i-1) + fitness_value_new(i);
end
end
for i=s+1:pop_size
r = rand * fitness_table(pop_size);
first = 1;
last = pop_size;
mid = round((last+first)/2);
idx = -1;
while (first <= last) && (idx == -1)
if r > fitness_table(mid)
first = mid;
elseif r < fitness_table(mid)
last = mid;
else
idx = mid;
break;
end
mid = round((last+first)/2);
if (last - first) == 1
idx = last;
break;
end
end
pop_new(i,:)=pop(idx,:);
end
end
MATLAB遗传算法之 无放回式随机余数选择算子（remainder stochastic sampling with replacement）效率不高，以过期！
%选择操作：
%无放回式随机余数选择
function selection()
global pop;
global fitness_value;
[n,m]=size(pop);
num_expectation = zeros(n,1);
for i=1:1:n
num_expectation(i)=n*fitness_value(i)/sum(fitness_value);%计算群体中每个个体在下一代群体中的生存期望数目。
end
fitness_value_new = zeros(n,1);%存放新的适应度举矩阵
num_expectation_int = floor(num_expectation);
%取生存期望数目整数部分作为对应个体在下一代群体中的生存数目。sum(num_expectation_int)可以确定出下一代群体的个数
for i=1:1:n
fitness_value_new(i)=fitness_value(i) - num_expectation_int(i)*sum(num_expectation_int)/n;
%确定下一代中还未确定的的 n - sum(num_expectation_int)个个体的新适应度
end
pop_new = zeros(n,m);%用于存放新的个体。
a=zeros(1,m);%用于标识群体中一行，即一个个体。
s=0;%计数，用于标记新的群体中个数。
%以下方法为按照生存期望选择的个体。
for i=1:1:n
if num_expectation_int(i)~=0
for j=1:1:num_expectation_int(i)
if j==1
a=pop(i,:);
else
a=[a;pop(i,:)];
end
end
pop_new(s+1:s+num_expectation_int(i),:)=a;
s=s+num_expectation_int(i);
end
end
%按照生成期望选择个体结束
%以下为按照轮盘赌选择剩余的个体。
%Tips:当V为全零向量时，该算法无效，将随机选择个体;否则算法将从重要性指标不为0的个体中选择。
%index 选的m个体的位置索引
m = n - s ;%计算剩下的需要轮盘赌的个数
if(m ~= 0)%当为0时不进行轮盘赌选择操作
n=size(fitness_value_new);%待挑选的个体数
if max(fitness_value_new)==0&min(fitness_value_new)==0%如果V是全零向量，随机选择
index=ceil(rand(1,m)*n);
else
%将适应度为0的个体驱逐出待选择范围
temindex=find(fitness_value_new~=0);
n=length(temindex);%待挑选的个体数目降低
fitness_value_new=fitness_value_new(temindex);
index=zeros(1,m);
%[V,I]=sort(V,'descend');
fitness_value_new=cumsum(fitness_value_new)/sum(fitness_value_new);
pp=rand(1,m);
for i=1:m,
for j=1:n,
if pp(i)<fitness_value_new(j)
index(i)=j;
break
end
end
end
index=temindex(index);%最终index中的值为原fitness_value_new数组中的位置。
end
for i = 1:m
pop_new = [pop_new;pop(index(i),:)];
end
end
for i=1:m
for j=1:chromo_size
pop(i,j) = pop_new(i,j);
end
end