码迷,mamicode.com
首页 > 其他好文 > 详细

matlab练习程序(演化策略ES)

时间:2014-10-26 14:18:26      阅读:348      评论:0      收藏:0      [点我收藏+]

标签:blog   http   io   os   ar   for   sp   on   2014   

% 算法步骤如下:
% 1.设定种群个体数和需要迭代的次数。
% 2.选择父代中的个体按照公式:
% z1=sqrt(-2*ln(u1))*sin(2*pi*u2)*m,
% z2=sqrt(-2*ln(u1))*cos(2*pi*u2)*m
% 这里u1,u2都是随机值,m是控制因子,演化次数越多m,
% m越小,父代通过与z1,z2相加得到后代。
% 3.计算后代的适应性。
% 4.选择后代中最优的适应性作为全局最优适应性。

bubuko.com,布布扣
%% matlab 代码开始
clc;
clear all;
close all;
%% 绘出所求函数图形
[x,y] = meshgrid(-100:100,-100:100);
sigma = 50;
img = (1/(2*pi*sigma^2))*exp(-(x.^2 + y.^2)/(2*sigma^2));
mesh(img);
axis([-100,100,-100,100,2e-5,8e-5]);
hold on
%% 初始化种群
pop_size = 20;%种群规模
pop = struct([]);
iterative = 100;
for i = 1:pop_size
pop(i).x = -100 + 200*rand();
pop(i).y = -100 + 200*rand();
pop(i).fit = compute_fit(pop(i));
end
pop_best = pop(1);

%% 演化开始
p1 = plot3(pop_best.x +100,pop_best.y+100,pop_best.fit,‘k.‘,‘markersize‘,30);
for i = 1:iterative
set(p1, ‘xdata‘,pop_best.x+100,‘ydata‘,pop_best.y+100,‘zdata‘,pop_best.fit)
drawnow
pause(.1)
[pop,pop_best] = select_and_combin(pop,pop_best,pop_size,i,iterative);
end

function [new_pop,pop_best] = select_and_combin(pop,pop_best,pop_size,k,iterative)
mul = (iterative - k) / iterative;
new_pop = pop;
% z1=sqrt(-2*ln(u1))*sin(2*pi*u2)*m,
% z2=sqrt(-2*ln(u1))*cos(2*pi*u2)*m
for i = 1:pop_size
%随机数
u1 = rand();
u2 = rand();
z1 = sqrt(-2*log(u1))*sin(2*pi*u2)*mul;
z2 = sqrt(-2*log(u1))*cos(2*pi*u2)*mul;
%变异
new_pop(i).x = pop(i).x + z1;
new_pop(i).y = pop(i).y + z2;
%计算适应值
new_pop(i).fit = compute_fit(new_pop(i));
if new_pop(i).fit < pop(i).fit
new_pop(i) = pop(i);
end
if pop_best.fit < new_pop(i).fit
pop_best = pop(i);
end
end
end

 

function re = compute_fit(pop)
x = pop.x;
y = pop.y;
sigma = 50;
if x < -100 || x > 100 || y < -100 || y > 100
re = 0;
else
re =(1/(2*pi*sigma^2))*exp(-(x^2 + y ^2)/(2 * sigma^2));
end

 http://www.cnblogs.com/tiandsp/p/3165388.html

matlab练习程序(演化策略ES)

标签:blog   http   io   os   ar   for   sp   on   2014   

原文地址:http://www.cnblogs.com/Kermit-Li/p/4051915.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!