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

qpggup

clc
clear all
close all
%% 画出函数图
/ g2 C$ B4 j1 }- s: ifigure(1);
lbx=-2;ubx=2;
lby=-2;uby=2;
4 ?$ e! R- @, a% h# F8 v3 A* zezmesh('y*sin(2*pi*x)+x*cos(2*pi*y)',[1bx,ubx,1by,uby],50);
hold on;
7 f& H" [) J9 e, p7 q%% 定义遗传算法参数
NIND=40;
+ S5 G" f( H" ~5 n6 E3 t1 }MAXGEN=20;
PRECI=20;
GGAP=0.95;
: L6 c" J! `3 j- U* ^: y3 ^px=0.7;
pm=0.01;
trace=zeros(3,MAXGEN);
FieldD=[PRECI PRECI;lbx 1by;ubx uby;1 1;0 0;1 1;1 1];
: g9 [, Q) r9 g) A; M) a/ c1 o; qChrom=crtbp(NIND,PRECI*2);
%% 优化
gen=0;
4 U$ A! Y: o* j0 h2 j6 cX=bs2rv(Chrom,FieldD);
X=XY(:,1);Y=XY(:,2);
ObjV=Y.*sin(2*pi*X)+X.*cos(2*pi*Y);
while gen<MAXGEN
' g0 Q& `1 G  H8 M7 ]- H% zFitnV=ranking(-ObjV);
SelCh=select('sus',Chrom,FitnV,GGAP);
SelCh=recombin('xovsp',SelCh,px);
SelCh=mut(SelCh,pm);
9 H* o/ Y: m/ V# e" nX=bs2rv(SelCh,FieldD);
X=XY(:,1);Y=XY(:,2);
ObjVSel=Y.*sin(2*pi*X)+X.*cos(2*pi*Y);
' ^  j3 ?- R2 {3 r[Chrom,ObjV]=reins(Chrom,SelCh,1,1,ObjV,ObjVSel);
X=bs2rv(Chrom,FieldD);
gen=gen+1;
%获取每代的最优解及其序号，Y为最优解，I为个体的序号
* G! C9 Y3 ^' M( m[Y,I]=min(ObjV);
- z, ~. f; t1 i) w6 k. y! }trace(1,2,gen)=XY(I,;
trace(3,gen)=Y;
" _0 R2 _* D$ N, T# Hend
plot3(trace(1,,trace(2,,trace(3,:),'bo');
0 r* q# s: n% h1 Fgrid on;
plot3(XY(:,1),XY(:,2),ObjV,'bo');
hold off
%% 画进化图
3 N8 ^& @0 A+ u5 q' O2 I  tfigure(2);
# l% b0 [5 P, N1 u3 u3 dplot(1:MAXGEN,trace(3,:));
grid on
xlabel('遗传代数')
  B/ C6 B4 S9 v. ]ylabel('解的变化')
title('进化过程')
bestZ=trace(3,end);
bestY=trace(2,end);
bestX=trace(1,end);
( w  u  j( C& b" Z3 ~fprintf(['最优解:\nX=',num2str(bestX),'\nY=',num2str(bestY),'\nZ=',num2str(bestZ),'\n'])
# q# i7 g$ _9 U% A: C/ W5 Z0 I

relchhiclty

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