Labfans是一个针对大学生、工程师和科研工作者的技术社区。 | 论坛首页 | 联系我们(Contact Us) |
![]() |
![]() |
#1 |
初级会员
注册日期: 2012-02-03
住址: 天津南开
年龄: 35
帖子: 3
声望力: 0 ![]() |
![]()
本人新手加菜鸟一枚,在网上找的关于洛伦兹程序,但是编辑的时候就不行了,出错了,但是又不知道在什么地方出错了,请大侠们指点一下
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]); |
![]() |
![]() |