Reference implementation 'impulse invariance transform'

Octave currently lacks a functions to perform the
transformation of a continuous time system to a
discrete time system and back, using the 'impulse
invariance transform'. The impulse invariance
transform is a key technique for digital filter and
dsp design. See also
http://www.mathworks.com/access/helpdesk/help/toolbox/signal/impinvar.html.

If anyone is interested, I have written a reference
implementation for this function that is compatible
with Octave. The implementation covers both continuous
to discrete and the reverse transform. It also handles
multiple order poles correctly (see
http://ieeexplore.ieee.org/search/srchabstract.jsp?arnumber=536691&isnumber=11401&punumber=78&k2dockey=536691@ieeejrns).

Let me know if there is any interest, and I will
gladly release under GPL.

Rudy.

__________________________________________________
Do You Yahoo!?
Tired of spam?  Yahoo! Mail has the best spam protection around 
http://mail.yahoo.com