Т. е. берётся

, где

— функция Хевисайда, потом сворачивается с какой-то суммой дельт, и потом пропускается через низкочастотный фильтр, что эквивалентно свёртке с

. Видимо, если немного поиграть с преобразованием Фурье или каким-то другим интегральным преобразованием, задача станет попроще (преобразование Фурье переводит свёртку в умножение). Правда, 40 — не очень удобное число для использования БПФ, но это, вроде, как-то обходится.
(Честно говорю, что не разбираюсь в сигналах.)