TA的每日心情 | 开心 2022-1-21 15:21 |
|---|
签到天数: 1 天 [LV.1]初来乍到
|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
本帖最后由 Heaven_1 于 2022-9-27 13:15 编辑
3 @5 V" k- _" D5 l. ^: k" @/ C: J
# T% Y( q! k: P+ vn=64;
$ s3 u- X) L$ |+ {' wkey=0.512;' K9 p5 W9 s' e9 s! Y
an=linspace(3.56,4,176);
+ y3 c6 B* ]! m$ {. ~( D' {hold on;box on;axis([min(an),max(an),-0.5,1.5]);7 \5 \( z6 ]4 h; Y$ b' I
N=n^2;
" B! U, K+ ]+ S8 Q: I8 m7 hxn=zeros(1,N);
/ W. s! ]0 i' u1 efor a1=an;! W* ?% q9 U# p+ I9 K% N2 y1 w
a2=an;/ @) U" j& W4 |! X& `8 P5 x
y=key;5 L, d( n( ~9 i7 Z; S3 d
x=key;
* e/ q% U1 @* x& A# [1 Gfor k=1:16;; p0 E0 z; C7 n
y=a1*y*(1-y);* j1 V7 [+ h9 h* W% [
x=(a2+(4-a2)*y).*x.*(1-x);%产生公式
% H, N9 i1 }8 p% L; ^end;- ` W7 \. p7 v+ ?
for k=1:N;
) D$ j. l; a" ~5 Vy=a1*y*(1-y);4 w9 b- M0 u) ^$ S
x=(a2+(4-a2)*y).*x.*(1-x);7 {3 g C2 u) V, l& }+ w
xn(k,1)=x;
5 ]/ f! O7 }( G; I' x7 [# Gb(k,1)=x;%一维矩阵记录迭代结果
2 f8 `+ p- d+ v0 P5 g3 Oend;9 C" W# I# H' M5 N* G
plot(a*ones(1,N),xn,'k.','markersize',1);
T8 c( B/ V: r, D8 Rend;
1 c4 G& n& e3 R9 z! Nfigure;
- l$ Q8 `" ~+ F3 q' s8 {imhist(b)
2 J0 p2 ^: g. k6 |
4 s7 {2 w- z& S# b, L, Y& j红色部分有问题6 i: i6 H/ X- O, k \2 O1 X; ~. \
6 l x% ~9 k( N. k& [
. w2 U$ I1 E/ ^/ E; \# p) U |
|