LORENZ吸引子

yunfenger

如下MATLAB程序。。。。MATLAB6.5
5 F- X, P5 e0 _2 M. V' Y4 L( i%a,b,c 为系统的三个变量，用ODE45数值解算法
. a# ?' R. ?& H; I6 S  Eclear
' Y# X) Y! S! ?9 ^% ~4 Oglobal a b c
3 I) \2 {( c1 Q6 r+ ^! ea=10;b=8/3;c=22;
- y, Z6 Z$ t* P' C[T,X]=ode45('Lorenz',[0 100],[10;4;10]);
9 `: M! K# W6 f$ d  R3 C' Dsubplot(2,2,1)
3 l( z0 \: |& ^9 {  y8 @plot(X(:,1),X(:,2))
# v3 m  \( Y3 n* ]5 b6 b3 ?2 l7 Wsubplot(2,2,2)
$ ~% O+ _% u/ f$ Iplot(X(:,1),X(:,3))
subplot(2,2,3)
plot(X(:,2),X(:,3))
subplot(2,2,4)
plot3(X(:,1),X(:,2),X(:,3))
0 u( P3 }9 U* G- fview(-20,60);
+ P/ V- J4 o# O) _% Txlabel('x');ylabel('y');zlabel('z');


, L0 h" }+ M: h
9 l( c" g. m- a0 y7 v% Pfunction dx=Lorenz(t,x)
5 L$ e3 }0 l0 t1 b' R6 l3 Lglobal a b c
dx=[a*(x(2)-x(1));c*x(1)-x(1)*x(3)-x(2);x(1)*x(2)-b*x(3)];
9 j! r- p" r) l9 t
2 @7 p. ?+ t) s# E
画出来的图形附件！根据不同的参数画出来的图
7 s, E4 Q6 f- j' O' M+ o0 }

zhao6884269

试试运气啦~~~~~~~~~~~

lilijie_shangha

来报道!!!!!!!!!!!

wy371tt1

鉴定完毕！  

fengzhisheng106

强人，佩服死了。呵呵，不错啊

想做就做

强人，佩服死了。呵呵，不错啊

547896321

我来了~~~~~~~~~ 闪人~~~~~~~~~~~~~~~~  

mybaby

哦～～

surfw

试试运气啦~~~~~~~~~~~

xiaoai

强人，佩服死了。呵呵，不错啊