LORENZ吸引子

yunfenger

如下MATLAB程序。。。。MATLAB6.5
( T  k6 x/ v% B7 C+ X8 G2 M+ w%a,b,c 为系统的三个变量，用ODE45数值解算法
9 @/ K! y) T: {0 e& dclear
' W8 l% W2 ~" g6 O( o# [, R% _global a b c
% b5 k; y5 v. i/ Ya=10;b=8/3;c=22;
6 B/ _, i" b. h) W5 F  t) _8 V3 d[T,X]=ode45('Lorenz',[0 100],[10;4;10]);
9 O% ]7 Y6 k. F8 |subplot(2,2,1)
plot(X(:,1),X(:,2))
subplot(2,2,2)
plot(X(:,1),X(:,3))
subplot(2,2,3)
3 _: t/ d: ^6 K) Z: _; [plot(X(:,2),X(:,3))
; f5 {( p. Z" e1 S9 s4 Rsubplot(2,2,4)
; H1 F1 j: t2 I2 l4 a3 Fplot3(X(:,1),X(:,2),X(:,3))
7 T' r# e  a" C! a7 _view(-20,60);
7 ~1 J  }! l- zxlabel('x');ylabel('y');zlabel('z');
% K2 C" v/ P7 M3 m9 H


( Y- j9 k5 ?: u. J, lfunction dx=Lorenz(t,x)
global 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)];

! z6 h( W" s  u/ [
画出来的图形附件！根据不同的参数画出来的图

; W/ z+ {, w( |9 T5 u# X

zhao6884269

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

lilijie_shangha

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

wy371tt1

鉴定完毕！  

fengzhisheng106

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

想做就做

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

547896321

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

mybaby

哦～～

surfw

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

xiaoai

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