TA的每日心情 | 开心 2022-1-21 15:21 |
|---|
签到天数: 1 天 [LV.1]初来乍到
|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
本帖最后由 Heaven_1 于 2022-9-27 13:15 编辑 8 |' Y& b4 a" \. }* o
8 W$ A6 V9 Z" I9 t+ tn=64;
; C7 D+ n7 i& i- Kkey=0.512;
! ]! A) x" K1 H( v0 g/ man=linspace(3.56,4,176); B* ~9 ^7 @ c. B0 D2 Q6 w$ X
hold on;box on;axis([min(an),max(an),-0.5,1.5]);" V' E. n2 d( R+ E, g
N=n^2;0 R8 v( T! H _! \/ P0 k
xn=zeros(1,N);- J+ C6 n9 @* h3 p9 G2 o% u* q
for a1=an;- I' y8 ]6 ?' q3 S6 \, q j
a2=an;3 F8 n* m4 f$ u* V" U5 H
y=key;
8 k# _" P) ?9 h( G' ]& `1 a5 \0 ax=key;
' N, t9 D6 a" Yfor k=1:16;% E& K( ?3 p: Z2 a6 F
y=a1*y*(1-y);
. h; z7 R$ D: @" J4 G2 x J$ x# Rx=(a2+(4-a2)*y).*x.*(1-x);%产生公式
& P g# a+ o) m. z1 a( pend;
C/ F# z# c& N- E7 U" Jfor k=1:N;
! Z* ~- C2 E3 m7 p- {/ by=a1*y*(1-y);; [4 m) b0 u6 f3 {0 ?
x=(a2+(4-a2)*y).*x.*(1-x);5 h# p) n }0 [6 v" [% g" Q# n
xn(k,1)=x;+ e. `. y+ A9 w& U! _; `7 R
b(k,1)=x;%一维矩阵记录迭代结果( J$ P7 h0 p' E% V3 b5 |# U3 x
end;& x/ ?7 a& R P7 y
plot(a*ones(1,N),xn,'k.','markersize',1);. G( v6 H! `+ P5 e4 T& h
end;, d1 q0 d2 H6 U3 `* ?" l5 }
figure;% T! P( Z4 }8 L( d* U; |
imhist(b)% K$ S2 o) F# {( z0 E. M8 s) O$ f
% W6 M: |1 O! o- D红色部分有问题/ [( V* g0 I1 i) N9 |- X; W( K1 U; `
' X* J2 w4 z5 y
6 U5 j4 B. f( x- A
|
|