LORENZ吸引子

yunfenger

如下MATLAB程序。。。。MATLAB6.5
: U* m9 T5 s3 L. }' |%a,b,c 为系统的三个变量，用ODE45数值解算法
0 l6 h! ?, C- y9 c& vclear
$ m( U& @7 o% k: V$ l9 R1 l& sglobal a b c
a=10;b=8/3;c=22;
[T,X]=ode45('Lorenz',[0 100],[10;4;10]);
; i+ c9 R9 t; q: ]( W- Ysubplot(2,2,1)
! f% ]$ |" K3 n( t) k/ E7 Gplot(X(:,1),X(:,2))
subplot(2,2,2)
plot(X(:,1),X(:,3))
  R2 X9 \: ^, w  D( m  s9 \subplot(2,2,3)
plot(X(:,2),X(:,3))
) \% L3 F2 c8 E' O. ~1 ~subplot(2,2,4)
; a9 p# o* B! W  ~" y* bplot3(X(:,1),X(:,2),X(:,3))
view(-20,60);
- {$ ^' l/ T1 m2 \8 Ixlabel('x');ylabel('y');zlabel('z');
! k6 P5 c: z0 E  _& g. t
; x4 Q% q& X" o: O% h! l2 ~) i

6 J' X0 U3 X2 w$ t% ]% h. B; Zfunction dx=Lorenz(t,x)
global a b c
: s4 m) r, K. [* d; ~+ X* Odx=[a*(x(2)-x(1));c*x(1)-x(1)*x(3)-x(2);x(1)*x(2)-b*x(3)];


画出来的图形附件！根据不同的参数画出来的图
1 B/ a6 F( I) \5 r* m

zhao6884269

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

lilijie_shangha

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

wy371tt1

鉴定完毕！  

fengzhisheng106

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

想做就做

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

547896321

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

mybaby

哦～～

surfw

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

xiaoai

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