马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
function dy=Mathieu(t,x)
dy=zeros(4,1);
dy(1)=x(2);
dy(2)=-x(1)+0.005*x(1)^3+1/8*(-x(2)+4*x(2)^3+sin(t));
dy(3)=x(4);
dy(4)=(1+0.015*x(1)^2)*x(3)+1/8*(-1+12*x(2)^2)*x(4);
tspan=linspace(1,1000,6000);
[t,y]=ode45('Mathieu',tspan,[0.3 -0.15 0.3 -0.15]);
plot(t(5000,10000),y(5000,10000))
??? Index exceeds matrix dimensions.
Error in ==> shijianbianli2 at 5
plot(t(5000,10000),y(5000,10000))
出现这个错误怎么改呢?请高手指教一下,我快郁闷死了,怎么改也不对
[ 本帖最后由 sigma665 于 2008-7-21 23:40 编辑 ] |