LORENZ吸引子

yunfenger

如下MATLAB程序。。。。MATLAB6.5
%a,b,c 为系统的三个变量，用ODE45数值解算法
clear
0 G4 l' }: |' k6 L+ ]5 _0 cglobal a b c
1 o2 g$ q) w" ga=10;b=8/3;c=22;
" A" v4 n6 @; P+ l" A[T,X]=ode45('Lorenz',[0 100],[10;4;10]);
& v+ m8 c) _* l1 j' Tsubplot(2,2,1)
plot(X(:,1),X(:,2))
2 u; q) j1 q5 S* esubplot(2,2,2)
plot(X(:,1),X(:,3))
- V$ z; Z$ M3 ]! \subplot(2,2,3)
6 j$ Z0 I5 P$ A# M+ _' bplot(X(:,2),X(:,3))
% l  M. n& I7 hsubplot(2,2,4)
. z5 m2 x- ]8 x: P9 ^plot3(X(:,1),X(:,2),X(:,3))
3 m4 H; A* }/ C' mview(-20,60);
xlabel('x');ylabel('y');zlabel('z');

9 K) ^+ ?) x5 y, x+ ?
3 k* p* b9 X" U% U- B4 N
function dx=Lorenz(t,x)
global a b c
7 O$ n/ j5 X  ]% Adx=[a*(x(2)-x(1));c*x(1)-x(1)*x(3)-x(2);x(1)*x(2)-b*x(3)];
, p  y3 G2 N- N! J( o5 s% u

, a& n! D  m7 i; O( R, U) a; N0 I% p# z画出来的图形附件！根据不同的参数画出来的图
8 H( O2 Y% P4 ~& L4 F( y0 R1 O
6 p; l% Y  R( h4 U$ G

zhao6884269

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

lilijie_shangha

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

wy371tt1

鉴定完毕！  

fengzhisheng106

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

想做就做

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

547896321

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

mybaby

哦～～

surfw

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

xiaoai

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