warden
2012-02-04, 11:19
本人新手加菜鸟一枚,在网上找的关于洛伦兹程序,但是编辑的时候就不行了,出错了,但是又不知道在什么地方出错了,请大侠们指点一下
function dy=Lorenz(t,y)
dy=zeros(3,1);
dy(1)=10*(-y(1)+y(2));
dy(2)=28*y(1)-y(2)-y(1)*y(3);
dy(3)=y(1)*y(2)-8*y(3)/3;
[t,y]=ode45('Lorenz',[0,30],[12,2,9]);
plot(t,y(:,1))
pause
plot(t,y(:,2))
pause
plot(t,y(:,3))
pause
plot3(y(:,1),y(:,2),y(:,3));
view([20,42]);
function dy=Lorenz(t,y)
dy=zeros(3,1);
dy(1)=10*(-y(1)+y(2));
dy(2)=28*y(1)-y(2)-y(1)*y(3);
dy(3)=y(1)*y(2)-8*y(3)/3;
[t,y]=ode45('Lorenz',[0,30],[12,2,9]);
plot(t,y(:,1))
pause
plot(t,y(:,2))
pause
plot(t,y(:,3))
pause
plot3(y(:,1),y(:,2),y(:,3));
view([20,42]);