24小时热门版块排行榜    

查看: 422  |  回复: 1

s20090281

银虫 (小有名气)

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

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代都保持相同的值不变,有时候值特别的大,若是收敛值也不理想,请教有什么地方不合适需要改进的地方,请指点,谢谢
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

s20090281

银虫 (小有名气)

请大家帮个忙,谢谢
2楼2011-07-01 15:18:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 s20090281 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 关于代码变化问题,想知道的进来 +6 且听虎啸 2026-08-07 9/450 2026-08-07 18:40 by naalan7001
[基金申请] 关于filecode +3 布布和一二 2026-08-07 6/300 2026-08-07 17:51 by 医学老男孩
[基金申请] 固定端突然变了,今天 +6 archvillain 2026-08-06 10/500 2026-08-07 16:03 by 医学老男孩
[基金申请] 听说今天filecode变了 +24 布布和一二 2026-08-06 47/2350 2026-08-07 16:02 by zhiyanjiang
[基金申请] filecode变化情况 +6 布布和一二 2026-08-07 22/1100 2026-08-07 14:45 by 且听虎啸
[论文投稿] 十年后又回来了,论文投稿求助 +4 哈哈114477 2026-08-01 4/200 2026-08-07 14:39 by jgy194592
[基金申请] 关于豆爷回答的JTJC与%2F数量 +4 yang182083 2026-08-06 6/300 2026-08-07 12:43 by 虫友是什么虫
[基金申请] 面上再次挂了,太难了,躺也躺不了,倦也卷不过,小学校之殇! +23 低垂的野花 2026-07-31 31/1550 2026-08-07 12:22 by speculiar
[基金申请] filecode +14 等待解的谜 2026-08-06 19/950 2026-08-07 12:20 by wlwhappy
[基金申请] filecode +8 布布和一二 2026-08-06 11/550 2026-08-06 20:41 by tangpu318
[基金申请] 求各位大神看下 100+6 hpkpkpkp 2026-08-05 33/1650 2026-08-06 14:49 by zhiyanjiang
[基金申请] 影响面上的因素 +8 布布和一二 2026-08-05 11/550 2026-08-06 10:41 by 宝贝虫子
[基金申请] 8月时间戳变的,举个手。玩一下,释放压力 +9 archvillain 2026-08-04 11/550 2026-08-05 20:06 by wlwhappy
[基金申请] 好消息?这个有何含义??? +8 Tide man 2026-08-05 10/500 2026-08-05 16:14 by xmuxiaoyu
[考博] 【2027博士申请】纳米药物递送方向 20+3 13586093586 2026-08-03 4/200 2026-08-05 09:59 by lfy8008
[基金申请] 有没有H口的?有收到消息的吗? +3 超级海虾 2026-08-04 3/150 2026-08-04 17:26 by 学教育滴
[基金申请] 纯娱乐,不喜欢勿喷 +7 Tide man 2026-08-04 10/500 2026-08-04 15:10 by loufangrui
[基金申请] 面上提前没消息,有中的吗 +14 archvillain 2026-08-02 18/900 2026-08-04 14:42 by archvillain
[基金申请] 娱乐 +4 Tide man 2026-08-03 4/200 2026-08-04 11:51 by wgch518
[基金申请] 什么时候能放榜呀? +3 Jacob678 2026-08-03 3/150 2026-08-03 16:14 by gltch
信息提示
请填处理意见