如何将常微分方程与脉冲式时变参数相结合？

Question

问题。

我正在模拟具有时变参数的常微分方程如下

p   = @(t) p0*( (T0 < t) & (t < T1) )
sol = ode15s(@(t,x) myode(t,x,p(t)),[t0 tf],x0)

其中p（t）是幅度p0和持续时间（T1-T0）的脉冲（如果T0

由于其自适应时间步长，积分器＆＃34;错过脉冲＆＃34;如果时间步长变得大于脉冲持续时间。

天真的解决方案。

一个天真的解决方案是将MaxStep约束到（T1-T0）/ 2，以确保积分器检测到脉冲。但是，这会限制MaxStep在不需要的时候。

更有效的解决方案？

我想知道是否有更有效的方法可以确保检测到脉冲。