TA的每日心情 | 开心 2022-1-21 15:21 |
|---|
签到天数: 1 天 [LV.1]初来乍到
|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
本帖最后由 Heaven_1 于 2022-9-27 13:15 编辑 6 b3 q6 @, }' Z5 z
* X9 W( V3 U. l% ?" An=64;( H7 U5 O6 B' x( |. }) U* W, d0 z
key=0.512;( u* I c, H6 R/ f, O" V
an=linspace(3.56,4,176);
6 Z! v6 y/ t/ A2 l6 q7 b, thold on;box on;axis([min(an),max(an),-0.5,1.5]);7 I, l( L2 a- h4 _9 j, |
N=n^2;/ O2 u2 y: D4 A/ Q6 M
xn=zeros(1,N);
8 b! X1 F2 P# L6 f( x+ Mfor a1=an;
" b( |' g0 s2 G! }6 C9 ^ a2=an;
' @! l s$ S! h. g6 w7 p4 Z y=key;+ x' |2 I; z1 z" ~$ O
x=key;: q# c9 Q' ~: k/ |; Q' G+ v4 i' R2 F
for k=1:16;5 V M: x; m, ]+ l7 P- {/ P
y=a1*y*(1-y);2 \9 | A$ W$ ^2 p9 @; f+ N
x=(a2+(4-a2)*y).*x.*(1-x);%产生公式
5 H0 e+ L* X: ^: y" Mend;* v+ J% i8 E! N3 {8 R( g
for k=1:N;
" m# S) `# O! E$ h; Ny=a1*y*(1-y);
. T( ^& Z, ?4 t0 P4 r) dx=(a2+(4-a2)*y).*x.*(1-x);- a% S. _4 o5 i
xn(k,1)=x;
; p2 L2 u% c X4 c/ P' E9 d f% X/ |b(k,1)=x;%一维矩阵记录迭代结果
8 W6 O) r8 d# d# K) \! j" E: [/ Pend;
5 U: \$ {9 |3 M \! Splot(a*ones(1,N),xn,'k.','markersize',1);
" l9 s8 D. }) e, C" E& xend;
! T2 d7 L+ t: [figure;, ~- q3 g/ ]) b6 ]+ e
imhist(b)* t, W, Z/ v6 u" g2 | l5 S6 w) \! m
0 o3 h& w7 n- L5 v6 l
红色部分有问题
3 |6 R; g' f9 h+ A) g. V& c- D
x& S& p8 H" P0 K" R
& L9 |; A! p. g3 v! U0 M |
|