0
$\begingroup$

I would like to take a time series of several spikes or narrow Gaussians and bandpass using an ift/ifft together with a frequency domain mask to create a series containing band-limited spikes. This is pretty simple but I have an additional desire - that the resulting band-limited spikes are causal - that is each band limited spike has a sharp onset (at the position of the original spike) followed by a tail dependent on the bandpass parameters.

The ultimate goal is a series that looks superficially like a superposition of a multitude of earthquakes.

Is there a simple way to do this?

Thanks very much for any guidance.

Improve this question
$\endgroup$
3
  • 2
    $\begingroup$ Why do you want to use FFT and not an ordinary IIR bandpass filter? IIR is always causal and much cheaper to implement. $\endgroup$
    – Hilmar
    Commented May 15, 2019 at 19:12
  • $\begingroup$ I am certainly open to that but I am not well educated in signal processing although I do have a little experience with ffts. Can you point to a simple reference? $\endgroup$
    – abby yorker
    Commented May 15, 2019 at 19:26
  • $\begingroup$ OK, inspired by the comment above, I looked around a bit and I came across the implementations of constructing and applying a filter, described here: stackoverflow.com/questions/13740348/… I ended up using lfilter rather than filtfilt (which is zero phase). $\endgroup$
    – abby yorker
    Commented May 15, 2019 at 20:21

0

Reset to default

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge you have read our privacy policy.

Browse other questions tagged or ask your own question.