61-免疫布谷鸟算法的函数寻优分析（MATLAB视频）

Halcom

免疫布谷鸟算法的函数寻优分析（MATLAB视频）
百度网盘链接：
链接：https://pan.baidu.com/s/1alG9Eqe7YF4ojHKayRl51g
具体链接在halcom.cn论坛，联系人QQ：3283892722
该论坛是一个学习交流平台，我会逐一的和大家分享学习。
欢迎大家录制视频，并提交给我，我来设置视频，你可在论坛进行打赏分享。
视频专用播放器：http://halcom.cn/forum.php?mod=viewthread&tid=258&extra=page%3D1
主程序如下：

1. % 基于免疫布谷鸟算法的函数寻优
   
2. clc,clear,close all
   
3. warning off
   
4. tic  % 计时开始
   
5. %% AIA算法参数
   
6. nvar = 2;   % 2个未知数
   
7. % 取值范围
   
8. popmin = [-1, -1];  % Lb
   
9. popmax =[ 1, 1];    % Ub
   
10. itermax = 30; % 迭代次数
    
11. sizepop = 30; % 种群
    
12. Ak = -0.01;   % 聚合适应度常数k
    
13. lamda = 0.9;  % 相似度常数[0.9，1]
    
14. %% CS算法参数
    
15. pa=0.25;        % 发现新巢的概率Discovery rate of alien eggs/solutions
    
16. %% 初始化种群
    
17. for i=1:sizepop
    
18.     pop(i,:) = code(popmin,popmax);    % 实数编码
    
19.     fitness(i) = fun(pop(i,:));        % 目标函数值
    
20. end
    
21. [bestfitness,index] = min(fitness); % 极小值
    
22. zbest = pop(index,:);        % 最佳个体
    
23. fitnesszbest = bestfitness;  % 全局最优
    
24. 
    
25. %% 迭代寻优
    
26. for i=1:itermax
    
27.     % 选择算子
    
28.     pop = select_AIA(sizepop,pop,fitness,popmin,popmax);
    
29.     % 布谷鸟的位置更新
    
30.     new_nest = get_cuckoos(pop,zbest,popmin,popmax); % Levy flights 粒子位置更新
    
31.     [fitnessgbest,bestgbest,pop,fitness]=get_best_nest(pop, fitness, new_nest);
    
32.     % 比较适应度值
    
33.     if fitnessgbest<fitnesszbest
    
34.         fitnesszbest = fitnessgbest;
    
35.         zbest = bestgbest;
    
36.     end
    
37.    
    
38.     % 计算相似度
    
39.     for j=1:sizepop
    
40.         % 交叉算子
    
41.         pop1 = cross_AIA2(sizepop,pop,popmin,popmax);
    
42.         AN = 1./( 1+norm( zbest-pop1, 2 ) );
    
43.         if(AN>lamda)
    
44.             Ci = 1;
    
45.         else
    
46.             Ci = 0;
    
47.         end
    
48.         fitness1 = fun( pop1 ) .* exp( Ak*Ci );   % 计算适应度值
    
49.         if fitness1<fitness(j)
    
50.             fitness(j) = fitness1;
    
51.             pop(j,:) = pop1;
    
52.         end
    
53.         
    
54.         % 变异
    
55.         pop2 = mutation_AIA2(sizepop,pop,popmin,popmax);
    
56.         AN = 1./( 1+norm( zbest-pop2, 2 ) );
    
57.         if(AN>lamda)
    
58.             Ci = 1;
    
59.         else
    
60.             Ci = 0;
    
61.         end
    
62.         fitness2 = fun( pop2 ) .* exp( Ak*Ci );   % 计算适应度值
    
63.         if fitness2<fitness(j)
    
64.             fitness(j) = fitness2;
    
65.             pop(j,:) = pop2;
    
66.         end
    
67.         
    
68.         % 全局比较
    
69.         if fitness(j)<fitnesszbest
    
70.             fitnesszbest = fitness(j);
    
71.             zbest = pop(j,:);
    
72.         end
    
73.         
    
74.         % 寻找新的鸟巢--新的解
    
75.         new_nest=empty_nests(pop, popmin, popmax, pa) ;
    
76.         [fitnessgbest,bestgbest,pop,fitness]=get_best_nest(pop, fitness, new_nest);
    
77.         % 比较适应度值
    
78.         if fitnessgbest<fitnesszbest
    
79.             fitnesszbest = fitnessgbest;
    
80.             zbest = bestgbest;
    
81.         end
    
82.         
    
83.     end
    
84.    
    
85.     fitness_iter(i) = fitnesszbest;
    
86. end
    
87. toc;  % 结束开始
    
88. figure,
    
89. plot(fitness_iter,'ro-','linewidth',2)
    
90. % loglog(fitness_iter,'ro-','linewidth',2)
    
91. grid on
    
92. axis tight
    
93. disp(['最优解   ', num2str(zbest)] )

复制代码

achi

有没有完整的代码？

mzch

有没有完整的代码？

13011185011

代码好像不完整，那个38块的教程是62个算法一起的价格，还是一个算法的价格？

雪花球团子

请问完整的代码在哪里呢？