如下MATLAB程序。。。。MATLAB6.5 ' ~9 Y1 b/ L( A%a,b,c 为系统的三个变量,用ODE45数值解算法* F4 X. i. n/ I- w: ~! k& l: v5 r& k) i
clear% m6 ~- `+ Q- u8 n0 F/ S: f7 _# `; `
global a b c; o. E, R7 N5 { Z' P( X+ N
a=10;b=8/3;c=22; ( z( T- F: U: _& |" c% Q[T,X]=ode45('Lorenz',[0 100],[10;4;10]); ) |6 L7 m' s4 K3 G+ lsubplot(2,2,1) " z* Z$ }! W. D7 O6 yplot(X(:,1),X(:,2)) % A) R+ ]/ M' c1 _# N nsubplot(2,2,2)3 R4 p; Z+ w9 R
plot(X(:,1),X(:,3))& u* l, z: h: O8 n
subplot(2,2,3) " d- x2 ~; l0 ^' F) @5 ?) }plot(X(:,2),X(:,3)) . n9 T! Z" I4 U# ^- Ssubplot(2,2,4)2 P- _8 d8 G$ u$ }
plot3(X(:,1),X(:,2),X(:,3))0 {$ Y* ]" \" H t+ `: P
view(-20,60);" L' _5 J7 D5 q6 v6 g
xlabel('x');ylabel('y');zlabel('z');7 L: d* o7 Y' z. w v; d
4 A9 B9 z$ ?& f" h% ` ~$ n# R: F3 W/ t* C2 _
8 m7 d6 Q1 D" {2 Rfunction dx=Lorenz(t,x), Q: a& b5 N* [. W
global a b c5 i r$ Z- J7 c
dx=[a*(x(2)-x(1));c*x(1)-x(1)*x(3)-x(2);x(1)*x(2)-b*x(3)];0 X- h7 m/ R1 d
0 G5 n6 I2 |8 c0 U( R- z9 r! }$ K4 m
画出来的图形附件!根据不同的参数画出来的图& Z% n& \' J$ d `- F- U' q
6 _9 L- P3 \; o