绘制目标函数图像
目标函数的图像绘制为(未考虑约束条件)
画法1代码
%有约束的三维函数图像绘制
x1 = -10:.1:10;
x2 = -6:.1:6;
[x1,x2] = meshgrid(x1,x2);
f = exp(x1).*(4*x1.^2+2*x2.^2+4*x1.*x2+2*x2+1);
%f(x1.*x2-x1-x2 > -1.5 | x1.*x2 < -10) = 100;%这里可以设为nan,设为100是为了便于观看图像
mesh(x1,x2,f)
axis([-10 10 -6 6 0 1.6e+7])
画法2代码:用蠢笨的两级for循环
clc
clear all
close all
x=-10:.1:10;
y=-6:.1:6;
for j=1:length(x)
for i=1:length(y)
g1(i,j)=1.5+x(j)*y(i)-x(j)-y(i);
g2(i,j)=-x(j)*y(i);
% if(g1(i,j)>0|g2(i,j)>10)
% f(i,j)=100;
% else
f(i,j)=exp(x(j))*(4*x(j)^2+2*y(i)^2+4*x(j)*y(i)+2*y(i)+1);
%end
end
end
[x,y]=meshgrid(x,y);%x,y每一行都是-20到20的500个数
mesh(x,y,f);
axis([-10 10 -6 6 0 1.6e+7])
画法3代码(就是不明白为什么这里不注释子函数里面描述约束的代码画出的图也是这样,不知道是这个图没考虑约束条件还是考虑了,晕死了)
clc
clear all
close all
x=-10:.1:10;
y=-6:.1:6;
[x,y]=meshgrid(x,y);
z=func1(x,y);
mesh(x,y,z);
axis([-10 10 -6 6 0 1.6e+7])
%子函数,适应度函数,同时也是目标函数
function f=func1(x1,x2);
g1=1.5+x1.*x2-x1-x2;
g2=-x1.*x2;
if(g1>0|g2>10)
f=100;
else
f=exp(x1).*(4*x1.^2+2*x2.^2+4*x1.*x2+2*x2+1);
end
用接力遗传算法求解(用ga工具箱)
命令行输出
Optimization terminated: maximum number of generations exceeded.
Bestx =
-8.5582 1.1071
BestFval =
0.0501
优化代码:
% 主程序:采用遗传算法接力进化
% 即把上次进化结束后得到的最终种群作为本次输入的初始种群
clc
clear all
close all
% 画目标函数图像
x=-10:.1:10;
y=-6:.1:6;
[x,y]=meshgrid(x,y);
z=func1(x,y);
mesh(x,y,z);
axis([-10 10 -6 6 0 1.6e+7])
T=100;%进化代数
optionsOrigin=gaoptimset('Generations',T/2);%返回一个结构体,有35个字段,Generations设为50
[x,fval,reason,output,finnal_pop]=ga(@func2,2,optionsOrigin);
%进行第二次接力进化
options1=gaoptimset('Generations',T/2,'InitialPopulation',finnal_pop,'PlotFcns',@gaplotbestf);
[x,fval,reason,output,finnal_pop]=ga(@func2,2,options1);
Bestx=x
BestFval=fval
axis([1 50 0 60])