24小时热门版块排行榜    

查看: 792  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[博后之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 Aj1rhIDL5ixY 2026-09-14 3/150 2026-09-15 20:59 by zNcEhTcH7ihq
[考博] 上海工程技术大学激光智能制造课题组招收博士研究生 +6 两三岁ss 2026-09-14 6/300 2026-09-15 20:17 by 只吃肉葡萄
[博后之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 17:59 by CgyNCDVNhGVg
[考研] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 s3fFTmArrBt6 2026-09-14 5/250 2026-09-15 16:59 by CgyNCDVNhGVg
[博后之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 s3fFTmArrBt6 2026-09-14 5/250 2026-09-15 16:50 by CgyNCDVNhGVg
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 s3fFTmArrBt6 2026-09-14 5/250 2026-09-15 16:35 by CgyNCDVNhGVg
[公派出国] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 s3fFTmArrBt6 2026-09-14 4/200 2026-09-15 16:26 by CgyNCDVNhGVg
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-13 4/200 2026-09-15 16:15 by CgyNCDVNhGVg
[考研] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +6 3n8v2C8RimXI 2026-09-13 7/350 2026-09-15 14:49 by VMDqgnlKimtZ
[考博] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +6 QUjhNVAcOSff 2026-09-13 6/300 2026-09-15 14:25 by VMDqgnlKimtZ
[博后之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +6 QUjhNVAcOSff 2026-09-13 7/350 2026-09-15 14:13 by VMDqgnlKimtZ
[公派出国] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +7 QUjhNVAcOSff 2026-09-13 9/450 2026-09-15 14:10 by VMDqgnlKimtZ
[考研] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +5 QUjhNVAcOSff 2026-09-13 5/250 2026-09-15 14:02 by VMDqgnlKimtZ
[博后之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 s3fFTmArrBt6 2026-09-14 3/150 2026-09-15 07:34 by C79jjtjAKjEn
[教师之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +4 s3fFTmArrBt6 2026-09-13 4/200 2026-09-15 07:02 by 5BDX0d0WFp7t
[教师之家] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 s3fFTmArrBt6 2026-09-13 4/200 2026-09-15 06:50 by 5BDX0d0WFp7t
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 s3fFTmArrBt6 2026-09-13 5/250 2026-09-15 06:49 by 5BDX0d0WFp7t
[公派出国] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 s3fFTmArrBt6 2026-09-13 5/250 2026-09-15 06:38 by 5BDX0d0WFp7t
[找工作] 浙江师范大学是怎么坑我的 +6 recruit123 2026-09-09 6/300 2026-09-14 12:00 by 萨仁其其格
[考博] 27秋季材料博士申请求助 10+3 Cahal1 2026-09-11 4/200 2026-09-14 08:50 by 北京莱茵编辑
信息提示
请填处理意见