求遗传算法的工具箱的initializega函数

fshmily0711

哪位仁兄有遗传算法的工具箱的initializega函数，能否共享一下？谢谢！！

[ 本帖最后由 eight 于 2007-9-1 16:58 编辑 ]

happy

1. function [pop] = initializega(num, bounds, evalFN,evalOps,options)
   
2. % function [pop]=initializega(populationSize, variableBounds,evalFN,
   
3. %                           evalOps,options)
   
4. %    initializega creates a matrix of random numbers with
   
5. %    a number of rows equal to the populationSize and a number
   
6. %    columns equal to the number of rows in bounds plus 1 for
   
7. %    the f(x) value which is found by applying the evalFN.
   
8. %    This is used by the ga to create the population if it
   
9. %    is not supplied.
   
10. %
    
11. % pop            - the initial, evaluated, random population
    
12. % populatoinSize - the size of the population, i.e. the number to create
    
13. % variableBounds - a matrix which contains the bounds of each variable, i.e.
    
14. %                  [var1_high var1_low; var2_high var2_low; ....]
    
15. % evalFN         - the evaluation fn, usually the name of the .m file for
    
16. %                  evaluation
    
17. % evalOps        - any options to be passed to the eval function defaults []
    
18. % options        - options to the initialize function, ie.
    
19. %                  [type prec] where eps is the epsilon value
    
20. %                  and the second option is 1 for float and 0 for binary,
    
21. %                  prec is the precision of the variables defaults [1e-6 1]
    
22. 
    
23. % Binary and Real-Valued Simulation Evolution for Matlab GAOT V2
    
24. % Copyright (C) 1998 C.R. Houck, J.A. Joines, M.G. Kay
    
25. %
    
26. % C.R. Houck, J.Joines, and M.Kay. A genetic algorithm for function
    
27. % optimization: A Matlab implementation. ACM Transactions on Mathmatical
    
28. % Software, Submitted 1996.
    
29. %
    
30. % This program is free software; you can redistribute it and/or modify
    
31. % it under the terms of the GNU General Public License as published by
    
32. % the Free Software Foundation; either version 1, or (at your option)
    
33. % any later version.
    
34. %
    
35. % This program is distributed in the hope that it will be useful,
    
36. % but WITHOUT ANY WARRANTY; without even the implied warranty of
    
37. % MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
    
38. % GNU General Public License for more details. A copy of the GNU
    
39. % General Public License can be obtained from the
    
40. % Free Software Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA.
    
41. 
    
42. if nargin<5
    
43.   options=[1e-6 1];
    
44. end
    
45. if nargin<4
    
46.   evalOps=[];
    
47. end
    
48. 
    
49. if any(evalFN<48) %Not a .m file
    
50.   if options(2)==1 %Float GA
    
51.     estr=['x=pop(i,1); pop(i,xZomeLength)=', evalFN ';'];  
    
52.   else %Binary GA
    
53.     estr=['x=b2f(pop(i,:),bounds,bits); pop(i,xZomeLength)=', evalFN ';'];
    
54.   end
    
55. else %A .m file
    
56.   if options(2)==1 %Float GA
    
57.     estr=['[ pop(i,:) pop(i,xZomeLength)]=' evalFN '(pop(i,:),[0 evalOps]);'];
    
58.   else %Binary GA
    
59.     estr=['x=b2f(pop(i,:),bounds,bits);[x v]=' evalFN ...
    
60.         '(x,[0 evalOps]); pop(i,:)=[f2b(x,bounds,bits) v];'];  
    
61.     end
    
62. end
    
63. 
    
64. 
    
65. numVars     = size(bounds,1);                 %Number of variables
    
66. rng         = (bounds(:,2)-bounds(:,1))'; %The variable ranges'
    
67. 
    
68. if options(2)==1 %Float GA
    
69.   xZomeLength = numVars+1;                 %Length of string is numVar + fit
    
70.   pop         = zeros(num,xZomeLength);         %Allocate the new population
    
71.   pop(:,1:numVars)=(ones(num,1)*rng).*(rand(num,numVars))+...
    
72.     (ones(num,1)*bounds(:,1)');
    
73. else %Binary GA
    
74.   bits=calcbits(bounds,options(1));
    
75.   xZomeLength = sum(bits)+1;                 %Length of string is numVar + fit
    
76.   pop = round(rand(num,sum(bits)+1));
    
77. end
    
78. 
    
79. for i=1:num
    
80.   eval(estr);
    
81. end

复制代码

mcx8305

sswwss

非常感谢！