LORENZ吸引子

yunfenger

如下MATLAB程序。。。。MATLAB6.5
0 l  {8 K/ P+ m/ }% c%a,b,c 为系统的三个变量，用ODE45数值解算法
# D* O& x3 s0 m+ S" nclear
4 i0 k: d! e) V1 w1 j8 G" Aglobal a b c
  g0 V% e' O" f: [# Ha=10;b=8/3;c=22;
[T,X]=ode45('Lorenz',[0 100],[10;4;10]);
subplot(2,2,1)
- m1 V" u5 w5 o: e. @  y4 r% \plot(X(:,1),X(:,2))
. Y& y  U, k) j4 P2 J% C' _subplot(2,2,2)
plot(X(:,1),X(:,3))
$ q& J9 \* `8 z, U, u$ q5 ]subplot(2,2,3)
plot(X(:,2),X(:,3))
subplot(2,2,4)
- a; ?- Y7 R. y( Y& q4 I* y5 Xplot3(X(:,1),X(:,2),X(:,3))
view(-20,60);
# C6 z5 Q, @% X& X9 nxlabel('x');ylabel('y');zlabel('z');

/ w% r  k/ X6 j  ?8 r

3 ^$ m3 d! Q: @0 ]: Rfunction dx=Lorenz(t,x)
/ [( a' d0 e/ n! O* T( Gglobal 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)];
) T& X4 V+ k$ H+ @! a; w
% R9 J' }! _* x  W* M% \
, B6 z- B7 B9 s- J* R9 \画出来的图形附件！根据不同的参数画出来的图

5 S2 e6 K2 `3 c& U( r

zhao6884269

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

lilijie_shangha

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

wy371tt1

鉴定完毕！  

fengzhisheng106

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

想做就做

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

547896321

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

mybaby

哦～～

surfw

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

xiaoai

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