LORENZ吸引子

yunfenger

如下MATLAB程序。。。。MATLAB6.5
%a,b,c 为系统的三个变量，用ODE45数值解算法
9 n# x& I3 f7 ]; O& {clear
5 }& ^- A8 X: bglobal a b c
a=10;b=8/3;c=22;
, d) H, Z7 A; Z# Y& `' [8 o[T,X]=ode45('Lorenz',[0 100],[10;4;10]);
subplot(2,2,1)
plot(X(:,1),X(:,2))
* g: A* X$ g' i# k: w. o1 Ssubplot(2,2,2)
plot(X(:,1),X(:,3))
- i' `4 n) o1 Qsubplot(2,2,3)
plot(X(:,2),X(:,3))
6 A6 E3 X2 q" k' K7 N4 Ksubplot(2,2,4)
plot3(X(:,1),X(:,2),X(:,3))
view(-20,60);
& d: J$ _0 {: e1 m: q6 w0 ]xlabel('x');ylabel('y');zlabel('z');

7 R5 n! X0 x, k; R3 B- \! L  R: p

. x3 g+ G/ k) r& W7 Q. `function dx=Lorenz(t,x)
: v0 k) B# e$ I$ o/ e9 qglobal a b c
' I0 h- F; F7 ]- ~# N, ddx=[a*(x(2)-x(1));c*x(1)-x(1)*x(3)-x(2);x(1)*x(2)-b*x(3)];
0 I+ [; o- Z5 @6 j  l# O9 ^
3 q  t+ D/ R1 }; i& Y: o
$ o7 X1 D. _2 ?画出来的图形附件！根据不同的参数画出来的图
. P2 N" `) {- {0 m7 ?* {
( ^6 H! G5 v* V1 I' H; i

zhao6884269

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

lilijie_shangha

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

wy371tt1

鉴定完毕！  

fengzhisheng106

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

想做就做

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

547896321

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

mybaby

哦～～

surfw

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

xiaoai

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