46改进的遗传算法-函数寻优

Halcom

改进的遗传算法-函数寻优：
百度网盘链接：
链接：http://pan.baidu.com/s/1jIkQs8u

具体链接在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. tic
   
4. popmin = [25]/100;  % Lb
   
5. popmax =[100]/100;  % Ub
   
6. pc = 0.75;  % 0-1之间
   
7. pm = 0.25;   % 0-1之间
   
8. nvar = 1;   % 1个未知数
   
9. itermax = 100; % 迭代次数
   
10. sizepop = 100; % 种群
    
11. % 初始化种群
    
12. for i=1:sizepop
    
13.     pop(i,:) = code(popmin,popmax);    % 实数编码
    
14.     fitness(i) = fun(pop(i,:));        % 目标函数值
    
15. end
    
16. [bestfitness,index] = min(fitness); % 极小值
    
17. zbest = pop(index,:);        % 最佳个体
    
18. fitnessgbest = fitness;      % 个体最优
    
19. gbest = pop;                 % 最优种群
    
20. fitnesszbest = bestfitness;  % 全局最优
    
21. %% 迭代寻优
    
22. for i=1:itermax
    
23.     % 选择算子
    
24.     pop = select(sizepop,pop,fitness,popmin,popmax);
    
25.     % 交叉算子
    
26.     pop = cross(sizepop,pc,pop,popmin,popmax);
    
27.     % 变异
    
28.     pop = mut(sizepop,pm,pop,popmin,popmax);
    
29.    
    
30.     for j=1:sizepop
    
31.         fitness(j) = fun( pop(j,:) );   % 计算适应度值
    
32.         if fitness(j)<fitnesszbest
    
33.             fitnesszbest = fitness(j);
    
34.             zbest = pop(j,:);
    
35.         end
    
36.     end
    
37.     fitness_iter(i) = fitnesszbest;
    
38. end
    
39. toc
    
40. disp(['最优解:',num2str(zbest*100)])
    
41. disp(['最优解目标值：', num2str(fitnesszbest)] )
    
42. figure,
    
43. plot(fitness_iter,'ro-','linewidth',2)

复制代码

选择算子：

1. function pop = select(sizepop,pop,fitness,popmin,popmax)
   
2. [a, b] = sort(fitness,'descend');
   
3. c = floor(sizepop/4);
   
4. % 替换了最差的个体，保留了精英种群
   
5. for i=1:c
   
6.     pop(b(i),:) = popmin + rand.*(popmax-popmin);
   
7. end
   

复制代码

交叉：

1. function pop = cross(sizepop,pc,pop,popmin,popmax)
   
2. for j=1:sizepop*2
   
3.     if rand>pc
   
4.         a = [0,0];
   
5.         while min(a)==0
   
6.             a = floor( rand(1,2)*sizepop );
   
7.         end
   
8.         k = rand;
   
9.         pop1 =  pop(a(1),:) ;
   
10.         pop2 =  pop(a(2),:) ;
    
11.         pop12 = (1-k)*pop1+k*pop2;
    
12.         pop21 = (1-k)*pop2+k*pop1;
    
13.         
    
14.         flag1 = test( pop12 ,popmin,popmax);
    
15.         if flag1==1   % 不可行
    
16.             pop(a(1),:) = pop12;
    
17.         end
    
18.         flag2 = test( pop21 ,popmin,popmax);
    
19.         if flag2==1   % 不可行
    
20.             pop(a(2),:) = pop21;
    
21.         end
    
22.         
    
23.     end
    
24. end

复制代码

变异：

1. function pop = mut(sizepop,pm,pop,popmin,popmax)
   
2. for j = 1:sizepop*2
   
3.     if rand>pm
   
4.         a=0;
   
5.         while a==0
   
6.             a = floor(rand*sizepop);
   
7.         end
   
8.         b = pop(a,:);
   
9.         if rand>0.5
   
10.             b = pop(a,:) + rand.* (popmax-pop(a,:));
    
11.         else
    
12.             b = pop(a,:) + rand.* (popmin-pop(a,:));
    
13.         end
    
14.         
    
15.         flag2 = test( b,popmin,popmax );
    
16.         if flag2==1   % 不可行
    
17.             pop(a,:) = b;
    
18.         end
    
19.         
    
20.     end
    
21. end

复制代码

编码：

1. function pop = code(popmin,popmax)
   
2. flag =1;
   
3. while flag==1
   
4.     pop = popmin + rand.*(popmax-popmin);  % 实数编码
   
5.     % 检验染色体的可行性
   
6.     flag = test(pop,popmin,popmax);
   
7. end
   

复制代码

染色体检验：

1. function flag = test(pop,popmin,popmax)
   
2. 
   
3. % 检验染色体的可行性
   
4. % if pop(1)<popmin(1) || pop(1)>popmax(1) || pop(2)<popmin(2) || pop(2)>popmax(2)
   
5. if pop<popmin || pop>popmax
   
6.     flag = 1;  % 不可行
   
7. else
   
8.     flag = 0;  % 可行
   
9. end

复制代码

目标函数：

1. function y = fun(x1)
   
2. x1 = x1*100;
   
3. 
   
4. Q = 1000;
   
5. fw = 0.6;
   
6. Rg = 2;
   
7. Qr = 38000;
   
8. nL = 0.8;
   
9. nb = 0.6;
   
10. Rd = 1.5;
    
11. T1 = 25;
    
12. P1 = 0.15;
    
13. L = 10000;
    
14. d = 0.2;
    
15. co = 4.2;
    
16. cw = 2;
    
17. % T2 = x1;
    
18. % P2 = x2;
    
19. % FrL = Rg*0.8/nL/Qr *( co*(1-fw) +cw*fw) *Q*(x1 - T1);
    
20. % Fdf = Rd*(x2-P1)*Q*(1-fw)/3.6/nb;
    
21. v = 6910.13-462.16382*x1+12.45893*x1^2-0.16682*x1^3+0.00111*x1^4-0.000002910778*x1^5;
    
22. Hf = 0.0246*(Q*(1-fw))^1.75*(v^1.75)*L/(d^4.75);
    
23. x2 = Hf/100 + 0.3;
    
24. 
    
25. FrL = Rg*0.8/nL/Qr *( co*(1-fw) +cw*fw) *Q*(x1 - T1);
    
26. Fdf = Rd*(x2-P1)*Q*(1-fw)/3.6/nb;
    
27. y = FrL + Fdf;

复制代码

Lshuai

好好学习，天天向上。

somedaytpp

:)好好学习，天天向上

why2018470

好好学习，天天向上