TA的每日心情 | 开心 2022-1-21 15:21 |
|---|
签到天数: 1 天 [LV.1]初来乍到
|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
本帖最后由 Heaven_1 于 2022-9-27 13:15 编辑 8 s& n& j. ^1 y% z- C
" F) @( U, k% j/ \0 f
n=64;! N( S* T( G0 T _9 v. O
key=0.512;
8 k& h2 o! B2 Zan=linspace(3.56,4,176);" T6 L9 F- {0 @( }$ l
hold on;box on;axis([min(an),max(an),-0.5,1.5]);' x4 ~$ q: O! G
N=n^2;
# }: i5 B3 J1 q# h4 Fxn=zeros(1,N);
4 U( T9 A7 R7 I: e0 u) Z2 Ofor a1=an;1 l! C+ h+ S Q0 q# G- _
a2=an;( k6 h; \3 z2 c# f0 \
y=key;, a" Q, E& ?1 ?; Q U# u3 U
x=key;% ]2 ~- `8 L9 T* q; ?
for k=1:16;
( s$ o9 {7 f, q) @; w& Xy=a1*y*(1-y);
+ D$ ] e1 y N3 ^4 S8 sx=(a2+(4-a2)*y).*x.*(1-x);%产生公式
- O9 a/ J) {& ]) C8 [end;
) Q n2 a- n$ r) E6 Z4 D) q8 o. wfor k=1:N;
/ ]+ \& a8 @: ` K6 {, V/ Dy=a1*y*(1-y);
( ]& u/ |8 C. ~' r/ ^x=(a2+(4-a2)*y).*x.*(1-x);, H. Z2 m6 a0 p- y2 g
xn(k,1)=x;
+ b2 w9 B) } X% f5 v) jb(k,1)=x;%一维矩阵记录迭代结果
9 e) i% _5 z9 T1 `4 x# Nend;
, `% s1 U( }/ v) n, C. Tplot(a*ones(1,N),xn,'k.','markersize',1);
4 r# T' I+ o* Mend;6 o7 v' t5 S9 U. k* G7 _
figure;1 ?, F& e+ T' k
imhist(b)* _. }3 n. z8 x
1 o+ U1 i: o" a) X% a+ C L5 ^# d红色部分有问题3 {2 _+ a! H* k0 [
( K1 O/ y' v* t5 \' c% Y- L
& x; V5 g2 W) m- F* v4 x |
|