十进制遗传算法完整程序.pptVIP

十进制遗传算法完整程序 %--------------------------------------------% % % % 主程序 % %38.8503 % %--------------------------------------------% function [pop]=main %-----------确定遗传算法基本参数--------------% nvars=2; %变量个数 psize=30; %种群大小 maxgen=1000; %最大进化代数 pc=0.9; %交叉概率 pm=0.1; %变异概率 %-------------初始化，产生初始种群-------------% [pop]=initialize(nvars,psize); %---------------计算个体的适应值---------------% [pop]=evalution(pop,nvars,psize); %----------------保持最优个体------------------% [pop]=keepbest(pop,nvars,psize); %---------------进行遗传算法操作---------------% gen=1; while gen=maxgen % epss=0.0001; % while (abs(pop(10,2))=epss) %------------------进行选择操作----------------% [pop]=select(pop,nvars,psize); %------------------进行交叉操作----------------% [pop]=crossover(pop,nvars,psize,pc); %------------------进行变异操作----------------% [pop]=mutate(pop,nvars,psize,pm,gen,maxgen); %----------计算新种群中个体的适应值------------% [pop]=evalution(pop,nvars,psize); %----------------引进精华模型------------------% [pop]=elitist(pop,nvars,psize); %-----------------更新进化代数-----------------% gen=gen+1; end %---------遗传算法优化结束，得到最优解---------% %-----画出k(x)实际值和计算结果的曲线比较图-----% %gen %--------------------------------------------% % % % 初始化 % % % %--------------------------------------------% function [pop]=initialize(nvars,psize) %--pop-种群，nvars-变量个数，psize-种群大小---% %------读入变量范围，注意数据文件的路径-------% load e:\abc.txt lbound=abc(:,1); %取矩阵的第一列为变量的下边届 ubound=abc(:,2); %取矩阵的第二列为变量的上边届 %----------------采用实数编码--------------% for i=1:psize for j=1:nvars