24小时热门版块排行榜

返回列表

s20090281

银虫 (小有名气)

应助: 0 (幼儿园)
金币: 460.9
散金: 177
红花: 1
帖子: 235
在线: 31.2小时
虫号: 1185039
注册: 2011-01-07
性别: MM
专业: 人工智能与知识工程

[求助] 遗传算法求助指点需要改进的地方

function [out]=my_function(input)
clear all;
clc;
size=50;
codel=4;
G=100;
T1=clock;
s=sprintf('程序正在运行中,请稍等......');
disp(s);
popu=2*rand(size,codel)-1;
for kg=1:G
time(kg)=kg;
for i=1:size
      xi=popu(i,;
      x1=xi(1);
      x2=xi(2);
      x3=xi(3);
      x4=xi(4);
      a=load('C:\testresult\控制.txt');
      b=load('C:\testresult\反馈.txt');
      c=0;
      for k=3:360
         f=[b(k)-x1*b(k-1)-x2*b(k-2)-x3*a(k-1)-x4*a(k-2)]^2;
         c=c+f;
         Msum(i,k)=c;
         Mfit(i,k)=1/c;
      end
      fitness=Mfit(:,k);
end
[Oderfi,Indexfi]=sort(fitness);
   Bestfi=Oderfi(size);
   BestS=popu(Indexfi(size),;
   bfi(kg)=Bestfi
   fi_sum=sum(fitness);
   fi_size=(Oderfi/fi_sum)*size;
   fi_s=floor(fi_size);
   r=size-sum(fi_s);
   for i=size:-1:size-r+1
      fi_s(Indexfi(i))=fi_s(Indexfi(i))+1; %Adding rest to equal Size
   end
   k=1;
   for i=size:-1:1
      for j=1:fi_s(i)
         newpopu(k,=popu(Indexfi(i),; %Selecting and Reproduce
         k=k+1;                         %k is used to reproduce
      end
   end
   popu=newpopu;
   pc=0.9;
   for i=1:2:size-1
      temp=rand;
      if pc>temp
         alfa=rand;
         newpopu(i,=alfa*popu(i+1,+(1-alfa)*popu(i,;
         newpopu(i+1,=alfa*popu(i,+(1-alfa)*popu(i+1,;
      end
   end
   A(size,

=BestS;
   popu=newpopu;

   pm=0.1;
   minx=[-1 -1 -1 -1];
   maxx=[1 1 1 1];
   Dif=maxx-minx;
   for i=1:size
      for j=1:codel
         if pm>rand(i,j)
               mean=(minx+maxx)/2
               newpopu(i,j)=mean(j)+Dif(j)*(rand-0.5);
         end
      end
   end
B(size,

=BestS;
popu=newpopu;
end
fmin=1/Bestfi
BestS
%观测终了代数的极小值（解）
figure(1)                      %画每代适应度最小值fmin、平均值
t=1:G;
plot(t,1./bfi,'r')
xlabel('t');
ylabel('每代适应度最小值fmin');
T2=clock;
%计算运行时间
elapsed_time=T2-T1;
if elapsed_time(6)<0
elapsed_time(6)=elapsed_time(6)+60; elapsed_time(5)=elapsed_time(5)-1;
end  %秒
if elapsed_time(5)<0
elapsed_time(5)=elapsed_time(5)+60;elapsed_time(4)=elapsed_time(4)-1;
end  %分
if elapsed_time(4)<0
elapsed_time(4)=elapsed_time(4)+24;elapsed_time(3)=elapsed_time(3)-1;
end %时
str2=sprintf('程序运行耗时 %d 小时 %d 分钟 %.f 秒',elapsed_time(4),elapsed_time(5),elapsed_time(6));
disp(str2);
end
这是我最近根据其他人改进的适合自己的遗传算法，运行速度是6、7分钟，存在以下问题，结果每次相差太大，有时候就是100代都保持相同的值不变，有时候值特别的大，若是收敛值也不理想，请教有什么地方不合适需要改进的地方，请指点，谢谢

回复此楼