LORENZ吸引子

yunfenger

如下MATLAB程序。。。。MATLAB6.5
%a,b,c 为系统的三个变量，用ODE45数值解算法
clear
# D- d, b( I$ J6 n+ V  Mglobal a b c
8 a6 v! K4 S: \. Y/ I" Qa=10;b=8/3;c=22;
" d$ j- O+ ^" k" Z% F[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))
subplot(2,2,3)
plot(X(:,2),X(:,3))
subplot(2,2,4)
9 B( A. I) S* p! o3 A. ~! ^! Uplot3(X(:,1),X(:,2),X(:,3))
view(-20,60);
xlabel('x');ylabel('y');zlabel('z');
2 a: R  `0 Q0 l7 |( J, e


function 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)];
  D. Q( \- }% i
/ i/ j1 W/ h& I+ d
画出来的图形附件！根据不同的参数画出来的图

天马绝尘

试试运气   

hngy

Lorentz混沌这么费体力

灼￡日

程序报错地方多

枫露之茗

鉴定完毕！
* l3 O2 d) k4 f7 O+ b: V$ p3 f

weixinmaths

guohuawin

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

ycy74_2005

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

sigh

强烈支持。楼主万岁

阿里wudi

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

xiaoai

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