代码:
clc
clear all
close all
% 定义遗传算法参数
lb=0.01,0.01,0.01;ub=1.2,0.5,pi;
NIND=100; %个体数目
MAXGEN=100; %最大遗传代数
PRECI=42; %变量的二进制位数
GGAP=0.95; %代沟
px=0.8; %交叉概率
pm=0.1; %变异概率
trace=zeros(4,MAXGEN); %寻优结果的初始值(原值为2)
FieldD=[PRECI;lb;ub;1;0;1;1]; %区域描述器
Chrom=crtbp(NIND,PRECI); %初始种群
% 优化
gen=0; %代计数器
x=bs2rv(Chrom,FieldD);
y=bs2rv(Chrom,FieldD);
theta=bs2rv(Chrom,FieldD); %计算初始种群的十进制转换
ObjV=-(0.2*(1.3*((x+0.5*cos(theta)).^2-x.^2)+1.4*(log(x+0.5*cos(theta))-log(x))+1.3*((y+0.5*sin(theta)).^2-y.^2)+1.4*(0.5*sin(theta)))); %计算目标函数值
while gen<MAXGEN
FitnV=ranking(ObjV); %分配适应度值
SelCh=select('sus',Chrom,FitnV,GGAP); %选择
SelCh=recombin('xovsp',SelCh,px); %重组
SelCh=mut(SelCh,pm); %变异
x=bs2rv(SelCh,FieldD); %子代个体的十进制转换
y=bs2rv(SelCh,FieldD); %子代个体的十进制转换
theta=bs2rv(SelCh,FieldD); %子代个体的十进制转换
ObjVSel=-(0.2*(1.3*((x+0.5*cos(theta)).^2-x.^2)+1.4*(log(x+0.5*cos(theta))-log(x))+1.3*((y+0.5*sin(theta)).^2-y.^2)+1.4*(0.5*sin(theta)))); %计算子代的目标函数值
[Chrom,ObjV]=reins(Chrom,SelCh,1,1,ObjV,ObjVSel); %重插入子代到父代,得到新种群
x=bs2rv(Chrom,FieldD);
y=bs2rv(Chrom,FieldD);
theta=bs2rv(Chrom,FieldD)
gen=gen+1; %代计数器增加
%获取每代的最优解及其序号,Y为最优解,I为个体的序号
[Q,I]=max(ObjV);
trace(1:2,gen)=y(I); %记下每代的最优值
trace(2:3,gen)=x(I); %记下每代的最优值
trace(3:4,gen)=theta(I); %记下每代的最优值
trace(2,gen)=Q; %记下每代的最优值
end
plot3(trace(1,:),trace(2,:),trace(3,:),'bo'); %画出每代的最优点
grid on;
figure(4);
plot(y,ObjV,'b*');hold on;
figure(5);
plot(x,ObjV,'r*');hold on; %画出最后一代的种群
figure(6);
plot(theta,ObjV,'black*');hold on; %画出最后一代的种群
% 画进化图
figure(2);
plot(1:MAXGEN,trace(2,:));
grid on
xlabel('遗传代数')
ylabel('解的变化')
title('进化过程')
bestQ=trace(2,end);
bestx=trace(1,end);
besty=trace(1,end);
besttheta=trace(1,end);
fprintf(['最优解:\nx=',num2str(bestx),'\ny=',num2str(besty),'\ntheta=',num2str(besttheta),'\nQ=',num2str(bestQ),'\n'])
截图: