找回密码
 注册
关于网站域名变更的通知
查看: 609|回复: 1
打印 上一主题 下一主题

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

[复制链接]

该用户从未签到

跳转到指定楼层
1#
发表于 2019-10-14 09:30 | 只看该作者 回帖奖励 |正序浏览 |阅读模式

EDA365欢迎您登录!

您需要 登录 才可以下载或查看,没有帐号?注册

x
本帖最后由 lupkpu9 于 2019-10-14 18:55 编辑
. {: o% \1 x9 Y4 O
5 ~/ O( f6 Y& L9 ~; Xclc: b0 N4 D' M1 i
clear all/ `4 z* G) o. y3 h2 g6 R! J
close all/ o) V3 _0 h0 M- A
%% 画出函数图+ ?4 i1 O4 Q/ {# F
figure(1);
) k% u; e7 k/ @lbx=-2;ubx=2;
# Q- ?' t! k- P: Z/ A0 @0 `lby=-2;uby=2;0 S7 n% ?4 B+ c: [) ~$ ]0 M
ezmesh('y*sin(2*pi*x)+x*cos(2*pi*y)',[1bx,ubx,1by,uby],50);4 C8 x& v. W# ]! b+ j0 c, A" _
hold on;( n) k8 ~0 k+ f7 s
%% 定义遗传算法参数
: a( o! C" T- {) }9 A: r, s$ WNIND=40;6 D: m& Z) f' [5 @0 K  L! J: Z
MAXGEN=20;1 e+ O& C( x! a1 V% ]
PRECI=20;
7 G' z8 v( `1 m/ `& sGGAP=0.95;
8 y3 z- B5 N2 o1 [3 P3 [px=0.7;+ D+ G8 }* E) r. U2 m/ ^( O! |
pm=0.01;
9 }  _* S* i! x/ s3 Ttrace=zeros(3,MAXGEN);6 n! d( x1 _1 S
FieldD=[PRECI PRECI;lbx 1by;ubx uby;1 1;0 0;1 1;1 1];' P+ v' Z" G6 K
Chrom=crtbp(NIND,PRECI*2);
! p! K) z9 t# W( ^%% 优化
# B* b2 {# ^  b  v) |5 ~. ggen=0;
' g) ?! _2 [. Q$ r! d! |8 H, n/ B  |X=bs2rv(Chrom,FieldD);
9 Q" K7 e5 S5 I& d8 Z  nX=XY(:,1);Y=XY(:,2);
; M$ Y9 S+ \/ ]( i& YObjV=Y.*sin(2*pi*X)+X.*cos(2*pi*Y);2 w0 S  ^* Q& W
while gen<MAXGEN
9 u) J& [" p; C' E- d: b' w* {1 \FitnV=ranking(-ObjV);; Q) p1 T3 b& \
SelCh=select('sus',Chrom,FitnV,GGAP);8 L5 Q5 s& y7 C+ A- {
SelCh=recombin('xovsp',SelCh,px);
5 z; o9 I; U1 _- l5 j  t/ B' f8 vSelCh=mut(SelCh,pm);
( |6 A, m! m* _6 }- e& a; {X=bs2rv(SelCh,FieldD);
9 P; x# K! V, F8 H" M* tX=XY(:,1);Y=XY(:,2);
; U# w# a7 i* M7 i9 C- M' ~ObjVSel=Y.*sin(2*pi*X)+X.*cos(2*pi*Y);9 z$ `6 B2 ?& x
[Chrom,ObjV]=reins(Chrom,SelCh,1,1,ObjV,ObjVSel);
- ~# \/ K; z5 N2 w8 [2 {X=bs2rv(Chrom,FieldD);
4 `6 [4 ~) b3 k9 s2 I# Hgen=gen+1;
/ a  a/ k+ t) ]- ]% w5 z%获取每代的最优解及其序号,Y为最优解,I为个体的序号
+ R' i& E: o! ~8 S: T[Y,I]=min(ObjV);
/ T- q5 B& |9 L. }7 {trace(1,2,gen)=XY(I,: );2 U# t1 P6 w8 H
trace(3,gen)=Y;
/ Z4 f4 U) x2 `, B" i7 V, n/ Aend
& r1 w( y  Y/ S& l% c% O  H, hplot3(trace(1,: ),trace(2,: ),trace(3,: ),'bo');
8 e' P6 U9 O2 J# e% Zgrid on;
1 o, F1 t& Z1 V0 M4 d3 s2 Xplot3(XY(:,1),XY(:,2),ObjV,'bo');
! z5 H5 u: F+ h$ U) i0 }& |hold off  P" m7 R1 S/ J! l1 e" l6 t2 M
%% 画进化图
4 {" z& S" h/ g2 p. _figure(2);
* {# j2 X5 T3 j$ y% a: \plot(1:MAXGEN,trace(3,: ));
3 s/ A. ?( \; P' }grid on& z! D2 H8 \4 U4 K$ l
xlabel('遗传代数')
; D' s; h* K* Q) c/ N) R3 U1 K. d" Z- Cylabel('解的变化')$ V5 _! h% G7 x; w
title('进化过程')
# d4 E) D& f8 l3 V3 t% ~bestZ=trace(3,end);; X1 s1 V! z/ R. S; k9 j& P9 _3 a8 m
bestY=trace(2,end);2 `( m1 {2 S4 Q+ s  O: P9 k+ C+ h
bestX=trace(1,end);
7 r1 Y) a* ~) ^' xfprintf(['最优解:\nX=',num2str(bestX),'\nY=',num2str(bestY),'\nZ=',num2str(bestZ),'\n'])

4 T* m* h) `' {3 j
' n  y4 P3 J% s" T

7 b7 ~* N0 F7 @, v2 v) \+ N$ H/ e6 D6 b' p9 g

' T8 }: g% }' K" z! u" z5 L3 |! N: g

9 O/ j/ T2 P% P9 Q1 r8 e( R: z! n8 s1 b3 m$ A# ~
2 _, E! V3 j! D3 g

该用户从未签到

2#
发表于 2019-10-14 18:52 | 只看该作者
1bx和1by里面的数字1改成小写字母l(L的小写),这样才跟6和7行的相同啊
您需要登录后才可以回帖 登录 | 注册

本版积分规则

关闭

推荐内容上一条 /1 下一条

EDA365公众号

关于我们|手机版|EDA365电子论坛网 ( 粤ICP备18020198号-1 )

GMT+8, 2025-11-23 18:39 , Processed in 0.156250 second(s), 24 queries , Gzip On.

深圳市墨知创新科技有限公司

地址:深圳市南山区科技生态园2栋A座805 电话:19926409050

快速回复 返回顶部 返回列表