LORENZ吸引子

yunfenger

如下MATLAB程序。。。。MATLAB6.5
%a,b,c 为系统的三个变量，用ODE45数值解算法
clear
3 M: O: Q+ }) o1 b+ Y& Q+ D- @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))
. P0 Y# L8 X7 V% P. @/ q$ m2 @2 U- j! Wsubplot(2,2,2)
& t; J7 U1 F+ D2 n! F. l2 `+ J" |6 ~* jplot(X(:,1),X(:,3))
0 [8 D- W. q% {+ ~! ^; Z6 }8 T- [subplot(2,2,3)
plot(X(:,2),X(:,3))
subplot(2,2,4)
5 p% ]7 \1 `6 e$ ?9 ?plot3(X(:,1),X(:,2),X(:,3))
view(-20,60);
xlabel('x');ylabel('y');zlabel('z');



function dx=Lorenz(t,x)
global a b c
$ L9 x9 Q5 i. l8 Sdx=[a*(x(2)-x(1));c*x(1)-x(1)*x(3)-x(2);x(1)*x(2)-b*x(3)];
$ ^; |8 S  ]. R9 m

$ M$ X- E1 V' N5 y1 y3 `2 V画出来的图形附件！根据不同的参数画出来的图
" a* P  r; x# o1 X8 J3 V
/ C1 F. }( ~. t3 x) T

天马绝尘

试试运气   

hngy

Lorentz混沌这么费体力

灼￡日

程序报错地方多

枫露之茗

鉴定完毕！

weixinmaths

guohuawin

声明一下：本人看贴和回贴的规则，好贴必看，精华贴必回。

ycy74_2005

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

sigh

强烈支持。楼主万岁

阿里wudi

楼主，你写得实在是太好了。我惟一能做的，就只有把这个帖子顶上去这件事了

xiaoai

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