LORENZ吸引子

yunfenger

如下MATLAB程序。。。。MATLAB6.5
%a,b,c 为系统的三个变量，用ODE45数值解算法
clear
global a b c
: X" l8 c$ |: i1 Y: qa=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))
subplot(2,2,3)
plot(X(:,2),X(:,3))
# }- ?! m" z4 |3 Qsubplot(2,2,4)
' j4 ?; m- w- D: ?$ w7 x+ @plot3(X(:,1),X(:,2),X(:,3))
5 x/ \6 E8 m  r) F! _view(-20,60);
: S( Z. m9 ~. n1 h- K: }xlabel('x');ylabel('y');zlabel('z');


, }( c! ?! ^3 m6 X3 f1 f
) ?0 X+ R3 ^) g0 E4 o2 M, n( j% ?function dx=Lorenz(t,x)
; c) f% ?: ~! ^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)];
# ]( r0 a% a+ g4 B* T

/ }* W5 q/ J$ s4 F- }画出来的图形附件！根据不同的参数画出来的图

zhao6884269

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

lilijie_shangha

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

wy371tt1

鉴定完毕！  

fengzhisheng106

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

想做就做

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

547896321

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

mybaby

哦～～

surfw

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

xiaoai

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