跪求版主上传多目标DE算法Selectzbest函数与uniqueRep函数

lancelotli

跪求版主上传多目标DE算法Selectzbest函数与uniqueRep函数

lancelotli

Halcom 发表于 2017-8-13 20:57
有啊。自己看参考链接啊

参考链接中这两个函数也没做实现

Halcom

lancelotli 发表于 2017-8-13 21:56
参考链接中这两个函数也没做实现

这两个函数，可要可不要，不要也行！Selectzbest函数是从rep中随机选择一个，uniquerep是对结果进行非重复结果提取。

Halcom

1. % 选择最优解
   
2. function zbest=Selectzbest(rep,beta)
   
3.     % 所有的非支配解节点索引数
   
4.     GridIndex=[rep.GridIndex];
   
5.     % 唯一出现的数值，从小到大排序
   
6.     unique_GridIndex=unique(GridIndex);
   
7.     % 初始化，唯一出现该数值的次数
   
8.     N=zeros(size(unique_GridIndex));
   
9.     for k=1:numel(unique_GridIndex)
   
10.         N(k)=numel(find(GridIndex==unique_GridIndex(k)));
    
11.     end
    
12.     % 选择概率
    
13.     P=exp(-beta*N);
    
14.     P=P/sum(P);
    
15.     % 轮赌法随机选择一个节点
    
16.     rand_Index=RandWheelSelection(P);
    
17.     % 获取对应节点的节点索引数
    
18.     rand_GridIndex=unique_GridIndex(rand_Index);
    
19.     % 在GridIndex中寻找rand_GridIndex所在的索引值
    
20.     index=find(GridIndex==rand_GridIndex);
    
21.     % 在【1到numel(index)】之间随机生成一个整数
    
22.     integer=randi([1 numel(index)]);
    
23.     % zbest
    
24.     zbest=rep(index(integer));
    
25. end

复制代码

Halcom

1. function rep = uniqueRep(rep)
   
2. % global nVar sizepop nObj
   
3. fitnessFA = [];
   
4. fitnessFA = [rep.fitness];  % 适应度值
   
5. % fitnessFA = [reshape( fitnessFA, [nObj, sizepop] )]';  % 适应度值
   
6. flag = [];
   
7. for i=1:size(fitnessFA,2)-1
   
8.     for j=i+1:size(fitnessFA,2)
   
9.         if fitnessFA(1,j)==fitnessFA(1,i) && fitnessFA(2,i)==fitnessFA(2,j)
   
10.             flag = [flag,j];
    
11.         end
    
12.     end
    
13. end
    
14. flag = unique(flag);
    
15. rep(flag) = [];

复制代码

lancelotli

Halcom 发表于 2017-9-24 14:02

感谢余少！！！！！