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

lupkpu9

: L$ a9 q5 B. f* l0 u* `: w, b) l
clc
5 k  Z' a' D  T7 Oclear all
close all
%% 画出函数图
( b0 t. H0 F  w; N% _% }$ V9 [figure(1);
lbx=-2;ubx=2;
# a- |% i4 M' J1 c, jlby=-2;uby=2;
; j# O& d' L* H2 n- `ezmesh('y*sin(2*pi*x)+x*cos(2*pi*y)',[1bx,ubx,1by,uby],50);
hold on;
7 |5 Z" }4 ~) c%% 定义遗传算法参数
NIND=40;
MAXGEN=20;
5 |" d! ~8 J% Z. a+ [/ z. O- cPRECI=20;
GGAP=0.95;
px=0.7;
pm=0.01;
' s$ Y* P& i3 I4 [trace=zeros(3,MAXGEN);
FieldD=[PRECI PRECI;lbx 1by;ubx uby;1 1;0 0;1 1;1 1];
Chrom=crtbp(NIND,PRECI*2);
1 d: h  V, b2 r3 S! m, J%% 优化
" M; P6 o+ X7 P. u9 K, `gen=0;
6 y" K8 j2 W- l' NX=bs2rv(Chrom,FieldD);
) {! f- w7 ]6 s5 PX=XY(:,1);Y=XY(:,2);
6 E+ M% Q* ^; s  M( qObjV=Y.*sin(2*pi*X)+X.*cos(2*pi*Y);
% }% N. l5 u7 J4 y- Dwhile gen<MAXGEN
( s" Z4 H3 ?0 d- K, m- n2 o  nFitnV=ranking(-ObjV);
: ~+ t- ]$ E& j) C2 D+ F, m" DSelCh=select('sus',Chrom,FitnV,GGAP);
SelCh=recombin('xovsp',SelCh,px);
SelCh=mut(SelCh,pm);
X=bs2rv(SelCh,FieldD);
- i) \- w* f; B( g: ^X=XY(:,1);Y=XY(:,2);
ObjVSel=Y.*sin(2*pi*X)+X.*cos(2*pi*Y);
[Chrom,ObjV]=reins(Chrom,SelCh,1,1,ObjV,ObjVSel);
X=bs2rv(Chrom,FieldD);
gen=gen+1;
: n" X$ F- z& l# _9 z+ R%获取每代的最优解及其序号，Y为最优解，I为个体的序号
3 C- c; U$ G0 V8 O3 Z[Y,I]=min(ObjV);
trace(1,2,gen)=XY(I,: );
trace(3,gen)=Y;
" g$ M! G$ @0 V  Y+ p+ hend
plot3(trace(1,: ),trace(2,: ),trace(3,: ),'bo');
grid on;
plot3(XY(:,1),XY(:,2),ObjV,'bo');
, n: r2 v7 v& K+ x/ lhold off
%% 画进化图
7 s- I. b6 x4 U5 |$ Y3 D/ ifigure(2);
plot(1:MAXGEN,trace(3,: ));
grid on
xlabel('遗传代数')
ylabel('解的变化')
/ P0 N4 O6 a1 D; y7 Z. Utitle('进化过程')
5 I! O' `# G) j5 l9 ?0 }' @bestZ=trace(3,end);
9 W$ ^2 v  r# ^! R6 F# tbestY=trace(2,end);
bestX=trace(1,end);
fprintf(['最优解:\nX=',num2str(bestX),'\nY=',num2str(bestY),'\nZ=',num2str(bestZ),'\n'])
' [% _9 {) x; U' x7 d1 p0 T

! z; E) l; _, D; ]2 j
! Q$ |9 ]* U  z

- L! J* {% A  N/ e3 H

yxlk

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