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

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

[复制链接]

该用户从未签到

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

EDA365欢迎您登录!

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

x
clc% L  c8 ^1 k# s/ i$ K' X
clear all- K% F6 C1 `% T7 z5 n2 I
close all
2 O6 q8 Z" \$ _- x7 r%% 画出函数图$ e2 X) q, Q+ D. R& K6 q
figure(1);
: l: u. `" Q  \$ j8 Q8 slbx=-2;ubx=2; 1 f. u2 L9 y9 m0 @, v, @
lby=-2;uby=2;
1 Z" c% K: @) Y) _  jezmesh('y*sin(2*pi*x)+x*cos(2*pi*y)',[1bx,ubx,1by,uby],50);  Z  P. {# B8 J4 {, \. i8 O' a& U# w
hold on;
' f) l* s3 q/ B- ]! Q) l% Q%% 定义遗传算法参数  z; B& v! l1 O/ t
NIND=40;) V: U9 Y  h( ^- K( r, @; G1 e
MAXGEN=20;
+ T8 {9 z& g& aPRECI=20;
; Y6 [0 X) y' `+ L/ V+ sGGAP=0.95;2 h3 ^  j1 V! h' `. X" ?* u3 P4 Z
px=0.7;
, _2 i' Z- q; m7 |4 Apm=0.01;
" Y, T* n7 ?/ u7 M# V' Ltrace=zeros(3,MAXGEN);7 Y2 a  O* U. i6 N/ h: X
FieldD=[PRECI PRECI;lbx 1by;ubx uby;1 1;0 0;1 1;1 1];
5 y8 t, D' o$ o: {! j, fChrom=crtbp(NIND,PRECI*2);5 I- }" `& G* V# E. n" x
%% 优化' k% F4 D7 H# h  ]% J& R
gen=0;# e: R& Y. c! V. E1 Q7 y' X& M
X=bs2rv(Chrom,FieldD);
8 [+ z" d+ @6 J1 j2 ]& RX=XY(:,1);Y=XY(:,2);
5 _8 D3 |# l9 k1 S7 O( W# q# E. n4 kObjV=Y.*sin(2*pi*X)+X.*cos(2*pi*Y);
3 W# c& Q4 A( W) f+ P  S* @& wwhile gen<MAXGEN* G- }; n: `  G6 V0 T# P# Y
FitnV=ranking(-ObjV);* S, B6 `' {4 ?: U1 x% P  V
SelCh=select('sus',Chrom,FitnV,GGAP);  [4 H' u, V! o: j  A* u
SelCh=recombin('xovsp',SelCh,px);
- ~0 }5 G  B& \! r  ^6 S  |SelCh=mut(SelCh,pm);! A  z& _9 `; K  H! C+ U8 e) b9 G
X=bs2rv(SelCh,FieldD);
9 U$ M0 e; r8 m+ C7 J. _8 ]X=XY(:,1);Y=XY(:,2);
/ e4 V) f. {2 e' Z; Z6 \3 l- q6 BObjVSel=Y.*sin(2*pi*X)+X.*cos(2*pi*Y);+ X3 b  X, x8 l. g6 T* y( @
[Chrom,ObjV]=reins(Chrom,SelCh,1,1,ObjV,ObjVSel);
2 _# H+ F: t1 `X=bs2rv(Chrom,FieldD);. {. ~, ~2 D  O0 z
gen=gen+1;
8 G0 f, M" a: q$ h%获取每代的最优解及其序号,Y为最优解,I为个体的序号) Q* K$ t' k7 ]5 S; ^2 H: v3 n
[Y,I]=min(ObjV);5 T+ H4 L. C( G" L& N
trace(1,2,gen)=XY(I,;
7 B- @9 P2 b8 Q, j' s( d7 u6 x" P& t1 Mtrace(3,gen)=Y;" v8 A( a# {) v, W0 D" D
end# d! Y( h/ }. l+ p
plot3(trace(1,,trace(2,,trace(3,:),'bo');
: @6 M, g, L% s6 @3 t1 _6 tgrid on;
: E: h- a+ |) B/ U8 v' K1 E, a3 A( dplot3(XY(:,1),XY(:,2),ObjV,'bo');& z/ L, P$ \) d4 [, `' `! j8 Z
hold off4 D# o$ M0 M7 p* M% I3 z, i, t
%% 画进化图5 D: ~1 \) u2 ]6 g. C
figure(2);
) l; ^  j0 r) N0 [  T% Jplot(1:MAXGEN,trace(3,:));
5 b# g* q, }& u5 E( Hgrid on( K7 y0 x- F& p# J
xlabel('遗传代数')
, ~" Y* d" h3 B7 D4 h/ W' Rylabel('解的变化')
4 G* o5 u" n3 r+ v6 Otitle('进化过程')
0 F0 o7 ?0 i2 C% K; |# P0 FbestZ=trace(3,end);* V# T; M& r# V1 c& H" l
bestY=trace(2,end);, v) R4 E! g0 f. S
bestX=trace(1,end);
% T- Y2 B/ Y$ c2 U; Z+ [9 y3 Qfprintf(['最优解:\nX=',num2str(bestX),'\nY=',num2str(bestY),'\nZ=',num2str(bestZ),'\n'])
" _' [- Q8 ]4 d, v6 N* r

该用户从未签到

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

本版积分规则

关闭

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

EDA365公众号

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

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

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

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

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