TA的每日心情 | 开心 2022-1-21 15:21 |
---|
签到天数: 1 天 [LV.1]初来乍到
|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
本帖最后由 Heaven_1 于 2022-9-27 13:15 编辑
+ [4 Q* l6 W; N0 b$ R4 [, X/ N8 s$ g* h. Z. M( e/ y! p$ f7 \) s# N
n=64;* A, q) Z0 z' a9 T$ e
key=0.512;" ~6 v. u5 K' \; L
an=linspace(3.56,4,176);3 d L9 s" x6 \) N3 y
hold on;box on;axis([min(an),max(an),-0.5,1.5]);
7 l/ T3 D6 m, M8 i# q5 l3 I. zN=n^2;# b: X: S& S9 d6 ^0 _9 ?/ V/ H# ]$ H
xn=zeros(1,N);( B7 t5 ~7 `- O
for a1=an;0 F+ x( V# o) `+ ?/ |
a2=an;7 Z. A) F, x7 F' r" Y6 D
y=key;
1 e4 D8 t+ ~7 O' @9 v! qx=key;
- s1 r8 M1 s* }: F9 ]# yfor k=1:16;
/ P7 h7 n7 p" dy=a1*y*(1-y);$ N+ f8 o/ U( ?# ?) r0 N
x=(a2+(4-a2)*y).*x.*(1-x);%产生公式& z" D4 ?, @- I4 l% M
end;& l+ \0 t8 G- S) s0 l' K
for k=1:N;
% y% t0 X* b% A. k- Y' s; Z# ~8 E ^5 Uy=a1*y*(1-y);
9 W7 X( [7 b! |" P. `' T" Jx=(a2+(4-a2)*y).*x.*(1-x);" B- N6 p3 O: O# w: p
xn(k,1)=x;& j" o2 o2 C& Y- {9 v
b(k,1)=x;%一维矩阵记录迭代结果
' b# o1 d3 ~+ j' ^( i" l5 B7 Nend;
1 H2 ~9 n/ J5 p6 Rplot(a*ones(1,N),xn,'k.','markersize',1);
& W: L: `7 L5 Iend; U/ b( D* X6 S3 t& H# J0 J
figure;
7 ^0 q0 |) P9 ]3 n4 n, fimhist(b)
- M6 Q; Z- g7 q( E4 `- R; {! Z) O% S4 e2 v& b G
红色部分有问题
( A" B( `5 r' r0 Y8 O
/ M9 [2 H- V. @) S. ~) t8 p
7 k5 g- E" P8 x1 g% E( c |
|