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

4 views (last 30 days)
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 Comments
Star Strider
Star Strider on 1 May 2014
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 on 8 Jul 2015
pleas send the formulae for the two equations 1) The mymode function 2) The pulse train Awaiting an early reply!

Sign in to comment.

Answers (0)

Categories

Find more on Mathematics in Help Center and File Exchange

Community Treasure Hunt

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

Start Hunting!