TA的每日心情 | 开心 2022-1-21 15:21 |
|---|
签到天数: 1 天 [LV.1]初来乍到
|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
本帖最后由 Heaven_1 于 2022-9-27 13:15 编辑
& l6 X+ r1 o9 m$ H9 Q/ r% ]" w0 J5 I) d' U
n=64;
; Y! S+ C! v$ r; Gkey=0.512;7 B3 p O% k4 ~5 q$ d* w# t
an=linspace(3.56,4,176);+ F/ G/ e5 A; c2 u* I. t0 s" R
hold on;box on;axis([min(an),max(an),-0.5,1.5]);5 y T5 F1 q7 K( c2 Q4 a: G) v
N=n^2;4 a W9 ?. E7 |$ v. |
xn=zeros(1,N);3 E" I3 n) V0 x
for a1=an;
( j6 G* h- K# c) s5 y a2=an;
' P( ?& F& w$ ` B1 V( R y=key;7 V% _ u9 x' r) o' a$ j8 @# k
x=key;
: c& I5 t2 y: L% ~for k=1:16;
( O% S$ a5 U$ d8 A2 u" u, sy=a1*y*(1-y);
4 a R& B' t& R; j- [x=(a2+(4-a2)*y).*x.*(1-x);%产生公式; F5 _* h: f0 g1 n& q
end;1 {: a: J" N1 V! \" k; W) v
for k=1:N;% s- l) e3 r0 x# N3 h
y=a1*y*(1-y);) \( n7 R" u5 r" A
x=(a2+(4-a2)*y).*x.*(1-x);
V1 o: b. w% E) x8 c: Y. }xn(k,1)=x;, `' O6 e6 L, K9 @0 I& r
b(k,1)=x;%一维矩阵记录迭代结果/ `- N' ~9 d+ O, m* T9 J) w& X
end;
0 w+ X3 J# O+ l& V/ fplot(a*ones(1,N),xn,'k.','markersize',1);
/ t- g7 K6 p' K& tend;5 }7 b9 f. {* ~) \9 n- m' V2 c& I) V
figure;
3 [( P2 k1 J; n/ q- B# |# Timhist(b)' _ r: r1 ^) \/ o
8 B) t% R* ~5 O- K红色部分有问题+ h) L% q+ l3 j3 U6 u. B& Q7 P* L0 M
2 w9 [" s( T9 ~. Q* {1 N7 J y( y% u* A, @5 h
|
|