TA的每日心情 | 开心 2022-1-21 15:21 |
|---|
签到天数: 1 天 [LV.1]初来乍到
|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
本帖最后由 Heaven_1 于 2022-9-27 13:15 编辑 9 ^3 M2 Z1 E; m1 P2 o
0 {8 a: ]5 F& x( H1 K: ~
n=64;
" J% q" A3 @- X) b2 _0 Lkey=0.512;# V2 g- M9 m& Z$ m9 \$ w
an=linspace(3.56,4,176);" |0 q; U! m1 ~1 m4 x9 {# [7 K
hold on;box on;axis([min(an),max(an),-0.5,1.5]);2 B4 B& n5 a- y9 j
N=n^2;8 w; B' n! i' P8 a9 S: v
xn=zeros(1,N);
" F' M* S& T" K+ j. qfor a1=an;
; V4 I5 a- b* {% l* C; `' N$ t a2=an;
2 q6 `2 k( L3 }% N, ? y=key;
+ k A& |& T8 i6 o# ^' \" Xx=key;
0 [3 u1 Z0 @9 A/ ^for k=1:16;, q8 [5 f3 b6 ] R$ a% ^
y=a1*y*(1-y);
# M7 D* W# Z4 F% X6 i. w& _1 Bx=(a2+(4-a2)*y).*x.*(1-x);%产生公式
2 I# [8 {2 ]& e& [end;5 V, }1 C* ^% F
for k=1:N;- L7 P: `; G* I" c* J* `2 _. L' q
y=a1*y*(1-y);
7 _6 e/ S$ e* D2 `) }1 Hx=(a2+(4-a2)*y).*x.*(1-x);
# Q0 \3 G$ x. S4 d nxn(k,1)=x;
2 S& w+ Y( b; n0 R& bb(k,1)=x;%一维矩阵记录迭代结果8 b( s: z1 t8 `2 q6 l9 B
end;& b0 p+ L$ H5 c% S: f) v
plot(a*ones(1,N),xn,'k.','markersize',1);
- p: F. v- o ?end;, F3 q$ b; r# `4 P
figure;, W! A! b" Z! z6 T# b
imhist(b)1 j A2 Z* ?$ P
2 _; [+ ^( f: T) `7 s( P) Z) H
红色部分有问题
" ~2 d% N3 F9 c0 A" r$ L( u5 E
. \8 W) w% w4 V+ \2 G7 t* S9 G1 g: V* x
|
|