LORENZ吸引子

yunfenger

如下MATLAB程序。。。。MATLAB6.5
9 I) H+ u8 X1 z" q) C" k%a,b,c 为系统的三个变量，用ODE45数值解算法
7 {+ x% W9 S' u  T: u: _clear
5 X" ]) S. i, G! D* L3 S1 yglobal a b c
- v* p1 z* g( G% O& Ya=10;b=8/3;c=22;
[T,X]=ode45('Lorenz',[0 100],[10;4;10]);
: K8 o# c" V/ o+ e( P% psubplot(2,2,1)
6 G1 F8 U7 T* [plot(X(:,1),X(:,2))
subplot(2,2,2)
plot(X(:,1),X(:,3))
subplot(2,2,3)
9 Z7 h3 q% V4 P' X8 Q( J5 C' b, nplot(X(:,2),X(:,3))
8 l5 ]0 g1 k+ u7 Csubplot(2,2,4)
7 t9 a! ~' }/ B' _8 G" }2 r, ?plot3(X(:,1),X(:,2),X(:,3))
view(-20,60);
xlabel('x');ylabel('y');zlabel('z');


  B; W6 G9 [& @" ^7 N& V
6 [# v/ Y) F) u: R/ ?7 h( w7 Ffunction dx=Lorenz(t,x)
global a b c
1 j' h: x8 {0 o* F0 g6 idx=[a*(x(2)-x(1));c*x(1)-x(1)*x(3)-x(2);x(1)*x(2)-b*x(3)];

1 O$ M" d5 r4 p3 @% D$ Q) I
画出来的图形附件！根据不同的参数画出来的图

+ {* u4 @" o* o; A) K$ D. ?2 B/ s

zhao6884269

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

lilijie_shangha

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

wy371tt1

鉴定完毕！  

fengzhisheng106

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

想做就做

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

547896321

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

mybaby

哦～～

surfw

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

xiaoai

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