LORENZ吸引子

yunfenger

如下MATLAB程序。。。。MATLAB6.5
%a,b,c 为系统的三个变量，用ODE45数值解算法
clear
global a b c
7 l$ O, R1 ^. H+ ~a=10;b=8/3;c=22;
[T,X]=ode45('Lorenz',[0 100],[10;4;10]);
5 z: d& Z/ I# w  q! r7 Jsubplot(2,2,1)
plot(X(:,1),X(:,2))
subplot(2,2,2)
  p3 @- C; l* Dplot(X(:,1),X(:,3))
subplot(2,2,3)
& q% J+ Z% V# f  V. @- ^plot(X(:,2),X(:,3))
# E5 U! l2 P: Wsubplot(2,2,4)
( w/ r2 A2 w: D! ]plot3(X(:,1),X(:,2),X(:,3))
  m8 S+ L& P, Oview(-20,60);
. Z7 c. a) K, I: c" v; U% Gxlabel('x');ylabel('y');zlabel('z');
/ S$ X) D6 L4 d& v! F  L- d
1 `. n1 @+ _; V, ]9 S

function dx=Lorenz(t,x)
global a b c
: \. ]8 j% Q1 c8 O, f* edx=[a*(x(2)-x(1));c*x(1)-x(1)*x(3)-x(2);x(1)*x(2)-b*x(3)];


0 S! ^1 X3 o$ `! E2 I. f画出来的图形附件！根据不同的参数画出来的图
, H: H3 N/ @4 p
2 z5 c/ H, R( B7 g

zhao6884269

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

lilijie_shangha

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

wy371tt1

鉴定完毕！  

fengzhisheng106

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

想做就做

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

547896321

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

mybaby

哦～～

surfw

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

xiaoai

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