Question

matlab ode45为什么显示错误

dxy=@(t,x)[-0.0106*x(1)+5400*(t>=0&t<1)+600*(t>=2&t<3)+13000*(t>=5&t<6)-0.0544*x(2)];
[t,xy]=ode45(dxy,[0:36],[0,21500]);
subplot(211),polt(t,xy(:,1),'r*',t,xy(:,2),'gD')
xlabel('时间t'),ylabel('人数'),legend('美军','日军')
subplot(212),plot(xy(:,1),xy(:,2))
xlabel('美军人数x'),ylabel('日军人数y')
错误使用 odearguments (line 95)
@(T,X)[-0.0106*X(1)+5400*(T>=0&T<1)+600*(T>=2&T<3)+13000*(T>=5&T<6)-0.0544*X(2)] 返回的向量的长度为 1，但初始条件向量的长度为
2。@(T,X)[-0.0106*X(1)+5400*(T>=0&T<1)+600*(T>=2&T<3)+13000*(T>=5&T<6)-0.0544*X(2)] 返回的向量和初始条件向量的元素数目必须相同。

Answer 1

错误的原因，有

1、自定义函数的格式不对。应该这样来写

dxy=@(t,x)[(-0.0106*x(1)+5400).*(t>=0&t<1)+600.*(t>=2&t<3)+13000.*(t>=5&t<6);-0.0544*x(2)];

2、polt(t,xy(:,1),'r*',t,xy(:,2),'gD')这个命令，书写格式也存在问题。应该这样来写

plot(t,xy(:,1),'r*',t,xy(:,2),'g-')

更改上述错误后，运行可得到如下图形。