LORENZ吸引子

yunfenger

如下MATLAB程序。。。。MATLAB6.5
8 O; V% O1 s1 y- i+ d) A%a,b,c 为系统的三个变量，用ODE45数值解算法
clear
8 a# |6 W7 ?+ [1 X* Nglobal a b c
7 X7 C2 C* y* t8 p( ?a=10;b=8/3;c=22;
[T,X]=ode45('Lorenz',[0 100],[10;4;10]);
$ A$ ]) T- V8 O" f7 zsubplot(2,2,1)
plot(X(:,1),X(:,2))
( m5 F. i# E) b( z& Gsubplot(2,2,2)
  ^) q! C: ]: c$ B8 Z4 ]" i  B) L9 w8 zplot(X(:,1),X(:,3))
+ x" i0 |6 @, c7 F: d8 ?" A) Fsubplot(2,2,3)
+ N3 V( z) X6 J3 v5 yplot(X(:,2),X(:,3))
5 l$ K: ^3 v$ b/ a1 Q- W2 hsubplot(2,2,4)
plot3(X(:,1),X(:,2),X(:,3))
0 x6 [- s& a; J: k2 Gview(-20,60);
" P0 `( {8 t- ~+ Q: e4 m8 D4 Hxlabel('x');ylabel('y');zlabel('z');

6 \, G7 w4 O& v! t

6 j: V2 O) m4 Z0 ^* B# Gfunction dx=Lorenz(t,x)
global a b c
0 {) A  S: \: t3 D0 ~; kdx=[a*(x(2)-x(1));c*x(1)-x(1)*x(3)-x(2);x(1)*x(2)-b*x(3)];
$ ~- |, \. @- x7 J1 r: `& ~5 `
# ?) k8 A! j4 J. r) _
画出来的图形附件！根据不同的参数画出来的图
/ G; ^! \% `; T6 P  \
! n/ k/ Q- ?1 h" V0 U% ~5 A6 A

zhao6884269

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

lilijie_shangha

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

wy371tt1

鉴定完毕！  

fengzhisheng106

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

想做就做

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

547896321

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

mybaby

哦～～

surfw

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

xiaoai

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