作为Matlab的初学者,我正在用test.m探索ode45函数,其中包含要求解的ode方程。
在编辑器中,我通过两个选项分别调用此函数。在一个案例中,我定义了选项1
t=[0 50];
y0=[0 2];
[t,y]=ode45(@(t,y)test(t,y),t,y0);
在选项2中,我要求一个结构输出
t=[0 50];
y0=[0 2];
sol=ode45(@(t,y)test(t,y),t,y0);
然而,与选项1相比,结构选项2得到的时间步长更小,因此我的图是"0";更粗糙";。
我找不到增加步骤数量的方法来完善选项2的解决方案……有什么想法吗?
将deval
函数与sol
和您选择的时间点一起使用。例如,
tSpan = 0:0.01:50;
y = deval(sol, tSpan);
plot(tSpan,y)
请参阅文档中的评估和扩展解决方案结构