12万有引力定律算法函数优化--视频分享

Halcom

12万有引力定律算法函数优化--视频分享
链接：http://pan.baidu.com/s/1o8HIehW 密码：0bnz

具体链接在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. % 万有引力定律算法
   
5. % GSA 参数
   
6. G0 = 100;
   
7. a = 20;
   
8. Vmin = -1;
   
9. Vmax = 1;
   
10. maxiter = 200;  % 迭代次数
    
11. sizepop = 20;  % 种群数量
    
12. popmin1 = -1;  popmax1 = 1; % x1
    
13. popmin2 = -1;  popmax2 = 1; % x2
    
14. % 初始化种群
    
15. for i=1:sizepop
    
16.     x1 = popmin1 + (popmax1-popmin1)*rand;
    
17.     x2 = popmin2 + (popmax2-popmin2)*rand;
    
18.     pop(i,1) = x1;
    
19.     pop(i,2) = x2;
    
20.     fitness(i) = fun([x1,x2]);
    
21. end
    
22. % 记录一组最优值
    
23. [bestfitness,bestindex]=min(fitness);
    
24. zbest=pop(bestindex,:);   % 全局最佳
    
25. gbest=pop;                % 个体最佳
    
26. fitnessgbest=fitness;     % 个体最佳适应度值
    
27. fitnesszbest=bestfitness; % 全局最佳适应度值
    
28. V = [0,0];
    
29. % 迭代寻优
    
30. for i=1:maxiter
    
31.    
    
32.     Gt = G0*exp(-a*i/maxiter);        % 万有引力定律常数
    
33.    
    
34.     [mfitness,index] = sort(fitness); % 从小到大排序
    
35.     for j=1:sizepop
    
36.         mt(j) = (fitness(j)-mfitness(end))./(mfitness(1)-mfitness(end)+eps);
    
37.     end
    
38.    
    
39.     for j=1:sizepop
    
40.         Mt = abs(mt(j))./sum(abs(mt))+eps;
    
41.         Rij = norm( pop(j,:)-gbest(j,:) ) +eps;
    
42.         Fij = rand(1,2).*( Gt*Mt*Mt/(Rij+eps) .*(pop(j,:)-gbest(j,:)+eps ) );
    
43.         Va = Fij./Mt;
    
44.         
    
45.         % 速度更新
    
46.         V = rand(1,2).*V - Va;
    
47.         % V--x1
    
48.         if V(1,1)>Vmax
    
49.             V(1,1)=Vmax;
    
50.         end
    
51.         if V(1,1)<Vmin
    
52.             V(1,1)=Vmin;
    
53.         end
    
54.         % V--x2
    
55.         if V(1,2)>Vmax
    
56.             V(1,2)=Vmax;
    
57.         end
    
58.         if V(1,2)<Vmin
    
59.             V(1,2)=Vmin;
    
60.         end
    
61.         
    
62.         % 位置更新
    
63.         pop(j,:) = pop(j,:) + V;
    
64.         % x1  越界限制
    
65.         if pop(j,1)>popmax1
    
66.             pop(j,1)=popmax1;
    
67.         end
    
68.         if pop(j,1)<popmin1
    
69.             pop(j,1)=popmin1;
    
70.         end
    
71.         % x2  越界限制
    
72.         if pop(j,2)>popmax2
    
73.             pop(j,2)=popmax2;
    
74.         end
    
75.         if pop(j,2)<popmin2
    
76.             pop(j,2)=popmin2;
    
77.         end
    
78.         
    
79.         % 适应度更新
    
80.         fitness(j) = fun(pop(j,:));
    
81.         
    
82.         % 比较  个体间比较
    
83.         if fitness(j)<fitnessgbest(j)
    
84.             fitnessgbest(j) = fitness(j);
    
85.             gbest(j,:) = pop(j,:);
    
86.         end
    
87.         if fitness(j)<bestfitness
    
88.             bestfitness = fitness(j);
    
89.             zbest =  pop(j,:);
    
90.         end
    
91.     end
    
92.     fitness_iter(i) = bestfitness;
    
93. end
    
94. disp('最优解')
    
95. disp(zbest)
    
96. fprintf('\n')
    
97. 
    
98. figure('color',[1,1,1])
    
99. plot(fitness_iter,'ro-','linewidth',2)
    
100. 
     
101. figure('color',[1,1,1])
     
102. loglog(fitness_iter,'ro-','linewidth',2)
     
103. axis tight
     
104. grid on

复制代码

lelouch

多谢楼主多谢多谢

huwenping

没注意到还有这个算法，学习。

wal1985

感谢群主，谢谢

mzch

多谢楼主多谢多谢

jiaonong1

谢谢楼主，非常好的讲解

zheng

我想知道密码

HXCS

学习一下代码实现，谢谢分享。

小回收

万有引力算法，想学习一下

lancelotli

谢谢余少，辛苦了。强力顶