main.m

用4个种群来优化函数

clc
clear

%该函数为多种群遗传算法的主函数

%参数定义
maxgen=50;                         %进化代数，即迭代次数
sizepop=12;                        %种群规模
%AimFunc=StrAimFunc;                 %目标函数(本算法以目标函数值为适应度)
fselect='roulette';                 %染色体的选择方法，您可以选择：锦标赛法－ 'tournament';轮盘赌法－'roulette'
fcode='float';                       %编码方法，您可以选择：浮点法－'float';grey法则－－'grey';二进制法－'binary' 
pcross=[0.7];                       %交叉概率选择，0和1之间
fcross='float';                   %交叉方法选择，您可以选择： 浮点交叉－'float';单点交叉－'simple';均匀交叉－'uniform'
pmutation=[0.3];                    %变异概率选择，0和1之间
fmutation='float';                 %变异方法选择，您可以选择：浮点法－'float';单点法－'simple';

lenchrom=[1 1];          %每个变量的字串长度，如果是浮点变量，则长度都为1
bound=[-2.048,2.048;-2.048 2.048];

individuals=struct('fitness',zeros(1,sizepop), 'chrom',[]);  %将种群信息定义为一个结构体
avgfitness=[];                      %每一代种群的平均适应度
bestfitness=[];                     %每一代种群的最佳适应度
bestchrom=[];                       %适应度最好的染色体

%用结构体定义四个种群
struct individuals individuals1;
struct individuals individuals2;
struct individuals individuals3;
struct individuals individuals4;
struct individuals individuals5;

%生成三个初始种群
[individuals1,avgfitness1,bestfitness1,bestchrom1,trace1] = PopInit(maxgen,sizepop,fselect,fcode,pcross,fcross,pmutation,fmutation,lenchrom,bound);
[individuals2,avgfitness2,bestfitness2,bestchrom2,trace2] = PopInit(maxgen,sizepop,fselect,fcode,pcross,fcross,pmutation,fmutation,lenchrom,bound);
[individuals3,avgfitness3,bestfitness3,bestchrom3,trace3] = PopInit(maxgen,sizepop,fselect,fcode,pcross,fcross,pmutation,fmutation,lenchrom,bound);
% %从三个种群中提取4个最佳个体组成第四个数组
[individuals4,avgfitness4,bestfitness4,bestchrom4,trace4] = PopMerge31(individuals1,individuals2,individuals3,sizepop,lenchrom,fcode,bound);

N = 300;
% 
for i = 1:N
     i
%    
     
     [individuals1,avgfitness1,bestfitness1,bestchrom1] = iterative(individuals1,maxgen,sizepop,fselect,fcode,pcross,fcross,pmutation,fmutation,lenchrom,bound,bestfitness1,bestchrom1);
     [individuals2,avgfitness2,bestfitness2,bestchrom2] = iterative(individuals2,maxgen,sizepop,fselect,fcode,pcross,fcross,pmutation,fmutation,lenchrom,bound,bestfitness2,bestchrom2);
     [individuals3,avgfitness3,bestfitness3,bestchrom3] = iterative(individuals3,maxgen,sizepop,fselect,fcode,pcross,fcross,pmutation,fmutation,lenchrom,bound,bestfitness3,bestchrom3);
     [individuals4,avgfitness4,bestfitness4,bestchrom4] = iterative(individuals4,maxgen,sizepop,fselect,fcode,pcross,fcross,pmutation,fmutation,lenchrom,bound,bestfitness4,bestchrom4);
    
     trace1=[trace1;avgfitness1 bestfitness1]; 
     trace2=[trace2;avgfitness2 bestfitness2];
     trace3=[trace3;avgfitness3 bestfitness3];
     trace4=[trace4;avgfitness4 bestfitness4];
    
    [individuals5,bestfitness5,bestchrom5] = PopMerge41(individuals1,individuals2,individuals3,individuals4,sizepop,lenchrom,fcode,bound,bestfitness4,bestchrom4);
    individuals4=individuals5;
     
end
    
    [individuals5,avgfitness5,bestfitness5,bestchrom5] = iterative(individuals4,maxgen,sizepop,fselect,fcode,pcross,fcross,pmutation,fmutation,lenchrom,bound,bestfitness5,bestchrom5);


hfig=findobj('Tag','trace1');
% See if it is open
if ishandle(hfig)
    figure(hfig);
else
    hfig=figure('Tag','trace1');
end
figure(hfig);
[r c]=size(trace1);
plot([1:r]',trace1(:,1),'r-',[1:r]',trace1(:,2),'b--');
title(['适应度曲线  ' '终止代数＝' num2str(N)]);
xlabel('进化代数');ylabel('适应度');
legend('平均适应度','最佳适应度');
disp('适应度                   变量');
x=Decode(lenchrom,bound,bestchrom,fcode);
% 窗口显示
disp([bestfitness x]);

%画出适应度曲线
hfig=findobj('Tag','trace2');
% See if it is open
if ishandle(hfig)
    figure(hfig);
else
    hfig=figure('Tag','trace2');
end
figure(hfig);
[r c]=size(trace2);
plot([1:r]',trace2(:,1),'r-',[1:r]',trace2(:,2),'b--');
title(['适应度曲线  ' '终止代数＝' num2str(N)]);
xlabel('进化代数');ylabel('适应度');
legend('平均适应度','最佳适应度');
disp('适应度                   变量');
x=Decode(lenchrom,bound,bestchrom,fcode);
% 窗口显示
disp([bestfitness x]);

%画出适应度曲线
hfig=findobj('Tag','trace3');
% See if it is open
if ishandle(hfig)
    figure(hfig);
else
    hfig=figure('Tag','trace3');
end
figure(hfig);
[r c]=size(trace3);
plot([1:r]',trace3(:,1),'r-',[1:r]',trace3(:,2),'b--');
title(['适应度曲线  ' '终止代数＝' num2str(N)]);
xlabel('进化代数');ylabel('适应度');
legend('平均适应度','最佳适应度');
disp('适应度                   变量');
x=Decode(lenchrom,bound,bestchrom,fcode);
% 窗口显示
disp([bestfitness x]);

%画出适应度曲线
hfig=findobj('Tag','trace4');
% See if it is open
if ishandle(hfig)
    figure(hfig);
else
    hfig=figure('Tag','trace4');
end
figure(hfig);
[r c]=size(trace4);
plot([1:r]',trace4(:,1),'r-',[1:r]',trace4(:,2),'b--');
title(['适应度曲线  ' '终止代数＝' num2str(N)]);
xlabel('进化代数');ylabel('适应度');
legend('平均适应度','最佳适应度');
disp('适应度                   变量');
x=Decode(lenchrom,bound,bestchrom,fcode);
% 窗口显示
disp([bestfitness x]);