This is a more general mathematical question, but I am asking this here since it may be a frequent problem and perhaps matlab implementations already exist...
I have a certain impulse response, and a vector of onset times. When I convolve these two, for closely spaced onset times, the amplitudes of the individual impulse responses just add up linearly, which is, however, not the case in the real world.
There are two nonlinear effects, which I would like to add to the simple convolution model:
1. There is an upper limit of signal amplitude, which cannot be exceeded by summing up more and more impulse responses.
2. During the repetition of many impulse responses, a signal "undershoot" is building up, which causes a very long slowly recovering negative signal, that becomes apparent after the positive impulse responses are finished.
I hope this was clear, I would appreciate any ideas, references or pseudo code. Thanks!