pulse_response() [1/2]

MAYAFLUX_API double MayaFlux::Kinesis::Discrete::pulse_response ( double  t, double  onset, double  duration, double  rise, double  fall, double  period = 0.0  )

noexcept

Level at a time of a leaky response to a rectangular pulse.

Parameters

t	Time
onset	Time the pulse begins
duration	Length of the pulse, held to the period when one is given
rise	Time constant while the pulse is on; zero or less jumps to one
fall	Time constant after the pulse ends; zero or less falls to zero at once
period	Repeat interval; zero or less means a single pulse

Returns

Level in [0, 1]

During the pulse the level climbs toward one as 1 - exp(-s / rise), with s the time since the pulse began. After it ends the level decays from where it stood as exp(-s / fall). Before the onset it is zero. With a period, the pattern repeats from rest every period, so a tail still decaying when the next pulse begins is cut off there.

Definition at line 7 of file Dynamics.cpp.

{
    double local = t - onset;
    if (local < 0.0)
        return 0.0;
 
    double pulse = std::max(duration, 0.0);
    if (period > 0.0) {
        local -= period * std::floor(local / period);
        pulse = std::min(pulse, period);
    }
 
    const auto peak_at = [rise](double s) {
        return rise > 0.0 ? 1.0 - std::exp(-s / rise) : 1.0;
    };
 
    if (local <= pulse)
        return peak_at(local);
 
    if (fall <= 0.0)
        return 0.0;
    return peak_at(pulse) * std::exp(-(local - pulse) / fall);
}

Referenced by pulse_response().

Here is the caller graph for this function: