TA的每日心情 | 开心 2022-1-21 15:21 |
---|
签到天数: 1 天 [LV.1]初来乍到
|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
本帖最后由 Heaven_1 于 2022-9-27 13:15 编辑
3 a2 ?' N& I; W: p. `) f$ U. ^/ o8 I' A( e3 c
n=64;; k8 `4 c. @. m# i) E+ O/ w9 S
key=0.512;4 p% F# z# q& E0 o3 E' K
an=linspace(3.56,4,176);
% X, q' Y" M1 U2 _+ g# G7 f+ Shold on;box on;axis([min(an),max(an),-0.5,1.5]);6 Y& Y8 S( V/ O3 {" h. j
N=n^2;
+ f# G3 h( c7 n5 W. Exn=zeros(1,N);
' }/ r) E5 u; K4 o* efor a1=an;
9 p2 X/ [" K0 P3 G5 H5 L a2=an;3 z) S t) y, h
y=key;
% N( ]4 N6 ^# b' l1 v& ax=key;/ m- [3 G! {8 T
for k=1:16;
+ Z1 ^ k$ m& J: M" _9 Cy=a1*y*(1-y);
" V( {5 G9 ~- D- Dx=(a2+(4-a2)*y).*x.*(1-x);%产生公式
5 u. R& y, i8 s# t D8 }0 Q( hend;
- \. P- _' m$ N1 f' @6 U+ T% W& Nfor k=1:N;
! {' n" q; J8 C0 W/ K% H, ]( Ly=a1*y*(1-y);( |- Z0 X7 \* H h! V! N7 R
x=(a2+(4-a2)*y).*x.*(1-x);
4 X& k5 `: C; F7 e% b( cxn(k,1)=x;8 L1 K5 l L" U( n1 u0 F
b(k,1)=x;%一维矩阵记录迭代结果
' G# t/ A7 F8 x: Mend;
/ n& k, f8 Q, ?' Y& Y3 Tplot(a*ones(1,N),xn,'k.','markersize',1);
2 H* B" M9 y- g+ n7 W* d& Q+ ^end;- n! V; }1 F+ Y8 h& V
figure;
( y# L3 B, W! X7 B2 \imhist(b)
$ i; p X A# u5 B; U4 Y2 y, z; w* P6 E/ I! K" g k/ S
红色部分有问题/ y( p* o$ V/ `7 |5 G7 }
! Y) B* e7 \3 ^' b: |9 y/ n$ Z. f8 p% W: w6 O) w+ H
|
|