LORENZ吸引子

yunfenger

如下MATLAB程序。。。。MATLAB6.5
%a,b,c 为系统的三个变量，用ODE45数值解算法
8 b: r) |( P7 vclear
global a b c
; f. s/ n& m; e, @1 z- Ka=10;b=8/3;c=22;
4 ]1 Z5 y, a! x9 T: j/ M( v[T,X]=ode45('Lorenz',[0 100],[10;4;10]);
* m0 f) I2 c' [% F) `" Xsubplot(2,2,1)
plot(X(:,1),X(:,2))
( w! @$ O& |3 v7 [( vsubplot(2,2,2)
plot(X(:,1),X(:,3))
: z+ n; A, e/ _' Q* Rsubplot(2,2,3)
plot(X(:,2),X(:,3))
# ^9 k* B. \2 ~* f) C( [9 wsubplot(2,2,4)
plot3(X(:,1),X(:,2),X(:,3))
view(-20,60);
xlabel('x');ylabel('y');zlabel('z');
* m' U  B' L8 \; N9 D1 u' o


function dx=Lorenz(t,x)
global a b c
4 _  O1 h3 X# Y. ydx=[a*(x(2)-x(1));c*x(1)-x(1)*x(3)-x(2);x(1)*x(2)-b*x(3)];

6 M- k2 K- m6 B6 M( D1 w
画出来的图形附件！根据不同的参数画出来的图

( m6 V3 C8 G, ]; \

zhao6884269

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

lilijie_shangha

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

wy371tt1

鉴定完毕！  

fengzhisheng106

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

想做就做

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

547896321

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

mybaby

哦～～

surfw

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

xiaoai

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