TA的每日心情 | 开心 2022-1-21 15:21 |
---|
签到天数: 1 天 [LV.1]初来乍到
|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
本帖最后由 Heaven_1 于 2022-9-27 13:15 编辑
( J: X8 z5 M9 L# T1 T
2 q: K# ?4 b, n! a7 q! nn=64;7 T+ ^' T4 c) n* ~1 T7 N
key=0.512;9 R) Q% Z. S5 O& G- j9 Q
an=linspace(3.56,4,176);
* B; K9 I, t6 Z) H, J jhold on;box on;axis([min(an),max(an),-0.5,1.5]);! y- o ] _, i" O
N=n^2;
8 ^, [4 x: m! F! [xn=zeros(1,N);
3 q4 B: q3 l& ^, _* q6 N {1 Tfor a1=an;
$ [. n7 u+ B& A& G0 P a2=an;
0 _; ~5 g4 y1 e. y y=key;
# s" p: d$ W% L5 G3 Rx=key;% {7 ?* ^( G8 A" m7 `
for k=1:16;7 @2 S8 H) B: o) j) a5 s1 Q y) k
y=a1*y*(1-y);
) Q# S. [9 q# r$ i) _: Q5 ex=(a2+(4-a2)*y).*x.*(1-x);%产生公式. G6 N: l8 Z% U9 u# n" S* r
end;2 p3 M8 j/ D* y8 F
for k=1:N;& X2 a; ~7 m! j3 z1 v
y=a1*y*(1-y);
) Q1 i5 ?# e- T. h' j# Nx=(a2+(4-a2)*y).*x.*(1-x);
6 p" e2 L2 n; E7 gxn(k,1)=x;! O2 u3 G" A* }1 r
b(k,1)=x;%一维矩阵记录迭代结果, b- N& X ]/ o% _) i# ?/ ^1 A
end;
J" q$ l% e; c8 b3 p1 ^, _7 f3 _plot(a*ones(1,N),xn,'k.','markersize',1);
: G4 W, }; T* T, mend;" Q( Z b3 i$ M% G O
figure;( L; g& X" }. l _5 G) b0 ~- A, M
imhist(b)/ n: E& b/ @# y6 a0 d
! r/ \ g: {* F
红色部分有问题
! S* A) R" s: C
6 d% D; l8 `1 o' n3 ^3 h* k: Z1 a0 T: J' n% k# H3 D. W6 V
|
|