[原创]一个遗传算法的小例子

qibbxxt

遗传算法matlab是有工具箱的，但是掌握其原理有助于更好的运用，下面是我编的一个遗传算法的小例子，对遗传算法有一定了解的初学者，应该会有帮助
% f(x1,x2)=21.5+x1.*sin(4*pi*x1)+x2.*sin(20*pi*x2);
% -3.0<=x1<=12.1
% 4.1<=x2<=5.8
% 求此函数的极大值

1. 
   
2. % f(x1,x2)=21.5+x1.*sin(4*pi*x1)+x2.*sin(20*pi*x2);
   
3. % -3.0<=x1<=12.1
   
4. % 4.1<=x2<=5.8
   
5. % 求此函数的极小值
   
6. % ga
   
7. function [x1,x2,ff,j]=Myga(iterN,err)
   
8. if nargin==0
   
9.     iterN=1000;err=0.8;
   
10. end
    
11. tic
    
12. x1a=-3.0;x1b=12.1;
    
13. x2a=4.1;x2b=5.8;
    
14. m1=18;m2=15;
    
15. M=100;
    
16. len=33;
    
17. Pc=0.25;Pm=0.01;
    
18. chrom=randint(M,len,[0,1]);
    
19. U=zeros(size(chrom));
    
20. hh= waitbar(0,'Please wait...');
    
21. rand('twister', sum(100*clock));
    
22. for j=1:iterN
    
23.     [evalU,valU]=compute(M,chrom,m1,m2,x1a,x1b,x2a,x2b);
    
24.     if max(evalU)-min(evalU)<err
    
25.         [cc,dd]=max(evalU);
    
26.         x1=valU(dd,1);
    
27.         x2=valU(dd,2);
    
28.         ff=max(evalU);
    
29.         return
    
30.     end
    
31.     %% 选择
    
32.     F=sum(evalU);
    
33.     Pk=evalU/F;
    
34.     Qk=cumsum(Pk);
    
35.     for i=1:M
    
36.         a=find(Qk>rand,1);
    
37.         U(i,:)=chrom(a,:);
    
38.     end
    
39.     %% 交叉
    
40.     randdata=rand(10,1);
    
41.     a=find(randdata<Pc);
    
42.     if length(a)>=2
    
43.         if mod(length(a),2)==1
    
44.             a(end)=[];
    
45.         end
    
46.         a1=a(1:2:end);
    
47.         a2=a(2:2:end);
    
48.         for i=1:length(a1)
    
49.             p=U(a1(i),:);
    
50.             q=U(a2(i),:);
    
51.             aa=randint(1,1,[1,len-1]);
    
52.             temp=p(aa+1:end);
    
53.             p(aa+1:end)=q(aa+1:end);
    
54.             q(aa+1:end)=temp;
    
55.             U(a1(i),:)=p;
    
56.             U(a2(i),:)=q;
    
57.         end
    
58.     end
    
59.     %% 变异
    
60.     aa=rand(M,len);
    
61.     [I,J]=find(aa<Pm);
    
62.     for i=1:length(I)
    
63.         U(I(i),J(i))=not(U(I(i),J(i)));
    
64.     end
    
65.     chrom=U;
    
66.     mcc=j/iterN;
    
67.     waitbar(mcc,hh,['now....',num2str(round(100*mcc)),'%']);
    
68. end
    
69. [cc,dd]=max(evalU);
    
70. x1=valU(dd,1);
    
71. x2=valU(dd,2);
    
72. ff=max(evalU);
    
73. close(hh);
    
74. toc
    
75. function [evalU,valU]=compute(M,chrom,m1,m2,x1a,x1b,x2a,x2b)
    
76. dat1=0:m1-1;
    
77. dat2=0:m2-1;
    
78. valU=zeros(M,2);
    
79. evalU=zeros(M,1);
    
80. for i=1:M
    
81.     valUt1=sum(2.^dat1.*chrom(i,1:m1));
    
82.     valU(i,1)=x1a+valUt1*(x1b-x1a)/(2^m1-1);
    
83.     valUt2=sum(2.^dat2.*chrom(i,m1+1:end));
    
84.     valU(i,2)=x2a+valUt2*(x2b-x2a)/(2^m2-1);
    
85.     evalU(i)=myfun(valU(i,1),valU(i,2));
    
86. end
    
87. function y=myfun(x1,x2)
    
88. y=21.5+x1.*sin(4*pi*x1)+x2.*sin(20*pi*x2);
    

复制代码

yhq20082012012

是挺好的，有多目标物资调度优化方面的吗