LORENZ吸引子

yunfenger

如下MATLAB程序。。。。MATLAB6.5
%a,b,c 为系统的三个变量，用ODE45数值解算法
clear
" p# D5 J: a4 a* ?+ {# [2 [9 cglobal a b c
: v5 c& h3 o  U2 b1 na=10;b=8/3;c=22;
! Q$ \7 X/ q8 R) h' ^- d' a, x[T,X]=ode45('Lorenz',[0 100],[10;4;10]);
subplot(2,2,1)
plot(X(:,1),X(:,2))
3 g- ~4 v) f) j0 L+ Fsubplot(2,2,2)
" r& k& `9 f/ x. W. gplot(X(:,1),X(:,3))
! r. Q. b( ]  j& ]1 c7 I3 I; isubplot(2,2,3)
plot(X(:,2),X(:,3))
subplot(2,2,4)
% J  A# f6 V  a( c/ D. L" mplot3(X(:,1),X(:,2),X(:,3))
  D, _5 w0 b$ [% s4 p8 _* tview(-20,60);
xlabel('x');ylabel('y');zlabel('z');



function dx=Lorenz(t,x)
$ p1 J( j/ a6 Bglobal a b c
- O0 K2 t/ g$ r: ldx=[a*(x(2)-x(1));c*x(1)-x(1)*x(3)-x(2);x(1)*x(2)-b*x(3)];
  J% J: L, I5 Z: ~: q# n; ^: O

' [1 N' S8 k+ U6 }; E5 a画出来的图形附件！根据不同的参数画出来的图

. ]' _% u1 d) F$ N3 \

zhao6884269

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

lilijie_shangha

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

wy371tt1

鉴定完毕！  

fengzhisheng106

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

想做就做

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

547896321

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

mybaby

哦～～

surfw

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

xiaoai

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