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

lupkpu9

, d5 E0 P% _9 ~2 |clc
clear all
/ X) \: V  q9 T  m( q6 n6 @4 n- i2 ~: xclose all
' M  X/ f6 P$ U2 q- f( a! c; y! q( P%% 画出函数图
figure(1);
& g: v! G* Y4 _# mlbx=-2;ubx=2;
lby=-2;uby=2;
* L9 ]' J! |$ {ezmesh('y*sin(2*pi*x)+x*cos(2*pi*y)',[1bx,ubx,1by,uby],50);
; T2 C) X) o! w% n. E; j# ihold on;
9 a9 I* R8 r3 o7 b, h3 w%% 定义遗传算法参数
& ~2 s  f( t% f* U! YNIND=40;
) {6 e# Z* X+ K* H; f* f$ A& a' {" `MAXGEN=20;
PRECI=20;
GGAP=0.95;
px=0.7;
pm=0.01;
1 {7 W5 ]9 M: t+ R* `2 L( ytrace=zeros(3,MAXGEN);
FieldD=[PRECI PRECI;lbx 1by;ubx uby;1 1;0 0;1 1;1 1];
% t/ o/ J  o5 Z, D2 M- s6 hChrom=crtbp(NIND,PRECI*2);
%% 优化
gen=0;
" I- {0 t& Y+ xX=bs2rv(Chrom,FieldD);
X=XY(:,1);Y=XY(:,2);
1 w& p$ l7 t* T' |- x% C  T; o7 |: VObjV=Y.*sin(2*pi*X)+X.*cos(2*pi*Y);
9 @8 k8 }( w$ G, v% ~while gen<MAXGEN
+ m! ^3 u+ J4 c& E+ RFitnV=ranking(-ObjV);
SelCh=select('sus',Chrom,FitnV,GGAP);
SelCh=recombin('xovsp',SelCh,px);
SelCh=mut(SelCh,pm);
' Q+ S! H7 P' c" A  V3 t1 \X=bs2rv(SelCh,FieldD);
X=XY(:,1);Y=XY(:,2);
ObjVSel=Y.*sin(2*pi*X)+X.*cos(2*pi*Y);
. k+ S  p7 w1 v+ ]5 M9 M4 y[Chrom,ObjV]=reins(Chrom,SelCh,1,1,ObjV,ObjVSel);
X=bs2rv(Chrom,FieldD);
gen=gen+1;
%获取每代的最优解及其序号，Y为最优解，I为个体的序号
' C; T) V3 K( C* U[Y,I]=min(ObjV);
trace(1,2,gen)=XY(I,: );
& J+ L- @  i  e4 G: f/ x% j; w' Ptrace(3,gen)=Y;
* a9 i7 ?' x. ^$ `0 y2 gend
9 |  N+ o6 I0 {plot3(trace(1,: ),trace(2,: ),trace(3,: ),'bo');
% E5 H' u, j) ?: l  W; Ogrid on;
) J4 ^3 G1 ?5 |* H  k* mplot3(XY(:,1),XY(:,2),ObjV,'bo');
hold off
' t+ ]; V( I5 a%% 画进化图
' H# m; g: R% ]figure(2);
plot(1:MAXGEN,trace(3,: ));
grid on
xlabel('遗传代数')
ylabel('解的变化')
8 W- S7 X0 P) G) p/ Y+ Y0 Ttitle('进化过程')
/ C) c/ H$ h/ y, qbestZ=trace(3,end);
2 G- Y/ S0 Q! j, dbestY=trace(2,end);
/ b' w  [! }% A2 T; ~9 i6 ?bestX=trace(1,end);
  l# f$ V8 G2 vfprintf(['最优解:\nX=',num2str(bestX),'\nY=',num2str(bestY),'\nZ=',num2str(bestZ),'\n'])
  @: u$ S+ ]3 Q3 m




# x4 b$ `* `7 S
: Z: s/ e4 K/ L* ]" g) D1 O4 n

yxlk

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