24小时热门版块排行榜    

查看: 782  |  回复: 2

witzwq

银虫 (小有名气)

[求助] 遗传算法matlab程序求解

figure(1);
fplot('variable.*sin(10*pi*variable)+2.0',[-1,2]);%画出函数曲线
%定义遗传算法参数
NIND=40;%个体数目(Number of individuals)
MAXGEN=25;%最大遗传代数(Maximum number.of generations)
PRECI=20;%变量的二进制位数(Precision of variables)
GGAP=0.9;%代沟(Generation gap)
trace=zeros(2,MAXGEN);%寻优结果的初始值
FieldD=[20;-1;2;1;0;1;1];%区域描述器(Build field descriptor)格雷编码,算术刻度,包含边界
Chrom=crtbp(NIND,PRECI);%初始种群
gen=0;%代计数器
variable=bs2rv(Chrom,FieldD);%计算初始种群的十进制转换
ObjV=variable.*sin(10*pi*variable)+2.0;%计算目标函数值
while gen     FitnV=ranking(-ObjV);%分配适应度值(Assign fitness values)
    SelCh=select('sus',Chrom,FitnV,GGAP);%选择
    SelCh=recombin('xovsp',SelCh,0.7);%重组
    SelCh=mut(SelCh);%变异
    variable=bs2rv(SelCh,FieldD);%子代个体的十进制转换
    ObjVSel=variable.*sin(10*pi*variable)+2.0;%计算子代的目标函数值
    [Chrom ObjV]=reins(Chrom,SelCh,1,1,ObjV,ObjVSel);%重插入子代的新种群
    gen=gen+1;%代计数器增加
    %输出最优解及其序号,并在目标函数图像中标出,Y为最优解,I为种群的序号
    [Y,I]=max(ObjV);
    hold on;
    plot(variable(I),Y,'bo');%b是blue,o是圆圈
    trace(1,gen)=max(ObjV);%遗传算法性能跟踪
    trace(2,gen)=sum(ObjV)/length(ObjV);
end
variable=bs2rv(Chrom,FieldD);%最优个体的十进制转换
hold on ;
grid on;%显示网格线
plot(variable',ObjV','b*');
figure(2);
plot(trace(1,');%plot函数可以接受字符串输入变量,这些字符串输入变量用来指定不同的颜色、线型和标记符号。:表示点线
hold on;
plot(trace(2,','-.');grid;%-.表示点划线
legend('解的变化','种群均值的变化');%图形加标注
   


麻烦问一下 染色体数目会变为36?
怎么样解决,使其一直为40?
回复此楼

» 收录本帖的淘帖专辑推荐

导航、制导与控制 matlab

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

web22

新虫 (初入文坛)

【答案】应助回帖

★ ★ ★ ★ ★
witzwq: 金币+5, 有帮助 2013-05-13 20:06:02
能不能将crtbp也贴上来,上面这个执行不了
2楼2013-05-02 23:39:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

web22

新虫 (初入文坛)

搞错了,是我的matlab没装遗传算法工具箱
3楼2013-05-02 23:52:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 学员Uaj2Wn 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 看来今天不会放榜了? +6 chengyan1220 2026-08-21 6/300 2026-08-21 14:47 by 20120902066
[基金申请] 我面上完蛋了 +7 且听虎啸 2026-08-20 8/400 2026-08-21 12:31 by 酷酷墨镜
[基金申请] 只有每年这种时候来逛逛小木虫 +23 yaoyewhu2008 2026-08-20 25/1250 2026-08-21 12:19 by Poppy1104
[基金申请] 时间戳又变了8-15 +15 archvillain 2026-08-15 27/1350 2026-08-21 11:32 by tim76
[基金申请] 今天基金会出结果吗?20260819 +15 kkkl_v 2026-08-19 16/800 2026-08-21 11:22 by 365372687
[基金申请] 科研孤儿太难了 +15 我4大白菜 2026-08-20 16/800 2026-08-21 11:07 by yming1319
[基金申请] 应该是下周三26日公布了吧? +4 哈哈蛤? 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[基金申请] 今天放榜没戏了吧 +9 yuleib84 2026-08-19 11/550 2026-08-21 10:06 by gltch
[基金申请] 放榜前的不淡定 40+4 snowwithsea 2026-08-19 13/650 2026-08-21 08:30 by 北京莱茵编辑
[基金申请] 估计是周四 +3 archvillain 2026-08-18 3/150 2026-08-21 01:48 by jnhyjjm
[基金申请] 今天系统多次维护,明天很可能放榜! +10 zju2000 2026-08-16 11/550 2026-08-20 20:02 by cl479861084
[基金申请] 时间戳变了,能看出什么问题? +18 基诺咪客 2026-08-17 23/1150 2026-08-20 17:19 by Godzela
[基金申请] 重要消息,中午系统在维护 +11 yuleib84 2026-08-18 12/600 2026-08-20 11:09 by xskun
[基金申请] 2027广东省杰青 +4 奶牛小黑 2026-08-15 10/500 2026-08-19 11:02 by wanfengnew
[基金申请] 朋友圈看到的 +6 wangzilk 2026-08-18 8/400 2026-08-19 10:55 by Haru815
[基金申请] 快农历七夕节了,轻松一下,男人悄悄话,女施主请不要进来。 +6 Tide man 2026-08-14 7/350 2026-08-19 09:56 by ZJTJZ
[基金申请] 明天放榜? +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
[基金申请] 93BebMhtakh前后11位开头都是大写 +4 且听虎啸 2026-08-17 5/250 2026-08-18 00:49 by 蔡棒棒菂
[精细化工] 招聘 金属平磨液,抛光液研发工程师 +3 小天0311 2026-08-14 3/150 2026-08-16 07:31 by H9PLUS
[基金申请] 各位道友,我要去昆明玩几天,回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
信息提示
请填处理意见