24小时热门版块排行榜

返回列表

当前主题已经存档。

angry

金虫 (小有名气)

应助: 2 (幼儿园)
金币: 1405.3
帖子: 206
在线: 77.1小时
虫号: 321021
注册: 2007-03-10
专业: 人工智能与知识工程

[交流] 【求助】谁帮忙看下我编写的进化规划算法(EP)，优化结果总是不对

在matlab中编写了基本的进化规划算法，用测试函数试了一下，但总是不正确，自己找不出毛病，哪位高手帮忙看下，谢过了！
%EP program  主程序
%initial population
clc;clear;
pop_size = 100;
dimension = 30;
MaxIteration = 2000;
xmax = 100; xmin = -100;
score = zeros(1,2*pop_size);
p = 10;
for k = 1:dimension
X(:,k)= unifrnd(xmin,xmax,pop_size,1);
end
ETA = 3*ones(pop_size,dimension);
% begin iteration
for iter = 1:MaxIteration
% calculate function value
for i = 1 : pop_size
      Funcvalue(i) = f1(X(i,: ));
end
min(Funcvalue)
%create offspring

for i = 1 : pop_size
%       rnd = rand(1);
      rnd = normrnd(0,1,1,1);
      rnd2 = normrnd(0,1,1,dimension);
      rnd3 = normrnd(0,1,1,dimension);
      exponent = rnd./(sqrt(2*dimension)) + rnd2./sqrt(2*sqrt(dimension));
      NewETA(i, : ) = ETA(i,: ).*exp(exponent);
      NewX(i,: ) = X(i,: ) + rnd3.*NewETA(i,: );
end
ETA = NewETA;
%calculate offspring's fitness
for i = 1 : pop_size
      NewFuncvalue(i) = f1(NewX(i,: ));
end
% put parent's and offspring's fitness into a variable
compFunValue = [Funcvalue, NewFuncvalue];
UnitX = [X;NewX];
compIndex = floor(rand(1,p).*2*pop_size);
for i = 1 : p
      while compIndex(i) == 0
         compIndex(i) =  floor(rand*2*pop_size);
      end
end
score = zeros(1,2*pop_size);
% condult the pair comparison
for i = 1 : 2*pop_size
      for j = 1 : p
         if compFunValue(i) < compFunValue(compIndex(j))
            score(i) = score(i) + 1;
         end
      end;
end
id(1:2*pop_size) = 1 : 2*pop_size;
[Y, ID] = BiDirectionBubble(score,id);  %sort the parent and offspring
for i = 1 : pop_size
      X(i,: ) = UnitX(ID(i),: );
end
end;

% f1.m 测试函数
function y=f1(x)
% This is sphere function
% x is a vector
d=length(x);
y=0;
for k=1:d
y=y+x(k)^2;
end

% sort function
function [Y3,ID]=BiDirectionBubble(x,id)
runtimes=size(x,2);
low=1;up=runtimes;
while up>low
t=low;
for i=low:up-1;
      if x(i)          temp=x(i);
         x(i)=x(i+1);
         x(i+1)=temp;
         idtemp=id(i);
         id(i)=id(i+1);
         id(i+1)=idtemp;
         t=i;
      end;
end;
up=t;
for i=up:-1:low+1;
      if x(i)>x(i-1);
         temp=x(i);
         x(i)=x(i-1);
         x(i-1)=temp;
         idtemp=id(i);
         id(i)=id(i-1);
         id(i-1)=idtemp;
         t=i;
      end;
end;
low=t;
end;
Y3=x;
ID=id;

回复此楼