TA的每日心情 | 开心 2022-1-21 15:21 |
---|
签到天数: 1 天 [LV.1]初来乍到
|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
本帖最后由 Heaven_1 于 2022-9-27 13:15 编辑
2 |1 f# X' L1 ^1 {- [
- t7 L3 J) u( `n=64;
# h% a% u6 A' m/ x. F2 _key=0.512;
) u& f& N$ o) D9 f* _) i+ p+ uan=linspace(3.56,4,176);
& c4 k: L7 g5 j9 \# {+ ?hold on;box on;axis([min(an),max(an),-0.5,1.5]);
6 Z: T+ ]& O2 _) n* w% rN=n^2;
) v5 W* F" C' n Sxn=zeros(1,N);/ M0 a( D- v! \& w1 Z
for a1=an;9 B0 e. R6 n1 O' B f% S; N) p
a2=an;
: u+ K4 ]* D; p4 n! g% P/ T5 t4 D& q. [ y=key;) A. a m% A& A! n+ P
x=key;
. c5 }# t) y! [! Zfor k=1:16;3 Z7 \! Y S1 V4 H8 f2 F. f% i
y=a1*y*(1-y);; _" t9 l. w0 R9 }5 p3 {5 u
x=(a2+(4-a2)*y).*x.*(1-x);%产生公式
j6 U2 T- ]$ w o; y- Xend;
' P: ^5 m; t6 G! Vfor k=1:N;. W9 _) f! X* i# k ^
y=a1*y*(1-y);: y+ [$ N( M# [0 F# y# s6 Q3 H
x=(a2+(4-a2)*y).*x.*(1-x);
$ R9 y- _/ d0 o8 }9 n* V F- B3 Lxn(k,1)=x;
, ]1 a. H" ^: c) Y: cb(k,1)=x;%一维矩阵记录迭代结果
& I( ]- ]& a6 n! v/ Kend;0 ~6 @3 M+ |5 |- k: b
plot(a*ones(1,N),xn,'k.','markersize',1);4 }$ b6 m9 r) x+ [6 O
end;$ n: S' w- z Y
figure;
, n: [. o$ s% U5 `imhist(b)
* z% B* T+ L4 j; I) Q, m1 S: ^
" J6 I4 f+ [; X* g红色部分有问题
+ s6 d( b. D: \8 m! ]7 ?0 v$ `% v% Z) j7 Y
w9 G# I+ N9 l4 L: d9 A |
|