TA的每日心情 | 开心 2022-1-21 15:21 |
---|
签到天数: 1 天 [LV.1]初来乍到
|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
本帖最后由 Heaven_1 于 2022-9-27 13:15 编辑 7 D1 G7 V. Q- b6 N
8 H# [2 Y7 @% v; e9 r% |0 W
n=64;4 T: W8 C- F4 c
key=0.512;6 r( J6 b1 u3 S0 }, L; Y
an=linspace(3.56,4,176);
) `9 Z- V. i# U2 G- [hold on;box on;axis([min(an),max(an),-0.5,1.5]);
8 K( f \- ]; q+ A/ Y- LN=n^2;) E8 K* p3 n+ H# \/ w! s- U
xn=zeros(1,N);
9 k/ ~4 q+ y0 L y3 jfor a1=an;
1 M& s& \- \, U. ?! f a2=an;
( G! K3 ]3 V; P6 _ y=key;
. V1 t1 ~* [7 z: A" rx=key;
# x2 q3 @% I( J5 p1 t! s) w: \for k=1:16;" _& z7 C) e; V/ c$ j
y=a1*y*(1-y);
& F; ~8 E+ D9 Y0 t3 r; t; fx=(a2+(4-a2)*y).*x.*(1-x);%产生公式7 I2 }4 t5 ?7 w9 B
end;
, F; N& c( m$ q+ D: I" F9 |for k=1:N;' Z9 o; m/ K% q+ u8 s6 d) B
y=a1*y*(1-y);1 u* b7 W+ W( b+ Z" [
x=(a2+(4-a2)*y).*x.*(1-x);9 ^1 y! l% p5 k. G2 _; w0 k
xn(k,1)=x;. C' D# W8 c6 w) W. h+ i+ L# \4 B
b(k,1)=x;%一维矩阵记录迭代结果* g# Y; Q* g- q/ F! p4 B- M* A
end;- M$ X- ?0 {2 x1 a- X5 V$ U
plot(a*ones(1,N),xn,'k.','markersize',1);
( s/ J. g6 m9 q' Dend;
- g1 b) p! N0 N t* s7 ofigure;6 x: v3 ?% ]+ r
imhist(b)" P8 a6 J7 c0 c8 x+ P& s
$ r; _! d) \: A) P ?9 U! X
红色部分有问题- g5 B; r$ j0 n1 N7 ]" s8 ` k
/ z/ @, [ P A" W3 X4 U
" }$ X! [) U4 M6 O8 Q
|
|