11人群搜索算法函数优化--视频分享

Halcom

11人群搜索算法函数优化--视频分享
百度网盘链接：链接：http://pan.baidu.com/s/1cMkmPw 密码：4muk

具体链接在halcom.cn论坛，联系人QQ：3283892722
该论坛是一个学习交流平台，我会逐一的和大家分享学习。
欢迎大家录制视频，并提交给我，我来设置视频，你可在论坛进行打赏分享。
视频专用播放器：http://halcom.cn/forum.php?mod=viewthread&tid=258&extra=page%3D1
具体的代码如下：

1. clc,clear,close all
   
2. warning off
   
3. format longG
   
4. %% SOA 参数
   
5. maxiter = 20;  % 迭代次数
   
6. sizepop = 10;  % 种群数量
   
7. Umax=0.9500;   % 最大隶属度值
   
8. Umin=0.0111;   % 最小隶属度值
   
9. popmin1 = -1;  popmax1 = 1; % x1
   
10. popmin2 = -1;  popmax2 = 1; % x2
    
11. fai1 = 0.5;
    
12. fai2 = 0.5;
    
13. w1 = 0.5;
    
14. % w2 = 0.5;
    
15. w2max = 0.7;
    
16. w2min = 0.2;
    
17. %% 初始化种群
    
18. for i=1:sizepop
    
19.     x1 = popmin1 + (popmax1-popmin1)*rand;
    
20.     x2 = popmin2 + (popmax2-popmin2)*rand;
    
21.     pop(i,1) = x1;
    
22.     pop(i,2) = x2;
    
23.     fitness(i) = fun([x1,x2]);
    
24. end
    
25. %% 记录一组最优值
    
26. [bestfitness,bestindex]=min(fitness); % 索引最小值
    
27. zbest=pop(bestindex,:);   % 全局最佳
    
28. gbest=pop;                % 个体最佳
    
29. fitnessgbest=fitness;     % 个体最佳适应度值
    
30. fitnesszbest=bestfitness; % 全局最佳适应度值
    
31. Frontpop = pop;
    
32. %% 迭代寻优
    
33. for i=1:maxiter
    
34.     for j=1:sizepop
    
35.        dego = gbest(j,:)-pop(j,:);
    
36.        dalt = zbest-pop(j,:);
    
37.        dpro = Frontpop(j,:) - pop(j,:);
    
38.             Frontpop = pop;  % 记录上一时刻的最优种群
    
39.        dij = sign( w1*dpro+fai1*dego+fai2*dalt );   % 确定方向
    
40.        % 确定步长
    
41.        [maxfitness,index] = sort(fitness,'descend');
    
42.        u = Umax - index(j)/sizepop *(Umax-Umin);
    
43.        u = u.*rand(1,2);          % 2个变量
    
44.        T = sqrt(-log(u) );
    
45.        xmin = pop(index(1),:);    % 最差的个体
    
46.        xmax = pop(index(end),:);  % 最优的个体
    
47.       
    
48.        w2 = w2max - j/sizepop*(w2max-w2min);
    
49.        delta = w2*abs( xmin-xmax );
    
50.       
    
51.        alpha = delta.*T;
    
52.       
    
53.        % 步长进行限制
    
54.        for k=1:length(alpha)
    
55.            if alpha(k)>0.6
    
56.                alpha(k)=0.6;
    
57.            elseif alpha(k)<-0.6
    
58.                alpha(k) = -0.6;
    
59.            end
    
60.        end
    
61.       
    
62.        % 位置更新
    
63.        pop(j,:) = pop(j,:) + dij.*alpha;        
    
64. 
    
65.         % x1  越界限制
    
66.         if pop(j,1)>popmax1
    
67.             pop(j,1)=popmax1;
    
68.         end
    
69.         if pop(j,1)<popmin1
    
70.             pop(j,1)=popmin1;
    
71.         end
    
72.         % x2  越界限制
    
73.         if pop(j,2)>popmax2
    
74.             pop(j,2)=popmax2;
    
75.         end
    
76.         if pop(j,2)<popmin2
    
77.             pop(j,2)=popmin2;
    
78.         end
    
79.         
    
80.         % 适应度更新
    
81.         fitness(j) = fun(pop(j,:));
    
82.         
    
83.         % 比较  个体间比较
    
84.         if fitness(j)<fitnessgbest(j)
    
85.             fitnessgbest(j) = fitness(j);
    
86.             gbest(j,:) = pop(j,:);
    
87.         end
    
88.         if fitness(j)<bestfitness
    
89.             bestfitness = fitness(j);
    
90.             zbest =  pop(j,:);
    
91.         end
    
92.         
    
93.     end
    
94.     fitness_iter(i) = bestfitness;
    
95.    
    
96. end
    
97. disp('最优解')
    
98. disp(zbest)
    
99. fprintf('\n')
    
100. 
     
101. figure('color',[1,1,1])
     
102. plot(fitness_iter,'ro-','linewidth',2)
     
103. 
     
104. figure('color',[1,1,1])
     
105. loglog(fitness_iter,'ro-','linewidth',2)
     
106. axis tight
     
107. grid on

复制代码

huwenping

人群搜索算法，以前未涉及。

Flora

好好好好好好，支持支持

醉生梦死2017

看看看看看看

kkarumi

感谢分享！

wal1985

非常好的代码，谢谢

东晶石材

人群搜索算法函数优化--视频分享 [修改]
高级模式

王宇666

我看看,人群搜索算法。

wht1129838255

谢谢分享！！！！！！

kidcad

多谢多谢学习一下