LORENZ吸引子

yunfenger

如下MATLAB程序。。。。MATLAB6.5
%a,b,c 为系统的三个变量，用ODE45数值解算法
# W( c9 y& O# v4 {8 fclear
global a b c
! \/ T4 G: t3 U' @a=10;b=8/3;c=22;
; P' C! i$ T1 Q  A[T,X]=ode45('Lorenz',[0 100],[10;4;10]);
& B/ y& b' Y% W7 qsubplot(2,2,1)
2 ~  U  r' R) Y7 J3 I/ U. B5 a- [6 g1 P$ pplot(X(:,1),X(:,2))
subplot(2,2,2)
1 d- ~) Q/ q  q: L6 T0 Y  \/ }plot(X(:,1),X(:,3))
subplot(2,2,3)
, m5 [3 x8 I* [- I$ Hplot(X(:,2),X(:,3))
' B# }6 z6 W& S" y. ?& z/ R+ Fsubplot(2,2,4)
/ M4 \8 f! E" t, w' P! I; nplot3(X(:,1),X(:,2),X(:,3))
, q" L  F1 g; Oview(-20,60);
  ~( O# E& M$ I; B+ P/ H: N- C6 Xxlabel('x');ylabel('y');zlabel('z');

1 z* }1 Z9 N7 {8 n2 u

* a; f! Z3 D7 i$ T9 Sfunction dx=Lorenz(t,x)
. k5 Z, [4 _9 D2 Bglobal 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)];


画出来的图形附件！根据不同的参数画出来的图
% f4 a- B8 M  i! L0 M

zhao6884269

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

lilijie_shangha

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

wy371tt1

鉴定完毕！  

fengzhisheng106

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

想做就做

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

547896321

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

mybaby

哦～～

surfw

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

xiaoai

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