洛仑兹系统对初始条件的敏感依赖性

gghhjj

没事做写着玩的，大家不要见笑

洛仑兹方程:

dx/dt=-σ(x-y);
dy/dt=rx-y-xz;
dz/dt=xy-bz;

程序：

1. function dy = rigid(t,y)
   
2. dy = zeros(3,1);    % a column vector
   
3. u=10;v=29;w=2.67;
   
4. 
   
5. dy(1)=-u*(y(1)-y(2));
   
6. dy(2)=v*y(1)-y(2)-y(1)*y(3);
   
7. dy(3)=y(1)*y(2)-w*y(3);

复制代码

1. options = odeset('RelTol',1e-4,'AbsTol',[1e-4 1e-4 1e-5]);
   
2. [T1,Y1] = ode45(@rigid,[0 20],[12 2 9],options);
   
3. [T2,Y2] = ode45(@rigid,[0 20],[12 2.01 9],options);
   
4. [T3,Y3] = ode45(@rigid,[0 20],[12 1.99 9],options);
   
5. plot(T1,Y1(:,1),'r',T2,Y2(:,1),'b',T3,Y3(:,1),'g')

复制代码

[ 本帖最后由 gghhjj 于 2006-8-26 02:37 编辑 ]

ieluxinhua

写的还行啊