LORENZ吸引子

yunfenger

如下MATLAB程序。。。。MATLAB6.5
4 J; ~- K! h; h%a,b,c 为系统的三个变量，用ODE45数值解算法
+ [: Q8 n0 |( R1 l9 a7 C1 W) fclear
& I& M: V. K3 G( y  d1 [global a b c
a=10;b=8/3;c=22;
[T,X]=ode45('Lorenz',[0 100],[10;4;10]);
subplot(2,2,1)
plot(X(:,1),X(:,2))
subplot(2,2,2)
plot(X(:,1),X(:,3))
4 U5 l. J# A# H) p! rsubplot(2,2,3)
plot(X(:,2),X(:,3))
subplot(2,2,4)
; N( j. b9 T1 D2 O% i! l7 eplot3(X(:,1),X(:,2),X(:,3))
view(-20,60);
0 O& }" ^$ a' \$ Nxlabel('x');ylabel('y');zlabel('z');
: C+ r) ~% H4 `2 G' f5 q! w  I
7 J# C- o: M$ {; D. z1 I

function dx=Lorenz(t,x)
+ E" x' `3 q# \4 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)];
: U6 i1 @0 d: z6 @! h0 T

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

/ L$ ~/ x7 Q* t$ T( a! q" ]

zhao6884269

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

lilijie_shangha

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

wy371tt1

鉴定完毕！  

fengzhisheng106

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

想做就做

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

547896321

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

mybaby

哦～～

surfw

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

xiaoai

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