洛仑兹系统对初始条件的敏感依赖性
没事做写着玩的，大家不要见笑洛仑兹方程:
dx/dt=-σ(x-y);
dy/dt=rx-y-xz;
dz/dt=xy-bz;
程序：
function dy = rigid(t,y)
dy = zeros(3,1); % a column vector
u=10;v=29;w=2.67;
dy(1)=-u*(y(1)-y(2));
dy(2)=v*y(1)-y(2)-y(1)*y(3);
dy(3)=y(1)*y(2)-w*y(3);
options = odeset('RelTol',1e-4,'AbsTol',);
= ode45(@rigid,,,options);
= ode45(@rigid,,,options);
= ode45(@rigid,,,options);
plot(T1,Y1(:,1),'r',T2,Y2(:,1),'b',T3,Y3(:,1),'g')
[ 本帖最后由 gghhjj 于 2006-8-26 02:37 编辑 ] 写的还行啊
页:
[1]