找回密码
 注册
关于网站域名变更的通知
查看: 1072|回复: 1
打印 上一主题 下一主题

文件:Untitled.m 行:8 列:40 不是所需的 MATLAB 表达式?

[复制链接]

该用户从未签到

跳转到指定楼层
1#
发表于 2019-8-16 14:35 | 只看该作者 回帖奖励 |倒序浏览 |阅读模式

EDA365欢迎您登录!

您需要 登录 才可以下载或查看,没有帐号?注册

x
clc$ K) G. j$ T- q7 `- g5 i4 r8 H
clear all
' d3 Z, x+ I- U8 Q! Eclose all
  Z  J" K9 D% z%% 画出函数图3 c) C6 \& x9 S( V4 Q6 }
figure(1);3 j0 `/ T+ W8 \2 J; w  }  A
lbx=-2;ubx=2;
) u) p3 ]" @; @; x) L! I# s+ h9 klby=-2;uby=2;6 a1 ]9 g( `. e! g& O; Z6 |
ezmesh('y*sin(2*pi*x)+x*cos(2*pi*y)',[1bx,ubx,1by,uby],50);5 z, h5 @+ i' D
hold on;
* c6 z# [4 \# v: Z  h0 ?( g%% 定义遗传算法参数
) R6 _3 B( V1 D% D- UNIND=40;: y9 v+ z, y# z4 D
MAXGEN=20;1 Y9 p8 ]8 L+ A8 q
PRECI=20;
, Y: e3 ?+ G! A7 H) X  QGGAP=0.95;$ c8 n4 A+ q) r! w* ^
px=0.7;3 n+ n8 U) A% q. E
pm=0.01;
- r# F: Y9 A! q+ ~1 H# l& j# Qtrace=zeros(3,MAXGEN);
) U! c# y; N# j) k- p# \6 v$ X3 MFieldD=[PRECI PRECI;lbx 1by;ubx uby;1 1;0 0;1 1;1 1];
0 ]+ ]  K+ _1 D8 _6 T6 Y. L( m) fChrom=crtbp(NIND,PRECI*2);
- t' X! y7 R; f%% 优化
! z. ]2 L& `* b3 Q4 v6 \$ Sgen=0;6 E0 X$ ~3 C' T1 v* x' U" x
X=bs2rv(Chrom,FieldD);: s+ h8 E$ \  n; ?: K
X=XY(:,1);Y=XY(:,2);
4 M% I  C) G2 H0 w6 `* Q( jObjV=Y.*sin(2*pi*X)+X.*cos(2*pi*Y);5 y5 f( ]3 k" D+ Y, g
while gen<MAXGEN
6 F" ]* l: W' {1 u0 JFitnV=ranking(-ObjV);) I; _' @' z; R9 a1 b; x  p( l
SelCh=select('sus',Chrom,FitnV,GGAP);0 f; ^8 Z/ E2 `2 v$ ]# P
SelCh=recombin('xovsp',SelCh,px);0 R: T" U# Z& N$ u" e% P
SelCh=mut(SelCh,pm);
8 J: c. O9 G7 q, S# UX=bs2rv(SelCh,FieldD);
- Q0 ^2 J8 I! L4 e- \X=XY(:,1);Y=XY(:,2);
+ t/ i( P9 U* ?5 zObjVSel=Y.*sin(2*pi*X)+X.*cos(2*pi*Y);2 e3 l, h& V# o0 j7 e
[Chrom,ObjV]=reins(Chrom,SelCh,1,1,ObjV,ObjVSel);
. j+ ^, J) Y3 M8 TX=bs2rv(Chrom,FieldD);6 k5 v- n3 w5 A) \3 u+ y. _
gen=gen+1;
3 O; B" b3 \  U  ^6 s# O! k%获取每代的最优解及其序号,Y为最优解,I为个体的序号/ |) N. A' K& e+ n5 _% K2 I" m
[Y,I]=min(ObjV);
7 Q" ]8 v4 Y1 ~trace(1,2,gen)=XY(I,;
0 Q+ J8 T  `3 btrace(3,gen)=Y;0 J# ?" S3 o+ L( L, T, }1 F  ^
end9 E" U6 W" H" `4 R
plot3(trace(1,,trace(2,,trace(3,:),'bo');
* o% U& Q* z) @" \: C! w) Ugrid on;
' I* n; ]# d3 Y1 m( `4 Iplot3(XY(:,1),XY(:,2),ObjV,'bo');% v" I. O1 O& F8 y
hold off" _: [3 i6 m7 i. m3 ?3 V! ^
%% 画进化图
- l: t; \8 w! C/ l# K0 E1 X# p3 {! afigure(2);: d4 X& G5 e3 {6 t, n+ B) s
plot(1:MAXGEN,trace(3,:));
+ T: g  C4 K2 G4 `grid on& N, o9 v! i7 e+ f  q, o
xlabel('遗传代数')1 a: d8 T* L% @5 h
ylabel('解的变化')4 L+ z. y. |6 G, r! E: C" W
title('进化过程')
; j* y6 S) D7 V) b* }/ G1 T$ `3 FbestZ=trace(3,end);# j8 ]7 m( b+ ]& J1 \$ A0 y
bestY=trace(2,end);" z2 H& G! M5 G! W
bestX=trace(1,end);
2 W4 t% X- ~$ r6 r2 ?fprintf(['最优解:\nX=',num2str(bestX),'\nY=',num2str(bestY),'\nZ=',num2str(bestZ),'\n'])
3 @4 Y* O7 ?; x

该用户从未签到

2#
发表于 2019-8-16 17:44 | 只看该作者
1bx和1by里面的数字1改成小写字母l(L的小写),这样才跟6和7行的相同啊
您需要登录后才可以回帖 登录 | 注册

本版积分规则

关闭

推荐内容上一条 /1 下一条

EDA365公众号

关于我们|手机版|EDA365电子论坛网 ( 粤ICP备18020198号-1 )

GMT+8, 2025-11-23 23:33 , Processed in 0.140625 second(s), 24 queries , Gzip On.

深圳市墨知创新科技有限公司

地址:深圳市南山区科技生态园2栋A座805 电话:19926409050

快速回复 返回顶部 返回列表