matlab bs2rv.m,Matlab智能算法之遗传算法（一）
Matlab智能算法之遗传算法(1)
以往写过的一篇文章了，旧了
1)Sheffield遗传算法工具箱的安装
我共享了下修改过文件名和后缀名的原版工具箱，地址为：http://pan.baidu.com/s/1inVKE
安装方法：
将整个文件夹复制到matlab安装文件夹中的toolbox文件夹
例如：C:\Program Files\MATLAB\R2013b\toolbox文件夹。
然后在Command Window里面输入：
str = ['C:\Program Files\MATLAB\R2013b\toolbox\genetic']
addpath(str)
----------------
mac:
str = ['/Applications/MATLAB.app/toolbox/genetic'];
addpath(str);
----------------
可能有些同学出现过这个问题，
Undefined function or method 'crtbp' for input arguments of type 'double'
或者是：
Cannot find an exact (case-sensitive) match for 'crtbp.m'
The closest match is C:\Program Files\MATLAB\R2012a\toolbox\gatbx\CRTBP.M
To change the file extension, cd to the file's folder, type:
movefile CRTBP.M CRTBP.m_bad; movefile CRTBP.m_bad CRTBP.m
and then cd back.
这是因为新旧Matlab版本对于M文件的文件名要求不尽相同，将其全部改为小写文件名和文件后缀名即可。
(2)应用实例(1)——单变量函数求最值
类似模板的东西，求解此函数最大值：
有如下代码：
%% Do Some Cleaning
clc
clear all;
close all;
%% Set the initial parameters
lb = 1; ub = 2;%x belongs to [1,2]
%% Plot
figure(1);
hold on;
ezplot('sin(10*pi*X)/X',[lb,ub]);
xlabel('x/X')
ylabel('y/Y')
%% Define the parameters of GA
NIND = 40;%size of the group
MAXGEN = 20;%max generations
PRECI = 20;%length of a individual
GGAP = 0.95;%gap
px = 0.7;%the possibility of cross production
pm = 0.1;%the possibility of mutation
trace = zeros(2,MAXGEN);%init value of algorithm寻优函数
FieldD = [PRECI;lb;ub;1;0;1;1];%区域描述器
Chrom = crtbp(NIND,PRECI);%creat random discrete group
%% Optimizations
gen = 0;%counter of generations
X = bs2rv(Chrom,FieldD);%bin to dec
ObjV = sin(10 * pi * X) ./ X;%cal the f(x)
while gen
SelCh= select('sus',Chrom,FitnV,GGAP);%selectSelCh= recombin('xovsp',SelCh,px);%recombineSelCh= mut(SelCh,pm);%mutate
X= bs2rv(SelCh,FieldD);%todec
ObjVSel= sin(10* pi * X)./ X;%cal the next-gen's target f(x)
[Chrom, ObjV]= reins(Chrom,SelCh, 1, 1, ObjV, ObjVSel);%reinsert them to father-gen
X= bs2rv(Chrom,FieldD);
gen= gen+ 1;%counter+=1%get every gen's answers and it's nums, Y stant for best f(x), I for
%nums;
[Y I]= min(ObjV);trace(1,gen)= X(I);trace(2,gen)= Y;end
%% Plot
plot(trace(1,:),trace(2,:),'bo');%plot every gen's answer
grid on;
plot(X,ObjV,'b*');
hold off
%% Plot the evolution
figure(2);
plot(1:MAXGEN,trace(2,:));
grid on
xlabel('count of generations')
ylabel('answer')
title('procedure')
bestY= trace(2,end)bestX= trace(1,end)
运行结果：
(3)应用实例(2)——双变量函数求最值
那么有如下代码：
%% Do Some Cleaning
clc
clear all;
close all;
%% Set the initial parameters
lbx = -2; ubx = 2;%x belongs to [-2,2]
lby = -2; uby = 2;%y belongs to [-2,2]
%% Plot
figure(1);
ezmesh('y*sin(2*pi*x) + x*cos(2*pi*y)',[lbx,ubx,lby,uby],50);
xlabel('x/X')
ylabel('y/Y')
hold on;
%% Define the parameters of GA
NIND = 40;%size of the group
MAXGEN = 20;%max generations
PRECI = 20;%length of a individual
GGAP = 0.95;%gap
px = 0.7;%the possibility of cross production
pm = 0.01;%the possibility of mutation
trace = zeros(3,MAXGEN);%init value of algorithm寻优函数
FieldD = [PRECI PRECI;lbx lby;ubx uby;1 1;0 0;1 1;1 1];%Field discriber
Chrom = crtbp(NIND,PRECI*2);%creat random discrete group
%% Optimizations
gen = 0;%counter of generations
XY = bs2rv(Chrom,FieldD);%bin to dec
X = XY(:,1);
Y = XY(:,2);
ObjV = Y.*sin(2*pi*X) + X.*cos(2*pi*Y);%cal the f(x,y)
while gen
SelCh= select('sus',Chrom,FitnV,GGAP);%selectSelCh= recombin('xovsp',SelCh,px);%recombineSelCh= mut(SelCh,pm);%mutate
XY= bs2rv(SelCh,FieldD);%binto dec
X= XY(:,1);Y= XY(:,2);ObjVSel= Y.*sin(2*pi*X)+ X.*cos(2*pi*Y);%cal the next-gen's target f(x)
[Chrom, ObjV]= reins(Chrom,SelCh, 1, 1, ObjV, ObjVSel);%reinsert them to father-gen
XY= bs2rv(Chrom,FieldD);gen= gen+ 1;%counter+=1%get every gen's answers and it's nums, Y stant for best f(x), I for
%nums;
[Y, I]= max(ObjV);trace(1:2,gen)= XY(I,:);trace(3,gen)= Y;end
%% Plot
plot3(trace(1,:),trace(2,:),trace(3,:),'bo');%plot every gen's answer
grid on;
plot3(XY(:,1),XY(:,2),ObjV,'b*');
hold off
%% Plot the evolution
figure(2);
plot(1:MAXGEN,trace(3,:));
grid on
xlabel('count of generations')
ylabel('answer')
title('procedure')
bestX= trace(1,end)bestY= trace(2,end)bestZ= trace(3,end)
运行结果：
(4)应用实例(3)——遗传算法接力优化(采用系统GAtool工具箱)
优化此函数：
第一个文件(主文件)：
%主程序：本程序采用遗传算法接力进化，
%将上次进化结束后得到的最终种群作为下次输入的初始种群
clc;
close all;
clear all;
%进化的代数
T=100;
optionsOrigin=gaoptimset('Generations',T/2);
[x,fval,reason,output,finnal_pop]=ga(@ff,2,optionsOrigin);
%进行第二次接力进化
options1=gaoptimset('Generations',T/2,'InitialPopulation',finnal_pop,...
'PlotFcns',@gaplotbestf);
[x,fval,reason,output,finnal_pop]=ga(@ff,2,options1);
Bestx=x
BestFval=fval
评价函数文件，与主文件一起保存，名字为ff.m
%子函数：适应度函数同时也是目标函数,函数存储名称为ch14_2f.m
function f=ff(x)
g1=1.5+x(1)*x(2)-x(1)-x(2);
g2=-x(1)*x(2);
if(g1>0||g2>10)
f=100;
else
f=exp(x(1))*(4*x(1)^2+2*x(2)^2+4*x(1)*x(2)+2*x(2)+1);
end
程序结果：