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

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

[复制链接]

该用户从未签到

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

EDA365欢迎您登录!

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

x
clc
1 ^! z) ^* o" R" q  d3 _clear all
0 P4 _/ O2 N9 w. R1 b1 m6 [5 p3 l' g% x, U: jclose all/ q" }, V1 ^3 M0 Y% W. u
%% 画出函数图$ A) y1 T) a5 Q+ G" o
figure(1);" Q+ b- Z. u! U7 f: F
lbx=-2;ubx=2;
2 U7 A4 Y5 r" m+ dlby=-2;uby=2;
, P' b6 S0 f6 n% \% F& jezmesh('y*sin(2*pi*x)+x*cos(2*pi*y)',[1bx,ubx,1by,uby],50);8 V: j  `8 [( T2 I* O- a
hold on;4 F% F4 [, Z; l' H8 t
%% 定义遗传算法参数
) {* M- X0 i7 I$ Y, VNIND=40;6 i* D. M$ h8 {# f
MAXGEN=20;/ G8 s' z7 H, U/ r
PRECI=20;
' p. n9 _' d5 \. P/ h* mGGAP=0.95;) |' S7 E1 N/ @8 T. g9 z7 C
px=0.7;
- a2 v5 c) {) }" b+ x9 c4 _pm=0.01;
, @$ b# I; o- y  y7 A5 Ztrace=zeros(3,MAXGEN);3 X0 n( _. s5 m! @+ ?/ ~
FieldD=[PRECI PRECI;lbx 1by;ubx uby;1 1;0 0;1 1;1 1];
4 d; L" o# w. V6 w  W6 I6 `Chrom=crtbp(NIND,PRECI*2);
5 M& o: S" _6 q%% 优化
# R7 t% m' b0 |0 V0 n' e" [gen=0;8 W7 W6 g- K9 v+ q* [# `
X=bs2rv(Chrom,FieldD);
4 A% Z0 n. R- EX=XY(:,1);Y=XY(:,2);
2 w" ]# H$ a0 t# X8 [ObjV=Y.*sin(2*pi*X)+X.*cos(2*pi*Y);7 }' X- G1 i* t# C" t7 o3 P  g, h
while gen<MAXGEN
5 F/ ]5 i% c1 EFitnV=ranking(-ObjV);
# @) J; U/ [$ |, a4 {6 zSelCh=select('sus',Chrom,FitnV,GGAP);
- u( i! C1 E, |" ^  ASelCh=recombin('xovsp',SelCh,px);6 o4 Z( i$ Q, v
SelCh=mut(SelCh,pm);  z9 a5 m$ U" G0 i
X=bs2rv(SelCh,FieldD);3 \" c+ f4 T0 `
X=XY(:,1);Y=XY(:,2);2 V* D3 I* m! }- Z$ c
ObjVSel=Y.*sin(2*pi*X)+X.*cos(2*pi*Y);
* n( E7 o7 E( a2 U[Chrom,ObjV]=reins(Chrom,SelCh,1,1,ObjV,ObjVSel);
" y: t3 H2 u; p  XX=bs2rv(Chrom,FieldD);) j) T$ T( u5 C0 }: s, X9 q$ E
gen=gen+1;
2 J% D, O% j' W. O/ m$ F& l& V%获取每代的最优解及其序号,Y为最优解,I为个体的序号* a7 d' x0 Y8 M3 b& r* }% U
[Y,I]=min(ObjV);
; A. R! W! f  Q4 M6 b% ^" I# o+ Qtrace(1,2,gen)=XY(I,;; R/ \( a8 h" r( ^& U* t3 \
trace(3,gen)=Y;1 g+ T% ?9 O' L6 J% k2 e3 C
end
" v9 U' Z( d" _! ^; F4 A# u. k) Fplot3(trace(1,,trace(2,,trace(3,:),'bo');
1 x  M. F; p/ }2 xgrid on;( R0 V& n" z* \! J5 T& r
plot3(XY(:,1),XY(:,2),ObjV,'bo');; f- _! t) P8 C7 j* o$ G& J
hold off
' h9 j/ j( K/ T( F) F9 V% V% i$ V%% 画进化图3 k7 F) B0 A  t0 w9 Q
figure(2);
. ~9 O0 R1 H2 ]3 A% \4 pplot(1:MAXGEN,trace(3,:));5 J$ k4 J  x& C: k  z1 z0 m) a
grid on# I) ~3 |8 }6 r7 v* B! u4 G
xlabel('遗传代数')
+ w6 d9 b" T% E, Cylabel('解的变化'); S8 m# s5 u: v5 ?
title('进化过程')+ s- j3 ?0 j( ?+ F8 ^
bestZ=trace(3,end);& s& @' [8 ~8 g! ~" B
bestY=trace(2,end);! T4 D$ n  K4 B5 O! m- W& K, G+ M1 W
bestX=trace(1,end);% A& ~5 ~& m! w& k. s
fprintf(['最优解:\nX=',num2str(bestX),'\nY=',num2str(bestY),'\nZ=',num2str(bestZ),'\n'])
0 m$ \6 q9 b8 ?/ k% x- _: x

该用户从未签到

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

本版积分规则

关闭

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

EDA365公众号

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

GMT+8, 2025-10-5 06:34 , Processed in 0.109375 second(s), 23 queries , Gzip On.

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

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

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