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

lupkpu9

clc
2 e' Z* r) o. V/ ?( yclear all
close all
%% 画出函数图
( I, h$ X9 Q3 O7 u8 Bfigure(1);
; B& P' c$ _3 X' o) e6 }7 l+ ^8 I3 Blbx=-2;ubx=2;
lby=-2;uby=2;
ezmesh('y*sin(2*pi*x)+x*cos(2*pi*y)',[1bx,ubx,1by,uby],50);
hold on;
%% 定义遗传算法参数
NIND=40;
MAXGEN=20;
! V7 A- `: U; |# P1 l0 ^1 zPRECI=20;
GGAP=0.95;
px=0.7;
pm=0.01;
- U# D9 O; F- _! @# strace=zeros(3,MAXGEN);
FieldD=[PRECI PRECI;lbx 1by;ubx uby;1 1;0 0;1 1;1 1];
9 w" y2 B% w$ e4 N7 TChrom=crtbp(NIND,PRECI*2);
%% 优化
! t! h* |4 o) `7 B0 A# j& \6 ?, Rgen=0;
X=bs2rv(Chrom,FieldD);
X=XY(:,1);Y=XY(:,2);
ObjV=Y.*sin(2*pi*X)+X.*cos(2*pi*Y);
2 j" V% E" O7 S2 M. @while gen<MAXGEN
FitnV=ranking(-ObjV);
+ c- s3 t$ |0 {4 i9 j1 gSelCh=select('sus',Chrom,FitnV,GGAP);
5 S9 j7 Q) z! W! W" e3 cSelCh=recombin('xovsp',SelCh,px);
SelCh=mut(SelCh,pm);
X=bs2rv(SelCh,FieldD);
+ o2 L7 w2 }: V/ A4 rX=XY(:,1);Y=XY(:,2);
* L7 p3 n* d1 {) KObjVSel=Y.*sin(2*pi*X)+X.*cos(2*pi*Y);
9 N1 L9 R3 I, n. O5 a% A1 h[Chrom,ObjV]=reins(Chrom,SelCh,1,1,ObjV,ObjVSel);
" x7 Q0 `( ]2 hX=bs2rv(Chrom,FieldD);
) d9 G$ Z- {0 \( O  qgen=gen+1;
5 ~! ]' k9 ]/ @%获取每代的最优解及其序号，Y为最优解，I为个体的序号
[Y,I]=min(ObjV);
trace(1,2,gen)=XY(I,: );
trace(3,gen)=Y;
end
  p) |1 J, C# b9 K" _plot3(trace(1,: ),trace(2,: ),trace(3,: ),'bo');
6 c5 L( ]9 ~  z# ]grid on;
plot3(XY(:,1),XY(:,2),ObjV,'bo');
# ?2 f' V" m" Xhold off
%% 画进化图
figure(2);
$ g: y) G$ G6 aplot(1:MAXGEN,trace(3,: ));
grid on
; O' g2 E" c  n( Mxlabel('遗传代数')
ylabel('解的变化')
5 K0 L9 T) p8 d+ ptitle('进化过程')
bestZ=trace(3,end);
bestY=trace(2,end);
1 f! c1 t7 y6 @bestX=trace(1,end);
fprintf(['最优解:\nX=',num2str(bestX),'\nY=',num2str(bestY),'\nZ=',num2str(bestZ),'\n'])


4 G8 M8 \$ T( S8 t5 _+ m3 A
9 M/ }  f: t  r

8 O) T4 d6 G! _6 p1 V( q8 P
7 ?' N; Z, y5 D4 i4 e; C: t
! `" t& m" l! X1 T
) c( K- o# j% e/ l8 W

yxlk

1bx和1by里面的数字1改成小写字母l（L的小写），这样才跟6和7行的相同啊