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

lupkpu9

% E  K- i! B- Q: `4 ~  Xclc
9 x: ]( i1 d+ l& l: R2 e9 D2 Jclear all
close all
%% 画出函数图
figure(1);
' y& \2 ]) ]  m2 W9 U4 elbx=-2;ubx=2;
lby=-2;uby=2;
ezmesh('y*sin(2*pi*x)+x*cos(2*pi*y)',[1bx,ubx,1by,uby],50);
  W" b+ Z2 J/ S2 Fhold on;
- N; u2 @/ h" g  R7 |' r# |%% 定义遗传算法参数
NIND=40;
MAXGEN=20;
PRECI=20;
( \. f4 M9 O4 b' p/ gGGAP=0.95;
( x1 i7 x8 |9 u5 C9 Y. b0 G$ Vpx=0.7;
pm=0.01;
. w1 y: {  S/ V3 }3 s: O. Ntrace=zeros(3,MAXGEN);
FieldD=[PRECI PRECI;lbx 1by;ubx uby;1 1;0 0;1 1;1 1];
Chrom=crtbp(NIND,PRECI*2);
# N/ p8 q; x( x& S" ?%% 优化
- ^+ x+ B9 P7 ^gen=0;
X=bs2rv(Chrom,FieldD);
" H: _$ B# r# B6 DX=XY(:,1);Y=XY(:,2);
  f8 L8 q" V8 N" m! SObjV=Y.*sin(2*pi*X)+X.*cos(2*pi*Y);
# S; O) M9 _9 G5 W2 Cwhile gen<MAXGEN
FitnV=ranking(-ObjV);
  h) Y% j0 S2 ^" F' aSelCh=select('sus',Chrom,FitnV,GGAP);
SelCh=recombin('xovsp',SelCh,px);
SelCh=mut(SelCh,pm);
8 u8 w8 O8 W2 `+ e/ YX=bs2rv(SelCh,FieldD);
X=XY(:,1);Y=XY(:,2);
. e5 R3 k5 P7 v0 `# L! _ObjVSel=Y.*sin(2*pi*X)+X.*cos(2*pi*Y);
, x  B0 m2 Q2 K* ~7 f5 e% o$ |5 R[Chrom,ObjV]=reins(Chrom,SelCh,1,1,ObjV,ObjVSel);
X=bs2rv(Chrom,FieldD);
gen=gen+1;
4 u0 {( ?4 z6 s& O%获取每代的最优解及其序号，Y为最优解，I为个体的序号
[Y,I]=min(ObjV);
! {3 R! c2 S' |- z- f" n! Ntrace(1,2,gen)=XY(I,: );
trace(3,gen)=Y;
8 w9 j9 d; O/ Q. Gend
plot3(trace(1,: ),trace(2,: ),trace(3,: ),'bo');
* s2 M- T' {- q, ?grid on;
plot3(XY(:,1),XY(:,2),ObjV,'bo');
% x/ [2 x* w% Ohold off
, n/ a2 `7 _1 ~  }1 I) d%% 画进化图
figure(2);
plot(1:MAXGEN,trace(3,: ));
# s2 E$ Y- c/ l; ygrid on
% `" \- x$ R5 H  G3 d- k, yxlabel('遗传代数')
ylabel('解的变化')
title('进化过程')
bestZ=trace(3,end);
bestY=trace(2,end);
bestX=trace(1,end);
fprintf(['最优解:\nX=',num2str(bestX),'\nY=',num2str(bestY),'\nZ=',num2str(bestZ),'\n'])
" Z' N3 f% X% L7 V# `
1 O+ P" s3 o2 r, n7 g4 e
0 {0 k* v& y; [/ a2 j
2 V0 b' M4 e7 m% _7 }. m' l
- Y! y. }6 c! N) J4 ]+ _6 [
3 s9 C5 R. \8 m& W
, d2 z! ?& Y& ?; a2 e5 Z

yxlk

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