24小时热门版块排行榜    

查看: 436  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 麻烦专家们看看评委们的意见(F口面上) +7 gdd2018 2026-08-28 12/600 2026-09-01 08:32 by 尼古拉斯小虫
[基金申请] 怎么看青基中了没有啊 +6 叶九微 2026-08-26 6/300 2026-08-31 23:54 by yudaoqian88
[基金申请] 投票:  有多少人是今天查系统知道结果的? +17 爱看书的可乐 2026-08-26 19/950 2026-08-31 21:30 by xiangy672
[基金申请] 哪位高人中了,把查询到的截图贴出来让我看看,让我长长见识 +6 yuleib84 2026-08-26 7/350 2026-08-31 19:46 by 鱼翔浅底1
[基金申请] 能否申诉? +7 echo8914667 2026-08-30 8/400 2026-08-31 17:00 by yihongxu
[文学芳草园] 梦想 +5 myrtle 2026-08-26 5/250 2026-08-31 14:40 by hahaboy
[基金申请] 中青基了要发朋友圈吗? +7 349506619 2026-08-28 7/350 2026-08-31 13:39 by 冼亮淀粉酶
[基金申请] 为什么到现在没收到通知? +5 tannykie 2026-08-29 5/250 2026-08-30 21:05 by purplejack
[基金申请] 有没有仍没收到信息的 +7 德尚中行 2026-08-27 8/400 2026-08-30 20:52 by purplejack
[基金申请] 我就是申请一个面上项目而已,这评审意见是按照杰青的条件评的吧? +6 gouxfjh 2026-08-28 11/550 2026-08-30 07:57 by gouxfjh
[基金申请] 2026年叶企孙基金 +4 bud_bud 2026-08-27 7/350 2026-08-29 07:23 by foolishmani
[基金申请] 为什么资助数各大高校都创新高,自己申请怎么就这么难 +12 Kittylucky 2026-08-27 13/650 2026-08-29 00:04 by superceng
[基金申请] 系统查不到 +11 董八千 2026-08-26 11/550 2026-08-28 18:06 by Leogzhya
[基金申请] 面上合作单位盖章 +5 ssyjh 2026-08-27 5/250 2026-08-27 20:50 by gdfollow
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜? 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[基金申请] 牛来!米来!面来! +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
[基金申请] 某些机构,以效率低为荣,以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
信息提示
请填处理意见