En realidad, la convolución en el filtrado lineal de señales de tiempo discreto en el dominio del tiempo (como SIR) no se calcula tal cual. Esta operación requiere conocer la señal de entrada para cualquier instante temporal, y eso es técnicamente imposible. Se trabaja siempre con convoluciones circulares de forma vectorial (algoritmos overlap-add, overlap-save y variantes) que paso a paso van reconstruyendo la salida dentro de una ventana de tiempo predefinida en el DSP, o bien, operando matricialmente con los coeficientes de sus transformadas DFT.
El tema de las convoluciones es algo bastante antiguo. Hoy en día pocos DSPs trabajan con ellas en el dominio del tiempo. Tanto DivX Mp3, Jpeg, como otros miles de formatos y algoritmos de DSPing aprovechan las las limitaciones de nuestra percepción, y características de la señal en los dominios transformados.
Si te interesa el tema, hay mucha bibliografía escrita y te puedo recomendar páginas o libros como biblias sobre el asunto.