Here's some code for calculating the gain and phase of a function in sinusoidal steady state. It works with both time-domain functions with sine/cosine terms and transient terms (By taking the Laplace transform and removing the sinuosidal term in the Laplace domain) and with transfer functions in the Laplace domain, assuming that t -> infinity.
It might make sense to simply roll these functions into some larger package, to avoid excessively fine grained packages.
Suggestion: Post your code to
http://sourceforge.net/apps/phpbb/maxima/viewforum.php?f=3&sid=c18d1b11a001b12002d3139e0c709503
Also, you could append your convolution function to this code.
Not a bug report. Changing the tracker to a feature request.
Dieter Kaiser