I obtain numerical discrete data of the form

$$ S_{raw}(\omega) = \sum_{j}w_{j} \delta(\omega-\omega_{j}) $$

to compare the result with experimental data the delta peaks need to be broadened according to

$$ S_{broad}(\omega) = \int d\omega^{\prime} K(\omega,\omega^{\prime}) S_{raw}(\omega^{\prime}) $$

where $$ K(\omega, \omega^{\prime})$$ is the broadening kernel i.e.

$$ K(\omega,\omega^{\prime}) = \frac{\theta(\omega\omega^{\prime})}{\sqrt{\pi}b|\omega|}\exp\left[-\left(\frac{1}{b}\ln\left|\frac{\omega}{\omega^{\prime}}\right| - \frac{b}{4}\right)^2\right]. $$

If $$K(\omega,\omega^{\prime}) = K(\omega-\omega^{\prime})$$ the integral equation above can be interpreted as a convolution and can be efficiently calculated by means of the convolution theorem and FFTs. But the example kernel $K$ does not have this property.

I ask me now if there is still a way to compute the broadened function by exploiting the convolution theorem and fast FFT algorithms ?


Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Browse other questions tagged or ask your own question.