How to integrate ordinary differential equations with pulse-like time-varying parameters?

2 次查看（过去 30 天）
Pierre S.
Pierre S. 2014-4-23
评论： RahulTandon 2015-7-8
Problem.
I am simulating an ordinary differential equation with time-varying parameters as follows
p = @(t) p0*( (T0 < t) & (t < T1) )
sol = ode15s(@(t,x) myode(t,x,p(t)),[t0 tf],x0)
where p(t) is a pulse of amplitude p0 and duration (T1-T0) (if T0 < T1).
Because of its adaptive time-step, the integrator "misses the pulse" if the time-step becomes larger that the pulse duration.
Naive solution.
A naive solution would be to constrain the MaxStep to (T1-T0)/2 to be sure that the pulse is detected by the integrator. However, this constrains the MaxStep at time where it is not really needed.
More efficient solution?
I am wondering if there is a more efficient way to do ensure that the pulse is detected.
  4 个评论
显示 2更早的评论
Star Strider
Star Strider 2014-5-1
Posting (or attaching — use the ‘paperclip’ icon) the code for myode would do for a start, as well as information on p0, T0, and T1. We really cannot suggest a solution to a problem we cannot experiment with ourselves.
RahulTandon
RahulTandon 2015-7-8
pleas send the formulae for the two equations 1) The mymode function 2) The pulse train Awaiting an early reply!

请先登录，再进行评论。

请先登录，再回答此问题。

回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Mathematics 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by