LORENZ吸引子

yunfenger

如下MATLAB程序。。。。MATLAB6.5
$ L, c0 W5 Z# P$ b& J%a,b,c 为系统的三个变量，用ODE45数值解算法
1 n! n; G, U; _0 o! L" @$ Q/ nclear
global a b c
- U% u4 q4 r6 e2 e+ T, wa=10;b=8/3;c=22;
& q7 o: Q! h- }) |[T,X]=ode45('Lorenz',[0 100],[10;4;10]);
subplot(2,2,1)
plot(X(:,1),X(:,2))
subplot(2,2,2)
/ [6 G' D, E% Jplot(X(:,1),X(:,3))
+ v6 h0 d& s8 O* Vsubplot(2,2,3)
- E* x# ?8 c8 lplot(X(:,2),X(:,3))
* ]: v$ s+ a! u, n% e6 g8 k% ~subplot(2,2,4)
plot3(X(:,1),X(:,2),X(:,3))
view(-20,60);
xlabel('x');ylabel('y');zlabel('z');


4 p* I8 f! z9 ?% r% R1 A
' c! j0 R- u4 P  u: |* y& Z! l( Z2 |* Wfunction 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)];


画出来的图形附件！根据不同的参数画出来的图

2 w' p, ?& n  ]

zhao6884269

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

lilijie_shangha

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

wy371tt1

鉴定完毕！  

fengzhisheng106

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

想做就做

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

547896321

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

mybaby

哦～～

surfw

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

xiaoai

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