LORENZ吸引子

yunfenger

如下MATLAB程序。。。。MATLAB6.5
* h! Y( b" j1 V6 D%a,b,c 为系统的三个变量，用ODE45数值解算法
clear
' V" B( {! B- A" {5 p" |; C* Qglobal a b c
' Q1 P5 ^4 U4 _4 ^1 ra=10;b=8/3;c=22;
[T,X]=ode45('Lorenz',[0 100],[10;4;10]);
subplot(2,2,1)
plot(X(:,1),X(:,2))
1 _9 U  d$ C5 x$ `subplot(2,2,2)
# D$ S0 ~- s* Q5 ]plot(X(:,1),X(:,3))
9 A6 y& n3 _5 Ysubplot(2,2,3)
plot(X(:,2),X(:,3))
subplot(2,2,4)
plot3(X(:,1),X(:,2),X(:,3))
( E+ |8 ^/ N' gview(-20,60);
0 Y9 L9 y) N. B7 u" a3 M8 ]xlabel('x');ylabel('y');zlabel('z');
. C: i5 c+ S2 [; X! `


function dx=Lorenz(t,x)
3 G4 u8 ]6 W* L1 Z0 |+ V# dglobal 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)];

9 y! N- ~$ [$ y5 E, y
4 j; a$ M6 r$ |$ ^' p' N. G画出来的图形附件！根据不同的参数画出来的图

6 b" T5 W8 I: w( b0 A! t

zhao6884269

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

lilijie_shangha

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

wy371tt1

鉴定完毕！  

fengzhisheng106

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

想做就做

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

547896321

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

mybaby

哦～～

surfw

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

xiaoai

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